| 查看: 874 | 回复: 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 |
回复此楼» 猜你喜欢
085700资源与环境308求调剂
已经有7人回复
-
265求调剂
已经有5人回复
-
279分求调剂 一志愿211
已经有3人回复
-
求调剂
已经有3人回复
-
一志愿吉林大学材料学硕321求调剂
已经有12人回复
-
286分人工智能专业请求调剂愿意跨考！
已经有3人回复
-
329求调剂
已经有5人回复
-
申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复
已经有4人回复
-
材料学硕318求调剂
已经有5人回复
-
一志愿中国海洋大学，生物学，301分，求调剂
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
