ex10_14.m

来自「Programs for the book Advanced Engineeri」· M 代码 · 共 14 行

%EX10_14.M Compute the convolution of 
%  h(n)=(0.9)^n and u(n) so y(n)=h(n)*u(n)
%
clear
N=9;                       % N+1 terma
n=0:N;
u=ones(1,length(n));       % Input
for I=1:N+1
   h(I)=(0.9)^(I-1);       % Impulse response
   y1(I)=10*(1-(0.9)^(I)); % Response y 
end
y=conv(u,h)                % Lenght is 2N-1
% Comparison
error=y(1:N+1)-y1