| 查看: 2577 | 回复: 13 | ||
[求助]
Matlab用龙格库塔法求二阶非线性方程
|
||
|
求助解上式方程，用龙格库塔方法，用[t，y]=ode45(‘fname’,tspan,y0)，如果u（t）是常数或正弦变化等，我会编程序，但是u（t）是一个随机的数，与t有关的，我编了个小程序来表示随机数，但是不知道怎么放到M文件里面。 |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求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
