数字信号处理上机实验代码.pdfVIP

数字信号处理上机实险代码

文件名：tstem.m(实验一、二需要)

程序：

functiontstem(xn,yn)

%时域序列绘图数

%xn:被绘图的信号数据序列，yn:绘图信号的纵坐标名称(字符串)

n=O:length(xn)-l;

stem(n,xn,.);

xlabelCn^ylabelCyn);

axis([n(end)min(xn)1.2*max(xn)]);

///

文件名：tplotm(实验一、四需要)

程序：

functiontplot(xn,T,yn)

%时域序列连续曲线绘图数

%xn:信号数据序列,yn:绘图信号的纵坐标名称(字符串)

%T为采样间隔

n=O;length(xn)-l;t=n*T;

plot(txn);

z

xlabel(t/s);ylabel(yn);

axis([t(end)min(xn)1.2*max(xn)]);

/,,

文件名：myplot.m(实验一、四需要)

计算时域离散系统损耗数并绘制曲线图。function

myplot(BA)

z

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

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

[H,W]=freqz(BA1)

m=abs(H);

plot(W/pi,2*logl(m/max(m)));gridon;

xlabel(\omega/\pi);ylabel(幅度(dB))

axis([,l,-8,5]);title(损耗数曲线’)；

文件名：mstem.m(实验一、三需要)

程序：

functionmstem(Xk)

%mstem(Xk)绘制频域采样序列向量Xk的幅频特性图

M=length(Xk);

k=:M-l;wk=2*k/M;%产生M点DFT对应的采样点频率(关于

pi归一化值)

stem(wkabs(Xk).);boxon;%绘制M点DFT的幅频特性图

//

xlabelCw/\pi);ylabelC幅度)

axis([,2,/1.2*max(abs(Xk))]);

文件名：mpplot.m(实验一需要)

%(2)mpplot;计算时域离散系统损耗数和相频特性数,并绘

制曲线图。

functionmpplot(B,ARs)

z

%mpplot(BARs)

z,

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

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

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

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

ifnargin3ymin=-8;elseymin二-Rs-2;end;%确定幅频日线

纵坐标最小值

[H,W]=[H,W]=freqz(BA1)

m=abs(H);

subplot(22l);

//

plot(W/pi2*logl(m/max(m)));gridon;

/

xlabel(\omega/\pi);ylabel(幅度(dB))

axis([,l,ymin,5]