ltitr.m

经典通信系统仿真书籍《通信系统与 MATLAB (Proakis)》源代码

%
%LTITR	Linear time-invariant time response kernel.
%	X = LTITR(A,B,U) calculates the time response of the
%	system:
%		 x[n+1] = Ax[n] + Bu[n]
%
%	to input sequence U. Matrix U must have as many columns as
%	there are inputs, u. Each row of U corresponds to a new
%	time point. LTITR returns a matrix X with as many columns
%	as there are outputs y, and with LENGTH(U) rows.
%	LTITR(A,B,U,X0) can be used if initial conditions exist.
%	Here is what it implements, in high speed:
%
%		for i=1:n
%			x(:,i) = x0;
%    			x0 = a * x0 + b * u(:,i);
%		end
%		x = x.';

%	Copyright (c) 1986-93 by the MathWorks, Inc.

% built-in function