第2章 matlab.docx

第2章 matlab

其次章 MATLAB绘图与插值拟合
§ 二维作图
1. 数值图
: t=0::2*pi; plotyy(t, sin(t),t,0.,'LineWidth')) %
[x,y]=meshgrid(-3::3,-3::3);%生成平面网格 z=x.^2-y.^2;
mesh(x,y,z)
[x,y]=meshgrid(-3::3,-2::2); % z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y); axis([-3,3,-2,2,-,])%设定坐标系范围 mesh(x,y,z)


§ 车灯光源投影区域的绘制〔CUMCM 2002 A案例〕










% p=;x=;
for y1=-::
y0=(-::)'*ones(1,73); z0=ones(73,1)*(-::); x0=(y0.^2+z0.^2)/(2*p);
xn=(p^3+4*x0*2*p.*x0+p*(-4*y1*y0+3*2*p*x0))./(2*(p^2+2*p*x0)); yn=(2*p*x0.*y0+p^2*(-y1+y0)+y1*(y0.^2-z0.^2))./(p^2+2*p*x0);
54


zn=(p^2+2*p*x0+2*y1*y0).*z0./(p^2+2*p*x0); y=y0+(yn-y0).*(x-x0)./(xn-x0); z=z0+(zn-z0).*(x-x0)./(xn-x0); plot(y,z,'b.') xlabel('y') ylabel('z') hold on end




§ 动画的绘制
% clear










a=[-8/3 0 0;0 -10 10;0 28 -1]; y=[35 -10 -7]'; h=;
p=plot3(y(1),y(2),y(3),'.','EraseMode','none','MarkerSize',5); %将擦除模式设置为none axis([0 50 -25 25 -25 25]) hold on for i=1:4000 a(1,3)=y(2); a(3,1)=-y(2); ydot=a*y; y=y+h*ydot;
set(p,'XData',y(1),'YData',y(2),'ZData',y(3))%设定图形目标的性质
55


drawnow%填充未完成的图形事务 i=i+1; end % clear
axis([0,2*pi,-1,1]) x=0::2*pi; plot(x,0) hold off hold on for j=1:30
t1=(j-1)*2*pi/30; t2=j*2*pi/30; t=t1::t2; plot(t,0,'.r') plot(t,sin(t),'.') end % clear
axis equal%创立一