Difference in values of output by manual and code

2 views (last 30 days)
I am trying to solve a difference equation with initial conditions using 'filter' function. The manual calculation of y(n) and the one procuced by 'fliter' seems to have some difference. What is the error.
The code is shown below:
disp('Solution of a difference equation:')
% The example difference equation considered is
% y(n) = 0.5y(n-1) + x(n)
num_coef = [ 1 ]; % coefficient of x
den_coef = [1 -0.5]; % coefficient of y
n = 0:4; % Considering five samples
x = 6 * ones(1,5) % x(n) = 6u(n)
init_cond = [1]; % y(-1) = 1
y = filter(num_coef, den_coef, x, init_cond) % Solution returned by 'filter'
yM = 12 - (11/2)*(1/2).^n % Solution obtained manually
plot(y);hold on; % PLots for comparison
plot(yM,'r')
  2 Comments
David Wilson
David Wilson on 29 Apr 2019
I've always had trouble with the intial condition which is where I think your problem lies. Note that Zi (in the help) is the initial delay, not the initial value of the y.
prashanth A
prashanth A on 30 Apr 2019
Do I have to delay the signal? How did you solve it?

Sign in to comment.

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!