q706.m

来自「MATLAB及在电子信息课程中的应用,电子工业出版社出版 陈怀琛 吴大正 高西全」· M 代码 · 共 43 行

%《MATLAB及在电子信息课程中的应用》第七章例7.6程序q706
% 序列卷积运算
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
% 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,'.');  
set(gcf,'color','w')                    % 设置图形背景色为白色
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)])