567haha
2007-05-20, 14:39
过两天就要交毕业设计了，哪位大虾能帮我完成啊！
我有程序，但不会用Matlab.谁能帮我完成下面一段程序的图和初始条件的设置啊！
clear all;
[T,Y]=ode45('chua',[0,400],[0.1,0.1,0.1]);
%subplot(2,2,1);
%figure;
plot3(Y(:,1),Y(:,2),Y(:,3),'-');
%xlabel('X1');
%ylabel('X2');
%zlabel('x3');
%title('lorenz system with t=0,680,initial 1 2 3,r=28.0,w=28,ampi=sqrt(5.4)' );
%subplot(2,2,2);
%hold on
figure;
plot(T,Y(:,1),'-');
%xlabel('T')
%ylabel('X1');
%subplot(2,2,3);
figure;
plot(T,Y(:,2),'-');
%xlabel('T');
%ylabel('X2')
%subplot(2,2,4);
figure;
plot(T,Y(:,3),'-');
%xlabel('T');
%ylabel('X3');
%hold off
figure;
plot(Y(:,1),Y(:,2),'-');
figure;
plot(Y(:,1),Y(:,3),'-');
figure;
plot(Y(:,2),Y(:,3),'-');
还有一段程序
function dy=chua(t,y)
%clear all
dy=zeros(3,1);
alfa=10.0;%9.2156;%coficient
beta=14.87;%15.994; %cofficient
a=-1.27;%4905;%cofficient
b=-0.65;%75735;
bp=1.0;
dy(1)=alfa*(y(2)-y(1)-(b*y(1)+0.5*(a-b)*(abs(y(1)+bp)-abs(y(1)-bp))));
dy(2)=y(1)-y(2)+y(3);
dy(3)=-beta*y(2);
我有程序，但不会用Matlab.谁能帮我完成下面一段程序的图和初始条件的设置啊！
clear all;
[T,Y]=ode45('chua',[0,400],[0.1,0.1,0.1]);
%subplot(2,2,1);
%figure;
plot3(Y(:,1),Y(:,2),Y(:,3),'-');
%xlabel('X1');
%ylabel('X2');
%zlabel('x3');
%title('lorenz system with t=0,680,initial 1 2 3,r=28.0,w=28,ampi=sqrt(5.4)' );
%subplot(2,2,2);
%hold on
figure;
plot(T,Y(:,1),'-');
%xlabel('T')
%ylabel('X1');
%subplot(2,2,3);
figure;
plot(T,Y(:,2),'-');
%xlabel('T');
%ylabel('X2')
%subplot(2,2,4);
figure;
plot(T,Y(:,3),'-');
%xlabel('T');
%ylabel('X3');
%hold off
figure;
plot(Y(:,1),Y(:,2),'-');
figure;
plot(Y(:,1),Y(:,3),'-');
figure;
plot(Y(:,2),Y(:,3),'-');
还有一段程序
function dy=chua(t,y)
%clear all
dy=zeros(3,1);
alfa=10.0;%9.2156;%coficient
beta=14.87;%15.994; %cofficient
a=-1.27;%4905;%cofficient
b=-0.65;%75735;
bp=1.0;
dy(1)=alfa*(y(2)-y(1)-(b*y(1)+0.5*(a-b)*(abs(y(1)+bp)-abs(y(1)-bp))));
dy(2)=y(1)-y(2)+y(3);
dy(3)=-beta*y(2);