HOW to draw FRF graph and mode shape

I have experimental data of impact hammer for FRF graph and mode shape

AND i have to find natural frequency

that date contain Amplitude and phase about Hz of each 55th places

but i don`t know how to draw FRF(frequency response function) graph and mode shape like under picture

can I receive help? (**attached data!!)

[KSSV]

Attach data..

[seong won jeong]

here is excel data!

[KSSV]

Do you have any other information? Or only this file?

[seong won jeong]

i have only this file !

it's impact hammer test!

problom saids draw using amplitude and phase for frf and mode shape and think Reciprocity characteristic

sorry for poor explain...

[Muhammad Affan Arif]

I have data files from an impact hammer test. The experiment used 1 input (hammer voltage) and there were 6 accelerometers. I also know the sampling frequency. The experiments were conducted 7 times. So, I need to average the input and response for minimal error.

1. How to get the data in the form of amplitude and phase?
2. How to get the corresponding frequencies? Since the data only shows the input values in 1 column and acceleration values from 6 accelerometers in 6 seperate columns.
3. How to plot FRF against frequencies?

Answer 1

[KSSV]

data = xlsread("11bun_11.xlsx") ;
f = data(:,1:3:end) ; f = data(:,1) ; 

amp = data(:,2:3:end) ; 
ph = data(:,3:3:end) ; 

f10 = f(21) ; 
amp10 = amp(21,:) ; 
ph10 = ph(21,:) ; 

x = linspace(0,30,11) ;
y = linspace(0,150,5) ; 

surf(x,y,reshape(amp10,5,11))

[seong won jeong]

wowow... thank you very much!

but shamefully.. problem is make 55mesh and Draw the value of the amplitude of the natural frequency and the sign of the Phase on a point in Mesh. could you draw that?

[KSSV]

data = xlsread("11bun_11.xlsx") ; 

f = data(:,1:3:end) ; f = data(:,1) ; 

amp = data(:,2:3:end) ; 
ph = data(:,3:3:end) ; 

f10 = f(21) ; 
amp10 = amp(21,:) ; 
ph10 = ph(21,:) ; 

x = linspace(0,30,11) ;
y = linspace(0,150,5) ; 
[X,Y] = meshgrid(x,y) ;
Z = reshape(amp10,5,11) ; 

m = 100 ; n= 200 ; 
xi =  linspace(0,30,m) ; 
yi = linspace(0,150,n) ; 
[Xi,Yi] = meshgrid(xi,yi) ;
Zi = interp2(X,Y,Z,Xi,Yi) ;

surf(Xi,Yi,Zi)
shading interp 

[seong won jeong]

thank you very much!!! Are you an angel?

thank you very much!!!