matlab信号处理经典例程.doc

实验一 离散时间信号时域和频域分析 信号相加并作图 数学描述： MATLAB实现： 设 x1=[1 0.7 0.4 0.1 0]; x2=[0.1 0.3 0.5 0.7 0.9 1] 参考程序： n1=1:5; x1=[1 0.7 0.4 0.1 0]; n2=3:8; x2=[0.1 0.3 0.5 0.7 0.9 1]; n=1:8; x1=[x1 zeros(1,8-length(n1))]; x2=[zeros(1,8-length(n2)) x2]; x=x1+ x2; Subplot(3,1,1); Stem(n,x1); Subplot(3,1,2); Stem(n,x2); Subplot(3,1,3); Stem(n,x); 信号相乘并作图 数学描述： MATLAB实现: x1，x2取值与上题相同。参考程序： n1=1:5; x10=[1 0.7 0.4 0.1 0]; n2=3:8; x20=[0.1 0.3 0.5 0.7 0.9 1]; n=1:8; x1=[x10 zeros(1,8-length(n1))]; x2=[zeros(1,8-length(n2)) x20]; x=x1.*x2; Subplot(3,1,1); Stem(n,x1); Subplot(3,1,2); Stem(n,x2); Subplot(3,1,3); Stem(n,x); 计算卷积 用MATLAB计算序列a=[-2 0 1 –1 3]和b=[1 2 0 -1]的离散卷积。 首先用手工计算，然后用MATLAB编程计算，作图验证手工计算结果。 参考程序： a=[-2 0 1 -1 3]; b=[1 2 0 -1]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c); xlabel(n); ylabel(幅度); 已知系统函数 ，作零极点图，并求z反变换。 参考程序： b=[0,1,2,0,1]; a=[1,-1.5,0.5,0,0]; [z,p,k]=tf2zp(b,a) zplane(b,a) clear b=[0,1,2,0,1]; a=[1,-1.5,0.5,0,0]; [r,p,k]=residue(b,a) r = 8 -13 6 2 p = 1.0000 0.5000 0 0 k = [ ] 5、对信号利用采样间隔0.1s和0.3s采样，在同一个窗口作出该信号及采样信号及其频谱，观察有什么变化。参考程序： t=0:0.01:2*pi; x1=2*sin(t+0.2); n=1:60; nT1=n*0.1; y1=2*sin(nT1+0.2); nT2=n*0.3; y2=2*sin(nT2+0.2); subplot(4,1,1);stem(y1); subplot(4,1,2);stem(fft(abs(y1))); subplot(4,1,3);stem(y2) subplot(4,1,4);stem(fft(abs(y2))) 6、已知一长度为16的有限长序列x(n)=sin(0.25πn)，计算其16点和28点DFT 参考程序： n=0:15; x=sin(0.25*pi*n); dft_16=fft(x); dft_128=fft(x,128); subplot(3,1,1); stem(x); subplot(3,1,2); stem(abs(dft_16)); xlabel(Sample Num); ylabel(Magnitude); title(16点DFT); subplot stem(abs(dft_128)); title(128点DFT); 、将C盘Window文件夹下Media子文件夹中的start.wav文件采样频率重设为8000Hz，PCM编码方式，然后读取其右声道数据作为研究对象，并对其增加不同SNR分别为0dB，30dB的噪声，在同一个窗口画出波形并对比。 参考程序：a=wavread(c:\windows\media\start.wav); left=a(:,1); y1=awgn(left,0,measured); y2=awgn(left,-30,measured); subplot(3,1,1);plot(left); subplot(3,1,2);plot(y1); subplot(3,1,3);plot(y2); 实验二 滤波器设计及应用 1、分别用butterworth和Che