Une page de Wikiversité.
L'équation de la courbe de l'harmonique peut se ramener à la forme :
y
=
y
0
+
S
∗
s
i
n
(
ω
∗
t
)
+
C
∗
(
1
−
c
o
s
(
ω
∗
t
)
)
{\displaystyle y=y_{0}+S*sin(\omega *t)+C*(1-cos(\omega *t))}
Avec les propriétés suivantes :
- La somme des carrés des différences des
y
i
{\displaystyle y_{i}}
aux
y
(
i
)
{\displaystyle y(i)}
est minimale pour un
ω
{\displaystyle \omega }
donné ( propriété de base de la régression )
Propriété P1 :
∑
i
=
−
k
i
=
+
k
(
y
i
−
y
(
i
)
)
2
{\displaystyle \sum _{i=-k}^{i=+k}(y_{i}-y(i))^{2}}
minimale
Parmi toutes ces harmoniques obtenues pour tous les
ω
{\displaystyle \omega }
, une pour chaque, certaines sont plus probables que d'autres. Exemples :
_ celles où
y
m
a
x
−
y
m
i
n
{\displaystyle y_{max}-y_{min}}
est minimale ou inférieure à une valeur donnée
ϵ
{\displaystyle \epsilon }
. (Condition C1)
_ celles où la valeur absolue de la différence entre
y
m
a
x
−
y
m
i
n
{\displaystyle y_{max}-y_{min}}
et {la moyenne des k plus grands y – la moyenne de k plus petits y} hors
y
0
{\displaystyle y_{0}}
est minimale ou inférieure à une valeur donnée
ϵ
{\displaystyle \epsilon }
(Condition C2).
_ celles où la puissance de y minimale :
(
C
2
+
S
2
)
ω
{\displaystyle (C^{2}+S^{2})\omega }
est minimale ou inférieure à une valeur donnée
ϵ
{\displaystyle \epsilon }
(Condition C3).
_ ou d'autres conditions. On quitte le domaines des maths pures pour les maths appliquées.
P1 : Se traduit par le système :
(
∑
i
=
−
k
i
=
+
k
(
y
i
−
y
(
i
)
)
2
)
S
′
=
0
{\displaystyle (\sum _{i=-k}^{i=+k}(y_{i}-y(i))^{2})'_{S}=0}
et
(
∑
i
=
−
k
i
=
+
k
(
y
i
−
y
(
i
)
)
2
)
C
′
=
0
{\displaystyle (\sum _{i=-k}^{i=+k}(y_{i}-y(i))^{2})'_{C}=0}
{
(
∑
i
=
−
k
i
=
+
k
(
y
i
−
y
(
i
)
)
2
)
S
′
=
0
(
∑
i
=
−
k
i
=
+
k
(
y
i
−
y
(
i
)
)
2
)
C
′
=
0
{\displaystyle {\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}}}
{
(
∑
i
=
−
k
i
=
+
k
(
y
(
i
)
−
y
0
−
S
∗
s
i
n
(
ω
i
)
−
C
∗
(
1
−
c
o
s
(
ω
i
)
)
)
2
)
s
′
=
0
(
∑
i
=
−
k
i
=
+
k
(
y
(
i
)
−
y
0
−
S
∗
s
i
n
(
ω
i
)
−
C
∗
(
1
−
c
o
s
(
ω
i
)
)
)
2
)
c
′
=
0
{\displaystyle {\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}}}
{
∑
i
=
−
k
i
=
+
k
s
i
n
(
ω
i
)
(
y
(
i
)
−
y
0
−
S
∗
s
i
n
(
ω
i
)
−
C
∗
(
1
−
c
o
s
(
ω
i
)
)
=
0
(
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
ω
i
)
)
(
y
(
i
)
−
y
0
−
S
∗
s
i
n
(
ω
i
)
−
C
∗
(
1
−
c
o
s
(
ω
i
)
)
=
0
{\displaystyle {\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}}}
{
∑
i
=
−
k
i
=
+
k
s
i
n
(
ω
i
)
(
y
(
i
)
−
y
0
)
=
S
s
i
n
(
ω
i
)
2
+
C
s
i
n
(
ω
i
)
(
1
−
c
o
s
(
ω
i
)
)
(
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
ω
i
)
)
(
y
(
i
)
−
y
0
)
=
S
s
i
n
(
ω
i
)
(
1
−
c
o
s
(
ω
i
)
)
+
C
(
1
−
c
o
s
(
ω
i
)
)
2
{\displaystyle {\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 :
{
∑
i
=
−
k
i
=
+
k
s
i
n
(
ω
i
)
y
(
i
)
=
S
s
i
n
(
ω
i
)
2
(
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
ω
i
)
)
y
(
i
)
=
C
(
1
−
c
o
s
(
ω
i
)
)
2
{\displaystyle {\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
=
∑
i
=
−
k
i
=
+
k
y
(
i
)
s
i
n
(
ω
i
)
∑
i
=
−
k
i
=
+
k
s
i
n
2
(
ω
i
)
{\displaystyle S={\frac {\sum _{i=-k}^{i=+k}y(i)sin(\omega i)}{\sum _{i=-k}^{i=+k}sin^{2}(\omega i)}}}
et
C
=
∑
i
=
−
k
i
=
+
k
y
(
i
)
(
1
−
c
o
s
(
ω
i
)
)
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
ω
i
)
)
2
{\displaystyle 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 ( phi déphasage ) et des maxi-mini [ modifier | modifier le wikicode ]
y
=
y
0
+
S
∗
s
i
n
(
ω
∗
t
)
+
C
∗
(
1
−
c
o
s
(
ω
∗
t
)
)
{\displaystyle y=y_{0}+S*sin(\omega *t)+C*(1-cos(\omega *t))}
y
′
(
t
)
=
ω
S
c
o
s
(
ω
∗
t
)
+
ω
C
s
i
n
(
ω
∗
t
)
)
{\displaystyle y'(t)=\omega Scos(\omega *t)+\omega Csin(\omega *t))}
Pour
t
M
{\displaystyle t_{M}}
et
t
m
{\displaystyle t_{m}}
on obtient les extrèmum et
y
t
′
(
t
M
;
t
m
)
=
0
{\displaystyle y'_{t}(t_{M};t_{m})=0}
D'où
ω
∗
t
m
,
m
{\displaystyle \omega *t_{m,m}}
par
t
a
n
(
ω
∗
t
M
,
m
)
=
−
S
C
{\displaystyle tan(\omega *t_{M,m})=-{\frac {S}{C}}}
ω
t
M
=
a
t
a
n
(
−
S
C
)
{\displaystyle \omega t_{M}=atan(-{\frac {S}{C}})}
ω
t
m
=
a
t
a
n
(
−
S
C
)
+
π
{\displaystyle \omega t_{m}=atan(-{\frac {S}{C}})+\pi }
Aux extrémum, en substituant
ω
∗
t
M
,
m
{\displaystyle \omega *t_{M,m}}
, on obtient
y
M
,
y
m
{\displaystyle y_{M},y_{m}}
soit :
y
(
w
)
=
y
0
+
S
∗
s
i
n
(
a
t
a
n
(
−
S
C
)
)
+
C
∗
(
1
−
c
o
s
(
a
t
a
n
(
−
S
C
)
)
{\displaystyle 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
∗
s
i
n
(
a
t
a
n
(
−
S
C
)
)
−
C
c
o
s
(
a
t
a
n
(
−
S
C
)
)
{\displaystyle y(w)=S*sin(atan(-{\frac {S}{C}}))-Ccos(atan(-{\frac {S}{C}}))}
Un programme ou le tableur donneront la valeur de w optimale cherchée
ω
o
p
t
{\displaystyle \omega _{opt}}
telle que :
y
m
i
n
(
w
)
=
S
(
ω
)
∗
s
i
n
(
a
t
a
n
(
−
S
(
ω
)
C
(
ω
)
)
)
−
C
(
ω
)
c
o
s
(
a
t
a
n
(
−
S
(
ω
)
C
(
ω
)
)
)
{\displaystyle y_{min}(w)=S(\omega )*sin(atan(-{\frac {S(\omega )}{C(\omega )}}))-C(\omega )cos(atan(-{\frac {S(\omega )}{C(\omega )}}))}
y
(
w
)
=
S
(
ω
o
p
t
)
s
i
n
(
ω
o
p
t
∗
t
)
+
C
(
ω
o
p
t
)
c
o
s
(
ω
o
p
t
∗
t
)
)
{\displaystyle y(w)=S(\omega _{opt})sin(\omega _{opt}*t)+C(\omega _{opt})cos(\omega _{opt}*t))}
Les harmoniques simples s'écrivent alors :
y
=
y
0
−
C
+
S
2
+
C
2
s
i
n
(
w
t
+
a
t
a
n
(
−
S
C
)
{\displaystyle y=y_{0}-C+{\sqrt {S^{2}+C^{2}}}sin(wt+atan(-{\frac {S}{C}})}
Parmi celles-ci, celle qui a une amplitude minimale ; pour
ω
m
i
n
{\displaystyle \omega _{min}}
on a
S
(
ω
m
i
n
)
2
+
C
(
ω
m
i
n
)
2
{\displaystyle {\sqrt {S(\omega _{min})^{2}+C(\omega _{min})^{2}}}}
minimal.
Parmi aussi, celle de
S
2
+
C
2
e
f
f
{\displaystyle {\sqrt {S^{2}+C^{2}}}_{eff}}
efficace qui donnera
w
e
f
f
{\displaystyle w_{eff}}
en faisant
S
2
+
C
2
e
f
f
=
S
2
+
C
2
{\displaystyle {\sqrt {S^{2}+C^{2}}}_{eff}={\sqrt {S^{2}+C^{2}}}}
.A noter que la valeur efficace est la valeur moyenne au sens Riemanien.
S
2
+
C
2
e
f
f
=
S
(
ω
e
f
f
)
2
+
C
(
ω
e
f
f
)
2
=
1
2
π
∫
ω
=
1
ω
=
2
π
S
(
ω
)
2
+
C
(
ω
)
2
d
ω
{\displaystyle {\sqrt {S^{2}+C^{2}}}_{eff}={\sqrt {S(\omega _{eff})^{2}+C(\omega _{eff})^{2}}}={\frac {1}{2\pi }}\int _{\omega =1}^{\omega =2\pi }{\sqrt {S(\omega )^{2}+C(\omega )^{2}}}d\omega }
y
(
i
)
=
y
I
(
i
)
+
y
P
(
i
)
{\displaystyle y(i)=y_{I}(i)+y_{P}(i)}
avec
y
I
=
0
+
S
s
i
n
(
ω
t
)
{\displaystyle y_{I}=0+Ssin(\omega t)}
et
y
P
=
y
0
+
C
(
1
−
c
o
s
(
ω
t
)
)
{\displaystyle y_{P}=y_{0}+C(1-cos(\omega t))}
avec
y
I
(
i
)
=
y
(
i
)
−
y
(
−
i
)
2
{\displaystyle y_{I}(i)={\frac {y(i)-y(-i)}{2}}}
et
y
P
(
i
)
=
y
(
i
)
+
y
(
−
i
)
2
{\displaystyle 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 :
{
(
∑
i
=
−
k
i
=
+
k
(
y
I
i
−
y
(
i
)
)
2
)
S
′
=
0
(
∑
i
=
−
k
i
=
+
k
(
y
P
i
−
y
(
i
)
)
2
)
C
′
=
0
{\displaystyle {\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}}}
{
∑
i
=
−
k
i
=
+
k
s
i
n
(
ω
i
)
(
y
I
i
−
y
(
i
)
)
=
0
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
ω
i
)
)
(
y
P
i
−
y
(
i
)
)
=
0
{\displaystyle {\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}}}
{
∑
i
=
−
k
i
=
+
k
s
i
n
(
ω
i
)
(
y
I
i
−
S
s
i
n
(
ω
t
)
)
=
0
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
ω
t
)
)
(
y
P
i
−
y
0
+
C
(
1
−
c
o
s
(
ω
t
)
)
=
0
{\displaystyle {\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
=
∑
i
=
−
k
i
=
+
k
y
I
i
s
i
n
(
ω
i
)
∑
i
=
−
k
i
=
+
k
s
i
n
2
(
ω
i
)
{\displaystyle S={\frac {\sum _{i=-k}^{i=+k}y_{Ii}sin(\omega i)}{\sum _{i=-k}^{i=+k}sin^{2}(\omega i)}}}
et
C
=
∑
i
=
−
k
i
=
+
k
y
P
i
(
1
−
c
o
s
(
ω
i
)
)
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
ω
i
)
)
2
{\displaystyle C={\frac {\sum _{i=-k}^{i=+k}y_{Pi}(1-cos(\omega i))}{\sum _{i=-k}^{i=+k}(1-cos(\omega i))^{2}}}}
Or :
∑
i
=
−
k
i
=
+
k
y
P
i
s
i
n
(
ω
i
)
∑
i
=
−
k
i
=
+
k
s
i
n
2
(
ω
i
)
=
0
{\displaystyle {\frac {\sum _{i=-k}^{i=+k}y_{Pi}sin(\omega i)}{\sum _{i=-k}^{i=+k}sin^{2}(\omega i)}}=0}
Et :
∑
i
=
−
k
i
=
+
k
y
I
i
(
1
−
c
o
s
(
ω
i
)
)
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
ω
i
)
)
2
=
0
{\displaystyle {\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
=
∑
i
=
−
k
i
=
+
k
y
(
i
)
s
i
n
(
ω
i
)
∑
i
=
−
k
i
=
+
k
s
i
n
2
(
ω
i
)
{\displaystyle S={\frac {\sum _{i=-k}^{i=+k}y(i)sin(\omega i)}{\sum _{i=-k}^{i=+k}sin^{2}(\omega i)}}}
et
C
=
∑
i
=
−
k
i
=
+
k
y
(
i
)
(
1
−
c
o
s
(
ω
i
)
)
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
ω
i
)
)
2
{\displaystyle 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
ω
{\displaystyle \omega }
ω
t
M
=
a
t
a
n
(
−
S
C
)
{\displaystyle \omega t_{M}=atan(-{\frac {S}{C}})}
ω
t
m
=
a
t
a
n
(
−
S
C
)
+
π
{\displaystyle \omega t_{m}=atan(-{\frac {S}{C}})+\pi }
Avec la condition C1 plus haut , il s'agira de minimiser la quantité :
y
(
w
)
=
S
∗
s
i
n
(
a
t
a
n
(
−
S
C
)
)
+
C
∗
(
1
−
c
o
s
(
a
t
a
n
(
−
S
C
)
)
)
{\displaystyle 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
ω
o
p
t
{\displaystyle \omega _{opt}}
.
y
=
y
0
+
S
(
ω
o
p
t
)
s
i
n
(
ω
o
p
t
∗
t
)
+
C
(
ω
o
p
t
)
∗
(
1
−
c
o
s
(
ω
o
p
t
∗
t
)
)
{\displaystyle 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
{\displaystyle t=0,5+k}
et
t
=
0
,
5
−
k
{\displaystyle t=0,5-k}
sans
y
0
{\displaystyle y_{0}}
couples
(
x
−
2
k
−
1
,
y
+
2
k
+
1
)
{\displaystyle (x_{-2k-1},y_{+2k+1})}
D'où S et C :
S
=
∑
i
=
−
2
k
−
1
i
=
+
2
k
+
1
y
I
i
s
i
n
(
ω
i
)
∑
i
=
−
2
k
−
1
i
=
+
2
k
+
1
s
i
n
2
(
ω
i
)
{\displaystyle 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
=
∑
i
=
−
2
k
+
1
i
=
+
2
k
+
1
y
P
i
(
1
−
c
o
s
(
ω
i
)
)
∑
i
=
−
2
k
−
1
i
=
+
2
k
+
1
(
1
−
c
o
s
(
ω
i
)
)
2
{\displaystyle 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}}}}
ω
t
M
=
a
t
a
n
(
−
S
C
)
{\displaystyle \omega t_{M}=atan(-{\frac {S}{C}})}
ω
t
m
=
a
t
a
n
(
−
S
C
)
+
π
{\displaystyle \omega t_{m}=atan(-{\frac {S}{C}})+\pi }
Aux extrémum, en substituant
ω
∗
t
M
,
m
{\displaystyle \omega *t_{M,m}}
, on obtient
y
M
,
y
m
{\displaystyle y_{M},y_{m}}
soit :
y
(
w
)
=
S
∗
s
i
n
(
a
t
a
n
(
−
S
C
)
)
+
C
∗
(
1
−
c
o
s
(
a
t
a
n
(
−
S
C
)
)
)
{\displaystyle y(w)=S*sin(atan(-{\frac {S}{C}}))+C*(1-cos(atan(-{\frac {S}{C}})))}
à minimiser pour trouver
w
o
p
t
{\displaystyle w_{o}pt}
y
(
2
t
+
0
,
5
)
=
S
(
ω
o
p
t
)
s
i
n
(
ω
o
p
t
∗
(
2
t
+
0
,
5
)
)
+
C
(
ω
o
p
t
)
∗
(
1
−
c
o
s
(
ω
o
p
t
∗
(
2
t
+
0
,
5
)
)
)
{\displaystyle 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.
Soit à décomposer un échantillonnage y* en une somme sin et cos de périodes n de 3 à N.
Comme précédemment, pour 2k+1 couples (i,y) avec i au pas de 1, on écrit y sous la forme :
y
∗
(
i
)
=
y
0
+
∑
n
=
3
n
=
N
S
n
∗
s
i
n
(
2
π
n
∗
i
)
+
∑
n
=
3
n
=
N
C
n
(
1
−
c
o
s
(
2
π
n
∗
i
)
)
{\displaystyle y^{*}(i)=y_{0}+\sum _{n=3}^{n=N}S_{n}*sin({\frac {2\pi }{n}}*i)+\sum _{n=3}^{n=N}C_{n}(1-cos({\frac {2\pi }{n}}*i))}
Procédure :
partager y* en partie paire ( cos et terme constant ) et impaire ( sin ) ;
écrire les conditions de régression pour chaque partie, comme quoi les sommes des carrés des différences sont minimales et que donc les dérivées par rapport à chaque S et chaque C sont nulles. Cela permet d'obtenir deux systèmes de n équations linéaires, un en S et un en C, résolvables en C et S.
On obtient:
S
n
=
|
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
1
∗
i
)
s
i
n
(
2
π
1
∗
i
)
⋯
∑
i
=
−
k
i
=
+
k
y
i
∗
s
i
n
(
2
π
∗
i
1
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
1
∗
i
)
s
i
n
(
2
π
N
∗
i
)
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
1
∗
i
)
s
i
n
(
2
π
2
∗
i
)
⋯
∑
i
=
−
k
i
=
+
k
y
i
∗
s
i
n
(
2
π
∗
i
2
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
2
∗
i
)
s
i
n
(
2
π
N
∗
i
)
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
1
∗
i
)
s
i
n
(
2
π
n
∗
i
)
⋯
∑
i
=
−
k
i
=
+
k
y
i
∗
s
i
n
(
2
π
∗
i
n
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
n
∗
i
)
s
i
n
(
2
π
N
∗
i
)
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
1
∗
i
)
s
i
n
(
2
π
N
∗
i
)
⋯
∑
i
=
−
k
i
=
+
k
y
i
∗
s
i
n
(
2
π
∗
i
N
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
N
∗
i
)
s
i
n
(
2
π
N
∗
i
)
|
|
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
1
∗
i
)
s
i
n
(
2
π
1
∗
i
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
∗
i
1
)
s
i
n
(
2
π
∗
i
1
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
1
∗
i
)
s
i
n
(
2
π
N
∗
i
)
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
1
∗
i
)
s
i
n
(
2
π
2
∗
i
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
∗
i
2
)
s
i
n
(
2
π
∗
i
2
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
2
∗
i
)
s
i
n
(
2
π
N
∗
i
)
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
1
∗
i
)
s
i
n
(
2
π
n
∗
i
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
∗
i
n
)
s
i
n
(
f
r
a
c
2
π
∗
i
n
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
n
∗
i
)
s
i
n
(
2
π
N
∗
i
)
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
1
∗
i
)
s
i
n
(
2
π
N
∗
i
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
∗
i
N
)
s
i
n
(
2
π
∗
i
N
)
⋯
∑
i
=
−
k
i
=
+
k
s
i
n
(
2
π
N
∗
i
)
s
i
n
(
2
π
N
∗
i
)
|
{\displaystyle S_{n}={\frac {\begin{vmatrix}\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{1}}*i)sin({\frac {2\pi }{1}}*i)&{\cdots }&\sum _{i=-k}^{i=+k}y_{i}^{*}sin({\frac {2\pi *i}{1}})&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{1}}*i)sin({\frac {2\pi }{N}}*i)\\\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{1}}*i)sin({\frac {2\pi }{2}}*i)&{\cdots }&\sum _{i=-k}^{i=+k}y_{i}^{*}sin({\frac {2\pi *i}{2}})&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{2}}*i)sin({\frac {2\pi }{N}}*i)\\{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }\\\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{1}}*i)sin({\frac {2\pi }{n}}*i)&{\cdots }&\sum _{i=-k}^{i=+k}y_{i}^{*}sin({\frac {2\pi *i}{n}})&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{n}}*i)sin({\frac {2\pi }{N}}*i)\\{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }\\\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{1}}*i)sin({\frac {2\pi }{N}}*i)&{\cdots }&\sum _{i=-k}^{i=+k}y_{i}^{*}sin({\frac {2\pi *i}{N}})&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{N}}*i)sin({\frac {2\pi }{N}}*i)\end{vmatrix}}{\begin{vmatrix}\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{1}}*i)sin({\frac {2\pi }{1}}*i)&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi *i}{1}})sin({\frac {2\pi *i}{1}})&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{1}}*i)sin({\frac {2\pi }{N}}*i)\\\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{1}}*i)sin({\frac {2\pi }{2}}*i)&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi *i}{2}})sin({\frac {2\pi *i}{2}})&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{2}}*i)sin({\frac {2\pi }{N}}*i)\\{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }\\\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{1}}*i)sin({\frac {2\pi }{n}}*i)&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi *i}{n}})sin(frac{2\pi *i}{n})&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{n}}*i)sin({\frac {2\pi }{N}}*i)\\{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }\\\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{1}}*i)sin({\frac {2\pi }{N}}*i)&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi *i}{N}})sin({\frac {2\pi *i}{N}})&{\cdots }&\sum _{i=-k}^{i=+k}sin({\frac {2\pi }{N}}*i)sin({\frac {2\pi }{N}}*i)\\\end{vmatrix}}}}
C
n
==
|
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
⋯
∑
i
=
−
k
i
=
+
k
y
i
∗
(
1
−
c
o
s
(
2
π
∗
i
1
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
(
1
−
c
o
s
(
2
π
2
∗
i
)
)
⋯
∑
i
=
−
k
i
=
+
k
y
i
∗
(
1
−
c
o
s
(
2
π
∗
i
2
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
2
∗
i
)
)
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
(
1
−
c
o
s
(
2
π
n
∗
i
)
)
⋯
∑
i
=
−
k
i
=
+
k
y
i
∗
(
1
−
c
o
s
(
2
π
∗
i
n
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
n
∗
i
)
)
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
⋯
∑
i
=
−
k
i
=
+
k
y
i
∗
(
1
−
c
o
s
(
2
π
∗
i
N
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
|
|
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
∗
i
1
)
)
(
1
−
c
o
s
(
2
π
∗
i
1
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
(
1
−
c
o
s
(
2
π
2
∗
i
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
∗
i
2
)
)
(
1
−
c
o
s
(
2
π
∗
i
2
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
2
∗
i
)
)
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
(
1
−
c
o
s
(
2
π
n
∗
i
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
∗
i
n
)
)
(
1
−
c
o
s
(
2
π
∗
i
n
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
n
∗
i
)
)
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
1
∗
i
)
)
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
∗
i
N
)
)
(
1
−
c
o
s
(
2
π
∗
i
N
)
)
⋯
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
(
1
−
c
o
s
(
2
π
N
∗
i
)
)
|
{\displaystyle C_{n}=={\frac {\begin{vmatrix}\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{1}}*i))(1-cos({\frac {2\pi }{1}}*i))&{\cdots }&\sum _{i=-k}^{i=+k}y_{i}^{*}(1-cos({\frac {2\pi *i}{1}}))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{1}}*i))(1-cos({\frac {2\pi }{N}}*i))\\\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{1}}*i))(1-cos({\frac {2\pi }{2}}*i))&{\cdots }&\sum _{i=-k}^{i=+k}y_{i}^{*}(1-cos({\frac {2\pi *i}{2}}))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{2}}*i))(1-cos({\frac {2\pi }{N}}*i))\\{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }\\\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{1}}*i))(1-cos({\frac {2\pi }{n}}*i))&{\cdots }&\sum _{i=-k}^{i=+k}y_{i}^{*}(1-cos({\frac {2\pi *i}{n}}))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{n}}*i))(1-cos({\frac {2\pi }{N}}*i))\\{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }\\\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{1}}*i))(1-cos({\frac {2\pi }{N}}*i))&{\cdots }&\sum _{i=-k}^{i=+k}y_{i}^{*}(1-cos({\frac {2\pi *i}{N}}))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{N}}*i))(1-cos({\frac {2\pi }{N}}*i))\end{vmatrix}}{\begin{vmatrix}\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{1}}*i))(1-cos({\frac {2\pi }{1}}*i))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi *i}{1}}))(1-cos({\frac {2\pi *i}{1}}))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{1}}*i))(1-cos({\frac {2\pi }{N}}*i))\\\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{1}}*i))(1-cos({\frac {2\pi }{2}}*i))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi *i}{2}}))(1-cos({\frac {2\pi *i}{2}}))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{2}}*i))(1-cos({\frac {2\pi }{N}}*i))\\{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }\\\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{1}}*i))(1-cos({\frac {2\pi }{n}}*i))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi *i}{n}}))(1-cos({\frac {2\pi *i}{n}}))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{n}}*i))(1-cos({\frac {2\pi }{N}}*i))\\{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }&{\cdots }{\cdots }\\\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{1}}*i))(1-cos({\frac {2\pi }{N}}*i))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi *i}{N}}))(1-cos({\frac {2\pi *i}{N}}))&{\cdots }&\sum _{i=-k}^{i=+k}(1-cos({\frac {2\pi }{N}}*i))(1-cos({\frac {2\pi }{N}}*i))\\\end{vmatrix}}}}
___________
Alors que la Transformée de Fourier décompose uniquement en sin et cos de périodes < 2k, intervalle de définition, cette méthode travaille sur les périodes inférieures ET supérieures à l'intervalle de définition 2k ! On peut donc commencer à envisager des extrapolations possibles.
Il faut balayer toutes les valeurs pour trouver le spectre. Une méthode déjà présentée dans un Laboratoire de Patrick Bréjon, désormais disparu permet d'accéder directement à ce spectre selon les valeurs décroissantes de n , croissantes de omega
Autre approche : harmonique double de pulsations 2pi/n et 2pi/n+k puis 2pi/n et k2pi/n [ modifier | modifier le wikicode ]
Suite à revoir afin de respecter les données du titre
y
=
y
0
+
S
∗
s
i
n
(
ω
∗
t
)
+
S
h
∗
s
i
n
(
ω
∗
h
t
)
+
C
∗
(
1
−
c
o
s
(
ω
∗
t
)
)
+
C
h
∗
(
1
−
c
o
s
(
ω
∗
h
t
)
)
{\displaystyle 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
)
{\displaystyle y(i)=y_{I}(i)+y_{P}(i)}
avec
y
I
=
0
+
S
s
i
n
(
ω
t
)
+
S
h
s
i
n
(
h
ω
t
)
{\displaystyle y_{I}=0+Ssin(\omega t)+S_{h}sin(h\omega t)}
et
y
P
=
y
0
+
C
(
1
−
c
o
s
(
ω
t
)
)
+
C
(
1
−
c
o
s
(
ω
t
)
)
{\displaystyle y_{P}=y_{0}+C(1-cos(\omega t))+C(1-cos(\omega t))}
avec
y
I
(
i
)
=
y
(
i
)
−
y
(
−
i
)
2
{\displaystyle y_{I}(i)={\frac {y(i)-y(-i)}{2}}}
et
y
P
(
i
)
=
y
(
i
)
+
y
(
−
i
)
2
{\displaystyle 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 :
{
(
∑
i
=
−
k
i
=
+
k
(
y
I
i
−
y
(
i
)
)
2
)
S
′
=
0
(
∑
i
=
−
k
i
=
+
k
(
y
I
i
−
y
(
i
)
)
2
)
S
h
′
=
0
{\displaystyle {\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}}}
{
(
∑
i
=
−
k
i
=
+
k
(
y
P
i
−
y
(
i
)
)
2
)
C
′
=
0
(
∑
i
=
−
k
i
=
+
k
(
y
P
i
−
y
(
i
)
)
2
)
C
h
′
=
0
{\displaystyle {\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 :
{
(
∑
i
=
−
k
i
=
+
k
s
i
n
(
ω
∗
i
)
(
y
I
i
−
y
(
i
)
)
=
0
(
∑
i
=
−
k
i
=
+
k
s
i
n
(
h
∗
ω
∗
i
)
(
y
I
i
−
y
(
i
)
)
=
0
{\displaystyle {\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}}}
{
(
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
ω
∗
i
)
)
(
y
P
i
−
y
(
i
)
)
=
0
(
∑
i
=
−
k
i
=
+
k
(
1
−
c
o
s
(
h
∗
ω
∗
i
)
)
(
y
P
i
−
y
(
i
)
)
=
0
{\displaystyle {\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}}}