数学模型个实验四.pptVIP

数学模型实验（四） 矩阵特征值和常微分方程求解 一.矩阵的特征值和特征向量 [V,D]=eig（A） A:n?n二维数组，待求特征值和特征向量的矩阵。 V:n?n二维数组，列向量为特征向量 D:对角矩阵，矩阵A的特征值 例1 a=[1,-3,2;4,4,-1;6,3,5];[v,d]=eig(a) 例2 b=[0,1,1;1,0,0;0,1,0];[m,n]=eig(b) 常微分方程求解 1. 二、三阶龙格库塔法求常微分方程的数值解 函数语句 [tout,yout]=ode23(ypfun,[to,tfinal],yo): %调用ode23.m文件解微分方程 参数说明： ypfun:输入参数，被积分的函数名。 to：输入参数，t的初始值 tfinal：输入参数，t的终值。 yo：输入参数，y的初始值。 2. 其他的求解语句 ode45 ode113 ode15s ode 23s ode23t ode23tb 其他的输出语句 odeplot odeprint odephas2 odephas3 例3. 求解 x’=rx(1-x/k)， 0 ≤ t ≤ 20, x(0)=1 建立M文件 logic. M function y=logic(t,x);r=0.3;k=20; y=r.*x.*(1-x./k); 操作 x0=1; [t,x]=ode23(‘logic’,[0,20],x0); plot(t,x) 例4.求解 y’’+(y2-1)y’-y=0, y(0)=0.25, y’(0)=0 化为一阶方程组 y1’=y1(1-y22)-y2, y2’=y1, 建立M文件 rhf . M function fyy=rhf(t,x) fyy=[y(1).*(1-y(2).^2)-y(2);y(1)]; 操作 y0=[0.25,0]’; [t,y]=ode23(‘rhf’,[0,0.25],y0); plot(t,y) * *