第10章实验程序(验证版)

《第10章实验程序(验证版)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、《数字信号处理(第三版)》第10章实验程序注意：每次实验前，请新建以下8个m文件并依实验需要(已注明)移动至matlab当前工作目录下，便于主程序调用。各程序均已验证，请直接复制即可。文件名:tstem.m(实验一、二需要)程序：functiontstem(xn,yn)%时域序列绘图函数%xn:被绘图的信号数据序列，yn:绘图信号的纵坐标名称(字符串)n=0:length(xn)-l;stemjn^n/.');xlabelCn'hylabelCyn1);axis([0/n(end)/min(xn)/1.2*max(xn)]);文

2、件名：tplot.m(实验一、四需要)程序：functiontplot(xnXyn)%时域序列连续曲线绘图函数%xn:信号数据序列,yn:绘图信号的纵坐标名称(字符串)%T为采样间隔n=O;length(xn)-l;t=n*T;plot(t,xn);xlabel('t/s');ylabel(yn);axis([0/t(end)/min(xn)/1.2*max(xn)]);文件名：myplot.m(实验一、四需要)程序：%(l)myplot;计算时域离散系统损耗函数并绘制曲线图。functionmyplot(B,A)%B为系统函数

3、分子多项式系数向量%A为系统函数分母多项式系数向量[H/W]=freqz(B,A,1000)m=abs(H);plot(W/pi/20*logl0(m/max(m)));gridon;xlabel('omega/pi');ylabel('幅度(dB)Jaxis([0几・80,5]);title('损耗函数曲线');文件名：mstem.m(实验一、三需要)程序：functionmstem(Xk)%mstem(Xk)绘制频域采样序列向量Xk的幅频特性图M=length(Xk);k=0:M-l;wk=2*k/M;%产生M点DFT对

4、应的釆样点频率(关于pi归一化值)stem(wk,abs(Xk)/.');boxon;%绘制M点DFT的幅频特性图xlabel('w/pi');ylabeir幅度')；axis([0/2/0,1.2*max(abs(Xk))]);文件名：mpplot.m(实验一需要)程序：%(2)mpplot;计算时域离散系统损耗函数和相频特性函数，并绘制曲线图。functionmpplot(B,A,Rs)%mpplot(B,AzRs)%时域离散系统损耗函数和相频特性绘图%B为系统函数分子多项式系数向量%A为系统函数分母多项式系数向量%Rs为

5、滤波器阻带最小衰减，省略则幅频曲线最小值取一80dBifnargin<3ymin=-80;elseymin=-Rs-20;end;%确定幅频曲线纵坐标最小值[H,W]=[H/W]=freqz(B/A/1000)m=abs(H);subplot(2,2,l);plot(W/pi/20*logl0(m/max(m)));gridon;xlabelCXomega/Xpi^ylabelC幅度(dB)Jaxis([0,l,ymin,5]);title(损耗函数曲线')；subplot(2,2,3);plot(W/pbp/pi);xlabG

6、l('omega/pi');ylabel('相位/pi');gridon;titleC(b)相频特性曲线J;文件名：mfftplot.m(实验一需要)程序：functionmfftplot(xn,N)%mfftplot(xn,N)计算序列向量xn的N点fft并绘制其幅频特性曲线Xk=fft(xn,N);%计算信号xn的频谱的N点采样%二工以下为绘图部分==k=0:N-l;wk=2*k/N;m=abs(Xk);mm=max(m);plot(wk,m/mm);gridon;xlabel('omega/pi');ylabel

7、(‘幅度(dB)');axis([0,2,0J.2]);titled幅度特性曲线J;文件名：mstg.m(实验四需要)程序：functionst=mstg%产生信号序列向量st,并显示st的时域波形和频谱%st=mstg返回三路调幅信号相加形成的混合信号，长度N=1600N=1600;%N为信号st的长度Fs=10000;T=l/Fs;Tp=N*T;%采样频率Fs=10kHz,Tp为采样时间t=O:T:(N-l)*T;k=O:N-l;f=k/Tp;fcl=Fs/10;%第1路调幅信号的载波频率fcl=1000Hzfml=fcl/

8、10;%第1路调幅信号的调制信号频率fml=100Hzfc2=Fs/20;%第2路调幅信号的载波频率fc2=500Hzfm2=fc?/10;%第2路调幅信号的调制信号频率fm2=50Hzfc3=FsA0；%第3路调幅信号的载波频率fc3=250Hzfm3=fc3