From: Susan on
I have a TF which has a numerator larger than my denominator. I used the following code to determine the remainder and quotient. In the end I want a plot of the remainder divided by the original denominator with a y-shift by the quotient.

v=[0.1741 137.8 511.8 888.8 662.1 166.9 0 0];%numerator of original TF
u=[0.002721 2.147 33.96 75.43 89.24 0 0 0];%denominator of original TF
[q,r] = deconv(v,u);
%q = 63.9838 %showing the quotient
%r = [0 .4 -1661.1 -3937.5 -5047.8 0166.9 0 0] %remainder

t=1:1:1000;
Duty_sys=tf(r, u)
[Y,T] = step(Duty_sys,t)
plot(T,Y+q);%I'm shifting the plot up by the value of the quotient.

I'm not getting the plot my advisor expected. My function is a ramp, but my advisor is saying I should be getting an "s-shaped" curve. Is this code the right way to divide a numerator by a denominator to get a remainder and quotient?