为什么得到的功率谱在刚开始时候值特别大？

mzbys

clear;
Fs=1024;          %采样频率
N=1000;           %采样点数
n=0:N-1;
t=n/Fs;     %采样时间序列s
nfft=1998;
f =Fs/2*linspace(0,1,nfft/2+1);
size(f)

F1=load('C:\Users\zhang\Documents\F1.txt');%有限元计算结果1000
F2=load('C:\Users\zhang\Documents\F2.txt');%有限元计算结果1000
F3=load('C:\Users\zhang\Documents\F3.txt');%有限元计算结果1000

window=blackman(100 ); %blackman窗
noverlap=20; %数据无重叠
range='onesided'; %频率间隔为[0 Fs/2]，只计算一半的频率
[Pxx,f]=pwelch(F1,window,noverlap,nfft,Fs,range)；
[Pxx1,f]=pwelch(F2,window,noverlap,nfft,Fs,range);
[Pxx2,f]=pwelch(F3,window,noverlap,nfft,Fs,range);
figure(1)
plot(f,Pxx);
figure(2)
plot(f,Pxx1);
figure(3)
plot(f,Pxx2);


结果如图：

四维

是不是没去掉均值？？

ZH----过客

之所以最开始的值非常大是因为有直流分量的存在！！！经过功率之后会更大，所以很明显！！！

vibmaster

average一下，然后原图相减averge

mzbys

ZH----过客 发表于 2016-7-4 09:57
之所以最开始的值非常大是因为有直流分量的存在！！！经过功率之后会更大，所以很明显！！！

不懂啊  我想是不是可以取一下对数，然后这个刚开始这些较大的数可以变得小一点。

ZH----过客

mzbys 发表于 2016-7-4 11:44
不懂啊  我想是不是可以取一下对数，然后这个刚开始这些较大的数可以变得小一点。

不是这个意思，你要想将最开始的数据变小一点需要先精确化，然后将直流分量减去，这样就可以降低最开始的直流分量值了。

敷衍会致命

1. clear,clc
   
2. close all
   
3. 
   
4. x = 0:1/1024:1-1/1024;
   
5. y = sin(2*pi*10*x);
   
6. z = fft(y);
   
7. figure
   
8. stem(-512:511,abs(fftshift(z))/1024);
   
9. title('无直流')
   
10. 
    
11. 
    
12. y0= y+ 0.1;
    
13. z2 = fft(y0);
    
14. figure
    
15. stem(-512:511,abs(fftshift(z2))/1024);
    
16. title('加直流')
    
17. 
    
18. y1=y-mean(y);
    
19. z3=fft(y1)
    
20. figure
    
21. stem(-512:511,abs(fftshift(z3))/1024);
    
22. title('去均值')

复制代码

大家帮忙看看，最后的去均值是这样处理吗？

ZH----过客

敷衍会致命 发表于 2016-7-4 13:31
大家帮忙看看，最后的去均值是这样处理吗？

理论上如果是整周期函数是对的，但是如果是半周期或者其它的一些杂乱信号，这样处理不会得到很好的效果，最开始的值还是会很大！！！

敷衍会致命

ZH----过客 发表于 2016-7-5 10:55
理论上如果是整周期函数是对的，但是如果是半周期或者其它的一些杂乱信号，这样处理不会得到很好的效果， ...

我就想知道去均值那块队不敌

旧言虐心

敷衍会致命 发表于 2016-7-5 13:10
我就想知道去均值那块队不敌

你去除均值那块感觉怪怪的。。。。