正在搞这个
下载看看
楼主你的问题解决了没有？我也遇到你的问题了。。。。
楼主太厉害了！
zhlei566 发表于 2014-5-11 10:10
常规设计，齿数比一般不为整数，这样的好处是齿轮在运转过程中每一对齿都有相互啮合的机会，进而磨损较为 ...
楼主回答的真好
看各位的讨论学到不少东西哈大家加油
原来不是我一个人啊哈哈，也不知道楼主解决了么？还希望不惜赐教啊。那我近期打算做newmark试试。有消息分享给大家。
到现在我的程序换了不同的参数运行了一下午和一晚上仍旧没有任何反映，真的是我的计算方法选择错了么？ode45跑着一个，ode15s跑着一个。。。
Edinburgh 发表于 2015-11-10 09:15
个人理解楼主的问题其实不是很明确
长时间运行不知道是否能够得到结果
很可能是由于模型的收敛性差造成 ...
看你们的讨论学到不少东西
谢谢分享
ode是专门用于解微分方程的功能函数，他有ode23,ode45,ode23s等等，采用的是Runge-Kutta算法。ode45表示采用四阶，五阶runge-kutta单步算法,截断误差为(Δx)^3。解决的是Nonstiff(非刚性)的常微分方程.是解决数值解问题的首选方法,若长时间没结果,应该就是刚性的,换用ode23来解.其他几个也是类似的用法 使用方法 = ode45(odefun,tspan,y0) odefun 是函数句柄,可以是函数文件名,匿名函数句柄或内联函数名
tspan 是区间
或者一系列散点
y0 是初始值向量
T 返回列向量的时间点
Y 返回对应T的求解列向量 = ode45(odefun,tspan,y0,options)
options 是求解参数设置,可以用odeset在计算前设定误差,输出参数,事件等 =ode45(odefun,tspan,y0,options) 每组(t,Y)之产生称为事件函数。每次均会检查是否函数等于零。并决定是否在零时终止运算。这可以在函数中之特性上设定。例如以events 或@events产生一函数。 =events(t,y)其中，value(i)为函数之值，isterminal(i)=1时运算在等于零时停止，=0时继续；direction(i)=0时所有零时均需计算（默认值）， +1在事件函数增加时等于零， -1在事件函数减少时等于零等状况。此外，TE, YE, IE则分别为事件发生之时间，事件发生时之答案及事件函数消失时之指针i。
sol =ode45(odefun,,y0...)
sol 结构体输出结果 应用举例 1 求解一阶常微分方程
程序:
odefun=@(t,y) (y+3*t)/t^2; %定义函数
tspan=; %求解区间
y0=-2; %初值
=ode45(odefun,tspan,y0);
plot(t,y) %作图
title('t^2y''=y+3t,y(1)=-2,1<t<4')
legend('t^2y''=y+3t') xlabel('t')
ylabel('y') % 精确解
% dsolve('t^2*Dy=y+3*t','y(1)=-2')
% ans =
% (3*Ei(1) - 2*exp(1))/exp(1/t) - (3*Ei(1/t))/exp(1/t)
2 求解高阶常微分方程
关键是将高阶转为一阶,odefun的书写.
F(y,y',y''...y(n-1),t)=0用变量替换,y1=y,y2=y'...注意odefun方程定义为列向量
dxdy=
程序:
function Testode45
tspan=; %求解区间
y0=; %初值
=ode45(@odefun,tspan,y0);
plot(t,x(:,1),'-o',t,x(:,2),'-*')
legend('y1','y2')
title('y'' ''=-t*y + e^t*y'' +3sin2t')
xlabel('t') ylabel('y')
function y=odefun(t,x)
y=zeros(2,1); % 列向量
y(1)=x(2);
y(2)=-t*x(1)+exp(t)*x(2)+3*sin(2*t);
end
end
为什么我换参数就不收敛了啊
计算的时间长是程序的问题还是硬件的问题
学习学习{:{39}:}
程序里也没有循环运算速度怎么会慢呢？