24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

value_1987

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 17
帖子: 25
在线: 5.1小时
虫号: 1417825
注册: 2011-09-26
专业: 工程热物理相关交叉领域

[求助] Matlab用龙格库塔法求二阶非线性方程

$X''=bX'+cx+u(t)$
求助解上式方程，用龙格库塔方法，用[t，y]=ode45(‘fname’,tspan,y0)，如果u（t）是常数或正弦变化等，我会编程序，但是u（t）是一个随机的数，与t有关的，我编了个小程序来表示随机数，但是不知道怎么放到M文件里面。

回复此楼

» 猜你喜欢

26年申博自荐-计算机视觉已经有4人回复
考博已经有4人回复
药化及相关博士的申请已经有3人回复
一篇MDPI论文改变了学习工作和生活已经有4人回复
一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇已经有4人回复
太白金星有点烦已经有3人回复
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业已经有3人回复
收到国自然专家邀请后几年才会有本子送过来评已经有3人回复
有没有快的中文核心比较快录用的，纳米材料光催化已经有4人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2013-03-03 11:22:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

value_1987

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 17
帖子: 25
在线: 5.1小时
虫号: 1417825
注册: 2011-09-26
专业: 工程热物理相关交叉领域

引用回帖:

6楼: Originally posted by 夕阳西下 at 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文件，就会产生不同随机数，应该要产生相同随机数才行啊