|
楼主 |
发表于 2012-6-13 15:17
|
显示全部楼层
我检查了好多遍了。。。可能当局者迷吧。。。这是我的程序,希望大家帮我纠纠错。
clear;
fs=1000000;
fid=fopen('2012-06-12-10-55-59-1.txt','r');
N=200000;
x=fscanf(fid,'%f',N);
fclose(fid);
figure(1);
subplot(3,1,1);
plot(x);
ylabel('幅值 A');
xlabel('样本序号 n');
nfft=length(x);
y1=abs(fft(x,nfft));
y=y1/N*2;
f=(0:nfft-1)*fs/nfft/1000;
subplot(3,1,2);
plot(f,y);
axis([0,300,min(y),max(y)]);
ylabel('幅值 A');
xlabel('频率 f/kHz');
%------------------------------------------------------
wpt=wpdec(x,5,'db10','shannon');
E = wenergy(wpt);
e=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16];
subplot(3,1,3);
bar(E(e),'group');
E(e) |
|