q706.m

这是matlab在电子信息课程中的应用

%《MATLAB在电子信息课程中的应用》第七章例7.6程序q706
% 序列卷积运算
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

% Nx：输入序列x(n)长度
% Nh：h(n)长度
% m：x(n)的移位样点数
clear;close all
Nx=20;Nh=10;m=5;
n=0:Nx-1;
x1=(0.9).^n;   %x1(n)
%产生x2(n)=x1(n-m)
x2=zeros(1,Nx+m);
for k=m+1:m+Nx
    x2(k)=x1(k-m);
end
nh=0:Nh-1;
h1=ones(1,Nh);   %h1(n)
h2=h1;             %h2(n)
y1=conv(x1,h1);  %y1(n)
y2=conv(x2,h2);  %y2(n)
Ny1=length(y1);  %y1(n)长度
Ny2=length(y2);  %y2(n)长度
subplot(3,2,1);stem(n,x1,'.');  
xlabel('n');ylabel('x1(n)');axis([0,Ny1,0,1]);
subplot(3,2,3);stem(nh,h1,'.');  
xlabel('n');ylabel('h1(n)');axis([0,Ny1,0,1]);
n=0:Ny1-1;
subplot(3,2,5);stem(n,y1,'.');  
xlabel('n');ylabel('y1(n)');axis([0,Ny1,0,max(y1)]);
n=0:length(x2)-1;
subplot(3,2,2);stem(n,x2,'.');
xlabel('n');ylabel('x2(n)');axis([0,Ny2,0,1]);
subplot(3,2,4);stem(nh,h2,'.');  
xlabel('n');ylabel('h2(n)');axis([0,Ny2,0,1]);
n=0:Ny2-1;
subplot(3,2,6);stem(n,y2,'.');
xlabel('n');ylabel('y2(n)');axis([0,Ny2,0,max(y2)])