2006-12-20

Scilab 複數運算

因為今天把工程計算機借給了某評審,導致我回家沒計算機用,事前沒想到會這樣...所以就打開Scilab來算啦!試了好久,做一下筆記。

表示虛數的 i, j,在 Scilab中用%i表示,首先要知道這個。

舉個例,我們要算某 Transfer Function T= R/[R + 1/(sC)]

令 R= 4.7KΩ, C= 0.01μF,並令 s= jω求 f= 1kHz的 Frenquency Response,這樣在Scilab中的算式就是:
T=4700./(4700+1./(%i.*2.*%pi.*1e3.*0.01e-6))

接著算Transfer Function的相移(複數平面的角度),用atan()函數來算:
atan(imag(T),real(T))

換算成度數:
atan(imag(T),real(T)).*360./%pi

atan()的第一個參數是虛數部份,第二個才是實數,要注意。real()跟imag()就是取實部跟虛部的意思。

No comments:

Post a Comment