24小时热门版块排行榜    

查看: 833  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[基金申请] fileCode有新解读? +9 Tide man 2026-08-08 17/850 2026-08-08 14:56 by archvillain
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见