c6ex7.m

来自「数字信号处理Matlab实现?中藕糯鞰atlab实现数字信号处理Matlab」· M 代码 · 共 13 行

%	lsim plot for Example 6-7
%
v0 = -2;					% Define IC value
t = 0:.05:40;			% Define time variable
num = [1 0];
den = [1 1];			% Define H(s) numerator and denominator polynomials
v = 5*cos(2*t)-v0;	% Define input in time domain; include IC
y_sim = lsim(num, den, v, t);	% Use lsim( ) to get output (without y = lsim( ), plot is done)
y_anal = (1-v0)*exp(-t)+4*(cos(2*t)-.5*sin(2*t));	% Analytical output for comparison
% Plot simulated and analytically computed output
subplot(2,1,1),plot(t, y_sim), xlabel('t'), ylabel('sim. y(t)')
subplot(2,1,2),plot(t, y_anal), xlabel('t'), ylabel('anal. y(t)')