《数字信号处理(第三版)》第10章实验程序(验证版)(西安电子科技大学出版_主编：高西全_丁玉美)综述.doc

《数字信号处理(第三版）》第10章实验程序 注意：每次实验前，请新建以下8个m文件并依实验需要（已注明）移动至matlab当前工作目录下，便于主程序调用。各程序均已验证，请直接复制即可。 文件名：tstem.m（实验一、二需要） 程序： function tstem(xn,yn) %时域序列绘图函数 %xn:被绘图的信号数据序列，yn:绘图信号的纵坐标名称（字符串） n=0:length(xn)-1; stem(n,xn,.); xlabel(n);ylabel(yn); axis([0,n(end),min(xn),1.2*max(xn)]); 文件名：tplot.m（实验一、四需要） 程序： function tplot(xn,T,yn) %时域序列连续曲线绘图函数 %xn:信号数据序列,yn:绘图信号的纵坐标名称（字符串） %T为采样间隔 n=0;length(xn)-1;t=n*T; plot(t,xn); xlabel(t/s);ylabel(yn); axis([0,t(end),min(xn),1.2*max(xn)]); 文件名：myplot.m（实验一、四需要） 程序： %(1)myplot;计算时域离散系统损耗函数并绘制曲线图。 function myplot(B,A) %B为系统函数分子多项式系数向量 %A为系统函数分母多项式系数向量 [H,W]=freqz(B,A,1000) m=abs(H); plot(W/pi,20*log10(m/max(m)));grid on; xlabel(\omega/\pi);ylabel(幅度(dB)) axis([0,1,-80,5]);title(损耗函数曲线); 文件名：mstem.m（实验一、三需要） 程序： function mstem(Xk) %mstem(Xk)绘制频域采样序列向量Xk的幅频特性图 M=length(Xk); k=0:M-1;wk=2*k/M;%产生M点DFT对应的采样点频率（关于pi归一化值） stem(wk,abs(Xk),.);box on;%绘制M点DFT的幅频特性图 xlabel(w/\pi);ylabel(幅度); axis([0,2,0,1.2*max(abs(Xk))]); 文件名：mpplot.m（实验一需要） 程序： %(2)mpplot;计算时域离散系统损耗函数和相频特性函数，并绘制曲线图。 function mpplot(B,A,Rs) %mpplot(B,A,Rs) %时域离散系统损耗函数和相频特性绘图 %B为系统函数分子多项式系数向量 %A为系统函数分母多项式系数向量 %Rs为滤波器阻带最小衰减，省略则幅频曲线最小值取－80dB if nargin3 ymin=-80;else ymin=-Rs-20;end;%确定幅频曲线纵坐标最小值 [H,W]=[H,W]=freqz(B,A,1000) m=abs(H); subplot(2,2,1); plot(W/pi,20*log10(m/max(m)));grid on; xlabel(\omega/\pi);ylabel(幅度(dB)) axis([0,1,ymin,5]);title(损耗函数曲线); subplot(2,2,3); plot(W/pi,p/pi); xlabel(\omega/\pi);ylabel(相位/\pi);grid on; title((b)相频特性曲线); 文件名：mfftplot.m（实验一需要） 程序： function mfftplot(xn,N) %mfftplot(xn,N)计算序列向量xn的N点fft并绘制其幅频特性曲线 Xk=fft(xn,N);%计算信号xn的频谱的N点采样 %===以下为绘图部分==== k=0:N-1;wk=2*k/N; m=abs(Xk);mm=max(m); plot(wk,m/mm);grid on; xlabel(\omega/\pi);ylabel(幅度(dB)); axis([0,2,0,1.2]); title(幅度特性曲线); 文件名：mstg.m（实验四需要） 程序： function st=mstg %产生信号序列向量st,并显示st的时域波形和频谱 %st=mstg 返回三路调幅信号相加形成的混合信号，长度N=1600 N=1600; %N为信号st的长度 Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz，Tp为采样时间 t=0:T:(N-1)*T;k=0:N-1;f=k/Tp; fc1=Fs/10; %第1路调幅信号的载波频率fc1=1000Hz fm1=fc1/10; %第1路调幅信号的调制信号频率