|
|
楼主 |
发表于 2011-6-14 16:51
|
显示全部楼层
本帖最后由 comoma2011 于 2011-6-14 19:50 编辑
回复 10 # ChaChing 的帖子
谢谢大家的指导。
关于
问题3：我也不知道。但是一直得不到好的频谱图。
问题4：大概知道一阶固有频率在3Hz左右。但是频谱图上不明显。这也是我困惑的地方。
问题5：我也试过更改不同的参数。我在下面贴上程序（改进使用论坛以前一个朋友的程序），还请进一步指教.
问题6：我认为不少，因为关心的桥的频率在15Hz以下.
程序：
clear;
clc;
clf;
close all;
A = importdata('Experimentaldata.txt');
x=A(:,1);
x=x';
y=A(:,3);
y=y';
% Display the shape of signal(Time Domains)
subplot(2,1,1);
plot(x,y)
xlabel('Time (s)');
ylabel('Acceleration');
title('Signal shape(time domains)');
grid on;
% FFT
y=y-mean(y);
Fs=length(x)/(max(x)-min(x));
N=length(y);
z=fft(y,N);
% Spectral analysis
f=(0:N-1)*Fs/N;
Mag=2*abs(z)/N;
Pyy=Mag.^2;
% Display the spectrum (frequency domains)
subplot(2,1,2)
% plot(f(1:N/2),Pyy(1:N/2),'r')
plot(f(1:N/2),Mag(1:N/2),'r')
%
% If replace Pyy with Mag, we will get the amplitude-frequency figure
% axis([min(f(1:N/2)) max(f(1:N/2)) 1.1*floor(min(Pyy(1:N/2))) 1.1*ceil(max(Pyy(1:N/2)))])
axis([min(f(1:N/2)) max(f(1:N/2)) 1.1*floor(min(Mag(1:N/2))) 1.1*ceil(max(Mag(1:N/2)))])
xlabel('Frequency (Hz)')
ylabel('Energy')
title('Spectrum (frequency domains)')
grid on;
% Xlim([0,10]);
% Get the frequency and period at the maximum energy value
[a b]=max(Pyy(1:N/2));
fprintf('\n Results obtained by FFT:\n')
fprintf(' FFT_f = %1.3f Hz\n',f(b))
fprintf(' FFT_T = %1.3f s\n',1/f(b))
|
|
[复制链接]
发表于 2011-6-9 17:29
提升卡
变色卡
显身卡
发表于 2011-6-12 12:16
发表于 2011-6-14 13:39
