信息处理与仿真实验(4)

5.海尔,联想，方正近四年销售的数据如下。

data=[3000,2334,3444,6663;2342,3422,2111,4453;3232,2334,5433,3500];

请做出如下图表。（参考bar，legend）

图3-5 销售比较图

a=[3000 2334 3444 6663;2342 3422 2111 4453;3232 2334 5433 3500]; bar(a);

label={'海尔','联想','方正'}; set(gca,'xticklabel',label);

title('各品牌销售比较图');

legend('2006','2007','2008','2009');

16

6.已知在三维空间里，x,y,z的变化规律如下，试用plot3来做出其图形。 z = 0:0.1:40; x = cos(z).^2; y = sin(z)+cos(z);

z = 0:0.1:40; x = cos(z).^2;

y = sin(z)+cos(z); plot3(x,y,z) grid on

axis square

7. 根据下面的方程绘制图形，y=sin(x)+x，然后再分别用多项式拟合函数和用拟合工

具进行拟合。(参考polyfit, polyval)

x=-pi:0.1:pi;

17

y=sin(x)+x;

p = polyfit(x,y,6); f = polyval(p,x);

plot(x,y,'o',x,f,'-'); axis([-1 4 0 5]);

8. 在0-450范围内，x,y均每隔50产生一个数，这样产生一个10*10网格区域，用randn函数

来产生一个10*10的随机矩阵（标准差是10，平均数是20），使用surf来作出这个三维的曲面图。然后使用interp2进行二次插值，插值函数使用三次杨条插值函数，在原来的第一个分隔区域里产生5个小区域，并最后画出图形。(参考meshgrid)

x=0:50:450; y=0:50:450;

[px,py]=meshgrid(x,y); z=10*randn(10,10)+20; subplot(2,1,1); surf(px,py,z); xi=0:10:450; yi=0:10:450;

[pxi,pyi]=meshgrid(xi,yi);

zi=interp2(x,y,z,xi,yi','spline'); subplot(2,1,2); surf(pxi,pyi,zi);

18

实验四 MATLAB函数编写和程序设计

一、实验目的

1．掌握函数的编写规则。 2．掌握函数的调用。

3．会用Matlab程序设计实现一些工程算法问题。 二、实验设备及条件

计算机一台（带有MATLAB6.0以上的软件环境）。 三、实验内容

1.生成谐波信号x=sin（100*pi*t）+sin（240*pi*t），绘制谐波信号的时域波形和傅里叶幅度谱。

t=0:0.00001:0.4;

x=sin(100*pi*t)+sin(240*pi*t); plot(t,x);

axis([0 0.4 -2 2]);

19

t=0:0.0001:0.4;

x=sin(100*pi*t)+sin(240*pi*t); y=fft(x); z=abs(y) plot(t,z);

axis([-0.1 0.5 0 10]);

2.生成方差为2的高斯白噪声，加入1中谐波信号，绘制含噪信号的时域波形和频谱

t=0:0.001:0.4;

x=sin(100*pi*t)+sin(240*pi*t); y=x+2*randn(size(x)); plot(t,y);

axis([0 0.4 -10 10]);

t=0:0.001:0.4;

x=sin(100*pi*t)+sin(240*pi*t); y=x+2*randn(size(x)); n=fft(y);

20