24小时热门版块排行榜    

查看: 736  |  回复: 2

liuhaizhen

木虫 (初入文坛)

[求助] 求MATLAB编程

dy(1)/dt=(a+b+c)*y(1)+g*y(2)+b*y(3)
dy(2)/dt=-(a+b)*y(3)
dy(3)/dt=(a-b-c)*y(3)-g*y(1)+2*g*y(2)*y(1)
其中a,b,c,g为常量
回复此楼

» 猜你喜欢

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

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

mathcfd

金虫 (著名写手)

【答案】应助回帖

★ ★
余泽成(金币+2): 谢谢参与应助! 2011-12-02 17:32:02
liuhaizhen(金币+20): 2011-12-07 10:07:34
程序如下:

global a b c g;
clc
% 常数赋值
a=1;
b=1;
c=1;
g=1;

tspan=[0,5]; % 积分区间
y0=[0;1;0];  % 初值向量
[t,y]=ode45('vdp1',tspan,y0);

figure(1)
plot(t,y(:,1),'-o');
hold on
figure(2)
plot(t,y(:,2),'-r');
figure(3)
plot(t,y(:,3),'r*');
Neverevergiveup!
2楼2011-12-01 14:44:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathcfd

金虫 (著名写手)

【答案】应助回帖


余泽成(金币+1): 鼓励应助! 2011-12-02 17:32:13
liuhaizhen(金币+10): 2011-12-07 10:08:36
另函数定义


function F=vdp1(t,y)
global a b c g;
F=[(a+b+c)*y(1)+g*y(2)+b*y(3);-(a+b)*y(3);(a-b-c)*y(3)-g*y(1)+2*g*y(2)*y(1)];
Neverevergiveup!
3楼2011-12-01 14:45:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liuhaizhen 的主题更新
信息提示
请填处理意见