资源描述：
《数字信号处理上机实验代码》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、WORD格式整理文件名：tstem.m（实验一、二需要）程序：functiontstem(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（实验一、四需要）程序：functiontplot(xn,T,yn)%时域序列连续曲线绘图函数%xn:信号数据序列
2、,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（实验一、四需要）学习参考资料分享WORD格式整理程序：%(1)myplot;计算时域离散系统损耗函数并绘制曲线图。functionmyplot(B,A)%B为系统函数分子多项式系数向量%A为系统函数分母多项式系数向量[H,W]=freqz(
3、B,A,1000)m=abs(H);plot(W/pi,20*log10(m/max(m)));gridon;xlabel('omega/pi');ylabel('幅度(dB)')axis([0,1,-80,5]);title('损耗函数曲线');文件名：mstem.m（实验一、三需要）程序：functionmstem(Xk)%mstem(Xk)绘制频域采样序列向量Xk的幅频特性图M=length(Xk);k=0:M-1;wk=2*k/M;%产生M点DFT对应的采样点频率（关于pi归一化值）stem(
4、wk,abs(Xk),'.');boxon;%绘制M点DFT的幅频特性图xlabel('w/pi');ylabel('幅度');axis([0,2,0,1.2*max(abs(Xk))]);文件名：mpplot.m（实验一需要）学习参考资料分享WORD格式整理程序：%(2)mpplot;计算时域离散系统损耗函数和相频特性函数，并绘制曲线图。functionmpplot(B,A,Rs)%mpplot(B,A,Rs)%时域离散系统损耗函数和相频特性绘图%B为系统函数分子多项式系数向量%A为系统函数分母多项式
5、系数向量%Rs为滤波器阻带最小衰减，省略则幅频曲线最小值取－80dBifnargin<3ymin=-80;elseymin=-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)));gridon;xlabel('omega/pi');ylabel('幅度(dB)')axis([0,1,ymin,5]);title('损耗函数曲线');subplot(2
6、,2,3);plot(W/pi,p/pi);xlabel('omega/pi');ylabel('相位/pi');gridon;title('(b)相频特性曲线');文件名：mfftplot.m（实验一需要）学习参考资料分享WORD格式整理程序：functionmfftplot(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
7、);mm=max(m);plot(wk,m/mm);gridon;xlabel('omega/pi');ylabel('幅度(dB)');axis([0,2,0,1.2]);title('幅度特性曲线');文件名：mstg.m（实验四需要）程序：functionst=mstg%产生信号序列向量st,并显示st的时域波形和频谱%st=mstg返回三路调幅信号相加形成的混合信号，长度N=1600N=1600;%N为信号st的长度Fs=10000;T=1/Fs;Tp=N*T;%采样频率Fs=10kHz，Tp
8、为采样时间t=0:T:(N-1)*T;k=0:N-1;f=k/Tp;fc1=Fs/10;%第1路调幅信号的载波频率fc1=1000Hz学习参考资料分享WORD格式整理fm1=fc1/10;%第1路调幅信号的调制信号频率fm1=100Hzfc2=Fs/20;%第2路调幅信号的载波频率fc2=500Hzfm2=fc2/10;%第2路调幅信号的调制信号频率fm2=50Hzfc3=Fs/40;%第3路调幅信号的载波频率fc3=250H
