Why can't I run my code? (tfest function)

Question

Wirote Tippimarnporn 2024-2-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2083358-why-can-t-i-run-my-code-tfest-function

回答： Star Strider 2024-2-18

data_thermal_pulse.txt

在 MATLAB Online 中打开

I can run this code below in Matlab Online. But when I run this code in Matlab application, it resulted in error.

M2 = tfest(data_training,sys_init_T1);

Error using tfest

Not enough input arguments.

How can I fix this problem? I'm so confused

x = readtable('data_thermal_pulse.txt');
time = x.Var1;
T1 = x.Var2-x.Var2(1);
input = x.Var4;
Ts = 0.07;
data = iddata(T1,input,Ts);
data_training = data(1:round(0.5*size(data,1)),:);
data_testing = data((round(0.5*size(data,1)))+1:end-1,:);
sys_T1 = idtf(NaN,[1 NaN],'ioDelay',NaN);
M1 = tfest(data_training,sys_T1);
sys_init_T1 = idtf(NaN,[1 NaN],'ioDelay',NaN);
sys_init_T1.Structure.num.Value = 0.0040;
sys_init_T1.Structure.num.Minimum = 0.0001;
sys_init_T1.Structure.den.Value = [1 0.01];
sys_init_T1.Structure.den.Minimum = [1 0.001];
sys_init_T1.Structure.ioDelay.Value = 1.96;
sys_init_T1.Structure.ioDelay.Minimum = 0.01;
sys_init_T1.Structure.ioDelay.Maximum = 2.04;
M2 = tfest(data_training,sys_init_T1);

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

请先登录，再进行评论。

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

Answer 1

Sulaymon Eshkabilov 2024-2-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2083358-why-can-t-i-run-my-code-tfest-function#answer_1411003

在 MATLAB Online 中打开

data_thermal_pulse.txt

Here it is running ok:

x = readtable('data_thermal_pulse.txt');
time = x.Var1;
T1 = x.Var2-x.Var2(1);
input = x.Var4;
Ts = 0.07;
data = iddata(T1,input,Ts);
data_training = data(1:round(0.5*size(data,1)),:);
data_testing = data((round(0.5*size(data,1)))+1:end-1,:);
sys_T1 = idtf(NaN,[1 NaN],'ioDelay',NaN);
M1 = tfest(data_training,sys_T1)
M1 =
  From input "u1" to output "y1":
                   0.003866
  exp(-1.96*s) * ------------
                 s + 0.009965
 
Continuous-time identified transfer function.

Parameterization:
   Number of poles: 1   Number of zeros: 0
   Number of free coefficients: 2
   Use "tfdata", "getpvec", "getcov" for parameters and their uncertainties.

Status:                                                   
Estimated using TFEST on time domain data "data_training".
Fit to estimation data: 90.13%                            
FPE: 0.0001156, MSE: 0.0001156                            
 
sys_init_T1 = idtf(NaN,[1 NaN],'ioDelay',NaN);
sys_init_T1.Structure.num.Value = 0.0040;
sys_init_T1.Structure.num.Minimum = 0.0001;
sys_init_T1.Structure.den.Value = [1 0.01];
sys_init_T1.Structure.den.Minimum = [1 0.001];
sys_init_T1.Structure.ioDelay.Value = 1.96;
sys_init_T1.Structure.ioDelay.Minimum = 0.01;
sys_init_T1.Structure.ioDelay.Maximum = 2.04;
M2 = tfest(data_training,sys_init_T1)
M2 =
  From input "u1" to output "y1":
                   0.003869
  exp(-2.04*s) * ------------
                 s + 0.009975
 
Continuous-time identified transfer function.

Parameterization:
   Number of poles: 1   Number of zeros: 0
   Number of free coefficients: 2
   Use "tfdata", "getpvec", "getcov" for parameters and their uncertainties.

Status:                                                   
Estimated using TFEST on time domain data "data_training".
Fit to estimation data: 90.16%                            
FPE: 0.0001151, MSE: 0.000115                             
 

2 个评论
显示无隐藏无

Wirote Tippimarnporn 2024-2-18

error.jpg

This is my result

升毕 2024-2-18

I also run this code without error on matlab2021b and 2018b. M1 and M2 is the same as @Sulaymon Eshkabilov

请先登录，再进行评论。

Answer 2

Star Strider 2024-2-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2083358-why-can-t-i-run-my-code-tfest-function#answer_1411093

在 MATLAB Online 中打开

Run this from a script or your Command Window:

which tfest -all
/MATLAB/toolbox/ident/ident/tfest.m

It should give you this identical result. If it shows anything else, that is likely the problem. The solution is to re-name any other scripts, variables, or functions that are not in the System Identification Toolbox to something that makes sense in the context you are using them in, and do not overshadow any MATLAB functions.

.