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

第10章实验程序

《数字信号处理(第三版）》

注意：每次实验前，请新建以下8个m文件并依实验需要（已注明）移动至matlab当前工作目录下，便于主程序调用。各程序均已验证，请直接复制即可。

1文件名：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)]);

2文件名：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)]);

3文件名：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('损耗函数曲线');

4文件名：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))]);

5文件名：mpplot.m（实验一需要）

程序：

%(2)mpplot;计算时域离散系统损耗函数和相频特性函数，并绘制曲线图。

function mpplot(B,A,Rs)

%mpplot(B,A,Rs)

%时域离散系统损耗函数和相频特性绘图

%B为系统函数分子多项式系数向量

%A为系统函数分母多项式系数向量

%Rs为滤波器阻带最小衰减，省略则幅频曲线最小值取－80dB

if nargin<3 ymin=-80;else ymin=-Rs-20;end;

%确定幅频曲线纵坐标最小值

[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)相频特性曲线');

6文件名：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('幅度特性曲线');

7文件名：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路调幅信号的调制信号频率fm1=100Hz

fc2=Fs/20; %第2路调幅信号的载波频率fc2=500Hz

fm2=fc2/10; %第2路调幅信号的调制信号频率fm2=50Hz

fc3=Fs/40;%第3路调幅信号的载波频率fc3=250Hz

fm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25Hz

xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t);

%产生第1路调幅信号

xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t);