请问我的位移频谱图为什么会是这样？
本帖最后由 wdhd 于 2016-9-19 11:26 编辑我编了个巴特沃兹低通滤波器，用于加速度和位移的低通滤波，然后进行FFT变换。为什么最后出来的位移频谱是这样的呢?在FFT之前，我已经消除了趋势项，加速度的直流分量已经很小，但位移的直流部分仍然最大，请问怎么回事呢?以下是我的程序：
>> Fi=30;Fh=50;
>> Wp=Fi*2/Fs;Ws=Fh*2/Fs;
>> =buttord(Wp,Ws,1,40)
n =
10
Wn =
0.2224
>> =butter(n,Wn);
>> =freqz(b,a);
>> figure(1)
>>
>> subplot(3,1,1)
>> plot(Fs*w/2/pi,abs(h)) %巴特沃兹幅频特性图
>> openfile %打开数据文件
>> %消除趋势项
>> n=length(gData(1,:))
n=4864
>> t=0:1/Fs:(n-1)/Fs;
>> v=polyfit(t,gData(1,:),2);
>> y=gData(1,:)-polyval(v,t);
>> y1=filter(b,a,y);%求滤完波后的数据
>> subplot(3,1,2)
>> plot(t,y1)%滤波后的时域图
>> %FFT变换
>> N=5120;
>> y1(n:N)=0;%n以后添加(N-n)个数据，使其满足2的k次幂关系
>> y2=fft(y1,N);
>> subplot(3,1,3)
>> plot((0:n/2-1)*Fs/n,abs(y2(1:n/2))*2/n)%位移幅频图
最后得到的三幅图见附图所示，其中的位移频谱中直流分量达到了0.01
还有一个问题就是，如果我不进行频谱校正，会影响我对各频率分量作用大小的判断吗?
谢谢!!!:@)
[ 本帖最后由 liuhanxiang 于 2009-9-3 12:46 编辑 ]
请各位路过留下言吧:handshake
页:
[1]