求教如何绘制功率谱图中的特征频率？
求教如何绘制功率谱图中的特征频率？% 本程序采用相关函数法计算给定信号的功率谱
fs=200;
fid = fopen('E:/xuexi/SBL/xiluodu/testfor90610027/mshuju/080724D10.txt');
a = fscanf(fid,'%g %g %g %g',);
fclose(fid);
n=size(a);
nfft=n(2);
num=n(1);
for i=1:num
xn=a(i,:);
cxn=xcorr(xn,'unbiased');
cxk=fft(cxn,nfft);
pxx=abs(cxk);
index=0:round(nfft/2-1);
k=index*fs/nfft;
plot_pxx=10*log10(pxx(index+1));%10*log10(pxx(index+1));
figure
plot(k,plot_pxx)
xlabel('Frequency/Hz');
ylabel('Power/W')
grid on
title('PSD')
end
for i=1:num
average(i)=mean(a(i,:));
vary(i)=(sum((a(i,:)-mean(a(i,:))).^2)/length(a(i,:)))^0.5 ;
end
average
vary
请问大家用何程序才能把图中主要频率值绘制出来？
拜求各位大侠了！！！
页:
[1]