Transfer function of the highly distorted system; impact hammer testing and frequency response of the system

Question

Daniel 2019-9-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/478806-transfer-function-of-the-highly-distorted-system-impact-hammer-testing-and-frequency-response-of-th

编辑： Daniel 2019-9-5

Hello, everyone!

For couple of weeks now, I am dealing with transfer function estimation. The system which I am currently examining is highly distorted at high frequencies due to its structural dynamics. I analysed its frequency response and it is getting really messy above 2,5 kHz.

I want to compensate this system with the use of the transfer function. I have some input and output data (attached you will find the graph illustrating input and output signals). It has been carried out by means of the impact hammer testing. To cut the long story short: input - hammer hit, output - response of the system (oscillations).

I tried to handle it with the System Identification Toolbox. What I got so far is:

Could you please provide me with any hints how to get this function into time domain? The aim of my work is to reproduce the input signal without the distortions.

Thank you very much for any feedback.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Raj 2019-9-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/478806-transfer-function-of-the-highly-distorted-system-impact-hammer-testing-and-frequency-response-of-th#answer_390337

Not sure why you took the trouble of posting so much details. From what I understand, the one line question is that you have got a transfer function in 'S' domain and you want to convert this to time domain. You can use the ilaplace command to do that.

Note: You should have symbolic math toolbox in order to use this.

7 个评论
显示 5更早的评论隐藏 5更早的评论

Daniel 2019-9-4

在 MATLAB Online 中打开

Thank you for your answer, Raj!

I was trying to do it with "ilaplace" command. However, what I get then is:

val =
 
-(16195*symsum((exp(t*root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k))*root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k)^2)/(2*(37500*root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k)^2 + 12249*root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k) + 48595000)), k, 1, 3))/2 - 6076750000000*symsum(exp(t*root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k))/(2*(37500*root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k)^2 + 12249*root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k) + 48595000)), k, 1, 3) - 70327500*symsum((exp(root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k)*t)*root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k))/(2*(12249*root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k) + 37500*root(s3^3 + (12249*s3^2)/25000 + (19438*s3)/5 + 177760000, s3, k)^2 + 48595000)), k, 1, 3)

And this is the part I am struggling with. Do you happen to know the reason?

Daniel 2019-9-4

在 MATLAB Online 中打开

Of course!

% Transfer function estimation
data = iddata(output_data, input_data, 0.00001); % Sampling rate was 100 kHz
sys = tfest(data, 3);
% Inverse Laplace transform
syms s
tf_s =  (-3.24*10^4*s^2 - 2.811*10^8*s - 2.423*10^13) / (s^3 + 4.879*10^4*s^2 + 3.866*10^8*s + 1.77*10^13)
tf_time = ilaplace(tf_s)

As a results, I get:

val =
 
- 281100000*symsum((exp(root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k)*t)*root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k))/(97580*root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k) + 3*root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k)^2 + 386600000), k, 1, 3) - 32400*symsum((exp(t*root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k))*root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k)^2)/(3*root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k)^2 + 97580*root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k) + 386600000), k, 1, 3) - 24230000000000*symsum(exp(t*root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k))/(3*root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k)^2 + 97580*root(s3^3 + 48790*s3^2 + 386600000*s3 + 17700000000000, s3, k) + 386600000), k, 1, 3)

As you can see it above, I decided to rewrite the function manually as I did not know how to get numerator and denumerator in the correct data type. Apart from that, please have a look at it.

Thank you for your help.

Raj 2019-9-5

编辑：Raj 2019-9-5

"dividing the output signal by the transfer function in time domain" - No that is not correct! Transfer function is output/input in 'S' domain not in time domain. In time domain 'convolution' will come into picture. See this for better understanding:

Daniel 2019-9-5

编辑：Daniel 2019-9-5

在 MATLAB Online 中打开

close all
clear all
clc
t = linspace(0,2,200000); % Sampling rate: 100 kHz; t - time [s]
t = t';
% Import the data and decrease the number of samples
input = importdata('C:\Users\Lechowicz\Desktop\input.txt', '');
output = importdata('C:\Users\Lechowicz\Desktop\output.txt', '');
% Plot the data
figure('Name', 'Original Data', 'NumberTitle', 'on');
subplot(2,1,1)
plot(t, input)
title('Input Signal');
xlabel('Time [s]');
ylabel('Force [N]');
subplot(2,1,2)
plot(t, output)
title('Output Signal');
xlabel('Time [s]');
ylabel('Force [N]');
% Estimate the transfer function and convert it to time domain
data = iddata(output, input, 0.00001); 
sys = tfest(data, 3, 2);
[num, den] = tfdata(sys, 'v');
syms s
tf_s = poly2sym((num),s)/poly2sym((den),s);
% Get the transfer function into time domain
time_domain_tf = ilaplace(tf_s);
time_domain_tf = vpa(time_domain_tf);
% Create matrix of time domain transfer function
for i = (linspace(0,2,200000))'
    tf_time_domain_matrix = 10.753107173906815697416094377728.*exp(-3.327484297190415992653272529398.*i) - 840.46923468767926446235542117068.*exp(-246.19890346395301270080891184263.*i).*cos(19126.216965822094847808546864243.*i) - 14175.620532280572917636521940787.*exp(-246.19890346395301270080891184263.*i).*sin(19126.216965822094847808546864243.*i);
end
% I am almost sure I messed up here
output_corrected = conv(tf_time_domain_matrix, output);
% Plot the data
figure('Name', 'Reproduced Input Signal', 'NumberTitle', 'on');
plot(linspace(0,2,length(output_corrected)), output_corrected);
title('Reproduced Output Signal');
xlabel('Time [s]');
ylabel('Force [N]');

I give up. Could you please point out where I did mistake?

As far as I understand it, by correcting my output data I should get something that looks like my input data, i.e. single impact. I can't achieve it for nuts.

I attached my data as well.

All in all, it's getting better and better! Thank you very much for your help.

请先登录，再进行评论。

Transfer function of the highly distorted system; impact hammer testing and frequency response of the system

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

回答（1 个）

7 个评论
显示 5更早的评论隐藏 5更早的评论

另请参阅

类别

标签

产品

Community Treasure Hunt

Transfer function of the highly distorted system; impact hammer testing and frequency response of the system

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

回答（1 个）

7 个评论 显示 5更早的评论隐藏 5更早的评论

另请参阅

类别

标签

产品

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

7 个评论
显示 5更早的评论隐藏 5更早的评论