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’)
显示图像如下：
版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请联系我们举报，一经查实，本站将立刻删除。
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/231016.html原文链接：https://javaforall.net
