输入u(n)及脉冲过渡函数h(n)卷积4.asv

<matlab在电子信息中的应用>一书中个章节的配套源代码。

%《MATLAB在电子信息课程中的应用》第六章例6.4程序q604b
% 给出任意长度的输入u(n)及脉冲过渡函数h(n)卷积的演示
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

% 本程序为读者深入理解和教师讲解卷积过程提供了演示工具
clear,close all
u1=input('输入u数祖 u=（例如ones(1,5)）  ');
lu=length(u1);
h1=input('输入h数祖 h=（例如exp(-0.1*[1:10])）  ');
lh=length(h1);
sprintf('输入循环长度N大于等于%d',max(lu,lh))
N=input('输入循环长度N=（例如10）  ');
if N<max(lu,lh) disp('N太小'),end
dt=1;
ur=zeros(1,N);
for iu=1:lu ur(mod(iu,N)+1)=ur(mod(iu,N)+1)+u1(iu); end
hr=zeros(1,N);
for ih=1:lh hr(mod(ih,N)+1)=hr(mod(ih,N)+1)+h1(ih); end
hf=fliplr(hr);					% 将h的左右翻转,称为hf
for k=0:N-1
   p1=cirshft(hf,k,N);		% 使hf向右循环移位
   y1=ur.*p1*dt;					% 使输入和翻转移位的脉冲过渡函数逐项相乘，再乘dt
   yk=sum(y1);					% 相加，相当于积分
   y(mod(k,N)+1)=yk;					% 将结果放入数组y
   % 绘图，注意如何用axis命令把各子图的横坐标统一起来，使纵坐标随数据自动调整
   subplot(4,2,2);stem([0:N-1],ur)		% 用stairs是为了避免plot函数在突跳点形成的斜边
   axis([-5*dt,(5+N)*dt,min(ur),max(ur)]),hold on
   ylabel('u(t)')
   subplot(4,2,4);stem([0:N-1],p1)
   axis([-5*dt,(5+N)*dt,min(p1),max(p1)+eps])
   ylabel('h(k-t)')
   subplot(4,2,6);stem([0:N-1],y1)
   axis([-5*dt,(5+N)*dt,min(y1),max(y1)+eps])
   ylabel('s=u.*h(k-t)')
   subplot(4,2,8);stem(mod(k,N)*dt,yk,'.')	% 用stem函数表示每一次卷积积分求和的结果	
   axis([-5*dt,(5+N)*dt,floor(min(y)+eps),ceil(max(y+eps))]),hold on
   ylabel('y(k)=sum(s)*dt')
   subplot(4,2,1);stem(t1,u)		% 用stairs是为了避免plot函数在突跳点形成的斜边
   axis([-lt*dt,2*lt*dt,min(u),max(u)]),hold on
   ylabel('u(t)')
   subplot(4,2,3);stem(t1,p)
   axis([-lt*dt,2*lt*dt,min(p),max(p)])
   ylabel('h(k-t)')
   subplot(4,2,5);stem(t1,y1)
   axis([-lt*dt,2*lt*dt,min(y1),max(y1)+eps])
   ylabel('s=u.*h(k-t)')
   subplot(4,2,7);stem(k*dt,yk,'.')	% 用stem函数表示每一次卷积积分求和的结果	
   axis([-lt*dt,2*lt*dt,floor(min(y)+eps),ceil(max(y+eps))]),hold on
   ylabel('y(k)=sum(s)*dt')
    if k==round(0.8*N)  disp('暂停，按任意键继续'), pause
   else pause(0.4)
   end
end