数学实验第三讲.ppt

数学实验之三;微分方程解的形式;一阶微分方程:;MATLAB软件实现;例①;例②常系数的二阶微分方程;x=dsolve((Dx)^2+x^2=1,x(0)=0);输入：

[x,y]=dsolve(Dx=3*x+4*y,Dy=-4*x+3*y)

[x,y]=dsolve(Dx=3*x+4*y,Dy=-4*x+3*y,

x(0)=0,y(0)=1);数值解;1)向前欧拉公式:（y’=f(x,y)）

y(xn+1)?y(xn)+hf(xn,y(xn))(迭代式)

yn+1?yn+hf(xn,yn)(近似式)

特点：f(x,y)取值于区间[xn,xn+1]的左端点.;yn+1?yn+hf(xn+1,yn+1)

特点：①f(x,y)取值于区间[xn,xn+1]的右端点.

②非线性方程,称‘隐式公式’。;例1;2）向前欧拉法：

yn+1=yn+h(-yn+xn+1)

=(1-h)yn+hxn+h

3）向后欧拉法：

yn+1=yn+h(-yn+1+xn+1+1)

转化yn+1=(yn+hxn+1+h)/(1+h);x1(1)=0;y1(1)=1;y2(1)=1;h=0.1;(died.m)

fork=1:10

x1(k+1)=x1(k)+h;

y1(k+1)=(1-h)*y1(k)+h*x1(k)+h;

y2(k+1)=(y2(k)+h*x1(k+1)+h)/(1+h);

end

x1,y1,y2,（y1——向前欧拉解，y2——向后欧拉解）

x=0:0.1:1;

y=x+exp(-x)（解析解）

plot(x,y,x1,y1,k:,x1,y2,r--);x;（2）步长h=0.01的数值解比较表;图形显示;

;梯形公式;以例1为例，用改进欧拉公式编程计算，再与精确解的比较。;x;使用泰勒公式;;1）首先建立M-文件（weif.m）

functionf=weif(x,y)

f=-y+x+1;

2）求解：[x，y]=ode23(‘weif’,[0,1],1)

3)作图形:plot(x,y,‘r’);

4)与精确解进行比较

holdon

ezplot(‘x+exp(-x)’,[0,1]);1、在解n个未知函数的方程组时，x0和x均为n维向量，m-函数文件中的待解方程组应以x的分量形式写成.;注意1：;例如：;例2Vanderpol方程：;（1）编写M文件(文件名为vdpol.m):

functionyp=vdpol(t,y);

yp=[y(2);(1-y(1)^2)*y(2)-y(1)];;蓝色曲线

——y（1）；

（原方程解）

红色曲线

——y（2）；;实验内容;3、Rossler微分方程组：