24小时热门版块排行榜

>论坛更新日志 (3556)
>虫友互识 (328)
>文献求助 (314)
>导师招生 (188)
>休闲灌水 (130)
>基金申请 (107)
>硕博家园 (106)
>考博 (88)
>博后之家 (46)
>招聘信息布告栏 (39)
>考研 (36)
>论文投稿 (34)
>论文道贺祈福 (29)
>海外博后 (25)
>教师之家 (21)
>找工作 (18)

返回列表

【悬赏金币】回答本帖问题，作者rrf2008将赠送您 100 个金币

rrf2008

铁杆木虫 (正式写手)

应助: 0 (幼儿园)
金币: 7160.8
帖子: 378
在线: 48.2小时
虫号: 418845
注册: 2007-07-03
专业: 数学

[求助] 这个程序如何运行已有1人参与

1.微分方程函数文件
function dy=cdsf1(t,y)
dy=zeros(4,1);
dy(1)=-y(1)*0.9+(y(2)*y(1)+y(4)+y(3));
dy(2)=1-y(2)*0.2-y(1)^2;
dy(3)=-y(3)*1.5-y(1);
dy(4)=-y(4)*0.17-y(1)*y(2)*0.2;
2.和对比效果图
m=100;
t1=0;t2=0;
y1=0.5;y2=0.5;y3=0.5;y4=0.5;
y11=0.6;y22=0.6;y33=0.6;y44=0.6;
%x的图形对比
for n=1:m
t1=(n-1)*0.3;
t2=n*0.3;
[T1,Y1]=ode45('cdsf1',[t1,t2],[y1 y2 y3 y4]);
[T2,Y2]=ode45('cdsf1',[t1,t2],[y11 y22 y33 y44]);
[H1,L1]=size(Y1);[H2,L2]=size(Y2);
y1=Y1(H1,1);y2=Y1(H1,2);y3=Y1(H1,3);y4=Y1(H1,4);
y11=Y2(H2,1)-0.0111*(Y2(H2,1)-Y1(H1,1));
y22=Y2(H2,2)-0.0111*(Y2(H2,2)-Y1(H1,2));
y33=Y2(H2,3)-0.0111*(Y2(H2,3)-Y1(H1,3));
y44=Y2(H2,4)-0.0111*(Y2(H2,4)-Y1(H1,4));
plot(T1,Y1(:,1))
hold on
plot(T2,Y2(:,1),'r:')
hold on
end
请问上述程序如何在Matlab软件中运行出这个图形（如图）

回复此楼

» 猜你喜欢

AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有8人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有12人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有11人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有10人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有9人回复
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了已经有7人回复
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有9人回复
有院领导为了换新车，用横向课题经费买了俩车已经有10人回复

奥运火炬手2008

1楼 2022-07-01 19:25:14

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

zyh700600

金虫 (小有名气)

应助: 47 (小学生)
金币: 9368.2
红花: 10
帖子: 266
在线: 137.1小时
虫号: 56711
注册: 2005-01-16
专业: 金融工程

【答案】应助回帖

感谢参与，应助指数 +1

1. 将下面的代码保存为M文件：

CODE:

 
function rrf2008

close all;clear,clc

%  2.和对比效果图

m=100;

t1=0;t2=0;

y1=0.5;y2=0.5;y3=0.5;y4=0.5;

y11=0.6;y22=0.6;y33=0.6;y44=0.6;

% x的图形对比

for n=1:m

    t1=(n-1)*0.3;

    t2=n*0.3;

    [T1,Y1]=ode45(@cdsf1,[t1,t2],[y1 y2 y3 y4]);

    [T2,Y2]=ode45(@cdsf1,[t1,t2],[y11 y22 y33 y44]);

    [H1,L1]=size(Y1);[H2,L2]=size(Y2);

    y1=Y1(H1,1);y2=Y1(H1,2);y3=Y1(H1,3);y4=Y1(H1,4);

    y11=Y2(H2,1)-0.0111*(Y2(H2,1)-Y1(H1,1));

    y22=Y2(H2,2)-0.0111*(Y2(H2,2)-Y1(H1,2));

    y33=Y2(H2,3)-0.0111*(Y2(H2,3)-Y1(H1,3));

    y44=Y2(H2,4)-0.0111*(Y2(H2,4)-Y1(H1,4));

    plot(T1,Y1(:,1))

    hold on

    plot(T2,Y2(:,1),'r:')

    hold on

end

end

% 1.微分方程函数文件

function dy=cdsf1(t,y)

dy=zeros(4,1);

dy(1)=-y(1)*0.9+(y(2)*y(1)+y(4)+y(3));

dy(2)=1-y(2)*0.2-y(1)^2;

dy(3)=-y(3)*1.5-y(1);

dy(4)=-y(4)*0.17-y(1)*y(2)*0.2;

end