数字信号的处理上机实验代码文档格式.docx

1、n=0;t=n*T;plot（t,xn）;t/sylabel（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;omega/pi幅度（dB）axis（0,1,-80,5）;title（损耗函数曲线mstem.m（实验一、三需要）function

2、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的幅频特性图w/pi幅度axis（0,2,0,1.2*max（abs（Xk）;mpplot.m（实验一需要）%（2）mpplot;计算时域离散系统损耗函数和相频特性函数，并绘制曲线图。function mpplot（B,A,Rs）%mpplot（B,A,Rs）%时域离散系统损耗函数和相频特性绘图%Rs为滤波器阻带最小衰减，省略则幅频曲线最小值取80d

3、Bif nargin3 ymin=-80;else ymin=-Rs-20;end;%确定幅频曲线纵坐标最小值H,W=H,W=freqz（B,A,1000）subplot（2,2,1）;axis（0,1,ymin,5）;subplot（2,2,3）;plot（W/pi,p/pi）;相位/pi（b）相频特性曲线mfftplot.m（实验一需要）function mfftplot（xn,N）%mfftplot（xn,N）计算序列向量xn的N点fft并绘制其幅频特性曲线Xk=fft（xn,N）;%计算信号xn的频谱的N点采样%=以下为绘图部分=N-1;wk=2*k/N;m=abs（Xk）;mm=ma

4、x（m）;plot（wk,m/mm）;axis（0,2,0,1.2）;幅度特性曲线mstg.m（实验四需要）function st=mstg%产生信号序列向量st,并显示st的时域波形和频谱%st=mstg 返回三路调幅信号相加形成的混合信号，长度N=1600N=1600; %N为信号st的长度Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz，Tp为采样时间t=0:T:（N-1）*T;f=k/Tp;fc1=Fs/10; %第1路调幅信号的载波频率fc1=1000Hzfm1=fc1/10; %第1路调幅信号的调制信号频率fm1=100Hzfc2=Fs/20; %第2路调

5、幅信号的载波频率fc2=500Hzfm2=fc2/10; %第2路调幅信号的调制信号频率fm2=50Hzfc3=Fs/40; %第3路调幅信号的载波频率fc3=250Hzfm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25Hzxt1=cos（2*pi*fm1*t）.*cos（2*pi*fc1*t）; %产生第1路调幅信号xt2=cos（2*pi*fm2*t）.*cos（2*pi*fc2*t）; %产生第2路调幅信号xt3=cos（2*pi*fm3*t）.*cos（2*pi*fc3*t）; %产生第3路调幅信号st=xt1+xt2+xt3; %三路调幅信号相加fxt=fft（st,

6、N）; %计算信号st的频谱%=以下为绘图部分，绘制st的时域波形和幅频特性曲线=subplot（3,1,1）;plot（t,st）;grid;s（t）axis（0,Tp/8,min（st）,max（st）;（a） s（t）的波形subplot（3,1,2）;stem（f,abs（fxt）/max（abs（fxt）,（b） s（t）的频谱axis（0,Fs/5,0,1.2）;f/Hzxtg.m（实验五需要）function xt=xtg（N）%实验五信号x（t）产生,并显示信号的幅频特性曲线%xt=xtg（N） 产生一个长度为N,有加性高频噪声的单频调幅信号xt,采样频率Fs=1000Hz%载

7、波频率fc=Fs/10=100Hz,调制正弦波频率f0=fc/10=10Hz.Fs=1000;fc=Fs/10;f0=fc/10; %载波频率fc=Fs/10，单频调制信号频率为f0=Fc/10;mt=cos（2*pi*f0*t）; %产生单频正弦波调制信号mt，频率为f0ct=cos（2*pi*fc*t）; %产生载波正弦波信号ct，频率为fcxt=mt.*ct; %相乘产生单频调制信号xtnt=2*rand（1,N）-1; %产生随机噪声nt%设计高通滤波器hn,用于滤除噪声nt中的低频成分,生成高通噪声%=fp=150; fs=200;Rp=0.1;As=70; % 滤波器指标fb=fp

8、,fs;m=0,1; % 计算remezord函数所需参数f,m,devdev=10（-As/20）,（10（Rp/20）-1）/（10（Rp/20）+1）;n,fo,mo,W=remezord（fb,m,dev,Fs）; % 确定remez函数所需参数hn=remez（n,fo,mo,W）; % 调用remez函数进行设计,用于滤除噪声nt中的低频成分yt=filter（hn,1,10*nt）; %滤除随机噪声中低频成分，生成高通噪声ytxt=xt+yt; %噪声加信号fst=fft（xt,N）;plot（t,xt）;x（t）axis（0,Tp/5,min（xt）,max（xt）;（a） 信

9、号加噪声波形plot（f,abs（fst）/max（abs（fst）;（b） 信号加噪声的频谱axis（0,Fs/2,0,1.2）;10.1 系统响应及系统稳定性close all;clear all;clc;%容1：调用filter解差分方程，由系统对u（n）的响应判断稳定性A=1,-0.9;B=0.05,0.05; %系统差分方程系数向量B和Ax1n=1 1 1 1 1 1 1 1 zeros（1,50）; %产生信号x1（n）=R8（n）x2n=ones（1,128）; %产生信号x2（n）=u（n）hn=impz（B,A,58）; %求系统单位脉冲响应h（n）y=h（n）;tstem（

10、hn,y）; %调用函数tstem绘图（a） 系统单位脉冲响应h（n）box ony1n=filter（B,A,x1n）; %求系统对x1（n）的响应y1（n）subplot（2,2,2）;y1（n）tstem（y1n,y）;（b） 系统对R8（n）的响应y1（n）y2n=filter（B,A,x2n）; %求系统对x2（n）的响应y2（n）subplot（2,2,4）;y2（n）tstem（y2n,y）;（c） 系统对u（n）的响应y2（n）%容2：调用conv函数计算卷积x1n=1 1 1 1 1 1 1 1 ;h1n=ones（1,10） zeros（1,10）;h2n=1 2.5 2.5 1 zeros（1,10）;y21n=conv（h1n,x1n）;y22n=conv（h2n,x1n）;figure（2）h1（n）tstem（h1n,y）;（d） 系统单位脉冲响应h1（n）y21（n）tstem（y21n,y）;（e） h1（n）与R8（n）的卷积y21（n）h2（n）tstem（h2n,y）;（f） 系统单位脉冲响应h2（n）y22（n）tstem（y22n,y）;（g） h2（n）与R8（n）的卷积y22（n）%容3：谐振器分析un=ones（1,256）; %产生信号u（n）255;xsin=sin（