原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《数字信号处理(第三版）》第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路调幅信号的调制信号频率fm1=100Hz
fc2=Fs
文档评论（0）