MatLab解延时常微分方程组
在Mathematica中，解延时常微分方程组是很直接的，在MatLab中，则有
些编程方面的麻烦，也许是因为本人不熟悉MatLab的语法所致。以下是学习笔
记。
function ddex1%DDEX1 Example 1for DDE23.%This is asimple example of Wille'and Baker that illustrates the%straightforward formulation,computation,and plotting of the solution%of asystem of delay differential equations(DDEs).%%The differential equations%%y'_1(t)=y_1(t-1)%y'_2(t)=y_1(t-1)+y_2(t-
0.2)%y'_3(t)=y_2(t)%%are solved on[0,5]with history y_1(t)=1,y_2(t)=1,y_3(t)=1 for%t=0.%%The lags are specified as avector[1,0.2],the delay differential%equations are coded in the subfunction DDEX1DE,and the history is%evaluated by the function DDEX1HIST.Because the history is constant it%could be supplied as avector:%sol=dde23(@ddex1de,[1,0.2],ones(3,1),[0,5]);%%See also
DDE23,FUNCTION_HANDLE.
%Jacek Kierzenka,Lawrence F.Shampine and Skip Thompson%Copyright 1984-2004 The MathWorks,Inc.%$Revision:1.2.4.2$$Date:2005/06/21 19:
24:16$
sol=dde23(@ddex1de,[2,15,5],@ddex1hist,[0,250]);figure;
plot(sol.x,sol.y)title('Delay P53 Pulse.');xlabel('time t');
ylabel('solution y');
%--
function s=ddex1hist(t)%Constant history function for DDEX1.s=[1;
1;110];
%--
function dydt=ddex1de(t,y,Z)%Differential equations function for DDEX1.ylag1=Z(:,1);ylag2=Z(:,2);ylag3=Z(:,3);
dydt=[1+20*ylag1(2)-
200/(10+ylag1(2))*y(1)30/(10+ylag2(1))*y(3)^4/(250+y(3)^4)-0.5*y(2)*ylag3(1)-0.28*ylag2(1)^4/(250+ylag2(1)^4)*y(3)
本文档为【MatLab解延时常微分方程组】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483