数字信号处理及MALAB实现 李辉主编课本例题MATLAB实现.docVIP

第1章 例1- 离散时间信号产生。编写程序产生下列基本脉冲 （1）单位脉冲序列：起点ns，终点ne，在n0处有一单位脉冲（ns≤n0≤ne）。 （2）单位阶跃序列：起点ns，终点ne，在n0前为0，在n0处及以后为1（ns≤n0≤ne）。 （3）实数指数序列： （4）正弦序列： 程序如下： ns=0;ne=15;n0=5; n1=ns:ne;x1=[(n1-n0)==0]; %单位脉冲序列 n2=ns:ne;x2=[(n2-n0)=0]; %单位阶跃序列 n3=ns:ne;x3=(0.9).^n3; %实数指数序列 n4=ns:ne;x4=sin(n4); %正弦序列 subplot(2,2,1),stem(n1,x1); title(单位脉冲序列) subplot(2,2,2),stem(n2,x2); title(单位阶跃序列) subplot(2,2,3),stem(n3,x3); title(实数指数序列) subplot(2,2,4),stem(n4,x4); title(正弦序列) 运行结果如图1-所示。 图1- 基本脉冲序列 动态演示信号序列，卷积和的过程 程序如下： clear len_x1=20; len_x2=20; nx1=0:len_x1-1; %建立x1的时间向量 x1=0.6.^nx1; %建立x1序列 %len_x1=length(nx1); %取x1时间向量的长度 nx2=0:len_x2-1; %x2的时间向量 %len_x2=length(nx2); %取x2的时间向量的长度 x2=ones(1,len_x2); %建立x2序列 lmax=max(len_x2,len_x1); %求最长的序列 if len_x2len_x1 nx2=0;nx1=len_x2-len_x1; %若x2比x1长，对x1补nx1个0 elseif len_x2len_x1 nx1=0; nx2=len_x1-len_x2; %若x1比x2长，对x2补nx2个0 else nx2=0;len_x1=0; %若x1和x2同长，不补0 end %长者为补0长度的基础 lt=lmax; %先将x2补得与x1同长，再将两边补最大长度的0 u=[zeros(1,lt),x2,zeros(1,nx2),zeros(1,lt)]; t1=(-lt+1:2*lt); %先将x1补得与x2同长，再将左边补2倍最大长度的0 x1=[zeros(1,2*lt),x1,zeros(1,nx1)]; %将x1左右反折 hx1=fliplr(x1); N=length(hx1); y=zeros(1,3*lt); %将y存储单元初始化 for k=0:2*lt %动态演示绘图 p=[zeros(1,k),hx1(1:N-k)]; %使hx1向右循环移位 y1=u.*p; %使输入和翻转移位的脉冲过渡函数逐项相乘 yk=sum(y1); %相加 y(k+lt+1)=yk; %将结果放入数组y subplot(4,1,1);stem(t1,u); subplot(4,1,2);stem(t1,p); subplot(4,1,3);stem(t1,y1); subplot(4,1,4);stem(k,yk); %通过图形表示每一次卷积的结果 axis([-10,50,0,5]);hold on %在图形窗口上保留每一次运行的图形结果 pause(0.5);