数字信号处理matlab程序.docxVIP

数字信号处理实验报告实验1 常见离散信号的MATLAB产生和图形显示1.程序(1) 单位抽样序列 如果在时间轴上延迟了k个单位，得到即： 程序如下：N=20;k=9;x=zeros(1,N);x(10)=1;n=[k-9:N+k-10];stem(n,x);(2) 单位阶跃序列 程序如下：N=20;x=ones(1,N);x(1:3)=0;n=[-3:N-4];stem(n,x);title(单位阶跃序列);(3) 正弦序列程序如下：N=50;n=0:N-1;A=1;f=50;Fs=f*N;fai=0.5*pi;x=A*sin(2*pi*f*n/Fs+fai);stem(n,x); (4) 复正弦序列程序如下：N=50;n=0:N-1;w=2*pi/N;x=exp(j*w*n);plot(x,*);复指数序列或，它具有实部与虚部,是复正弦的数字域频率。对第一种表示形式，可以写成。如果用极坐标表示,则,若，则x(n)为衰减的复正弦,其实部和虚部分别为衰减振荡的正弦分量；若实部，则实部和虚部分别为增大的正弦分量；若，则实部和虚部分别为等幅振荡。本实验中使用的信号为时的情况，即等幅振荡信号。实验图示为极坐标下的复指数序列图。横坐标表示实部，纵坐标表示虚部。从图中明显看出，序列各点幅值相等，只有相角的周期性变化。(5) 指数序列程序如下：N=50;a=1.1;n=0:N-1;x=a.^n;plot(n,x,*);实验2 离散系统的差分方程、冲激响应和卷积分析实验目的：加深对离散系统的差分方程、单位脉冲响应和卷积分析方法的理解。实验原理：离散系统其输入、输出关系可用以下差分方程描述：?输入信号分解为单位脉冲序列，。记系统单位脉冲响应，则系统响应为如下的卷积计算式：?当时，h[n]是有限长度的（n：[0，M]），称系统为FIR系统；反之，称系统为IIR系统。在MATLAB中，可以用函数y=Filter(p,d,x)求解差分方程，也可以用函数y=Conv(x,h)计算卷积。IIR的冲激响应和阶跃响应：FIR的冲激响应和阶跃响应：1.程序N=20;p1=[1,-1];d1=[1,0.6,0.08];p2=[0,0.2,0.2,0.2,0.2,0.2];d2=1;delta=[1,zeros(1,N-1)];u=ones(1,N);labelx=0:N-1;h1=filter(p1,d1,delta);y1=filter(p1,d1,u);h2=filter(p2,d2,delta);y2=filter(p2,d2,u);subplot(2,2,1);stem(labelx,h1);title(系统1单位冲激响应,fontsize,15);subplot(2,2,2);stem(labelx,y1);title(系统1单位阶跃响应,fontsize,15);subplot(2,2,3);stem(labelx,h2);title(系统2单位冲激响应,fontsize,15);subplot(2,2,4);stem(labelx,y2);title(系统2单位阶跃响应,fontsize,15);2.程序计算结果计算结果精确到小数点后四位，只计算前20个点。1）系统1单位冲激响应：1 -1.6 0.88 -0.4 0.1696 -0.0698 0.0283 -0.0114 0.0046 -0.0018 0.0007 -0.0003 0.0001 0 0 0 0 0 0 0单位阶跃响应：1 -0.6 0.28 -0.12 0.0496 -0.0202 0.0081 -0.0033 0.0013 -0.0005 0.0002 -0.0001 0 0 0 0 0 0 0 02）系统2单位冲激响应：0 0.2 0.2 0.2 0.2 0.2 0 0 0 0 0 0 0 0 0 0 0 0 0 0单位阶跃响应：0 0.2 0.4 0.6 0.8 1 1 1 1 1 1 1 1 1 1 1 1 1 1 13.两系统单位冲激响应和单位阶跃响应图形4.理论计算结果求取单位冲激响应时，求取单位阶跃响应时，系统为因果系统，有1）对系统1：，单位冲激响应：单位阶跃响应：2）对系统2：单位冲激响应：单位阶跃响应：对比理论计算结果和计算机程序计算结果，两者保持一致。对不同的系统，因其系统函数不同和收敛域不同，单位冲激响应和单位阶跃响也不同。实验3 离散系统的频率响应分析和零