fct_generateur_calibrage.cpp

Aller à la documentation de ce fichier.

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fct_generateur_calibrage.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:52 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#pragma hdrstop
 
#include <math.h>
#include "fct_generateur_calibrage.h"
#pragma package(smart_init)
 
 
 
FCT_GENERATEUR_CALIBRAGE::FCT_GENERATEUR_CALIBRAGE(double dn,double eniarg,double courb,double lambdaarg,double nor):FCT_GENERATEUR_ECHANTILLONAGE_ESTIMATION(dn,4,nor),eni(eniarg),lambda(lambdaarg),courbure(courb)
{
}
 
FCT_GENERATEUR_CALIBRAGE::~FCT_GENERATEUR_CALIBRAGE()
{
}
 
void FCT_GENERATEUR_CALIBRAGE::enregistrer(char *mess)
{}
void FCT_GENERATEUR_CALIBRAGE::lire(char *mess)
{}
void FCT_GENERATEUR_CALIBRAGE::construit(int n1,int n2,MG_GROUPE_TOPOLOGIQUE* mggt)
{
    pas_echantillon=0;
    pas_estimation=0;
    echantillonner();
    xmin=-55*yentredeux;
    ymin=-55*yentredeux;
    zmin=-55*yentredeux;
    xmax=55*yentredeux;
    ymax=55*yentredeux;
    zmax=55*yentredeux;
    BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
    grille.initialiser(boite.get_xmin(),boite.get_ymin(),boite.get_zmin(),boite.get_xmax(),boite.get_ymax(),boite.get_zmax(),1,1,1);
    int nb_point=lstpoint.get_nb();
    for (int i=0;i<nb_point;i++)
        grille.inserer(lstpoint.get(i));
}
 
double FCT_GENERATEUR_CALIBRAGE::get_y(void)
{
    return yentredeux;
}
 
void FCT_GENERATEUR_CALIBRAGE::echantillonner(MG_GROUPE_TOPOLOGIQUE* mggt)
{
    double deltaen=0.25*eni;
    double rc;
    if (courbure<0) rc=deltaen/((dng-eni-deltaen)*zonetrans+deltaen);
    else rc=-0.05/courbure/zonetrans/(deltaen-dng+eni)*(sqrt((deltaen*(zonetrans-1)-zonetrans*(dng-eni))*(deltaen*(zonetrans-1)-zonetrans*(dng-eni))-40*courbure*deltaen*zonetrans*(deltaen-dng+eni))+deltaen*(zonetrans-1)-zonetrans*(dng-eni));
    double dc=2*lambda*(dng-eni)*rc;
    yentredeux=dc/2. ;
    for (int i=0;i<10;i++)
    {
        double x=0.,y,z=0.;
        y=(9-2*i)*dc/2.;
        double dist=(dng-eni)*lambda;
        BOITE_3D boite(x-dist,y-dist,z-dist,x+dist,y+dist,z+dist);
        DOUBLEN<1> ecart(&eni);
        CT_POINT<1>* pt=new CT_POINT<1>(x,y,z,ecart,boite,courbure);
        lstpoint.ajouter(pt);
    }
}
 
 
 
long fibo(long n)
{
    if (n<2) return 1;
    return fibo(n-1)+fibo(n-2);
}
 
 
 
double f(double x,long nb,double *xfonc,double *fonc,double eng,double eni,double lambda,double nor,double *fonc2)
{
    FCT_GENERATEUR_CALIBRAGE carte(eng,eni,x,lambda,nor);
    carte.construit();
    double y=carte.get_y();
    double res=0.0;
    for (int i=0;i<nb;i++)
    {
        double val=carte.estimer(xfonc[i],0.,0.);
        if (fonc2!=NULL) fonc2[i]=val;
        val=val-fonc[i];
        res=res+val*val;
        val=carte.estimer(xfonc[i],y,0.);
        if (fonc2!=NULL) fonc2[i+nb]=val;
        val=val-fonc[i];
        res=res+val*val;
    }
    return res/2./nb;
}
 
 
double calibre_courbure(long nb,double *xfonc,double *fonc,double binf,double bsup,double nor,double dx,double *fonc2)
{
    double eng=fonc[nb-1];
    double eni=fonc[0];
    double lambda=(xfonc[nb-1]-xfonc[0])/(eng-eni);
    double lon=dx;
    double eps=dx/1000.;
    double a=binf;
    double b=bsup;
    int n=2;
    do
    {
        n++;
    }
    while (fibo(n)<(b-a)/lon);
    double l=a+fibo(n-2)*1./fibo(n)*(b-a);
    double u=a+fibo(n-1)*1./fibo(n)*(b-a);
    double tetai=f(l,nb,xfonc,fonc,eng,eni,lambda,nor);
    double tetaii=f(u,nb,xfonc,fonc,eng,eni,lambda,nor);
    int k=1;
    do
    {
        if (tetai>tetaii)
        {
            a=l;
            l=u;
            u=a+fibo(n-k-1)*1./fibo(n-k)*(b-a);
        }
        else
        {
            b=u;
            u=l;
            l=a+fibo(n-k-2)*1./fibo(n-k)*(b-a);
        }
        tetai=f(l,nb,xfonc,fonc,eng,eni,lambda,nor);
        tetaii=f(u,nb,xfonc,fonc,eng,eni,lambda,nor);
        k++;
    }
    while (k<n-1);
    u=l+eps;
    tetai=f(l,nb,xfonc,fonc,eng,eni,lambda,nor);
    tetaii=f(u,nb,xfonc,fonc,eng,eni,lambda,nor);
    if (tetai>tetaii) a=l;
    else b=l;
    double res=0.5*(a+b);
    if (fonc2!=NULL)
        f(res,nb,xfonc,fonc,eng,eni,lambda,nor,fonc2);
    return 0.5*(a + b);
 
}