推 ntust661 : 放進程式跑看看阿~ 04/23 15:47
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 114.40.69.140
※ 文章網址: https://www.ptt.cc/bbs/Math/M.1492798647.A.AD1.html
想做出 s+1 / s+2的步階響應
先用MATLAB程式模擬 H = tf([1,1],[1,2]); step(H)
用反拉氏轉換也可以
但我想把它拆成時域的疊代形式,成功跑出一樣的圖,方法如下
(因為要打程式,所以微分不用f(t+h)-f(t)/h的定義,f(t+h)是未來的值,不可能知道)
As+B Y(s)
---- = ----
Cs+D U(s)
(As+B)U(s) = (Cs+D)Y(s) 做反拉氏轉換
u(t)-u(t-h) y(t)-y(t-h)
A( ---------- ) + B(impulse) = C( ----------- ) + Dy(t)
h h
(A/h)*u(t) - (A/h)*u(t-h) + (C/h)*y(t-h)
整理後為 y(t) = -----------------------------------------
D + C/h
y的初始都為0,h先設為1好了
U(s)的反拉氏是impulse,不知存不存在,那一項要計入嗎?
我這樣的推導是正確的嗎?
--