| 查看: 933 | 回复: 4 | |||
[交流]
MATLAB7.0出错了？ 已有2人参与
|
|
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.047t); 为什么MATLAB7.0程序运行到上面一行出错了？感觉也没什么错误啊 怎么改，谢谢了  ODE45_main文件: clc;clear tspan=[0,1]; y0=[0; 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.047t); 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'); 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.047t)+0.001124759; dy(2)=0.007229182*y(5)*y(1)-0.013867729*y(3)-0.005151943*sin(2.047t)+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.047t)-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.047t)-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.047t)+0.110594984; dy=[dy(1);dy(2);dy(3);dy(4);dy(5)]; |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有49人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有21人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有8人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有7人回复
-
filecode
已经有10人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MS7.0安装问题
已经有5人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- Matlab 7.0 基础教程 清华大学
已经有29人回复
- matlab一运行就出现的问题
已经有7人回复
- 运行matlab出现错误
已经有4人回复
- matlab程序出错调试
已经有7人回复
- Quantum-ESPRESSO4.3安装问题，求助！
已经有4人回复
- 问下这个问题怎么解决.....真心啊
已经有4人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- xcysden-1.5.53.tar.gz安装问题求助
已经有5人回复
- 求助 matlab交点问题
已经有19人回复
- 一笔画问题matlab实现
已经有5人回复
- matlab出错 求助高手解答
已经有5人回复
- 采用matlab中nlinfit拟合，出现问题了，不知道是哪错了，该怎么解决，求高手指点一下
已经有9人回复
- 【求助】matlab 中，几个m文件调用
已经有14人回复
- MATLAB安装后出现一堆java的问题
已经有12人回复
amoroso_xmc
新虫 (初入文坛)
- 应助: 16 (小学生)
- 金币: 342.5
- 散金: 5
- 帖子: 45
- 在线: 5.8小时
- 虫号: 2482315
- 注册: 2013-05-26
- 专业: 流体力学
2楼2015-04-05 20:41:14
3楼2015-04-06 07:06:14
|
谢谢回答了 2.047*t了，这一行还是出错啊 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); 还有：[t,y]=ode45('ODE45_fun',tspan,y0); 没有错误吧 Warning: Failure at t=5.501826e-001【【【【y0=换不同的初始值，Failure at t不同】】】】. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.776357e-015) at time t. the smallest value allowed怎么设置大一些？可以牺牲一些准确度 ODE45_main文件: clc;clear tspan=[0,180]; y0=[1e-5;1e-5;1e5;1e-5;1e-5]; [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); 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)]; |
4楼2015-04-06 14:43:44
5楼2015-04-06 19:16:41
