如何从Matlab

Question

我想找到使用冲击锤激发的测试样品的固有频率及其阻尼比，并通过加速度计测量响应。到目前为止，我通过傅里叶变换输入和输出数据得到了频率响应函数（Accelerance）。加速度传递函数具有对数标度。现在我陷入了下一步的困境。如何获得固有频率及其阻尼比？到目前为止的代码：

N =  125000;   %number of datapoints
fs = 12500;   %sampling frequency
ts = 1/fs;


%set the x-axis of graph
x_max = (N-1)*ts;
x_axis = 0:ts:x_max;

%data
Time = data(1,:);
Force = data(2,:);
Acceleration = data(3,:);

Input = abs(fft(Force));    %FFT input force
Output = abs(fft(Acceleration));     %FFT output acceleration

%Accelerance transfer function
Accelerance = Output./Input;

Accelerance_dB = 20*log10(Accelerance);   %Accelerance with logarithmic scale

我希望有人可以帮助我！