| 查看: 313 | 回复: 2 | |||
[交流]
如何设置使程序运行下去呢?
|
|
【 Failure at t=5.422028e-006. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.355253e-020) at time t.】针对这一句，可否人为的设置，使程序能够继续运行下去呢？ 比如减小step size？ 或增大integration tolerances？  我要t知道360，现在连0.1都不到啊！！！！！  步长太小.jpg |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有4人回复
-
基金委咋了？2026年的指南还没有出来？
已经有9人回复
-
基金申报
已经有5人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有5人回复
-
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
推荐一本书
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- origin 如何准确地设置text 的位置
已经有10人回复
- 请问下面这段c++小程序为什么不能运行
已经有6人回复
- 求助：win7 64位环境下，用什么软件运行fortran90?
已经有5人回复
- 求助，怎样去掉双重for循环加快程序运行
已经有10人回复
- 气相色谱分析时候升温程序如何设置
已经有23人回复
- 在linux下运行fortran有什么好处吗？
已经有7人回复
- 梯度荧光定量的程序如何设置
已经有9人回复
- linux下如何让一个程序运行n次？
已经有5人回复
- 有个matlab程序，但是运行不下来，求高手
已经有23人回复
- 岛津LC20AT如何设置样品分析完毕后让泵停止运行？
已经有8人回复
- MS能带计算出现“另一程序正在运行”怎么回事？？--求助啊，弄了一下午都没解决啊。。
已经有16人回复
- 如何在dos运行的程序下打开文件
已经有3人回复
- linux下如何设置脚本文件让matlab自动运行程序？
已经有4人回复
- linux下如何退出matlab界面而不关闭正在运行的程序
已经有7人回复
- ivf在windows下如何实现omp并行计算，vs2008应如何设置
已经有3人回复
- fortran程序运行黑框闪下就消失
已经有5人回复
- 【求助】如何将源代码在linux下运行？
已经有8人回复
- 【求助】弱问下Mathematica里怎么运行*.nb文件
已经有6人回复
|
程序 ODE45_main： clc;clear tspan=[0,180]; y0=[1e-5;0;0;0;0]; [t,y]=ode45('ODE45_fun',tspan,y0); [m,n]=size(y); for i=1:m y(i,6)=-178075801.6*y(i,4)*y(i,5)+490910.1296*y(i,3)*y (i,2)+1128974.05+7481.104004*sin(2.047*t(i)); %sin(2.047*t); % y(i,4)=-3*y(i,2)*y(i,3)-4*y(i,1)*y(i,3)+8*y(i,1)*y(i,3); end data=[t,y]; save ODE45_data.txt data -ascii subplot(2,3,1),plot(t,y(:,1)),title('y(1)') Xlabel('t');Ylabel('y'); subplot(2,3,2),plot(t,y(:,2)),title('y(2)') Xlabel('t');Ylabel('y'); subplot(2,3,3),plot(t,y(:,3)),title('y(3)') Xlabel('t');Ylabel('y'); subplot(2,3,4),plot(t,y(:,4)),title('y(4)') Xlabel('t');Ylabel('y'); subplot(2,3,5),plot(t,y(:,5)),title('y(5)') Xlabel('t');Ylabel('y'); subplot(2,3,6),plot(t,y(:,6)),title('y(6)') Xlabel('t');Ylabel('y'); % plot(t,y(:,1),'bo',t,y(:,2),'rx',t,y(:,3),'gv',t,y(:,4),'r-'); grid on ODE45_fun： function dy=ODE45_fun(t,y) dy(1)=-1.918298553*y(3)*y(4)-121.6697369*y(5)*y(2)+0.006472085*y(2)*y (2)+15.25250926*y(5)*y(5)-0.518363603*sin(2.047*t)+0.001124759; dy(2)=0.007229182*y(5)*y(1)-0.013867729*y(3)-0.005151943*sin (2.047*t)+33.43424564*y(4)*y(5)-0.092169794*y(3)*y(2)-0.698266828; dy(3)=72.10986245*y(2)+0.52129529*y(4)*y(1)+0.025471074*y(3)+921.886526*y(4)/y(1)- 0.47870471*y(4)+0.025471074*y(1)*t-0.38220722*cos(2.047*t)-4.62279911; dy(4)=-57.37263009*y(5)+0.001053501*y(3)*y(1)+20.28825016/y(1)-0.001053501*y (3)+0.064741605*y(4)+20.28825016*t-0.006201915*sin(2.047*t)-3651.904822; dy(5)=0.017284293*y(4)+0.00278644*y(1)*y(2)-0.551218454*y(2)*y(5)*y (5)+0.010839281*y(2)*y(2)*y(5)+0.020353114*cos(2.047*t)+0.110594984; dy=[dy(1);dy(2);dy(3);dy(4);dy(5)]; |
2楼2015-04-13 16:33:09
3楼2015-04-14 08:24:02