Utilisateur:Ereduverseau/Travaux personnels en attente de validation

Harmonique simple de régression unipulsatoire

Méthode

L'équation de la courbe de l'harmonique peut se ramener à la forme :

$y=y_{0}+S*sin(\omega *t)+C*(1-cos(\omega *t))$

Avec les propriétés suivantes :

- La somme des carrés des écarts des $y_{k}$ aux $y(k)$ est minimale pour un $\omega$ donné ( propriété de base de la régression )

Propriété P1 : $\sum _{i=-k}^{i=+k}(y_{i}-y(i))^{2}$ minimale

Parmi toutes ces harmoniques obtenues pour tous les $\omega$, une pour chaque, certaines sont plus probables que d'autres. Exemples :

_ celles où $y_{max}-y_{min}$ est minimale ou inférieure à une valeur donnée$\epsilon$. (Condition C1)

_ celles où la valeur absolue de la différence entre $y_{max}-y_{min}$ et {la moyenne des k plus grands y – la moyenne de k plus petits y} hors $y_{0}$ est minimale ou inférieure à une valeur donnée $\epsilon$(Condition C2).

_ celles où la puissance de y minimale  : $(C^{2}+S^{2})\omega$ est minimale ou inférieure à une valeur donnée $\epsilon$(Condition C3).

_ ou d'autres conditions. On quitte le domaines des maths pures pour les maths appliquées.

Application

2k+1 couples

Méthode directe

Calcul de S et C de régression

P1 : Se traduit par le système : $(\sum _{i=-k}^{i=+k}(y_{i}-y(i))^{2})'_{S}=0$ et $(\sum _{i=-k}^{i=+k}(y_{i}-y(i))^{2})'_{C}=0$

${\begin{cases}(\sum _{i=-k}^{i=+k}(y_{i}-y(i))^{2})'_{S}=0\\(\sum _{i=-k}^{i=+k}(y_{i}-y(i))^{2})'_{C}=0\\\end{cases}}$

${\begin{cases}(\sum _{i=-k}^{i=+k}(y(i)-y_{0}-S*sin(\omega i)-C*(1-cos(\omega i)))^{2})'_{s}=0\\(\sum _{i=-k}^{i=+k}(y(i)-y_{0}-S*sin(\omega i)-C*(1-cos(\omega i)))^{2})'_{c}=0\end{cases}}$

