Journal d'un Terrien

Web log de Serge Boisse

On line depuis 1992 !

Publicité

Si cette page vous a plu, Copiez son adresse et partagez-la !
http://sboisse.free.fr/science/maths/composition-fractionnaire-de-fonctions.php

composition-fractionnaire-de-fonctions
Metadata
Serge Boisse
Le 26/03/2023 à 11:03
web/MOC
oui
oui

Composition fractionnaire de fonctions mathématiques

Je vous propose ici une extension nouvelle de la composition de fonctions mathématiques (traditionnellement notée ). Si on note
f<2>(x)=f(f(x))=ff(x), et f<3>=fff ,etc., alors on va définir ici des choses comme f<1/2> ou f<2/3>, un  sujet d'étude qui, je crois, est entièrement nouveau !

On s'intéresse ici à des fonctions continues et au moins surjectives d'un domaine D sur un sous-ensemble ou une extension de D, par exemple d'un intervalle réel sur lui même. Par exemple si trois fonctions f, g, h sont (avec les bons domaines de définition, que je ne détaillerai pas ici) telles que :

g<2>=gg=f et
h<3>=hhh=f

alors on a g<2>=h<3> et on peut poser par convention 
g=h<3/2>De même si ff=g alors on peut poser par convention f=g<1/2>
Plus généralement si f<p>=g alors par définition f=g<1/p>

Par exemple il est facile de vérifier que (x4+2x2+2)<1/2>=x2+1

On a ainsi un moyen de composer "fractionnellement" des fonctions, ou de définir des puissances fractionnaires de composition. Naturellement lorsque "l'exposant" entre crochets <> est entier, on a la composition de fonction classique avec ses propriétés usuelles, telles que l'associativité et le fait que  f<1>(f(x))=f<0>= identité

notations

La notation f<1> me semble bien meilleure que le f1  des américains (sans crochets), qui est détestable à cause de la confusion possible avec 1/f. Ainsi on pourrait remplacer les discutablesarcsin et argsinh par sin<1> et sinh<1> sans risque d'erreur.

Il est également possible d'utiliser la notation alternative fn au lieu de f<n> Ainsi f(1/2)2=f

exemples

Bien que l'on ait x2<2>=(x2)2=x4, c'est le seul cas où f2=f<2>
Ainsi (x2+1)<2>=x4+2x2+2 alors que (x2+1)2=x4+2x2+1

Exemple (tiré d'une thèse de Feliks Burdecki, datant de 1926 !) :
Soit f(x)=ax2+c
il est facile de voir que f<2>(x)=a2x2+c(a+1) et que par conséquent pour p entier positif,
f<p>(x)=apx2+c/(1a)(1ap)
d'où l'on tire notamment que si 0<a<1 alors lorsque p la limite f<p>(x) vaut c/(1a)

il se trouve que cette formule pour f<p>(x) est aussi valable pour les p fractionnaires !
par exemple si g(x)=f<1/2>(x)=a1/2x2+c/(1a)(1+a1/2), alors g(g(x))=f(x)

Autre exemple : soit l'application "logistique" bien connue f(x)=4x(1x),  de l'intervalle [0,1] vers lui-même.
En posant x=sin2(y) , on a f(x)=f(sin2(y))=(2sin(y)cos(y))2=sin2(2y)
et plus généralement, pour p entier positif, f<p>(sin2(y))=sin2(2py)
mais cela ne marche pas  si p n'est pas entier...

problèmes

D'une manière générale, trouver f<1/2> connaissant f est un problème assez difficile, trouver f<1/3> ou f<1/p> est encore plus ardu !

La composition fractionnaire de fonctions pose des problèmes nouveaux. Par exemple on aimerait bien que, pour toute fonction f composable avec elle-même,  f<22/7> soit "proche" de  f<355/113> parce que 22/7 et 355/113 sont des réduites successives de π en fraction continue. Peux-t-on définir cette proximité, et cette conjecture est-elle vraie ?

Autre conjecture :
lorsque p, et pour toute fonction continue f de D sur D, limf<1/p>=identité ?

Enfin, peut-on définir des choses comme f<π> ou f<i> ? (Une composition imaginaire de fonctions !) Je serais intéressé par une expression de sin<π> ou e<i> !

Tous ces problèmes sont ouverts. Mathématiciens, à vous de jouer !

Dernière minute :

Je suis tombé sur un intéressant article de wiképédia, iterated function (en anglais). Il apparaît que l'idée d'itération fractionnaire n'est pas nouvelle (Zut, moi qui croyais l'avoir inventé...) et qu'elle est reliée à la notion mathématique de flot.

Par ailleurs plusieurs personnes ont cherché une solution f à f(f(x))=exp(x). On peut trouver un résumé de ces recherches ici.  Il apparaît que la solution n'est pas unique, mais que si l'on s'intéresse aux x réels positifs,  alors elle a une croissance plus rapide que tout polynôme tout en étant moins rapide que toute exponentielle. C'est donc une nouvelle classe de complexité (ou croissance). C'est intéressant.

