Pages

2011年5月4日 星期三

[Matlab] 積分

迅速記錄積分的作法:
以下為將y(t)對 t 軸積分。
itg(1) = 0;
for i=1:length(t)-1
   itg(i+1) = itg(i)+y(i)*ts;
end
%t = [0:ts:t0]
%for為什麼要到length(t)-1呢?畫圖理解如下:

itg2
上圖中,itg陣列索引值走到哪,長方形面積就積到哪。所以總共有(length(t) - 1)個長方形,for迴圈只要從1走到3。
而itg的索引值為0時,沒有積任何長方形,itg索引值為1時積了一個長方形,以此類推。

沒有留言:

張貼留言