在matlab内绘制函数图像
http://blog.youkuaiyun.com/superdont/archive/2009/12/07/4955303.aspx
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')
显示图像如下:
本文介绍了如何使用Matlab绘制两个不同的函数图像。首先通过定义函数和设置参数来创建图像,然后利用ode45和ode15s求解器生成数值解,并最终绘制出图像。这些例子展示了Matlab在解决常微分方程并可视化结果方面的强大功能。
9123

被折叠的 条评论
为什么被折叠?



