Use the Matlab command conv to calculate the response of the system to a unit step input, x[n]=u[n]. Consider . Show what you type into the Matlab command window. Also, submit a plot of the output. Be sure to label your axes.
Matlab Tip: Remember that you can easily get help on a command in Matlab by typing help with the command at the prompt:
EDU>>help conv
Matlab Tip: Also, don’t forget that in Matlab when you do certain operations on an array element-by-element, you must precede the operator with a period. For instance, in Matlab, if you define an array n as follows:
EDU>>n=[0:20];
then the first term of the unit-pulse response is written as 0.5.^n
Matlab Tip: Sometimes it is useful to know the size of arrays you have created. This can be done in several ways. For example, you can use the command size:
EDU>> size(n)
ans =
1 21
The size of the array is returned. In this case, the array n has a size of 1 row by 21 columns.
Alternatively, you can switch the display from the Current Directory to the Workspace. Then all current variables and arrays and their sizes are displayed. Just click on Workspace in the upper left corner as indicated below.
Use the Matlab function recur to calculate the response of the system to a unit step input, x[n]=u[n]. Again consider . Show all that you type into the Matlab command window. Submit a plot of the output with the axes labeled.
The Matlab function recur is not a function which comes with Matlab or any of the toolboxes you have installed. Instead, it is a that must be added to your Matlab directory.
The Matlab function recur can be found in a zip file at the Mathworks site at the following location:
http://www.mathworks.com/matlabcentral/fileexchange/loadFile.do?objectId=2148
Extract the file recur.m to the current directory shown at the top of the Matlab window.
On the other hand, you can copy and paste the script below into a text file like Notepad and save it as recur.m in the current directory shown at the top of the Matlab window.
function y = recur(a,b,n,x,x0,y0);
%
% y = recur(a,b,n,x,x0,y0)
% solves for y[n] from:
% y[n] + a1*y[n-1] + a2*y[n-2]… + an*y[n-N]
% = b0*x[n] + b1*x[n-1] + … + bm*x[n-M]
%
% a, b, n, x, x0 and y0 are vectors
% a = [a1 a2 … aN]
% b = [b0 b1 … bM]
% n contains the time values for which the solution will be computed
% y0 contains the initial conditions for y, in order,
% i.e., y0 = [y[n0-N], y[n0-N+1], …,y[n0-1]]
% where n0 represents the first element of n
% x0 contains the initial conditions on x, in order
% i.e., x0 = [x[n0-M],…,x[n0-1]]
% the output, y, has length(n)
%
N = length(a);
M = length(b)-1;
if length(y0) ~= N,
error(‘Lengths of a and y0 must match’)
end
if length(x0) ~= M,
error(‘Length of x0 must match length of b-1’)
end
y = [y0 zeros(1,length(n))];
x = [x0 x];
a1 = a(length(a):-1:1); % reverses the elements in a
b1 = b(length(b):-1:1);
for i=N+1:N+length(n),
y(i) = -a1*y(i-N:i-1)’ + b1*x(i-N:i-N+M)’;
end
y = y(N+1:N+length(n));
To get information on how to use the function, at the Matlab prompt, type:
EDU>>help recur
It might be helpful to change the index of the difference equation such that the highest index is n rather than n+2.