24小时热门版块排行榜    

查看: 718  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 310求调剂 +3 baibai1314 2026-03-16 3/150 2026-03-21 03:56 by JourneyLucky
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 279分求调剂 一志愿211 +11 chaojifeixia 2026-03-19 12/600 2026-03-21 01:49 by 星空星月
[考研] 354求调剂 +5 Tyoumou 2026-03-18 8/400 2026-03-21 00:35 by JourneyLucky
[考研] 317求调剂 +8 申子申申 2026-03-19 13/650 2026-03-21 00:09 by 刘国森
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
信息提示
请填处理意见