Trouver une fonction de transfer à partir de courbes

Salut à tous !!

Alors voilà : actuellement en stage je dois trouver la fonction de transfert suite à une essai en soufflerie. On mesure un flux de chaleur. Voilà où j'en suis :



Donc voici ci-dessus la commande qui a été programmé dans la soufflerie, c'est la courbe bleue qui est une succession d'échelon. C'est ce que demande la machine en entrée. Ensuite la réponse de la sonde machine en verte. En fait j'ai extrait quelques points témoins de la réponse originale qui contenant plus de 4000 points, or je veux tout faire sous Matlab et donc je doi avec le même nombre de données en entrée et en sortie. Et enfin la courbe rouge qui ne représente que les calculs qui ont amené à la courbe d'échelons, donc la rouge est inutile pour la suite.


Ensuite voici les commandes Matlab sur lesquelles je pars :

clear
clc
[numc] = xlsread('excelversmatlabconsigne2.xls');
[numr] = xlsread('excelversmatlabreponse.xls');
[tempsc] = xlsread('tempsconsigne2.xls');
[tempsr] = xlsread('tempsreponse.xls');

figure
plot(tempsc,numc,tempsr,numr)

Pnum=polyfit(tempsr,numr,2)
Pden=polyfit(tempsc,numc,2)
sys=tf(Pnum,Pden)
figure
bode(sys)

Voici ce que renvoie Matlab :
Pnum =

-0.2375 350.0361


Pden =

-0.2064 344.0407


sys =

0.2375 s - 350
--------------
0.2064 s - 344

Alors déjà je voudrais savoir pourquoi les signes sont inversés entre les polynômes "Pnum" et "Pden" et dans la fraction de "sys" ?? Et ensuite ce tracé change beaucoup si je modifie le degré des polynômes "Pnum" et "Pden" , à l'ordre 2 sur ce test. Donc quel ordre choisir ?

J'obtiens donc logiquement une fonction de transfert, mais qui n'est pas la même si j'utilsie la System Identification TOolbox avec une entrée type "time domain" et en entrée mes données "numc" qui sont les point de la consignes et le sortie "numr" qui sont les points de la courbe de réponse. Donc dois-je préférer cette toolbox ?



Et ensuite la deuxième image est le tracé de Bode avec la courbe d'échelon en entrée et la verte en sortie. Alors le boss m'a dit de transformer les échelons avec les transformées de Laplace pour la comparer à la réponse verte. Sauriez-vous faire sur Matlab ?

Merci d'avance ! je sais ça fait beaucoup de question mais au moins tout est posé !