| 查看: 855 | 回复: 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 |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有6人回复
-
求助文献
已经有3人回复
-
聘U V热熔胶研究人员
已经有9人回复
-
垃圾破二本职称评审标准
已经有10人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
-
博士申请都是内定的吗？
已经有14人回复
-
谈谈两天一夜的“延安行”
已经有13人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
-
论文投稿求助
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
关于用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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人