24小时热门版块排行榜    

查看: 642  |  回复: 11

ybkooo

至尊木虫 (著名写手)

^^

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

求出y

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

» 猜你喜欢

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

abcd
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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的回帖
普通回帖

ybkooo

至尊木虫 (著名写手)

^^

有没有大神
abcd
2楼2014-12-07 17:47:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

ybkooo

至尊木虫 (著名写手)

^^

微分方程组,2-3式
abcd
6楼2014-12-08 12:08:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

8楼2014-12-08 16:42:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^

怎么画出F-t的图像
abcd
9楼2014-12-08 16:56:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^

t=0:1/115:1;
>> p1=polyfit(t',Y(:,1),6);
dp=polyder(p1);
dy=polyval(dp,t');
x=sin(2*pi*t');
F=930*dy+5.4*(x-18.9);
>> plot(t,F)
>> plot(x,F)
abcd
10楼2014-12-08 20:01:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ybkooo 的主题更新
信息提示
请填处理意见