shili16.m

自己手敲并验证的matlab 6.0实例的代码 及运行程序所得的图像结果。此部分为前20个例子。为matlab的基础部分

function shili16
h0=figure('toolbar','none',...
    'position',[200 200 500 500],...
    'name','shili16');
A=[-8/3 0 0;0 -10 10; 0 28 -1];
%初始化矩阵A
y=[35 -10 -7]';
%初始化坐标
h=0.01;
P=plot3(y(1),y(2),y(3),'.',...
    'erasemode','none',...
    'markersize',5);
%在三维坐标中绘制点；
%设置erasemode属性值none
%设置坐标轴的显示范围
axis([0 50 -25 25 -25 25])
hold on 
%保留当前图形
for i=1:2000
    %建立循环
    P=findobj('xdata',y(1),...
        'ydata',y(2),...
        'zdata',y(3));
    %查找符合条件的图形的句柄
    A(1,3)=y(2);
    A(3,1)=-y(2);
    %修改矩阵A
    ydot=A*y;
    %劳伦兹方程表达式
    y=y+h*ydot;
    %建立新坐标
    set(P,'xdata',y(1),...
        'ydata',y(2),...
        'zdata',y(3))
    %修改句柄为P的图形的坐标属性
    drawnow
    %绘图
    i=i+1;
    %结束循环的条件
end