| 查看: 951 | 回复: 0 | ||
lijiawei2013金虫 (小有名气)
|
[求助]
龙格库塔计算结果问题
|
|
我在matlab用龙格库塔oed45计算下面微分方程组时，为什么随着求解区间的增大，它的计算结果为NaN？ function dx=fun1(t,x) k12=4.85*10^9;afa=pi/9;r1=0.48745;r2=0.1078,T1=2000;k1x=2*10^9;k1y=2*10^9 k2x=2.5*10^9;k2y=2.5*10^9;I1=152.16;I2=0.69;m1=1096.04;m2=101.7 dx=zeros(12,1) dx(1)=x(2) dx(2)=(k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11))*cos(afa)-k1x*x(1))/m1 dx(3)=x(4) dx(4)=(-k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11))*sin(afa)-k1y*x(3))/m1 dx(5)=x(6) dx(6)=(T1+r1*k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11)))/I1 dx(7)=x(8) dx(8)=(-k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11))*cos(afa)-k2x*x(7))/m2 dx(9)=x(10) dx(10)=(k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11))*sin(afa)-k2y*x(10))/m2 dx(11)=x(12) dx(12)=(-r2*k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11)))/I2 end |
回复此楼» 猜你喜欢
【长期招聘】复旦大学 AI 制药实验室 博士后招聘
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
关于用1stpro拟合数值的问题。
已经有4人回复
- 动力学方程组参数求解问题
已经有9人回复
- 求助！！对于调用龙格库塔法如何对结果判断反馈
已经有6人回复
- Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！
已经有4人回复
- 编程有点问题，求助
已经有5人回复
- 请问什么是分段龙格库塔的方法？需要用matlab软件采用这种方法求二阶微分方程的解
已经有3人回复
- 求助：matlab解17个未知数，17个方程。
已经有4人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 龙格-库塔（runger-kuta）求解微分方程
已经有3人回复
- Matlab用龙格库塔法求二阶非线性方程
已经有13人回复
- 龙格库塔法 matlab
已经有6人回复
- 4阶龙格库塔方法疑问
已经有6人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- matlab-常微分方程参数估计
已经有12人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 【求助】Matlab中利用四阶龙格-库塔法求解微分方程！！！！
已经有9人回复
- 【求助】请教二阶微分方程的特解（含初始条件）
已经有8人回复
- 【求助】四阶龙格-库塔方法！！！！
已经有17人回复
- 【求助】这个怪异公式是如何推导出来的？
已经有9人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
