Interpolation par un polynôme
%---------------------------------------------------------- % interpolation par un polynôme % utilisation de la fonction errorbar % d'aprés 'Apprendre et maitriser Matlab' %---------------------------------------------------------- echo on; %---------------------------------------------------------- % interpolation par un polynôme % fonction errorbar %--------------------------------------------------------- % ex : sigmoïde + bruit % 1 % ----------- % 1 + exp(-x) echo off; clear all; x=-5:0.1:5; % domaine n=length(x); y=1./(1+exp(-x)) + 0.05*randn(1,n); % courbe bruitée %cherche les coefficients du polynôme par les moindres carrés % qui satisfasse les données p=polyfit(x,y,5); % ordre 5 pour interpolation % évalue le polynôme sur le domaine poly=polyval(p,x); % calcul de l'erreur err=y-poly; errorbar(x,y,err); % trace la courbe et les erreurs posX1=n/2; posX2=11*n/16; text(x(posX1),y(posX1),'mesures'); text(x(posX2),poly(posX2),'modèle'); grid;
Dernière Modification : Lun 19 Mars 2007 13:28
Copyright © 1999-2010 Jean-Paul Molina Tous droits réservés.
Dernière Modification : Lun 19 Mars 2007 13:28
Copyright © 1999-2010 Jean-Paul Molina Tous droits réservés.