matlab程式：
1 直接法
>> clear
>> u = wgn(1,2000,0);
>> b = [1 0.5];
>> a = [1 -1.5 0.8];
>> xn = filter(b,a,u);
>> mk = [1000,800,600];
>> for i=1:3
n = mk(i);
xn = xn(1,1:n); %取x的1000個樣本點分析
nfft=1024; %取1024點fft運算
perw=abs(fft(xn,nfft)).^2/n; %按公式先計算x(n)的離散傅利葉變換，得x(k)，然後再取其幅值的平方，並除以n
t=0:round(nfft/2-1);
f=t*n/nfft;
perw_1=10*log10(perw(t+1));
figure;
plot(f,perw_1);
title(['直接法功率譜n=' num2str(n)]);
xlabel('頻率/hz');
ylabel('功率譜密度');
mean_perw = mean(perw_1);
var_perw= std(perw_1);
disp(['perw的均值為：' ,num2str(mean_perw)]);
disp(['perw的方差為：' ,num2str(var_perw)]);
end2間接法
>> clear
>> u = wgn(1,2000,0);
>> b = [1 0.5];
>> a = [1 -1.5 0.8];
>> xn = filter(b,a,u);
>> mk = [1000,800,600];
>> nfft=1024;
>> for i=1:3
m = mk(i);
xn = xn(1,1:m);
cxn=xcorr(xn,'unbiased'); %計算序列的自相關函式
cxk=fft(cxn,nfft);
pbt=abs(cxk);
index=0:round(nfft/2-1);
k=index*m/nfft;
figure;
plot(k,10*log10(pbt(index+1)));
title(['間接法（自相關函式法）功率譜','m=',num2str(m)]);
xlabel('頻率/hz');
ylabel('功率譜密度');
mean_pbtw = mean(10*log10(pbt(index+1)));
var_pbtw= std(10*log10(pbt(index+1)));
disp(['pbtw',num2str(i),'的均值為：' ,num2str(mean_pbtw)]);
disp(['pbtw',num2str(i),'的方差為：' ,num2str(var_pbtw)]);
end3 ar法
>> clear
>> u = wgn(1,2000,0);
>> b = [1 0.5];
>> a = [1 -1.5 0.8];
>> xn = filter(b,a,u);
>> n = 1000;
>> xn = xn(1,1:n);
>> nfft=1024;
>> erder = 5;
>> [pxx,w] = pyulear(xn,erder,nfft);
>> t=0:round(nfft/2-1);
>> f = t*n/nfft;
>> pxx_1 = 10*log10(pxx(t+1));
>> plot(f,pxx_1);
>> title('ar模型法功率譜');
>> xlabel('頻率/hz');
>> ylabel('功率譜密度');
>> mean_pxx = mean(pxx_1);
>> var_pxx= std(pxx_1);
>> disp(['perw的均值為：' ,num2str(mean_pxx)]);
>> disp(['perw的方差為：' ,num2str(var_pxx)]);
數字訊號處理之經典譜估計與現代譜估計
1 直接法 clc clear all u wgn 1,2000,0 產生高斯白雜訊訊號樣本點2000個 b 1 1 0.24 a 1 1.5 0.56 濾波器係數 xn filter b,a,u u通過濾波器的輸出xn n 1000 xn xn 1,1 n 取x的1000個樣本點分析 nfft 1...
數字訊號處理與數字訊號處理器
數字訊號處理與數字訊號處理器 數字訊號處理 digital signal processing,dsp 與數字訊號處理器 digital signal processor，dsp 兩者的縮寫都是dsp，但凡是電子資訊類專業的同學，都知道dsp是這個學科最重要的專業基礎課程之一，但實際上，很多朋友對這...
數字訊號處理之緒論
1 訊號 2 系統 3 訊號處理 訊號 資訊的物理表現形式 傳遞資訊的 函式 實質是函式 訊號的分類 連續時間訊號 離散時間訊號 數碼訊號 週期訊號 非週期訊號 確定訊號 隨機訊號 能量訊號 功率訊號 按自變數與函式值的取值形式不同分類 訊號處理是研究用系統對含有資訊的訊號進行處理 變換 以獲得人們...