非线性ode及在工程中的应用第一次作业更正版

（课程为天津大学研究生公选课，名为非线性ode及在工程中的应用）

第一次作业

1.画方向场

2.求曲线族

解答过程不知道标准步骤是什么，但是这道题可以直接看出来答案：

   

          

3.分析捕食者和食饵之间的关系（x为捕食者，y为食饵）

                 

    

使用数值方法求解，将时间步分成一个个小段，类似于差分，本道题目的解答以代码和例证的形式给出：

代码：

function E=euler111(f,g,a,b,xa,ya,M)

h=(b-a)/M;

T=zeros(1,M+1);

X=zeros(1,M+1);

Y=zeros(1,M+1);

T=a:h:b;

X(1)=xa;

Y(1)=ya;

for j=1:M

X(j+1)=X(j)+h*feval(f,T(j),X(j),Y(j));

Y(j+1)=Y(j)+h*feval(g,T(j),X(j),Y(j));

end

E=[T' X' Y'];

例证：

种群数3000食饵，120捕食者

参数：

f=@(t,x,y)2*x-0.02*x*y

g=@(t,x,y)0.0002*x*y-0.8*y

捕食者波动图像：

食饵波动图像：

4.地球半径为6437km g=9.8m/s^2，求逃逸速度v0

推导：动能全部转化为势能（对引力从r到无穷远积分即可）：

得到逃逸速度为11.2km/s

5. 求齐次方程的通解，系数矩阵A分别是：

(1) [-5, -10, -20; 5,5,10; 2,4,9]

特征值：[ 5, 2 + 1i, 2 - 1i]

复数空间的通解是

通过线性变化将解返回实数空间，可以得到微分方程组的通解为：

(2) [1, 2/3,-2/3; 0, 2/3, 1/3; 0, -1/3, 4/3]

解其特征值：[ 1, 1, 1]

通解：

（2）更正

两个方程均可以通过matlab程序求解，将本次作业的代码均上传至csdn审核，审核通过之后将会附上链接，如果作业或者代码有误，欢迎随时交流。