clc; 
clear all; 
close all; format short g


t=[0:0.001:10]'; %vector columna de tiempo
%parámetros del sistema
a1=4.456; a2=7.456;  a3=-0.089;  a4=12.37;
%sistema a identificar
y=a1*t+a2*sin(t)+a3*cos(t)+a4*tanh(t);
fi=[t, sin(t), cos(t), tanh(t)];
[~, theta] =mincuad(y,fi);
theta%resultado de la identificación