Je me suis toujours demandé si l'on ne peut pas faire un parallèle entre les les classes de croissance de fonctions et les ordinaux.  Par exemple si les polynômes P sont associés à l'ordinal P(ω), et si exp est associé à l'ordinal ωω, qu'y a-t-il entre deux ? Cette question me fascine.

Exploration plus en détail

Remarque

Dans ce qui suit on utilisera la notation fp(x) au lieu de f<p>(x)

Etant donné une fonction bijective d'une partie D de ℝ ou ℂ dans elle-même, et un entier n, on notera fon(x) la composée de f par elle même.

Par exemple fo2(x)=f(f(x)), et fo1(x) est la réciproque de f, généralement notée abusivement (mais pas dans ce texte) f1(x).😒

Ne pas confondre fo2(x), la composée de f par elle-même, et f2(x), le carré de f(x). Ainsi, fo22(x) est le carré de la composée, et f2o2(x) est la composée du carré.

On a bien sûr : fon(fom(x))=fo(n+m)(x)

On va explorer ici le cas où n n'est pas entier. Dans ce qui suit r sera donc un réel.

Exemple : l'application logistique

C'est un exemple très connu. Soit la fonction de [0..1] dans [0..1] définie par :
L(x)=4x(1x)

Y axisX axisExpression 1Expression 2

En posant x=sin2(y), on obtient

L(x)=4sin2(y)cos2(y)=sin2(2y)=sin2(2arcsin(x))

et donc

Lo2(x)=L(L(x))=L(sin2(2y))=sin2(4y)=sin2(4arcsin(x))

et

Lok(x)=L(sin2(2y))=sin2(4y)=sin2(2karcsin(x))

La surprise est que cela marche aussi pour k non entier !

Ainsi

Lo(1/2)(x)=sin2(2arcsin(x))

Pasted image 20220410132450.png

Mais, du fait que Cette application a un maximum en x0=sin2(π22) ≃ 0.802..., Lok(Lok(x)) n'est pas égal à L(x) si x>x0 :
Pasted image 20220410133828.png

dérivée fractionnaire (1)

Dans cette vidéo on définit la dérivée fractionnaire comme l'inverse de l'intégrale fractionnaire, elle-même définie à partir de la relation de Cauchy pour les intégrales multiples, en considérant l'intégrale comme une transformée que l'on appellera I :
si on pose If(x)=0xf(t)dt alors

Inf(x)=1(n1)!ax(xt)n1f(t)dt

a étant un paramètre qu'on choisit (non localité) on on devrait écrire Ianf(x)=
et donc pour n réel positif,

Inf(x)=1Γ(n)ax(xt)n1f(t)dt

pour n > 0
Et donc on définira

Dnf=dndxn(Innf)

(à nouveau on devrait écrire Danf(x)=)

Ce qui permet de définir un opérateur intégro-différentiel (differintegral operator) Jnf, qui vaut Dnf si n>0, f si n=0, et Inf si n<0

On a les propriétés suivantes :

Da(xn)=Γ(n+1)γ(n+1a)xna
D1/2(1)=1πx
Da(sin(x))=sin(x+aπ2)
Da(ekx)=kaekx

attention la "chain rule" et la "product rule" ne marchent plus...

Remarque

Nous pouvons également définir simplement la dérivée fractionnaire à l'aide de la transformation de Fourier. Comme la transformation de Fourier d'une fonction dérivée n fois est mise à l'échelle par la puissance n de la fréquence, nous pouvons remplacer n par une valeur réelle et utiliser la transformation de Fourier inverse pour obtenir le résultat.

Dérivées fractionnelles (2)

On peut être tenté de tracer le graphe des itérés successifs d'un réél x0 de départ, c'est à dire tracer le graphe d'une fonction g(x0,x) telle que g(xo,x)=fox(x0)

Cela conduit à se demander s'il peut exister une fonction f telle que , pour un x0 donné, et pour tout ϵ positif, f(x0+ϵ)=fo(ϵ+1)(x0)

Par exemple si ϵ=1, on a f(x0+1)=f(f(x0)), et bien sûr si ϵ=0, on retrouve f(x0)=f(x0)

Supposons qu'il existe une telle fonction telle que cela soit vrai pour x0 quelconque, alors

Publicité
Commentaires

Commentaires (0) :

Page :



Ajouter un commentaire (pas besoin de s'enregistrer)

Pseudo :
Message :


image de protection
En cliquant sur le bouton "Envoyer" vous acceptez les conditions suivantes : Ne pas poster de message injurieux, obscène ou contraire à la loi, ni de liens vers de tels sites. Respecter la "netiquette", ne pas usurper le pseudo d'une autre personne, respecter les posts faits par les autres. L'auteur du site se réserve le droit de supprimer un ou plusieurs posts à tout moment. Merci !
Ah oui : le bbcode et le html genre <br>, <a href=...>, <b>b etc. ne fonctionnent pas dans les commentaires. C'est voulu.
< Retour en haut de la page