${\begin{cases}\sum _{i=-k}^{i=+k}sin(\omega i)(y(i)-y_{0}-S*sin(\omega i)-C*(1-cos(\omega i))=0\\(\sum _{i=-k}^{i=+k}(1-cos(\omega i))(y(i)-y_{0}-S*sin(\omega i)-C*(1-cos(\omega i))=0\end{cases}}$

${\begin{cases}\sum _{i=-k}^{i=+k}sin(\omega i)(y(i)-y_{0})=Ssin(\omega i)^{2}+Csin(\omega i)(1-cos(\omega i))\\(\sum _{i=-k}^{i=+k}(1-cos(\omega i))(y(i)-y_{0})=Ssin(\omega i)(1-cos(\omega i))+C(1-cos(\omega i))^{2}\end{cases}}$

Et après simplification vu les parités :

${\begin{cases}\sum _{i=-k}^{i=+k}sin(\omega i)y(i)=Ssin(\omega i)^{2}\\(\sum _{i=-k}^{i=+k}(1-cos(\omega i))y(i)=C(1-cos(\omega i))^{2}\end{cases}}$

D'où S et C : $S={\frac {\sum _{i=-k}^{i=+k}y(i)sin(\omega i)}{\sum _{i=-k}^{i=+k}sin^{2}(\omega i)}}$ et $C={\frac {\sum _{i=-k}^{i=+k}y(i)(1-cos(\omega i))}{\sum _{i=-k}^{i=+k}(1-cos(\omega i))^{2}}}$

Calcul de wt optimal et des maxi-mini

$y=y_{0}+S*sin(\omega *t)+C*(1-cos(\omega *t))$

$y'(t)=\omega Scos(\omega *t)+\omega Csin(\omega *t))$

Pour $t_{M}$ et $t_{m}$ on obtient les extrèmum et $y'(t_{M};t_{m})=0$

D'où $\omega *t_{m,m}$ par $tan(\omega *t_{M,m})=-{\frac {S}{C}}$

$\omega t_{M}=atan(-{\frac {S}{C}})$

$\omega t_{m}=atan(-{\frac {S}{C}})+\pi$

Aux extrémum, en substituant $\omega *t_{M,m}$, on obtient $y_{M},y_{m}$ soit :

$y(w)=y_{0}+S*sin(atan(-{\frac {S}{C}}))+C*(1-cos(atan(-{\frac {S}{C}})))$

Avec la condition C1 plus haut , il s'agira de minimiser la quantité :

$y(w)=S*sin(atan(-{\frac {S}{C}}))+C*(1-cos(atan(-{\frac {S}{C}})))$

Un programme ou le tableur donneront la valeur de w optimale cherchée $\omega _{opt}$.

$y=y_{0}+S(\omega _{opt})sin(\omega _{opt}*t)+C(\omega _{opt})*(1-cos(\omega _{opt}*t))$

Méthode par décomposition paritaire

$y(i)=y_{I}(i)+y_{P}(i)$ avec $y_{I}=0+Ssin(\omega t)$ et $y_{P}=y_{0}+C(1-cos(\omega t))$

avec $y_{I}(i)={\frac {y(i)-y(-i)}{2}}$ et $y_{P}(i)={\frac {y(i)+y(-i)}{2}}$ par définition des parties impaires et paires d'une fonction ou d'un échantillonnage.

Chercher la régressions de y, ce peut être dans certains cas de superpositions de phénomènes, chercher séparément les régressions propres de y_I et y_P puis les ajouter.

On écrit les conditions de régression pour chaque partie et on les traduit :

${\begin{cases}(\sum _{i=-k}^{i=+k}(y_{Ii}-y(i))^{2})'_{S}=0\\(\sum _{i=-k}^{i=+k}(y_{Pi}-y(i))^{2})'_{C}=0\end{cases}}$

${\begin{cases}\sum _{i=-k}^{i=+k}sin(\omega i)(y_{I}i-y(i))=0\\\sum _{i=-k}^{i=+k}(1-cos(\omega i))(y_{P}i-y(i))=0\end{cases}}$

${\begin{cases}\sum _{i=-k}^{i=+k}sin(\omega i)(y_{I}i-Ssin(\omega t))=0\\\sum _{i=-k}^{i=+k}(1-cos(\omega t))(y_{P}i-y_{0}+C(1-cos(\omega t))=0\end{cases}}$

D'où S et C : $S={\frac {\sum _{i=-k}^{i=+k}y_{Ii}sin(\omega i)}{\sum _{i=-k}^{i=+k}sin^{2}(\omega i)}}$ et $C={\frac {\sum _{i=-k}^{i=+k}y_{Pi}(1-cos(\omega i))}{\sum _{i=-k}^{i=+k}(1-cos(\omega i))^{2}}}$

Or : ${\frac {\sum _{i=-k}^{i=+k}y_{Pi}sin(\omega i)}{\sum _{i=-k}^{i=+k}sin^{2}(\omega i)}}=0$

Et : ${\frac {\sum _{i=-k}^{i=+k}y_{Ii}(1-cos(\omega i))}{\sum _{i=-k}^{i=+k}(1-cos(\omega i))^{2}}}=0$

D'où en sommant terme à terme :

$S={\frac {\sum _{i=-k}^{i=+k}y(i)sin(\omega i)}{\sum _{i=-k}^{i=+k}sin^{2}(\omega i)}}$ et $C={\frac {\sum _{i=-k}^{i=+k}y(i)(1-cos(\omega i))}{\sum _{i=-k}^{i=+k}(1-cos(\omega i))^{2}}}$

Ce qui par ailleurs affranchit donc de la parité du nombre de couples ( 2k ou 2k+1 )

La suite des calculs 122112 est la même afin d'obtenir une valeur unique de $\omega$

$\omega t_{M}=atan(-{\frac {S}{C}})$

$\omega t_{m}=atan(-{\frac {S}{C}})+\pi$

Avec la condition C1 plus haut , il s'agira de minimiser la quantité :

$y(w)=S*sin(atan(-{\frac {S}{C}}))+C*(1-cos(atan(-{\frac {S}{C}})))$

Un programme ou le tableur donneront la valeur de w optimale cherchée $\omega _{opt}$.

$y=y_{0}+S(\omega _{opt})sin(\omega _{opt}*t)+C(\omega _{opt})*(1-cos(\omega _{opt}*t))$

Si le nombre de couples est pair $t=0,5+k$ et $t=0,5-k$ sans $y_{0}$

4k+2 couples

couples $(x_{-2k-1},y_{+2k+1})$

D'où S et C : $S={\frac {\sum _{i=-2k-1}^{i=+2k+1}y_{Ii}sin(\omega i)}{\sum _{i=-2k-1}^{i=+2k+1}sin^{2}(\omega i)}}$ et $C={\frac {\sum _{i=-2k+1}^{i=+2k+1}y_{Pi}(1-cos(\omega i))}{\sum _{i=-2k-1}^{i=+2k+1}(1-cos(\omega i))^{2}}}$

$\omega t_{M}=atan(-{\frac {S}{C}})$

$\omega t_{m}=atan(-{\frac {S}{C}})+\pi$

Aux extrémum, en substituant $\omega *t_{M,m}$, on obtient $y_{M},y_{m}$ soit :

$y(w)=S*sin(atan(-{\frac {S}{C}}))+C*(1-cos(atan(-{\frac {S}{C}})))$

à minimiser pour trouver $w_{o}pt$

$y(2t+0,5)=S(\omega _{opt})sin(\omega _{opt}*(2t+0,5))+C(\omega _{opt})*(1-cos(\omega _{opt}*(2t+0,5)))$

Par la suite on emploiera la méthode directe par partition paritaire des échantillonnages et des fonctions de régression.

Harmonique double de pulsations w et hw

2k+1 couples

$y=y_{0}+S*sin(\omega *t)+S_{h}*sin(\omega *ht)+C*(1-cos(\omega *t))+C_{h}*(1-cos(\omega *ht))$

$y(i)=y_{I}(i)+y_{P}(i)$ avec $y_{I}=0+Ssin(\omega t)+S_{h}sin(h\omega t)$ et $y_{P}=y_{0}+C(1-cos(\omega t))+C(1-cos(\omega t))$

avec $y_{I}(i)={\frac {y(i)-y(-i)}{2}}$ et $y_{P}(i)={\frac {y(i)+y(-i)}{2}}$ par définition des parties impaires et paires d'une fonction ou d'un échantillonnage.

Les conditions de régression sont :

${\begin{cases}(\sum _{i=-k}^{i=+k}(y_{Ii}-y(i))^{2})'_{S}=0\\(\sum _{i=-k}^{i=+k}(y_{Ii}-y(i))^{2})'_{S_{h}}=0\end{cases}}$

${\begin{cases}(\sum _{i=-k}^{i=+k}(y_{Pi}-y(i))^{2})'_{C}=0\\(\sum _{i=-k}^{i=+k}(y_{Pi}-y(i))^{2})'_{C_{h}}=0\end{cases}}$

Soit donc :

${\begin{cases}(\sum _{i=-k}^{i=+k}sin(\omega *i)(y_{Ii}-y(i))=0\\(\sum _{i=-k}^{i=+k}sin(h*\omega *i)(y_{Ii}-y(i))=0\end{cases}}$

${\begin{cases}(\sum _{i=-k}^{i=+k}(1-cos(\omega *i))(y_{Pi}-y(i))=0\\(\sum _{i=-k}^{i=+k}(1-cos(h*\omega *i))(y_{Pi}-y(i))=0\end{cases}}$