q2522.m

电路与系统专业的matlab程序。有很多经典的实例。

%《MATLAB在电子信息课程中的应用》第二章第五节演示程序二q252
% 三维绘图和屏幕控制
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

clf
disp('三维曲线'),echo on
pause, z=0:0.1:4*pi;x=cos(z);y=sin(z);plot3(x,y,z,'linewidth',2)
pause,text(0,0,0,'*')

disp('三维曲面')
pause, x=-8:0.5:8;y=x';
X=ones(size(y))*x; %33*33阵
Y=y*ones(size(x)); %33*33阵
pause,subplot(1,2,1),plot(X,Y),
pause,subplot(1,2,2),plot(Y,X),
pause,subplot(1,1,1), R=sqrt(X.^2+Y.*Y); z=sin(R)./R; mesh(z)
%直接生成网格点
[X,Y]=meshgrid(-8:0.5:8,-8:0.5:8);
pause,subplot(1,2,1),plot(X,Y),
pause,subplot(1,2,2),plot(Y,X),
pause,disp('去掉分母为0点') 
figure(1)
pause,R=sqrt(X.^2+Y.*Y)+eps; z=sin(R)./R; mesh(z)
pause,R=abs(X)+abs(Y)+eps; z1=sin(R)./R; figure(2),mesh(z1)
pause,title('R=abs(x)+abs(y) sin(R)/R')
%pause,text(10,30,1.1,'R=abs(X)+abs(Y),z=sin(R)./R'),view(20,0),
pause,disp('屏幕分割')
pause,disp('坐标设定')
pause,subplot(2,2,1),R=sqrt(X.^2+Y.*Y); z=sin(R)./R; meshc(z)
pause,title('meshc(z),shading flat'),shading flat
pause,subplot(2,2,2), R=sqrt(X.^2+Y.*Y)+eps; z=sin(R)./R;meshz(z)
pause,title('meshz(z),shading interp'),shading interp
pause,subplot(2,2,3), R=abs(X)+abs(Y)+eps; z1=sin(R)./R; surfc(z1)
pause,title('surfc(z1),shading flat'),shading flat,
%colormap(gray)
pause,subplot(2,2,4); surfl(z1),view(20,0)
pause,title('surfl(z1),view(20,0)')