Risposto

nlarx model initial conditions

You can prefix estimation data (both input and output signals) with nd zeros, where nd = maximum lag in the model. Initial condi...

nlarx model initial conditions

You can prefix estimation data (both input and output signals) with nd zeros, where nd = maximum lag in the model. Initial condi...

oltre 2 anni fa | 0

| accettato

Risposto

How to force tfest to estimate the process with "only real poles" ?

TFEST cannot guarantee real poles. If you can work with <=3 poles and <=1 zero, try PROCEST. This is a process model estimator...

How to force tfest to estimate the process with "only real poles" ?

TFEST cannot guarantee real poles. If you can work with <=3 poles and <=1 zero, try PROCEST. This is a process model estimator...

oltre 2 anni fa | 3

| accettato

Risposto

What is the difference between FRD and IDFRD

Both FRD and IDFRD are used to store Freqyency Response Data, that is, the complex frequency response vector (Mag.*exp(i*Phase))...

What is the difference between FRD and IDFRD

Both FRD and IDFRD are used to store Freqyency Response Data, that is, the complex frequency response vector (Mag.*exp(i*Phase))...

oltre 2 anni fa | 1

| accettato

Risposto

How to estimate the parameter in a customized transfer function

Grey-box identification is an option. You will need to write a function that takes K0 and a0 as inputs, and returns state-space ...

How to estimate the parameter in a customized transfer function

Grey-box identification is an option. You will need to write a function that takes K0 and a0 as inputs, and returns state-space ...

oltre 2 anni fa | 0

Risposto

Minimum input data resolution

Look up Nyquist Sampling Theorem. If you are sampling (hopefully with anti-aliasing) at 1Hz then you cannot theoretically captur...

Minimum input data resolution

Look up Nyquist Sampling Theorem. If you are sampling (hopefully with anti-aliasing) at 1Hz then you cannot theoretically captur...

oltre 2 anni fa | 0

| accettato

Risposto

Well-identified fitted process model does not behave like data on simulink

You are almost there. Convert the model into state-space form and use it for simulation. For initial conditions, you will need t...

Well-identified fitted process model does not behave like data on simulink

You are almost there. Convert the model into state-space form and use it for simulation. For initial conditions, you will need t...

oltre 2 anni fa | 0

| accettato

Risposto

How to calculate the transfer function for a 16 input system?

Try also TFEST. Although you might want to reduce the number of inputs by PCA or PLS analysis.

How to calculate the transfer function for a 16 input system?

Try also TFEST. Although you might want to reduce the number of inputs by PCA or PLS analysis.

oltre 2 anni fa | 0

Risposto

nlarx model compare and predict (horizon kept 1) fit totally differs

The difference between (finite-horizon) prediction and simulation is a fundamental concept, something you could read books/artic...

nlarx model compare and predict (horizon kept 1) fit totally differs

The difference between (finite-horizon) prediction and simulation is a fundamental concept, something you could read books/artic...

oltre 2 anni fa | 0

| accettato

Risposto

Can nlgreyest() estimate open-loop unstable models?

With greyest, either parameterize K matrix using the ODE function, or choose to esitmate it separately by using the "Disturbance...

Can nlgreyest() estimate open-loop unstable models?

With greyest, either parameterize K matrix using the ODE function, or choose to esitmate it separately by using the "Disturbance...

oltre 2 anni fa | 0

Risposto

System Identification of Closed Loop Data and Unstable Plant

The first reference: [1] System Identification — Theory For the User, Lennart Ljung, Section 13.4-13.5, 2nd ed, PTR Prentice Ha...

System Identification of Closed Loop Data and Unstable Plant

The first reference: [1] System Identification — Theory For the User, Lennart Ljung, Section 13.4-13.5, 2nd ed, PTR Prentice Ha...

oltre 2 anni fa | 0

| accettato

Risposto

Accessing the GUI function programatically.

I will repeat Aditya Baru's comment as an answer. The App now supports MATLAB code generation (creating a function from the tas...

Accessing the GUI function programatically.

I will repeat Aditya Baru's comment as an answer. The App now supports MATLAB code generation (creating a function from the tas...

oltre 2 anni fa | 1

| accettato

Risposto

How to set parameters of Recursive Polynomial Model Estimator in Simulink

The Recursive Polynomial Model Estimator supports single output estimations only.

How to set parameters of Recursive Polynomial Model Estimator in Simulink

The Recursive Polynomial Model Estimator supports single output estimations only.

oltre 2 anni fa | 0

| accettato

Risposto

System Identification - Frequency Domain

Make an attempt with stability enforced. opt = tfestOptions('EnforceStability', true); model=tfest(f_data,6,opt) Also, you ...

System Identification - Frequency Domain

Make an attempt with stability enforced. opt = tfestOptions('EnforceStability', true); model=tfest(f_data,6,opt) Also, you ...

oltre 2 anni fa | 2

| accettato

Risposto

System Identification Toolboox error dialog

What dataset are you using for validation? Does it suply the inputs and outputs that the model needs?

System Identification Toolboox error dialog

What dataset are you using for validation? Does it suply the inputs and outputs that the model needs?

oltre 2 anni fa | 0

Risposto

Does the order (index) of inputs and outputs matter in MIMO system identification?

Yes the order matters since within a given noise level, there are many models that can explain the data. Settings related to sea...

Does the order (index) of inputs and outputs matter in MIMO system identification?

Yes the order matters since within a given noise level, there are many models that can explain the data. Settings related to sea...

oltre 2 anni fa | 0

| accettato

Risposto

System Identification Toolbox: How can we modify the starting parameters for the armax-algorithm?

You can set the A, B, C values explicitly, as in estimatedPolymodel.A = ARCoeff Or, call the IDPOLY constructor with A, B, C p...

