Problems using 'tfest' for generating transfer function

Question

RITAM BASU 2022-11-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1843678-problems-using-tfest-for-generating-transfer-function

评论： Mathieu NOE 2022-11-28

采纳的回答： Mathieu NOE

transfer_function_trial_share.zip

在 MATLAB Online 中打开

Hello guys...

I am trying to find a transfer function for my signal, that matches with the moving average block of my signal. I used tfest to recreate the signal using the generated coefficients. However, the accuracy is very low at 2%. I am wondering how to increase it, or am I doing something wrong.

In another case, to verify whether the function works, I tried with tf = 1/(s+1), generated the signals, used 'tfest' to get the transfer function, I got 0.1404/(s+0.1998) , which is way-way off.

let me know what you think of this.

Have a nice day..

data_0= iddata(T_avg_mean,T_pedal,0.01);
>> tf_estimated = tfest(data_0,2)
tf_estimated =
 
  From input "u1" to output "y1":
     -0.0188 s + 0.005945
  --------------------------
  s^2 + 7.891e-06 s + 0.4232
  
  
 Continuous-time identified transfer function.
Parameterization:
   Number of poles: 2   Number of zeros: 1
   Number of free coefficients: 4
   Use "tfdata", "getpvec", "getcov" for parameters and their uncertainties.
Status:                                            
Estimated using TFEST on time domain data "data_0".
Fit to estimation data: 2.051%                     
FPE: 128.3, MSE: 128.3  

10 个评论
显示 8更早的评论隐藏 8更早的评论

RITAM BASU 2022-11-16

Heyy.. I hope you are doing well..Sorry, I got stuck with some other deadline. I needed some time to understand your solution.

Thanks for sharing the functions and plots. It is really helpful.

I looked into your code, and what i understand is that Txy is the variable for transfer function estimate. When i checked it is a vector of complex doubles (8193x1). Then using invfreqnz we can find the numerator and Denominator for an estimated transfer function model. the variables B and A contains the coefficients of num and den respectively.

Now, i tried this for tf_1 , i.e. a transfer function (1/(s+1)), I set orders NB=0 and NA=1. I got B=2.37e-08 and A = [1 -0.99]. also the bode plots did not match.

Q : Does it make a transfer function (0/s-1) approximatly?. in which order should i read this coefficients?

next, I set orders NB=1 and NA=2; i got B=[3.18e-08 0.0099] and A = [1 -0.99 5.62e-06]. This time the bode plots did match.

Q : Does it make a transfer function (0.01/s-1) approximatly?

I am trying to get B as 1 and A as [1 1] finally.

After that, I want to try this same model with the signal generated with the mean block.. what i see now is , it is closely matching it in BODE with NB=200 and NA=100.. around 4 or 5 as orders, bode does not match.

Q : do you think making a small order (max. 4th , 5th) with same/similar response (even if it does not match exactly in Bode plots) is possible?

(Currently , I solved the problem using trial-error. like for filtering purpose, a 200 sample averaging closely matches with a tf 1/s+1..But i am trying to find a better solution )

I hope you understand what i am trying to do. Maybe some suggestions to tweak these parameters would be very helpful.

Thanks a lot for your time and for sharing the codes.

Have a nice day ahead.

RITAM BASU 2022-11-28

I dont have any problems... The process is exactly as per my requirement...But the values were not satisfactory.. that is why i made further comments.. Thanks a lot for your help again..

Have a nice day...

Mathieu NOE 2022-11-28

My pleasure !

请先登录，再进行评论。

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

Answer 1

Mathieu NOE 2022-11-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1843678-problems-using-tfest-for-generating-transfer-function#answer_1102918

移动：Mathieu NOE 2022-11-17

hello again

here a bunch of files to post process your data

I don't have the ID toolbox but I have the Signal Processing Tbx ( I use invfreqz to get the num/den model from a frequency domain model (cf Bode plot)

the Bode plot and coherence plot is obtained with my function mytfe_and_coh.m (also attached)

coherence plot is also good to have to make sure the input signal (x) has energy distributed accross the entire frequency range (otherwise the results may be inaccurate in frequency domain that has been poorly excited);

I don't put all the codes here in the window otherwise it's gonna be very clutered

see the attachements and maybe that can help you

see identified models Bode plots are precisely matchind the tfe Bode plot

NB : I identified the models in the discrete domain, use d2cm to get back to continuous if needed

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

请先登录，再进行评论。

Problems using 'tfest' for generating transfer function

10 个评论
显示 8更早的评论隐藏 8更早的评论

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Problems using 'tfest' for generating transfer function

10 个评论 显示 8更早的评论隐藏 8更早的评论

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

10 个评论
显示 8更早的评论隐藏 8更早的评论

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