在matlab内绘制函数图像

1 函数图像1

详细情况可参考matlab的帮助文档

对于函数

建立如下，程序：

程序1：

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);

程序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),'.')

则显示图像：

2 函数图像2

函数如下：

建立文件，

文件1：

function dy = vdp1000(t,y)

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

dy(1) = y(2);

dy(2) = 1000*(1 - y(1)^2)*y(2) - y(1);

文件2：

[t,Y] = ode15s(@vdp1000,[0 3000],[2 0]);

% plot(T,Y(:,1),'-o')

% plot(t,Y);

plot(t,Y(:,1),'-o')

显示图像如下：