实验MATLAB在数字信号处理中应用.docxVIP

实验5：MATLAB在数字信号处理中的应用(4学时）　掌握MATLAB对离散时间信号与系统的时域和频域描述及分析【实验目的】掌握离散时间信号与系统的时域和频域分析方法【实验内容】练习书上实例，验证结果录制一段各人自己的语音信号,存成.wav文件,录制时间不要太长，5-10秒即可。利用MATLAB提供的wavread函数进行读取，而后对录制的信号x(t)进行采样；画出采样后语音信号的时域波形和频谱图。（提示，仿例题7.10,选择不同的采样频率fs,观察频谱）仿照教材7.1节对该语音信号进行各种时域处理（如截取，合成，添加随机噪声，周期延拓，移位），观察波形图，若可能可将其用MATLAB存成wav文件，听一下。（提示：函数randn可构造随机噪声序列；人的语音信号的范围20Hz-20KHz,选择合适的采样频率，在数字通信中，根据采样定理, 最小采样频率为语音信号最高频率的2倍。电话通令的频带是300-3400Hz，所以最低采样频率为6800Hz。根据CCITT的长途通信协议的标准，目前各国一般采用的声音采样频率为8kHz；所以，为了将来系统与其他通信设备交流的方便，最终采样频率一般确定为8000Hz)。用FFT函数对上述语音信号，进行频域分析，观察该语音信号的频谱范围.（仿例题7.19）.构造一个理想数字滤波器，要求保留1000HZ-2000HZ以内的语音信号，观察信号波形图。【实验报告】Q701x1=[0,1,2,3,4,3,2,1,0];ns1=-2; x2=[2,2,0,0,0,-2,-2]; ns2=2; nf1=ns1+length(x1)-1; nf2=ns2+length(x2)-1;ny= min(ns1,ns2):max(nf1,nf2); a1 = zeros(1,length(ny)); xa2 = xa1; xa1(find((ny=ns1)(ny=nf1)==1))=x1; xa2(find((ny=ns2)(ny=nf2)==1))=x2; ya = xa1 + xa2 yp = xa1.* xa2 subplot(4,1,1), stem(ny,xa1,.) subplot(4,1,2), stem(ny,xa2,.)line([ny(1),ny(end)],[0,0]) subplot(4,1,3), stem(ny,ya,.)line([ny(1),ny(end)],[0,0]) subplot(4,1,4), stem(ny,yp,.)line([ny(1),ny(end)],[0,0]) clear;close alln0=input(输入序列起点:n0=);N=input(输入序列长度:N=);n1=input(输入位移:n1=);n=n0:n1+N+5; u=[(n-n1)=0]; x1=[(n-n1)=0]-[(n-n1-N)=0]x2=[(n=n1)(n(N+n1))];x3=exp(j*n*pi/8).*x2;subplot(2,2,1);stem(n,x1,.);xlabel(n);ylabel(x1(n));axis([n0,max(n),0,1]);subplot(2,2,3);stem(n,x2,.);xlabel(n);ylabel(x2(n));axis([n0,max(n),0,1]);subplot(2,2,2);stem(n,real(x3),.);xlabel(n);ylabel(x3(n)的实部);line([n0,max(n)],[0,0]);axis([n0,max(n),-1,1]);subplot(2,2,4);stem(n,imag(x3),.);xlabel(n);ylabel(x3(n)的虚部);line([n0,max(n)],[0,0]);axis([n0,max(n),-1,1]);输入序列起点:n0=-6输入序列长度:N=15输入位移:n1=3x1 = Columns 1 through 13 0 0 0 0 0 0 0 0 0 1 1 1 1 Columns 14 through 26 1 1 1 1 1 1 1 1 1 1 1 0 0 Columns 27 through 30 0 0 0 0Q03c