24小时热门版块排行榜    

查看: 726  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂院校信息 +4 CX 330 2026-03-21 4/200 2026-03-21 23:48 by ms629
[考研] 328求调剂,英语六级551,有科研经历 +5 生物工程调剂 2026-03-17 9/450 2026-03-21 23:32 by zhujy1982
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 311求调剂 +13 冬十三 2026-03-15 14/700 2026-03-21 22:10 by peike
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 一志愿华中科技大学071000,求调剂 +3 沿岸有贝壳6 2026-03-21 3/150 2026-03-21 10:35 by 暮云清寒
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
信息提示
请填处理意见