program_4_3.m

这是国外一本经典教材附带的程序

% Program 4_3
% Determination of the input signal from
% given impulse response and output signal
% of same length N
h1 = input('Type in the impulse response samples = ');
y = input('Type in the output response samples = ');
N = length(y);
x = [y(1)/h1(1) zeros(1,N-1)];
for k = 2:N;
	x(k) = (y(k) - fliplr(h1(2:k)*x(1:k-1)'))/h1(1);
end
disp('Input Samples');disp(x);