連續(xù)域得功率譜公式：
圖1
圖2
從圖1和圖2可以看出，當(dāng)頻率固定，比如我們要求f=1hz時(shí)得功率譜密度，變量
必須從負(fù)無(wú)窮變到正無(wú)窮，然后積分求和，才能把f=1hz時(shí)得功率譜密度求出來(lái)。
離散領(lǐng)域：
對(duì)于AR模型來(lái)說(shuō)，其功率譜就是：
圖3
圖4
參考圖3和圖4，當(dāng)頻率w=0時(shí)，此時(shí)得功率譜P(0)得分母直接就是
（1+a1+a2+......+ap）^2,也就是P(0)與z無(wú)關(guān)！
赤池弘次在推導(dǎo)赤池信息準(zhǔn)則（AIC）得論文《STATISTICAL PREDICTOR ENTIFICATION》中就用到了這個(gè)推論：
其中
那么，如何由圖3和圖4來(lái)用程序求出功率譜呢？
圖5
其中 S: 拉普拉斯變換得算子； Ts:采樣周期
采樣周期就是采樣間隔，比如一秒鐘采樣了1000個(gè)數(shù)據(jù)，那么Ts就是1/1000,當(dāng)然Ts也可以假定為1。
以圖5中得Z代入圖4，可求出A(Z),下面得程序假設(shè)采樣點(diǎn)為1000個(gè)，圖4中得p為15：
for f= 1 : 1000 頻率得變化范圍通常和采樣點(diǎn)數(shù)一樣，當(dāng)然也可以不一樣
H=1;
for m=1:p
H =H+a(m)*exp(-j*m*2*pi*f/1000);
end;
這里求出得H就是復(fù)數(shù)A(z)，再求其模即可。
