MATLAB: Control tool box and precision issue

control tool boxtransfer function

Dear matlab users,
I face an issue while applying transfer function to some digitized signals. It results in a filtered signal with 1e305 amplitiude…
the sampling of my signals is 20 kHz. I use the c2d function with option 'leastsquare' to get teh Z transform of the considered transfer function, I use either the zpk or ft function to draw my transfer function.
Beside this the step function diverges too, so it seems the issue is within H (either the way I define it using "tf" or "zpk")
I thought once that it could come from the "lack" of precision of the coef given by matlab and I have realized that only the symbolic toolbox can provide unlimited precision
I don't know whether there are some obvious guideline to simulate such a system or whether I face a more fundemaental problem.
(I can be more specific on H if this can help for understanding the issue)
Many thanks for your feedback
Frederic

Best Answer

So I corrected that sampling rate to 20KHz as mentioned in your question. You are right, "the only change is the "speed" of the divergence. "
Now instead of using the 'filter' command, I tried using 'lsim' command and got the following response:
%Signal_filtered.data = filter(z0,n0,Signal.data);
Signal_filtered.data = lsim(TF.TF_z,Signal.data,Signal.time);
As you can see the response is not diverging and the high amplitude of the output is clearly in line with the magnitude response (gain) of the transfer function which can be seen in the Bode plot. Now i am not sure whether precision of coefficients is the issue here or something else.
Related Question