畫不出 step response，請指教

[Hsiao Chun Wei]

clear
clf
s = poly(0,'s');
G = syslin('c',(-6.67e10/(9.444e9*s+3.339e8*s^2+20228.333*s^3+s^4)));
T = (G)/(1+G);
t = 0:0.5:10;
y = csim('step',t,T);
plot(t,y),xgrid;

以上畫出來的圖都是一直線，請各位先進指導一下，謝謝

[Yung-Jang Lee]

分母部分 "9.444e9*s+3.339e8*s^2+20228.333*s^3+s^4 "
4 次方的係數和其他項差太大所致.

這種狀況太極端了, 實際上的物理模型應該不太有這種狀況.

2009/11/12 Hsiao Chun Wei <chunw...@gmail.com>:

[Hsiao Chun Wei]

感謝李老師的指導，
scilab:
clear
clf
kc= -6.13

s = poly(0,'s');

G = syslin('c',(-6.67e10*kc/(9.444e9*s

+3.339e8*s^2+20228.333*s^3+s^4)));
T = (G)/(1+G);
t = 0:0.5:10;
y = csim('step',t,T);
plot(t,y),xgrid;

matlab:
clear
clf
kc=-6.14
num=[-6.67e10*kc]
den=[1 20228.333 3.339e8 9.444e9 0]
G=tf(num,den)
T=feedback(G,1)
step(T)

相同的EQ在MATLAB中可以正確畫出，但在SCIALB則畫一直線，是SCILAB的code有問題，還是有其它問題，請再指點一下，感謝

[Yung-Jang Lee]

我沒使過 Matlab , 實在無法比較, 你是否能比對兩者的 T 函數是否相同 ？

在 Scilab 中為,
 --------------------------------
T

roots(T.den)

 
-->T
 T  =
 
                                                       2            3                
                    3.861D+21 + 1.365D+20s + 8.271D+15s + 4.089D+11s                 
    -------------------------------------------------------------------------------- 
                                       2            3            4            5      
    3.861D+21 + 2.257D+20s + 6.315D+18s + 1.119D+17s + 1.353D+13s + 1.077D+09s       
                   6   7                                                             
       + 40456.666s + s                                                              
 
-->
 
-->roots(T.den)
 ans  =
 
  - 28.33248               
  - 14.129061 + 32.046879i 
  - 14.129061 - 32.046879i 
  - 10100.037 + 15209.116i 
  - 10100.037 - 15209.116i 
  - 10100. + 15209.131i    
  - 10100. - 15209.131i    
 
--------------------------------------

其中分母部份產生很大的根, 可能是一長的原因. Matlab 的現象如何？

2009/11/18 Hsiao Chun Wei <chunw...@gmail.com>