通信系统建模与仿真06.pptVIP

5.5.4 图形的裁剪处理 例5-25 绘制三维曲面图，并进行插值着色处理，裁掉图中x和y都小于0部分。 程序如下： [x,y]=meshgrid(-5:0.1:5); z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2)/4); surf(x,y,z);shading interp; pause %程序暂停 i=find(x=0y=0); z1=z;z1(i)=NaN; surf(x,y,z1);shading interp; 为了展示裁剪效果，第一个曲面绘制完成后暂停，然后显示裁剪后的曲面。 5.6 图像处理与动画制作 5.6.1 图像处理 1．imread和imwrite函数 imread和imwrite函数分别用于将图像文件读入MATLAB工作空间，以及将图像数据和色图数据一起写入一定格式的图像文件。MATLAB支持多种图像文件格式，如.bmp、.jpg、.jpeg、.tif等。 2．image和imagesc函数 这两个函数用于图像显示。为了保证图像的显示效果，一般还应使用colormap函数设置图像色图。 例5-26 有一图像文件flower.jpg，在图形窗口显示该图像。 程序如下： [x,cmap]=imread(flower.jpg); %读取图像的数据阵和色图阵 image(x);colormap(cmap); axis image off %保持宽高比并取消坐标轴 5.6.2 动画制作 MATLAB提供getframe、moviein和movie函数进行动画制作。 1．getframe函数 getframe函数可截取一幅画面信息(称为动画中的一帧)，一幅画面信息形成一个很大的列向量。显然，保存n幅图面就需一个大矩阵。 2．moviein函数 moviein(n)函数用来建立一个足够大的n列矩阵。该矩阵用来保存n幅画面的数据，以备播放。之所以要事先建立一个大矩阵，是为了提高程序运行速度。 3．movie函数 movie(m,n)函数播放由矩阵m所定义的画面n次，缺省时播放一次。 例5-27 绘制了peaks函数曲面并且将它绕z轴旋转。 程序如下 [X,Y,Z]=peaks(30); surf(X,Y,Z) axis([-3,3,-3,3,-10,10]) axis off; shading interp; colormap(hot); m=moviein(20); %建立一个20列大矩阵 for i=1:20 view(-37.5+24*(i-1),30) %改变视点 m(:,i)=getframe; %将图形保存到m矩阵 end movie(m,2); %播放画面2次 例5-28 单摆运动的仿真 %制作动画 %挂摆横梁 plot([-0.2;0.2],[0;0],color,y,linestyle,-,... linewidth,10); %画初始位置的单摆 g=0.0098; %重力加速度,可以调节摆的摆速 l=1; theta0=pi/4; x0=l*sin(theta0); y0=(-1)*l*cos(theta0); axis([-0.75,0.75,-1.25,0]); axis(off); %不显示坐标轴 %创建摆锤 head=line(x0,y0,color,r,linestyle,.,... erasemode,xor,markersize,40); %创建摆杆 body=line([0;x0],[0;y0],color,b,linestyle,-,... erasemode,xor); %摆的运动 t=0; dt=0.01; while 1 t=t+dt; theta=theta0*cos(sqrt(g/l)*t); x=l*sin(theta); y=(-1)*l*cos(theta); set(head,xdata,x,ydata,y); set(body,xdata,[0;x],ydata,[0;y]); drawnow; end 例5.29 模拟布朗运动。 程序如下： n=30; %指定布朗运动的点数 s=0.02; %指定温度或速率 %产生n个随机点(x,y)，处于-0.5到0.5之间 x=rand(n,1)-0.5; y=rand(n,1)-0.5; h=plot(x,y,.); %绘制随机点 axis([-1 1 -1 1]) ; axis sq