3

我正在尝试在 matlab 中解决这个特殊问题

"考虑系统 dr/dt = r*(l - r^2) + mu*r cos (theta)

d(θ)/dt = 1

当 mu = o 时,在 r =1 处有一个稳定的极限环,表明只要 mu 足够小,对于 ,mu > 0 仍然存在闭合轨道。”

我能够绘制一个封闭的极限环,但无法为 x 与 y 图绘制矢量场。

这是我的代码..

function file->
function dy = tst(t,y)
%UNTITLED Summary of this function goes here
% Detailed explanation goes here
dy = zeros(2,1);
dy(1) = y(1)*(1-y(1)*y(1))+1.5*y(1)*cos(y(2));
dy(2) = 1;
end

script file->

[t,r] = ode45(@tst,[0, 40],[1,0]);
x = r(:,1).*cos(r(:,2)); 
y = r(:,1).*sin(r(:,2));
figure(1)
plot(x,y);

hold on

[R,T] = meshgrid(r(:,1),r(:,2));

X = R.*cos(T);

Y = R.*sin(T);

quiver(R,T,X,Y);

我的事情我做错了,绘制矢量场的正确方法是什么?问题出在脚本文件中。

4

1 回答 1

1

这个想法是采用参数的微分:

quiver(x(1:end-1),y(1:end-1),diff(x)./diff(t),diff(y)./diff(t))

在此处输入图像描述

于 2013-02-25T21:56:08.500 回答