close?all
clear
clc
%%%%%%%%%%%%%%%%%%%????輸入信號???%%%%%%%%%%%%%%%%%%%
N?=?512;fs=100;IP=80;;
f1?=?0.2*fs;f2?=?0.3*fs;
fan1?=?0;?fan2?=?0;
t?=?0:1./fs:（N-1）/fs;
x?=?sin（2*pi*f1*t?+?fan1）?+?sin（2*pi*f2*t?+?fan2）;
y?=?awgn（x10）;
plot（tx）;title（‘無噪聲原始信號‘）;
figure;plot（ty）;title（‘加噪輸入信號‘）;
%%%%%%%%%%%%%%%%%%%???周期圖法求功率譜?%%%%%%%%%%%%%%%%
y_z?=?fftshift（abs（fft（yN））.^2/N）;
df=fs./N;
f=（0:df:df*（N-1））-fs./2;
figure;plot（f10*log10（y_z））;title（‘周期圖法（dB）譜估計‘）;

%%%%%%%%%%%%%%%%%%%%%%%?p=0?????%%%%%%%%%%%%%%%%%%%%%%
ef=zeros（IPN）;
eb=zeros（IPN）;
ef（1:）=y;??????%?ef（0:）
eb（1:）=y;??????%?eb（0:）
row（1）=y*y‘/N;??%ρ=row
k（1）=0;?????????%k（0）
%%%%%%%%%%%%%計算a（pi）?power（p）?p=1?2?3?...IP%%%%%%%%%%?
a=zeros（IP+1IP+1）;
power（1）=row（1）;
for?p=2:IP+1
????numerator=0;
????denominator=0;
????for?n=p:N
????numerator=?numerator+?（-2）*ef（p-1n）*eb（p-1n-1）;
????denominator=denominator+（ef（p-1n））^2+（eb（p-1n-1））^2;
????end
????k（p）=numerator./（denominator+0.0001）;??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1732??2006-11-29?13:57??功率普估計\burg.m

?????文件???????1332??2006-11-10?09:21??功率普估計\levenson.m

?????文件????????361??2006-11-10?00:05??功率普估計\zhouqitu.m

?????目錄??????????0??2009-06-07?19:36??功率普估計

-----------?---------??----------?-----??----

?????????????????3425????????????????????4