画相位图【matlab源码】.doc

该【画相位图【matlab源码】 】是由【windurst】上传分享，文档一共【3】页，该文档可以免费在线阅读，需要了解更多关于【画相位图【matlab源码】 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。画相位图【matlab源码】
毕业论文,设计,
题目学院学院专业学生姓名学号年级级指导教师
毕业教务处制表毕业
??ÏàÎ?Í?
一、程序说明
本团队长期从事matlab编程与仿真工作,擅长各类毕业设计、数据处理、图表绘制、理论分析等,程序
代做、数据分析
具体信息联系
二、程序示例
functionvarargout=saxplaxltest31(varargin)
clear,clc
a=0;b=30;y0=[1029];hmax=;hmin=;e=10^(-4);
[x,y]=RKF1(a,b,y0,e,hmax,hmin);%n=length(x);
%n1=length(x);
%{
fprintf('xny1ny2ny3n\n');fori=1:n1
fprintf('%f%f%f\n',x(i),y(i,1),y(i,2),y(i,3));
end
%}
plot3(y(:,1),y(:,2),y(:,3))
view([2040])
functionu=f(x,y)
n=length(y);
u=zeros(1,n);
u(1)=10*(-y(1)+y(2));
u(2)=28*y(1)-y(2)-y(1)*y(3);
u(3)=y(1)*y(2)-8*y(3)/3;
function[x,y]=RKF1(a,b,y0,e,hmax,hmin)t=a;
x=a;
h=hmax;
i=1;
y=y0;
while(t<b)
k1=h*f(t,y0);
k2=h*f(t+h/4,y0+1/4*k1);
k3=h*f(t+h*3/8,y0+3/32*k1+9/32*k2);
k4=h*f(t+12/13*h,y0+1932/2197*k1-7200/2197*k2+7296/2197*k3);
k5=h*f(t+h,y0+439/216*k1-8*k2+3680/513*k3-854/4104*k4);
k6=h*f(t+h/2,y0-8/27*k1+2*k2-3544/2565*k3+1859/4104*k4-11/40*k5);
R=abs(1/360*k1-128/4275*k3-2197/75240*k4+1/50*k5+2/55*k6)/h;
RL=max(R);
delta=*(e/RL)^(1/4);
ifR<=e
t=t+h;
x=[x;t];
y0=y0+25/216*k1+1408/2565*k3+2197/4104*k4-1/5*k5;
y=[y;y0];
end
ifdelta<=
h=*h;
elseifdelta>=4
h=4*h;
else
h=delta*h;
end
ifh>=hmax
h=hmax;
end
ifh<hmin
error('Minimumhexceeded!!!');
end
end
%disp('Theprocedureiscomplete!');