Inverse Z Transform recursion code not working?

Question

Joel Okanta 2021-4-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/799876-inverse-z-transform-recursion-code-not-working

回答： Shadaab Siddiqie 2021-4-15

So I was working on this code that basically uses the power series to use recursion to do an inverse Z transform from a given function. ~

syms a b i n x z % Create symbolic variables
prompt = 'What is the value of n? ';
x_n = input(prompt)
x_n = (a*n - symsum(x*(n-i))*b*i/b);
x_n=simplify(x_n)       
g=iztrans(x_n,z,n);   %Inverse Z-transform for xn[z]
g=simplify(g)       
Fs = 100;           % Sampling frequency
T = (1 / Fs)*pi;    % Sampling period (time interval between samples)
L = 10*pi;          % Length of signal duration in seconds
t= (0:T:(L-T));     % Time interval set
g=subs(g,n,t);          %Substitute the values of t into the double value h
figure(1);                         % Open a new figure
stem(t,g)             % Plot the transfer function 
                                                    
hold on 
h2=sin(pi.*t/3);    
stem(t,(h2))
hold off
xlabel('t / s');                
ylabel('h(t)');