| 查看: 2858 | 回复: 13 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
Matlab用龙格库塔法求二阶非线性方程
|
|||
|
求助解上式方程，用龙格库塔方法，用[t，y]=ode45(‘fname’,tspan,y0)，如果u（t）是常数或正弦变化等，我会编程序，但是u（t）是一个随机的数，与t有关的，我编了个小程序来表示随机数，但是不知道怎么放到M文件里面。 |
回复此楼» 猜你喜欢
范进中举一文的中心思想
已经有9人回复
-
2026国自然函评费到账
已经有24人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有25人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求matlab 用四阶龙格-库塔法求解微分方程
已经有4人回复
- 龙格库塔法 matlab
已经有6人回复
- 隐式欧拉法求解一阶常微分方程
已经有7人回复
- 【求助】Matlab中利用四阶龙格-库塔法求解微分方程！！！！
已经有9人回复
|
function [y,z]=Runge_kutta(a,b,y0,z0,h) x=a:h:b; y(1)=y0; z(1)=z0; n=(b-a)/h+1; for i=2:n K(1,1)=f1(x(i-1),y(i-1),z(i-1)); K(2,1)=f2(x(i-1),y(i-1),z(i-1)); K(1,2)=f1(x(i-1)+h/2,y(i-1)+K(1,1)*h/2,z(i-1)+K(2,1)*h/2); K(2,2)=f2(x(i-1)+h/2,y(i-1)+K(1,1)*h/2,z(i-1)+K(2,1)*h/2); K(1,3)=f1(x(i-1)+h/2,y(i-1)+K(1,2)*h/2,z(i-1)+K(2,2)*h/2); K(2,3)=f2(x(i-1)+h/2,y(i-1)+K(1,2)*h/2,z(i-1)+K(2,2)*h/2); K(1,4)=f1(x(i-1)+h,y(i-1)+K(1,3)*h,z(i-1)+K(2,3)*h); K(2,4)=f2(x(i-1)+h,y(i-1)+K(1,3)*h,z(i-1)+K(2,3)*h); y(i)=y(i-1)+h/6*(K(1,1)+2*K(1,2)+2*K(1,3)+K(1,4)); z(i)=z(i-1)+h/6*(K(2,1)+2*K(2,2)+2*K(2,3)+K(2,4)); end y(2); plot(y,'r') 自己编写了随机数的M文件，但是每个K都是调用一次M文件，就会产生不同随机数，应该要产生相同随机数才行啊 |
8楼2013-03-04 15:43:55
mechaniclmw
木虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 2548
- 散金: 197
- 红花: 2
- 帖子: 292
- 在线: 205.8小时
- 虫号: 2001977
- 注册: 2012-09-15
- 性别: GG
- 专业: 动力学与控制
2楼2013-03-03 14:18:14
3楼2013-03-03 14:45:31
4楼2013-03-03 14:47:54
20