在matlab内绘制函数图像

小鱼儿 2022-05-16 08:14 300阅读 0赞

1 函数图像1

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

对于函数

[![clip\_image002][clip_image002]][clip_image002_clip_image002]

建立如下，程序：

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

则显示图像：

[![clip\_image004][clip_image004]][clip_image004_clip_image004]

2 函数图像2

函数如下：

[![clip\_image006][clip_image006]][clip_image006_clip_image006]

建立文件，

文件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')

显示图像如下：

[![clip\_image008][clip_image008]][clip_image008_clip_image008]

[clip_image002]: https://p-blog.csdn.net/images/p_blog_csdn_net/superdont/598588/o_clip_image002_thumb_633957813351718750.jpg
[clip_image002_clip_image002]: /images/20220516/008fcd368c47433b99159060586faebf.png
[clip_image004]: https://p-blog.csdn.net/images/p_blog_csdn_net/superdont/598588/o_clip_image004_thumb_633957813373437500.jpg
[clip_image004_clip_image004]: /images/20220516/52119ec789ad4c6f8250c8dfc75b120d.png
[clip_image006]: https://p-blog.csdn.net/images/p_blog_csdn_net/superdont/598588/o_clip_image006_thumb_633957813380468750.jpg
[clip_image006_clip_image006]: /images/20220516/c0c4672a8cf44b31b68ad114ca66914d.png
[clip_image008]: https://p-blog.csdn.net/images/p_blog_csdn_net/superdont/598588/o_clip_image008_thumb_633957813431250000.jpg
[clip_image008_clip_image008]: /images/20220516/2ecbcd4cd52644bfb3700756cc0b924e.png