24小时热门版块排行榜    

查看: 645  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ybkooo

至尊木虫 (著名写手)

^^

[求助] 求matlab程序,求解微分方程组 已有1人参与

求出y

求matlab程序,求解微分方程组
现象模型.png
求matlab程序,求解微分方程组-1
123.png
回复此楼

» 猜你喜欢

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

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

zlt11

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
ybkooo: 金币+50, ★★★★★最佳答案 2014-12-08 16:42:23
function df = odefunction(t,f,fp)
%隐式微分方程组
%   记 y(t)=f(1);z(t)=f(2);y'=fp(1);z'=fp(2)
x=@(t) 0.01*sin(2*pi*t);
dx=@(t) 0.01*2*pi*cos(2*pi*t);
alpha=963;c0=53;k0=14;c1=930;k1=5.4;gamma=200;beta=200;n=2;A=207;x0=18.9
df=[ fp(1)-1/(c0+c1)*(alpha*f(2)+c0*dx(t)+k0*(x(t)-f(1)))   ;
    fp(2)-(-gamma*abs(dx(t)-fp(1))*f(2)*abs(f(2)).^(n-1)-beta*(dx(t)-...
    fp(1))*abs(f(2)).^n+    A*(dx(t)-fp(1)));
]
end

[T,Y]=ode15i(@odefunction,[0,1],[0;0],[0;0])
plot(T,Y)%where Y(1) denotes y(t) and Y(2) z(t)
%here it is supposed that y(0)=z(0)=0;y'(0)=z'(0)=0 at the initial points which you might need to alter accordingly
abracadabra
7楼2014-12-08 14:58:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

mofakaka

新虫 (初入文坛)

不好意思,我也不会解你这个方程。但是以前解过其他的微分方程。
数值求解微分方程,第一步你也要先把方程离散化,而且你的方程式隐式的,你也要先确定y、z的差分格式,然后再求解。我知道的有什么龙格-库塔法,New Mark法之类的微分方程解法,你可以看看
3楼2014-12-07 22:27:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^

最近在求解这个方程,只是知道把一式带入2,3式可以求出y,z;在带去4式q求出F
用matlab的ode45可以吧
abcd
4楼2014-12-07 22:47:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^

求大神,之要有正解,可以多给金币的;在给50
abcd
5楼2014-12-07 22:50:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见