- Load the data into MATLAB: You can load the demodulated magnitude R, phase (theta), in-phase component X, and quadrature component Y values for both input signals into MATLAB. You can load the data using the load function or read it from a file using the appropriate file I/O functions.
- Calculate the FRF: To calculate the FRF, you need to divide the demodulated values of signal input 1 by the demodulated values of signal input 2.
- Convert the FRF to inertance and receptance: To visualize the plot in the inertance (acceleration) and receptance (displacement) format, you need to convert the FRF from complex to polar coordinates (magnitude and phase). Then you can use the following equations to convert the FRF to inertance and receptance: Plot the FRF in the inertance and receptance format: Finally, you can plot the FRF in the inertance and receptance format using the semilogx function to plot the magnitude and semilogy function to plot the phase:

# How to find the system FRF with available input and output FRF?

17 views (last 30 days)

Show older comments

Yaser Arafath Gulam Dhasthagir
on 13 May 2023

Commented: Yaser Arafath Gulam Dhasthagir
on 14 Jun 2023

Hello,

I did a frequency sweep from 10hz-15kHz through a lock in amplifier (LIA) with two inputs to the LIA. The LIA gives me the demodulated magnitude R, phase (theta), in-phase component X and quadrature component Y for both the input signal.

Now i would like to calculate the total FRF which is the output of the system divided by the input. In this case, demodulated values of signal input 1 is the output and demodulated values of signal input 2 is the input of the system. I would also like to visualise the plot in inertance (acceleration) and receptance (displacement) format.

Can anyone guide me on how to compute the total FRF? Because I was expecting the peak resonance to occur at -90 degree phase and the real part has to be zeo at that frequency. I attached a matlab file, where you can find the r,phase,x and y values. Here the demod 1 is the demodulated values of signal input 1 & demod 4 is the demodulated values of signal input 2.

##### 0 Comments

### Answers (1)

Gokul Nath S J
on 25 May 2023

Hi Yaser,

To compute the total FRF and visualize the plot in inertance (acceleration) and receptance (displacement) format, you can follow the steps below:

FRF = demod1 ./ demod4;

Inertance = -imag(FRF) / (2*pi*f)^2

Receptance = real(FRF) / (2*pi*f)

figure;

semilogx(freq, abs(Inertance));

xlabel('Frequency (Hz)');

ylabel('Inertance (m/s^2)');

figure;

semilogx(freq, abs(Receptance));

xlabel('Frequency (Hz)');

ylabel('Receptance (m/N)');

with regards,

Gokul Nath S J

### See Also

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!