资源描述:
《数据和函数可视化.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、数据和函数的可视化【例】用图形表示离散函数。n=0:12;%产生一组自变量数据y=1./abs(n-6);%计算相应点的函数值plot(n,y,'r*','MarkerSize',20)%用红花标出数据点gridon%画坐标方格Warning:Dividebyzero.图7.1.1-1离散函数的可视化【例】用图形表示连续调制波形。t1=(0:11)/11*pi;%<1>y1=sin(t1).*sin(9*t1);t2=(0:100)/100*pi;%<3>y2=sin(t2).*sin(9*t2);subplot(2,2,1),plot(t1,y1,'r.'),a
2、xis([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)')图7.1.2-1连续函数的图形表现方法【例】简单例题,比较方便的试验指令。t=(0:pi/50:2*pi)';k=0.4:0.
3、1:1;Y=cos(t)*k;plot(t,Y)图7.2.1-1plot指令基本操作演示【例】用图形表示连续调制波形及其包络线。t=(0:pi/100:pi)';%长度为101的时间采样列向量<1>y1=sin(t)*[1,-1];%包络线函数值,是(101x2)的矩阵<2>y2=sin(t).*sin(9*t);%长度为101的调制波列向量<3>t3=pi*(0:9)/9;%<4>y3=sin(t3).*sin(9*t3);plot(t,y1,'r:',t,y2,'b',t3,y3,'bo')%<5>axis([0,pi,-1,1])%控制轴的范围<6>图7.2
4、.1-2【例】用复数矩阵形式画利萨如(Lissajous)图形。(在模拟信号时代,Lissajous图形常用来测量信号的频率。)t=linspace(0,2*pi,80)';%<1>X=[cos(t),cos(2*t),cos(3*t)]+i*sin(t)*[1,1,1];%(80x3)的复数矩阵plot(X)%<3>axissquare%使坐标轴长度相同<4>legend('1','2','3')%图例图7.2.1-3Lissajous图【例】采用模型画一组椭圆。th=[0:pi/50:2*pi]';%长度为101的列向量a=[0.5:.5:4.5];%长度为9
5、的行向量X=cos(th)*a;%(101x9)的矩阵Y=sin(th)*sqrt(25-a.^2);%(101x9)的矩阵plot(X,Y),axis('equal'),xlabel('x'),ylabel('y')title('AsetofEllipses')图7.2.1-4一组椭圆【例】观察各种轴控制指令的影响。演示采用长轴为3.25,短轴为1.15的椭圆。注意:采用多子图表现时,图形形状不仅受“控制指令”影响,而且受整个图面“宽高比”及“子图数目”的影响。t=0:2*pi/99:2*pi;x=1.15*cos(t);y=3.25*sin(t);%y为长轴,
6、x为短轴subplot(2,3,1),plot(x,y),axisnormal,gridon,title('NormalandGridon')subplot(2,3,2),plot(x,y),axisequal,gridon,title('Equal')subplot(2,3,3),plot(x,y),axissquare,gridon,title('Square')subplot(2,3,4),plot(x,y),axisimage,boxoff,title('ImageandBoxoff')subplot(2,3,5),plot(x,y),axisimagef
7、ill,boxofftitle('ImageandFill')subplot(2,3,6),plot(x,y),axistight,boxoff,title('Tight')图7.2.3.1-1各种轴控制指令的不同影响【例】通过绘制二阶系统阶跃响应,clf;t=6*pi*(0:100)/100;y=1-exp(-0.3*t).*cos(0.7*t);tt=t(find(abs(y-1)>0.05));ts=max(tt);%<2>subplot(1,2,1),plot(t,y,'r-','LineWidth',3),gridon%<3>axis([0,6*pi,0
8、.6,ma