%Pisarenko算法估計

function?[AW]=pisar（xnsin_num）
load?testdat;
sin_num=3;
[N?t]=size（xn）
rxx=xcorr（xn‘biased‘）;
rxx=rxx（N:（2*sin_num）+N）;
?
%頻率估計
Rxx=toeplitz（rxx）;
ev=eig（Rxx）;
[S?i]=min（ev）;
[V?D]=eig（Rxx）;
a=V（:i）;
rts=roots（a）;
w_est=[];
for?i=1:2*sin_num
???w_est（i）=?abs（angle（rts（i）））;
end
W=（w_est/（2*pi））;
?
?
%幅度估計
mcos=[];
for?n=1:sin_num
???vcos=[];
???for?i=1:sin_num
??????vcos=[vcos?cos（n*w_est（i））];
???end
???mcos=[mcos;vcos];
end
rxx=rxx