资源描述:
《数字信号处理上机实验代码》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
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