MATLAB计算结果可视化..doc

第六章MATLAB计算结果可视化 6.1连续函数和离散函数的可视化 【例6-1】用图形表示离散函数。 n=0:12; %产生一组自变量数据 y=1./abs(n-6); %计算相应点的函数值 plot(n,y,r*,MarkerSize,20) %用红花标出数据点 grid on %画坐标方格 Warning: Divide by zero. 图 6-1 离散函数的可视化 【例6-2】用图形表示连续调制波形。 t1=(0:11)/11*pi; y1=sin(t1).*sin(9*t1); t2=(0:100)/100*pi; y2=sin(t2).*sin(9*t2); subplot(2,2,1),plot(t1,y1,r.),axis([0,pi,-1,1]),title(子图 (1)) subplot(2,2,2),plot(t2,y2,r.),axis([0,pi,-1,1]),title(子图 (2)) subplot(2,2,3),plot(t1,y1,t1,y1,r.) axis([0,pi,-1,1]),title(子图 (3)) subplot(2,2,4),plot(t2,y2) axis([0,pi,-1,1]),title(子图 (4)) 图 6-2 连续函数的图形表现方法 6.2二维曲线绘图的基本操作 plot的基本调用格式 【例6-3】用图形表示连续调制波形 及其包络线。 t=(0:pi/100:pi); %长度为101的时间采样列向量 y1=sin(t)*[1,-1]; %包络线函数值，是（101x2）的矩阵 y2=sin(t).*sin(9*t); %长度为101的调制波列向量 t3=pi*(0:9)/9; y3=sin(t3).*sin(9*t3);plot(t,y1,r:,t,y2,b,t3,y3,bo) axis([0,pi,-1,1]) %控制轴的范围 图6-3 6.2.2泛函绘图指令fplot 【例6-4】fplot 与一般绘图指令的绘图效果比较。 [x,y]=fplot(cos(tan(pi*x)),[-0.4,1.4],0.2e-3);n=length(x); subplot(1,2,1),plot(x,y) title(\fontsize{20}\fontname{隶书}泛函绘图指令效果) t=(-0.4:1.8/n:1.4); subplot(1,2,2),plot(t,cos(tan(pi*t))) title(\fontsize{20}\fontname{隶书}等分采样绘图) 图6-4 右图曲线中底部有明显绘图缺陷 注：min(diff(x))=(b-a)*tol; 6.2.3曲线的色彩、线型和数据点形 【例6-5】用图形演示平面上一个方块四个顶点在仿射投影（Affine Projection）下的位置、形状变化。 %平面上的四个点和它们构成的方块 p1=[-0.5,0,1];p2=[-0.5,1,1];p3=[0.5,1,1];p4=[0.5,0,1]; Sq=[p1,p2,p3,p4,p1]; %平移投影：沿x轴移动0.5 , 沿y轴移动1 。 dx=0.5;dy=1;T=[1,0,dx;0,1,dy;0,0,1]; %旋转投影：逆时针旋转30度。 th=pi/6;R=[cos(th),-sin(th),0;sin(th),cos(th),0;0,0,1]; %刻度投影：x方向放大到2倍，y方向放大到3倍。 alpha=2;beta=3;S=[alpha,0,0;0,beta,0;0,0,1]; E=eye(3,3);%为编程方便，设计一个单位阵。在它投影下，任何形状都不变。 TRS={E,T,R,S}; %用元胞数组存放四个变换矩阵 ss={r^,rd,rp,rh}; %用元胞数组存放数据点形设置值 tt={Original Square,Translation,Rotation,Scaling}; %用元胞数组存放四张子图的图名 for i=1:4 W=TRS{i}*Sq; %进行仿射变换 subplot(2,2,i) for k=1:4 plot(W(1,k),W(2,k),ss{k}); %用不同点形标志四个顶点 axis([-3,3,-1,5]),axis equal hold on %使以后图形画在当前子图上 end plo