System Identification Toolbox: How can we modify the starting parameters for the armax-algorithm?

You can set the A, B, C values explicitly, as in estimatedPolymodel.A = ARCoeff Or, call the IDPOLY constructor with A, B, C p...

oltre 2 anni fa | 0

| accettato

Risposto

How to identify a sytem by the System Identification Toolbox that is invertable ?

Tyically yes. If you are estimating state-space model, use "feedthough" name-value pair, as in ssest(Data, order, 'Feedthrough'...

How to identify a sytem by the System Identification Toolbox that is invertable ?

Tyically yes. If you are estimating state-space model, use "feedthough" name-value pair, as in ssest(Data, order, 'Feedthrough'...

oltre 2 anni fa | 0

Risposto

how to convert a xls to a data ensemble for import into diagnostic feature app

If the data is not too big to fit into MATLAB memory, I would suggest importing it into MATLAB first as a set of tables or timet...

how to convert a xls to a data ensemble for import into diagnostic feature app

If the data is not too big to fit into MATLAB memory, I would suggest importing it into MATLAB first as a set of tables or timet...

oltre 2 anni fa | 0

Risposto

How can I find initial states for simulation?

Initial states show the effect of the environment on the system. They are not a property of the system to be determined uniquely...

How can I find initial states for simulation?

Initial states show the effect of the environment on the system. They are not a property of the system to be determined uniquely...

oltre 2 anni fa | 1

| accettato

Risposto

System Identification toolbox: how to print estimated ARMA-coefficients for each iteration step in armax-algorithm

Use "full" display option, as in: opt=armaxOptions('Display','full'); estimatedPolymodel=armax(iddata(outputdata,inputdata,tsa...

System Identification toolbox: how to print estimated ARMA-coefficients for each iteration step in armax-algorithm

Use "full" display option, as in: opt=armaxOptions('Display','full'); estimatedPolymodel=armax(iddata(outputdata,inputdata,tsa...

oltre 2 anni fa | 1

| accettato

Risposto

how to use state space model?

If you want to reproduce the response of "predict" by (manual) simulation, you will need to generate the right prediction model ...

how to use state space model?

If you want to reproduce the response of "predict" by (manual) simulation, you will need to generate the right prediction model ...

oltre 2 anni fa | 1

Risposto

how to plot on the same bode plot a manual function plot() with function bode()?

You could try: G = frd(f1(w),w); % assuming f1(w) is a complex numeric vector bode(G,f2)

how to plot on the same bode plot a manual function plot() with function bode()?

You could try: G = frd(f1(w),w); % assuming f1(w) is a complex numeric vector bode(G,f2)

oltre 2 anni fa | 0

Risposto

Residual analysis of 100% fit model using system identification toolbox

With simulated data with no noise, it is difficult to read the residual results since there is no baseline noise floor. That is,...

Residual analysis of 100% fit model using system identification toolbox

With simulated data with no noise, it is difficult to read the residual results since there is no baseline noise floor. That is,...

oltre 2 anni fa | 0

| accettato

Risposto

Fitting complex function to measurement values

Let X be your data matrix. [~,I]=unique(X(:,1),'stable'); h=X(I,2).*exp(1i*X(I,3)/180*pi); w=X(I,1)*2*pi; G=idfrd(h,w,'Ts',0...

Fitting complex function to measurement values

Let X be your data matrix. [~,I]=unique(X(:,1),'stable'); h=X(I,2).*exp(1i*X(I,3)/180*pi); w=X(I,1)*2*pi; G=idfrd(h,w,'Ts',0...

oltre 2 anni fa | 0

Risposto

MISO system identification tool box step response

Yes, use LSIM with input: U = [u, u], where u is a step signal.

MISO system identification tool box step response

Yes, use LSIM with input: U = [u, u], where u is a step signal.

oltre 2 anni fa | 0

| accettato

Risposto

Merging MISO ARX models

What you need is horizontal concatenation, not MERGE which is about statistical merger of identical (same I/Os and model structu...

Merging MISO ARX models

What you need is horizontal concatenation, not MERGE which is about statistical merger of identical (same I/Os and model structu...

oltre 2 anni fa | 0

| accettato

Risposto

How to Improve the computation speed for 'ARMAX' function?

If you have access to Parallel Computing Toolbox, you could consider replacing the for-loop with a "parfor" loop. Other things t...

How to Improve the computation speed for 'ARMAX' function?

If you have access to Parallel Computing Toolbox, you could consider replacing the for-loop with a "parfor" loop. Other things t...

oltre 2 anni fa | 0

Risposto

How to estimate the inital state with an ssest output model

In short, you cannot. SSEST is a black-box identification function (unless you pass in a full initialized @idss model as input),...

How to estimate the inital state with an ssest output model

In short, you cannot. SSEST is a black-box identification function (unless you pass in a full initialized @idss model as input),...

oltre 2 anni fa | 0

| accettato

Risposto

Why cant I predict kstep ahead when adding System Identification models?

The issue is that algebra on identified models (plus, minus, series, parallel, feedback, inv etc) are not natively supported. Th...

Why cant I predict kstep ahead when adding System Identification models?

The issue is that algebra on identified models (plus, minus, series, parallel, feedback, inv etc) are not natively supported. Th...

oltre 2 anni fa | 0

| accettato

Risposto

System Identification fitness criterion (NRMSE vs NMSE)

A good model should have goodness of fit measure less than 1, indicating that the error (measured_data - model_reponse) is small...

System Identification fitness criterion (NRMSE vs NMSE)

A good model should have goodness of fit measure less than 1, indicating that the error (measured_data - model_reponse) is small...

circa 3 anni fa | 1

| accettato