饿狼追兔的可视化matlab,高阶常微分方程模型饿狼追兔问题数学建模实例

高阶常微分方程模型饿狼追兔问题数学建模实例

现有一只兔子、一匹狼，兔子位于狼的正西110米处，假设兔子与狼同时发现对方并一起起跑，兔子往正北70米处的巢穴跑，而狼在追兔子。已知兔子、狼是匀速跑且狼的速度是兔子的两倍。

要求：(1)建立狼的运动轨迹微分模型。

(2)画出兔子与狼的运动轨迹图形。

(3)用解析方法求解，问兔子能否安全回到巢穴？

(4)用数值方法求解，问兔子能否安全回到巢穴？

【注】常微分方程高阶初值问题的MATLAB库函数为：ode45。

语法为：[t,Y] =ode45(odefun,tspan,y0)

例如函数： function dy = rigid(t,y)

dy = zeros(3,1);    % a column vector

dy(1) = y(2) * y(3);

dy(2) = -y(1) * y(3);

dy(3) = -0.51 * y(1) * y(2);

设置选项：

options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);

求解得：

[t,Y] = ode45(@rigid,[0 12],[0 1 1],options);

画出解函数曲线图形：

plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')

3个回答

二阶偏微分方程的一般形式为

A*Uxx+2*B*Uxy+C*Uyy+D*Ux+E*Uy+F*U=0

其特征方程为

A*(dy)^2-2*B*dx*dy+C*(dx)^2=0

若在某域内B^2-A*C<0则在此域内称为椭圆形方程

若在某域内B^2-A*C=0则在此域内称为抛物形方程

若在某域内B^2-A*C>0则在此域内称为双曲形方程

其实主要是按特征方程的曲线类型分的

注：

Uxx表示U对x求二阶偏导，Uyy表示U对y求二阶偏导，Uxy表示对x求一阶偏导后再对y求一阶偏导，Ux表示U对x求一阶偏导，Uy表示U对y求一阶偏导

partial符号实在打不出来

◆◆

评论读取中....

请登录后再发表评论!

◆◆

修改失败，请稍后尝试

165

◆◆

评论读取中....

请登录后再发表评论!

◆◆

修改失败，请稍后尝试

精华知识