mzbys 发表于 2016-7-3 16:51

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

本帖最后由 mzbys 于 2016-7-3 17:07 编辑

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'; %频率间隔为，只计算一半的频率
=pwelch(F1,window,noverlap,nfft,Fs,range)；
=pwelch(F2,window,noverlap,nfft,Fs,range);
=pwelch(F3,window,noverlap,nfft,Fs,range);
figure(1)
plot(f,Pxx);
figure(2)
plot(f,Pxx1);
figure(3)
plot(f,Pxx2);


结果如图：



四维 发表于 2016-7-4 09:37

是不是没去掉均值？？

ZH----过客 发表于 2016-7-4 09:57

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

vibmaster 发表于 2016-7-4 10:04

average一下，然后原图相减averge

mzbys 发表于 2016-7-4 11:44

ZH----过客 发表于 2016-7-4 09:57
之所以最开始的值非常大是因为有直流分量的存在！！！经过功率之后会更大，所以很明显！！！
不懂啊我想是不是可以取一下对数，然后这个刚开始这些较大的数可以变得小一点。

ZH----过客 发表于 2016-7-4 12:49

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

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

敷衍会致命 发表于 2016-7-4 13:31

clear,clc
close all

x = 0:1/1024:1-1/1024;
y = sin(2*pi*10*x);
z = fft(y);
figure
stem(-512:511,abs(fftshift(z))/1024);
title('无直流')


y0= y+ 0.1;
z2 = fft(y0);
figure
stem(-512:511,abs(fftshift(z2))/1024);
title('加直流')

y1=y-mean(y);
z3=fft(y1)
figure
stem(-512:511,abs(fftshift(z3))/1024);
title('去均值')大家帮忙看看，最后的去均值是这样处理吗？

ZH----过客 发表于 2016-7-5 10:55

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

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

敷衍会致命 发表于 2016-7-5 13:10

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

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

旧言虐心 发表于 2016-7-8 08:46

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

你去除均值那块感觉怪怪的。。。。
页: [1]
查看完整版本: 为什么得到的功率谱在刚开始时候值特别大？