Pour réaliser un contrôle, on transforme le modèle FEMLAB en fonction MATLAB en le paramétrant.
function [fem,Tm]=radiateur(fem0,DT,Vent)
La boucle de simulation est une boucle en temps dans laquelle on met en route le ventilateur dès que la température de contrôle dépasse le seuil. On a tracé sur la figure ci-dessous l'évolution de la température pour différents cas de fonctionnement.
heures=3600; Tmax=4; Npas=20; DT=Tmax/(Npas-1); [FEM0,Tm0]=radiateur; Tmc=zeros(1,Npas); FEM=FEM0; Tmc(1)=Tm0; vent=0.0; for i=2:Npas [FEM,Tmc(i)]=radiateur(FEM,DT*heures,vent); if (Tmc(i)>=40) vent=1.0; else if (Tmc(i)<=38) vent=0.0; end end end; figure(1); plot([0:DT:Tmax],Tmc); title('Temperature de controle');xlabel('t');ylabel('T');