| 查看: 2486 | 回复: 13 | ||
[求助]
Matlab用龙格库塔法求二阶非线性方程
|
||
|
求助解上式方程，用龙格库塔方法，用[t，y]=ode45(‘fname’,tspan,y0)，如果u（t）是常数或正弦变化等，我会编程序，但是u（t）是一个随机的数，与t有关的，我编了个小程序来表示随机数，但是不知道怎么放到M文件里面。 |
回复此楼» 猜你喜欢
求个博导看看
已经有19人回复
-
论文终于录用啦！满足毕业条件了
已经有13人回复
-
2025年遐想
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求matlab 用四阶龙格-库塔法求解微分方程
已经有4人回复
- 龙格库塔法 matlab
已经有6人回复
- 隐式欧拉法求解一阶常微分方程
已经有7人回复
- 【求助】Matlab中利用四阶龙格-库塔法求解微分方程！！！！
已经有9人回复
5楼2013-03-03 16:37:39
夕阳西下
金虫 (著名写手)
- 应助: 138 (高中生)
- 贵宾: 0.041
- 金币: 9247.1
- 散金: 73
- 红花: 28
- 沙发: 1
- 帖子: 1268
- 在线: 1164.4小时
- 虫号: 220750
- 注册: 2006-03-20
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 物理
6楼2013-03-04 10:14:56
|
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
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
7楼2013-03-04 13:03:10
9楼2013-03-04 15:45:22
10楼2013-03-04 15:56:50