matlab在科学计算中的应用6+-+微分方程问题的解法 .pptVIP

例： 编写函数???： function dx=c7exdde(t,x,z) xlag1=z(:,1); %第一列表示提取 xlag2=z(:,2); dx=[1-3*x(1)-xlag1(2)-0.2*xlag2(1)^3-xlag2(1); x(3); 4*x(1)-2*x(2)-3*x(3)]; 历史数据函数： function S=c7exhist(t) S=zeros(3,1); 求解： lags=[1 0.5]; tx=dde23(c7exdde,lags,zeros(3,1),[0,10]); plot(tx.x,tx.y(2,:)) ％与ode45()等返回的x矩阵不一样，它是按行排列的。 6.4边值问题的计算机求解 6.4.1 边值问题的打靶算法 数学方法描述： 以二阶方程为例 编写函数： 线性的 function [t,y]=shooting(f1,f2,tspan,x0f,varargin) t0=tspan(1); tfinal=tspan(2); ga=x0f(1); gb=x0f(2); [t,y1]=ode45(f1,tspan,[1;0],varargin); [t,y2]=ode45(f1,tspan,[0;1],varargi