Fitting transfer function to amplitude phase data
7 次查看(过去 30 天)
显示 更早的评论
Hello,
I have amplitude-phase (Bode) data from my experiment, which I have obtained by sweeping a range of frequencies. My system has the delay term exp(-Tau*s). I am looking to fit transfer function models of the form exp(-Tau*s)*(N(s)/(D(s)) to my data with the lowest possible order. Any help would be greatly appreciated!
Thanks and Regards, Hari
0 个评论
回答(2 个)
Star Strider
2015-3-17
You are fitting a matrix (amplitude and phase) with your model, so the most appropriate function to use is the Optimization Toolbox function lsqcurvefit. Your objective function will have to produce two outputs (amplitude and phase) to match your data.
Arkadiy Turevskiy
2015-3-17
编辑:Arkadiy Turevskiy
2015-3-17
Actually the function specifically created for this task is tfest from System Identification Toolbox. Take a look at the doc page and look at the section "Examples". For the lowest possible order, specify the number of poles =1, and number of zeros =0. If you do not get a good fit, try to increase the number of poles and zeros.
tfest takes your data as an input argument. Data has to be in the idfrd format. You will need to convert your mag-phase data to complex-valued response with something like this:
>>response=mag.*exp(i*phase);
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Transfer Function Models 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!