Code covered by the BSD License  

Highlights from
Logiciels temps réel

Logiciels temps réel

by

 

Real-time software concerning the domains of the signal processing, feedback and regulation

[F,Yn]=matrice_F(x,y,t,d_num,d_den)
function [F,Yn]=matrice_F(x,y,t,d_num,d_den)

F=zeros(length(t),d_num+d_den);
yi=y;
F(:,d_num+d_den+1)=-y;
indice=d_num+d_den+1;
for nb_integ=1:1:d_den-1
    indice=indice-1;
    yi=cumtrapz(t,yi);
    F(:,indice)=-yi;
end;
Yn=cumtrapz(t,yi);
xi=cumtrapz(t,x);
for nb_integ=2:1:d_den-d_num
    xi=cumtrapz(t,xi);
end;
indice=indice-1;
F(:,indice)=xi;
for nb_integ=d_den-d_num+1:1:d_den
    indice=indice-1;
    xi=cumtrapz(t,xi);
    F(:,indice)=xi;
end;

Contact us