数字信号处理备用总程序Word格式文档下载.docx

1、,filled）单位阶跃序列3. 在MATLAB中利用数组运算符“.”来实现一个实指数序列。如：%函数zhishu（a,n1,n2）function y=zhishu（a,n1,n2）y=（a）.n;y=zhishu（0.3,0,50）;n=0:50;实指数序列4. 在MATLAB中用函数sin或cos产生正余弦序列，如：x=11*sin（0.3*pi*n+0.2*pi）+5*cos（0.3*pi*n）;stem（n,x,正余弦序列5. 已知，试显示在区间的波形。%函数yuxian（a,n1,n2）function y=yuxian（a,n1,n2）y=3*cos（2*pi/10*（n-a）;

2、y1=yuxian（0,0,20）y2=yuxian（3,0,20）y3=yuxian（-3,0,20）subplot（3,1,1）stem（n,y1,rx（n）subplot（3,1,2）stem（n,y2,gx（n-3）subplot（3,1,3）stem（n,y3,yx（n+3）6. 参加运算的两个序列维数不同，已知，求。%函数u（n0,n1,n2）function y=u（n0,n1,n2）n1=-4:6;n2=-5:8;x1=u（-2,-4,6）;x2=u（4,-5,8）;y1=0 x1 0 0;y2=x2 ;y=y1+y2;stem（n2,y,x（n）=x1（n）+x2（n）实验2

3、1. 利用interp1函数重构采样信号已知一个连续时间信号，取最高有限带宽频率，对x（t）进行等间隔采样，采样频率为fs ，要求：（1）画出原连续时间信号x（t）的波形；程序如下f0=1; %设置正弦波信号的频率T0=1/f0;t=0:0.01: 3*T0;xt=sin（2*pi*f0*t）+1/3*sin（6*pi*f0*t）; %产生正弦波信号fs=15*f0; %设置采样频率Ts=1/fs;3* T0/Ts; %采样点数xn=1/3*sin（2*pi*f0*n*Ts）; %产生采样信号subplot（2,1,1）plot（t,xt）; %绘制原连续信号正弦信号的采样subplot（2,

4、1,2）stem（n,xn）; %绘制采样信号仿真波形：（2）设三种情况下对连续信号分别进行采样，画出采样信号波形，并利用内插公式重建原信号，对结果进行分析。程序1：xn=sin（2*pi*f0*n*Ts）+1/3*sin（6*pi*f0*n*Ts）;tn=0:Ts:3*T0; %采样时间yt=interp1（tn,xn,t,spline %内插恢复信号采样信号的恢复fs=3fhplot（t,yt）; %绘制内插恢复信号仿真波形1：程序2：fs=10*f0;仿真波形2：程序3： T0=1/f0; subplot（3,1,1）plot（t,xt）原正弦信号fs=5*f0; subplot（3,1

5、,2）采样信号的恢复fs=fh仿真波形3：分析：当fs=2fh进行采样时，在进行恢复时可以原样恢复；当fs=0;h2=dstep（b,a,50）;subplot（2,1,2）,stem（h2）,title（dstep function4. 在MATLAB中利用filtic和filter子函数求解离散系统的单位脉冲响应（1） 若，输入为，计算系统的响应，并画出图形（）；pulse=1,zeros（1,25）;a=1,-0.9,0.5;25;h1=filter（b,a,pulse）;stem（n,h1）y（n）=0;n（2）画出以及图形，并进行比较。 %求X（k）subplot（4,2,1）, s

6、tem（k,abs（Xk）;subplot（4,2,2）, stem（k,angle（Xk）;*k2）; %求G（k）subplot（4,2,3）, stem（k2,abs（Gk）;|G（k）|subplot（4,2,4）, stem（k2,angle（Gk）;arg|G（k）|*k3）; %求Y（k）subplot（4,2,5）, stem（k3,abs（Yk）;|Y（k）|subplot（4,2,6）, stem（k3,angle（Yk）;arg|Y（k）|*k4）; %求H（k）subplot（4,2,7）, stem（k4,abs（Hk）;|H（k）|subplot（4,2,8）, s

7、tem（k4,angle（Hk）;arg|H（k）|比较分析： 序列的DFT是序列频谱的等间隔采样。G（k）与X（k）的频谱是相对应的，G（k）比X（k）频谱间隔小，谱线密；Y（k）是X（k）的重复，y（n）与x（n）相比，相当于改变了对信号的取样频率，H（k）与X（k）相比较，改变了谱间隔及谱线的比例。实验44.1已知离散时间系统的系统函数为求系统的零极点，画出零极点分布图，分析系统的因果稳定性。B=0.2 0.1 0.3 0.1 0.2;A=1 -1.1 1.5 -0.7 0.3;r1=roots（B） %求分子多项式的根，即系统的零点r2=roots（A） %求分母多项式的根，即系统的极

8、点figure（1）zplane（B,A）; %调用zplane函数画零极点图系统的零极点：零点：r1 = -0.5000 + 0.8660i -0.5000 - 0.8660i 0.2500 + 0.9682i 0.2500 - 0.9682i极点：r2 = 0.2367 + 0.8915i 0.2367 - 0.8915i 0.3133 + 0.5045i 0.3133 - 0.5045i系统的因果稳定性分析：因果稳定系统的充要条件是系统函数的极点都位于Z平面单位圆内部，不包括单位圆。由系统零极图可知该系统的全部极点都在单位圆内，所以该系统因果稳定。4.2已知离散时间系统的系统函数为画出系

9、统在频率范围内的幅频响应和相频响应H,w=freqz（B,A）; %求离散系统频响特性Hf=abs（H）; %取幅度Hx=angle（H）; %取相角plot（w,Hf）;幅频特性曲线plot（w,Hx）; %画相位谱相频特性曲线4.3已知系统的单位脉冲响应，利用freqz函数画出系统在B=8,4,2,1;A=1,0,0,0;N=1024;H,w=freqz（B,A,N,whole %求hn的离散时间傅里叶变换Hksubplot（2,1,1）,plot（w/pi,abs（H）;离散时间傅里叶变换后的幅度谱subplot（2,1,2）,plot（w/pi,angle（H）;离散时间傅里叶变换后的相位谱4.4已知离散时间系统的系统函数为频率范围内的幅频响应和相频响应图形，并与课本上例2-9进行对比。B=conv（1,-exp（pi/2*1i）,1,-exp（-pi/2*1i）;A=conv（1,-0.8*exp（pi/4*1i）,1,-0.8*exp（-pi/4*1i）;plot（w,abs（H）;幅频特性plot（w,angle（H）相频特性4.5已知系统单位脉冲响应为，如果输入为，利用圆周卷积定理求系统输出（1）用FFT实现：19;hn=cos（0.5*n