Transfer function extraction from frequency response

28 ビュー (過去 30 日間)
Ganesh Prasad
Ganesh Prasad 2024 年 2 月 20 日
コメント済み: Ganesh Prasad 2024 年 2 月 27 日
I am trying to get a step response to a system whose magnitude and impulse reponse data for different frequencies are known.
num= xlsread('C:\Users\PrasadNGanes\Desktop\Try4_TF.xlsx');
FrequencyVector= num(:,1);
mag = num(:,2);
phase = num(:,3);
phase = deg2rad(phase);
ResponseData = horzcat(mag, phase);
sys = frd(ResponseData,FrequencyVector);
[num,den] = bode(sys);
H = tf(num,den);
I am getting an error stating :
Error using DynamicSystem/bode
The "bode" command operates on a single model when used with output arguments.
How to resolve this error.


Star Strider
Star Strider 2024 年 2 月 20 日
You will need to use the System Identification Toolbox for this.
ResponseData = mag.*exp(1j*phase);
data = idfrd(ResponseData,FrequencyVector,SamplingInterval);
H = tfest(data)
And proceed from there.
See the documentation on idfrd for details.
  17 件のコメント
Star Strider
Star Strider 2024 年 2 月 27 日
My pleasure.
Apparently that has some system defined in it, however I was unable to determine what the system was, or how it was defined. If it is a plant or filter, you may be able to get the information you want from it.
Ganesh Prasad
Ganesh Prasad 2024 年 2 月 27 日
The system(A probe in my case) is defined in the form of Scattering parameter(S2P) matrix.
The theory i read in the internet is that impulse response is the IFFT of S21 vector(transmission coefficients) and the input is determined by the convolution of output with derived inverse impulse response.
However, that theory is not working in my case.


その他の回答 (1 件)

Aquatris 2024 年 2 月 20 日
First of all, bode command does not give you numerator and denominator. It gives you magnitude and phase of the system, or FRD in your case.
In order to get a step response with the 'step' function, you need to provide a transfer function or state space representation of a system to the step function.
That is why you first need to fit a model to your FRD. You can use the system identification toolbox if you have access to it. Alternatively, if you know the structure of your model, you can also formulate a simple optimization problem or manually try to fit a model by adjusting the parameters.
Once you fit a model to your FRD, then you can use the step function to obtain the step response.


Help Center および File ExchangeTime and Frequency Domain Analysis についてさらに検索




Community Treasure Hunt

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

Start Hunting!

Translated by