| 查看: 917 | 回复: 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 |
回复此楼» 猜你喜欢
不应该看fileCode
已经有6人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有7人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
关于Filecode分析方法
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
找到一些相关的精华帖子，希望有用哦~
关于用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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
