程序
x=[4,3,2,1,0]
N = length(x)
n = [0:N-1]
k = [0:N-1]
w = exp(-j*2*pi/N)
nk = n' * k
wnk = w.^(nk)
Xk = x * wnk
結(jié)果
x =
4 3 2 1 0
N =
5
n =
0 1 2 3 4
k =
0 1 2 3 4
w =
0.3090 - 0.9511i
nk =
0 0 0 0 0
0 1 2 3 4
0 2 4 6 8
0 3 6 912
0 4 81216
wnk =
1.0000 1.0000 1.0000 1.0000 1.0000
1.0000 0.3090 - 0.9511i-0.8090 - 0.5878i-0.8090 + 0.5878i 0.3090 + 0.9511i
1.0000-0.8090 - 0.5878i 0.3090 + 0.9511i 0.3090 - 0.9511i-0.8090 + 0.5878i
1.0000-0.8090 + 0.5878i 0.3090 - 0.9511i 0.3090 + 0.9511i-0.8090 - 0.5878i
1.0000 0.3090 + 0.9511i-0.8090 + 0.5878i-0.8090 - 0.5878i 0.3090 - 0.9511i
Xk =
10.0000 2.5000 - 3.4410i 2.5000 - 0.8123i 2.5000 + 0.8123i 2.5000 + 3.4410i不是啊，是數(shù)字信號處理里的運(yùn)算?。≡趺从?jì)算??？？？這就是按照DSP的計(jì)算步驟來做的。我們平時(shí)都是計(jì)算四點(diǎn)，旋轉(zhuǎn)因子是e（-j*pi/4）=-j。如果是五點(diǎn)，旋轉(zhuǎn)因子應(yīng)為e（-j*pi/5)。不過，加一個(gè)點(diǎn)，計(jì)算量就很大了，結(jié)果從上面可以看到，比較復(fù)雜。但是步驟都一樣，另外最好利用周期性化簡一下。e（-j*pi/5）=e（-j*pi/5+j*2*pi)=e（j*pi*9/5）。