Bibm@th

Forum de mathématiques - Bibm@th.net

Bienvenue dans les forums du site BibM@th, des forums où on dit Bonjour (Bonsoir), Merci, S'il vous plaît...

Vous n'êtes pas identifié(e).

#1 Entraide (supérieur) » Expansion asymptotique (Correlation canonique) » 12-09-2024 20:59:57

pentium mix
Réponses : 0

Bonsoir tout le monde.
S'il vous plait je suis entrain de lire un article sur la corrélation canonique sauf qu'il y'a un paragraphe qui m'embête beaucoup. Je n'y comprend absolument rien.


Le paragraphe:


Soit $\tilde{r}_1^2 \ge ... \ge \tilde{r}_p^2$ les racines de l'équation $|\tilde{M} - \tilde{r}^2_i S_{11}^{*-1}| = 0$ où $\tilde{M} = S_{12} S_{22,n}^{-1}S_{21}$,
et
$$
U =(u_{ij}) = \sqrt{n}(S - \Sigma) = \begin{pmatrix} V & Z \\ Z' & W \end{pmatrix}.
$$
Dans le contexte ci-dessus, nous fournissons le développement asymptotique de la distribution de $\tilde{y}_i = \sqrt{n}(\tilde{r}_i^2 - \rho_i^2)$. Nous exprimons le
développement asymptotique de $\tilde{M}$ comme suit :
\begin{equation*}
    \tilde{M} = \Omega + n^{-1/2}\tilde{M}^{(1)} + n^{-1}\tilde{M}^{(2)} + O_p(n^{-3/2}).
\end{equation*}
Ici, $\Omega = diag(\rho_1^2,...,\rho_p^2)$ et $\tilde{M}^{(k)} = \begin{pmatrix} \tilde{m}_{ij}^{(k)} \end{pmatrix}$ pour $k=1,2$.
Où on considère un échantillon
\begin{equation}\label{complete_data}
    \begin{pmatrix}
        x_{1,1} \\
        x_{2,1}
    \end{pmatrix},
    \begin{pmatrix}
        x_{1,2} \\
        x_{2,2}
    \end{pmatrix}, \ldots,
    \begin{pmatrix}
        x_{1,n} \\
        x_{2,n}
    \end{pmatrix},
\end{equation}
où $x_k = (x'_{1,k}, x'_{2,k})'$, $k = 1, \ldots, n$, sont des observations (i.i.d.) indépendantes et identiquement distribuées de $x$. Désignons par $\bar{x}$ et $S$ la moyenne de l'échantillon et la matrice de covariance correspondantes, exprimées comme suit :
\[
\bar{x} = (\bar{x}'_1, \bar{x}'_{2,n})' = \frac{1}{n} \sum_{j=1}^n x_j,
\]
\[
S = \begin{pmatrix}
    S_{11} & S_{12} \\
    S_{21} & S_{22,n}
\end{pmatrix} = \frac{1}{n} \sum_{j=1}^n (x_j - \bar{x})(x_j - \bar{x})'.
\]




Je ne comprend pas comment on fait pour ecrire l'expansion asymptotique de la matrice M

Merci d'avance

#2 Entraide (supérieur) » Convexité et Optimisation » 06-05-2024 12:31:13

pentium mix
Réponses : 0

Bonjour tout le monde .
S'il vous plaît je souhaite écrire l'algorithme du gradient conjugué avec animation video sur python.

J'ai commencé avec l'algorithme a pas fixe et  mon  problème c'est que je n'arrive pas a faire l'animation vidéo sur python.

Puis-je avoir de l'aide s'il vous plaît.

Merci

#3 Entraide (supérieur) » Factorisations LU » 01-01-2024 23:54:24

pentium mix
Réponses : 1

Bonne année à tous.
S'il vous plaît j'ai un problème sur la factorisation LU des matrices par blocs. Par exemple sur l'exercice https://www.cjoint.com/c/NAbvZGYihAy je ne sais comment montrer l'unicité de la factorisation pour la question 3 et pour le 4 je suis complètement bloqué.

Merci

#4 Re : Programmation » generate a pascal triangle well » 18-12-2023 05:47:23

Bonjour
Merci beaucoup
Voici une idée qui m'est venu en tète actuellement: utiliser des dataframes,

 import numpy as np
import pandas as pd
def TriangleDePascal(n):
    T=np.eye(n, dtype=int)
    for i in range(n):
        T[i,0]=1
    for i in range(2,n):
        for j in range (1,i):
            T[i,j]=T[i-1,j-1]+T[i-1,j]
    return T
     
           
n=int(input("n="))
df = pd.DataFrame(T)
df[df == 0] = "" # remplace toutes les valeurs 0 par du vide
print(df)
Fred a écrit :

Bonsoir,

  Je n'ai pas compris si tu voulais ou si tu ne voulais pas utiliser de liste, mais en tout cas, en partant d'un tableau, sauf à changer la façon d'afficher à la fin, tu auras toujours des zéros. Voici une solution qui fonctionne avec des listes.


def TriangleDePascal(n):
    T=[]
    for i in range(n):
        T.append([1])
    for i in range(2,n):
        for j in range (1,i-1):
            T[i].append(T[i-1][j-1]+T[i-1][j])
        T[i].append(1)
    for i in range(n):
        print(*T[i],'   ')
 

#5 Programmation » generate a pascal triangle well » 17-12-2023 21:36:06

pentium mix
Réponses : 2

bonsoir tout le monde, je souhaite générer le triangle de pascal tel qu'il est dans les livres mais ne sais comment faire pour ne plus afficher les 0 et puis je ne veux pas aussi travailler avec des listes

mon programme


import numpy as np
def TriangleDePascal(n):
    T=np.eye(n, dtype=int) ## matrice identité
    for i in range(n):
        T[i,0]=1        # remplir la première colonne
    for i in range(2,n):
        for j in range (1,i):
            T[i,j]=T[i-1,j-1]+T[i-1,j]
    for elt in T:
          print(*elt, sep="   ")

#6 Re : Programmation » Afficher le debut de la plus longue suite consecutive de zeros » 14-12-2023 05:38:04

Infiniment merci
Au final j'ai l'impression de n'avoir pas assez réfléchi pour cet exercice

Glozi a écrit :

Bonsoir,
Pour comprendre un algorithme, le plus concret c'est de prendre un exemple de tableau T, et de suivre avec crayon et papier les différentes étapes de la boucle. Ça permet de comprendre ce qui se passe.

avec une boucle for

Un algo pas du tout optimisé mais j'espère assez limpide.


def max_zero_sequence(T):
    """
    Entrée :
        T - une liste d'entiers/booléens
    Sortie :
        index_max, length_max
        où - length_max représente la longueur de la plus longue séquence de zéros consécutifs dans T (length_max=0 s'il n'y a aucun zéro dans T)
        - index_max représente l'indice de départ de cette séquence maximale (index_max=None s'il n'y a aucun zéro dans T)
    """
    index_max=None #contient l'indice du début de la plus longue séquence de zéros consécutifs découverte
    length_max=0 #la longueur de cette plus longue séquence découverte
   
    index_start_current = None #contient l'indice du début de la séquence de zéros consécutifs que nous sommes en train d'explorer
    length_current=0 #la longueur de cette séquence actuelle
   
    currently_in_zero_sequence = False #vaut True si et seulement si l'indice i précédent se situe dans une séquence de zéros.
   
    for i in range(len(T)): #boucle d'exploration
        if T[i]==0: ## si l'entrée qu'on regarde est un zéro alors :
            if currently_in_zero_sequence: #si on a déjà trouvé des zéros avant cette entrée
                length_current +=1 #alors on augmente juste la longueur de la séquence actuelle.
            else: #sinon c'est que ce zéro est le premier d'une séquence
                currently_in_zero_sequence = True #on démarre la séquence courante
                length_current = 1 #la longueur de la séquence actuelle vaut 1 pour le moment
                index_start_current = i #on dit que l'indice de départ de la séquence actuelle vaut i
        else: #sinon, si l'entrée qu'on regarde n'est pas un zéro :
            currently_in_zero_sequence = False # on dit qu'on n'est pas en train d'explorer une séquences de zéros
        if length_current > length_max: #si la longueur de la séquence courante est plus grande que la longueur max enregistrée
            length_max = length_current #on met à jour
            index_max = index_start_current #on met à jour
   
    return index_max, length_max
#Exemples :
print(find([0,0,1,0,1,0,0])) #(0,2)
print(find([1,0,1,0,1,0,0])) #(5,2)
print(find([0,0,0,0])) #(0,4)
print(find([1,1,1,1])) #(None, 0)
print(find([1,0,1])) #(1,1)
 

Bonne soirée

#7 Re : Programmation » Afficher le debut de la plus longue suite consecutive de zeros » 13-12-2023 19:28:13

oh la j'ai carement oublié de preciser que le tableau n'est constitué que de 0 et 1

Wiwaxia a écrit :

Bonjour,

Il faudra convenir de plusieurs choses:

a) la nature des éléments du tableau: il pourra s'agir d'entiers non signés au format Byte (par ex. {0. 1. 2}, ou éventuellement de booléens si l'on tient à allouer au tableau le minimum d'espace mémoire; les instructions dépendront du choix fait à ce niveau;
b) l'emplacement du premier zéro envisageable, par ex. la position (1); on aura par conséquent T(0) > 0 ;
c) l'origine ou la génération du tableau étudié, afin que cet objet soit directement observable : s'agit-il des décimales d'un nombre irrationnel, ou d'une séquence d'une suite modulaire arbitrairement choisie ? Mais peut-être as-tu une idée précise sur ce dernier point ...

Il faudra aussi deux variables de type entier pour mémoriser la plus grande longueur observée (Lmax) et la position du premier terme (Kini); toutes deux initialisées à zéro.

#8 Re : Programmation » Afficher le debut de la plus longue suite consecutive de zeros » 13-12-2023 19:27:01

mince je suis un peu perdu
je ne comprend pas comment vous proceder

Fred a écrit :

Re-

Une fonction comme la suivante devrait fonctionner :


def suite_tableau(T):
    index=-1
    longueur=-1
    i=0
    indexencours=0
    lencours=0
    l=0
    while (i<len(T)):
        if (T[i]==0):
            if (lencours==0):    # il y avait un 1 avant
                indexencours=i
            lencours+=1
            if (lencours>longueur):
                longueur=lencours
                index=indexencours
        else:
            lencours=0
        i+=1
    return (index,longueur)
 

#9 Programmation » Afficher le debut de la plus longue suite consecutive de zeros » 13-12-2023 15:04:38

pentium mix
Réponses : 11

bonjour tout le monde. j'aimerai écrire l’algorithme qui retourne la position i dans le tableau T telle que T (i) est le début de la plus longue suite consécutive de zéros.

malheureusement je ne sais quel fonction écrire au préalable pour aboutir. J'ai penser a la fonction qui retourne la position du premier zeros mais cela ne m'aide pas

Merci

#10 Entraide (supérieur) » Processus stochastique » 08-07-2023 19:02:36

pentium mix
Réponses : 2

Bonsoir a vous
S'il vous plaît depuis pratiquement une semaine je bloque sur une question.
S'il vous plaît dans quel cas un processus gaussien centré qui n'est pas un mouvement Brownien est un processus de Markov??
Merci

#11 Re : Entraide (supérieur) » Arithmétique vous plat » 18-04-2023 15:22:16

Glozi a écrit :

Bonjour,
Il faut connaitre (et utiliser) la definition de l'écriture d'un nombre dans une base.
Ex en base 10 : $354 = 3\times 10^2 +5\times 10^1+4\times 10^0$
Dans ta base $x$ on a donc par exemple $z= 1\times x^2 + 0 \times x + 1\times x^0$
Bonne journée

Merci bien

#12 Entraide (supérieur) » Arithmétique vous plat » 18-04-2023 09:45:43

pentium mix
Réponses : 2

Bonjour a tous
S'il vous plaît je bloque sur cet exercice depuis des jours.
https://www.cjoint.com/c/MDshR4zo8Vy

J'ai besoin d'aide s'il vous plaît
Merci

#13 Re : Entraide (supérieur) » Matrice d'une rotation » 31-03-2023 11:34:11

Michel Coste a écrit :

Puisque [tex]3a[/tex] est la trace de la matrice, [tex]a=\dfrac13(1+2\cos\theta)[/tex], où [tex]\theta[/tex] est l'angle de rotation. On porte cette valeur dans [tex]a^3-a^2+\dfrac4{27} \sin^2\varphi=0[/tex], et on obtient six valeurs possibles pour [tex]\theta[/tex] en fonction de [tex]\varphi[/tex]. C'est normal, vu les permutations sur [tex]a,b,c[/tex].


Waouhhh
Grand merci

Je n'avais pas compris ça


Merci beaucoup

#14 Re : Entraide (supérieur) » Matrice d'une rotation » 30-03-2023 15:59:20

pentium mix a écrit :
Michel Coste a écrit :

Avec plaisir.
Je me demande si pentium mix a compris mon indication pour le calcul de l'angle de rotation (et même s'il a trouvé l'axe de rotation).

Bien sur je comprend très bien
Sauf que dans mes recherches je suis tombé sur ceci
Du coup je suis un peu perdu
https://www.cjoint.com/c/MCEn0lIk1nc

La trace de la matrice d'une rotation est égale a 2cosa + 1 ou a est l'angle de la rotation
De la , j'ai l'angle de ma rotation
De plus le sous espace associé a la valeur propre 1 est l'axe de la rotation


C'est juste que je n'aboutit pas au résultat ci-dessus (pour l'angle seulement) et que nulle part dans mon raisonnement je n'utilise l'hypothèse k=......
Et, ça laisse croire que cette hypothèse ne sert a rien

Merci pour vos interventions.
J'ai beaucoup appris

#15 Re : Entraide (supérieur) » Matrice d'une rotation » 30-03-2023 15:54:02

Michel Coste a écrit :

Avec plaisir.
Je me demande si pentium mix a compris mon indication pour le calcul de l'angle de rotation (et même s'il a trouvé l'axe de rotation).

Bien sur je comprend très bien
Sauf que dans mes recherches je suis tombé sur ceci
Du coup je suis un peu perdu
https://www.cjoint.com/c/MCEn0lIk1nc

#16 Re : Entraide (supérieur) » Étude des suites recurentes » 29-03-2023 14:29:01

Ginger40 a écrit :
pentium mix a écrit :

Justement comment définir On et &n .

Je te conseille d'essayer de rédiger en LaTeX, ça aidera à la compréhension (même si j'ai compris que c'est $\theta_n$ et $\alpha_n$), ce n'est pas très compliqué et yoshi a fait un petit tuto pour débuter ici. (Petite astuce : pour connaître le code lié à une formule LaTeX écrite sur un post, on peut faire clique droit, Show Maths As, TeX Commands)

Sinon, comme je l'ai dit il faut utiliser $\arccos$ (ou $\cos^{-1}$, ce sont les mêmes fonctions) pour définir $\theta_n$.
Essaie d'appliquer d'abord formellement $\arccos$ à $c_n = \cos(\theta_n)$ pour trouver à quoi $\theta_n$ serait à priori égal.
Puis de manière rigoureuse (i.e. en regardant l'ensemble de définition et l'ensemble d'arrivée de $\arccos$), montre que tu peux définir de manière unique un $\theta_n$ vérifiant $c_n = \cos(\theta_n)$ et $\theta_n \in \left[0,\frac{\pi}{2}\right]$.

D'accord
Merci bien

#17 Re : Entraide (supérieur) » Matrice d'une rotation » 29-03-2023 12:57:13

Michel Coste a écrit :

OK.
Il n'est pas trop difficile de trouver l'axe de la rotation (droite propre de vecteur propre associé 1).
Pour l'angle de rotation [tex]\theta[/tex], on peut commencer par se souvenir de ce qu'est la trace d'une matrice de rotation d'angle [tex]\theta[/tex].

Quand je fait cela je n'utilise l'hypothèse donné nulle part

#18 Re : Entraide (supérieur) » Matrice d'une rotation » 29-03-2023 12:51:38

L'énoncé dis que a, b, c sont racine du polynôme. Ces racines sont réelles et c'est un polynôme de degré 3 sur R

plem06 a écrit :

Merci pour ce complément Michel, mais il manque encore quelques pièces à mon puzzle.
Ca m'énerve, je me dis que ça doit être juste devant mes yeux, mais je n'y arrive pas :

1)
a) Pour le déterminant, en fait c'est plutôt a3+b3+c3-3abc  (= 1)  non ?
selon première colonne : (-1)² a(a²-bc) + (-1)3 c(ab-c²)+(1)4b(b²-ac) ?
b) Si je ne me suis pas trompé ci-dessus (une vérif à la calculatrice semble me dire que non), comment trouve-t-on le a+b+c=1 alors ?
c) Merci pour mon oubli sur le fait que matrice orthogonale --> vecteurs orthonormés donc a2+b2+c2=1
Et merci aussi pour les relations racines/coeffs dont j'avais oublié certaines !

2) Pour la réciproque :
Quel est ce polynôme en fait ? Pas le polynome carac quand même ? Une rotation, c'est pas diagonalisable en général, donc  si c'était lui, il devrait pas avoir 3 racines de toute façon...
Tu me dis que les coeffs d'une rotation sont réels : soit, mais, sauf erreur, son polynome carac est pas scindé : en dimension 3, on n'a qu'une racine si elle n'est pas l'identité ou une symétrie axiale, non ?

Merci !

#19 Re : Entraide (supérieur) » Matrice d'une rotation » 27-03-2023 22:09:48

Michel Coste a écrit :

Bonjour,
Quelle est ta réponse à la première question ?

Pour que la matrice soit celle d'une rotation, il faut qu'elle soit orthogonale et de déterminant 1
C'est a dire que ab+ac+bc=0 , a^2+b^2+c^2=1 et a+b+c=1
Cela signifie que les coefficients a,b,c sont solution du polynôme donné avec K=bac
Ensuite pour terminer, ce trinôme doit avoir exactement 3 racines reelles
Ce qui permet, après études des variations du polynôme d'ecrire 0<=K<=4/27

Voilà

#20 Re : Entraide (supérieur) » Étude des suites recurentes » 27-03-2023 22:03:35

Justement comment définir On et &n .


Je ne sais pas comment les choisir pour que ça marche

#21 Entraide (supérieur) » Étude des suites recurentes » 27-03-2023 12:57:49

pentium mix
Réponses : 5

Bonjour
S'il vous plaît je ne sais vraiment pas comment faire
https://www.cjoint.com/c/MCBk4Dsgwh6 il n'ya pas autre indication
Merci

#22 Entraide (supérieur) » Matrice d'une rotation » 27-03-2023 12:37:22

pentium mix
Réponses : 20

Bonjour
S'il vous plait je n'arrive pas a déterminer les éléments caractéristiques de cette rotation
https://www.cjoint.com/c/MCBkHIpHdh6
J'ai déjà répondu a la première question je suis au niveau de la seconde
Merci

#23 Re : Entraide (supérieur) » Fonction bornée » 15-01-2023 13:50:11

Gui82 a écrit :

Bonjour,

En fait, il n'y a pas de [tex]- \infty[/tex] car on est dans [tex]\mathbb{R}^n[/tex] et on considère [tex]\|x\|[/tex] qui est positif.
Comme tu l'as dit, pour [tex]\displaystyle \varepsilon=1,\, \exists R>0[/tex] tel que [tex]\forall x \in \mathbb{R}^n,\, \|x\| \ge R \Longrightarrow |f(x)| \le 1[/tex]
Ensuite, sur [tex]\bar{B}(0,R)[/tex] tu peux voir ce qu'il se passe avec un argument simple de topologie.

Waouhhh

Merci bien.
Grand merci

#24 Re : Entraide (supérieur) » Fonction bornée » 15-01-2023 13:00:26

Fred a écrit :

Bonjour,

  Si tu veux traduire cela avec des epsilons, alors "tout simplement" :
$$\forall\varepsilon>0,\ \exists A>0,\ \|x\|\geq A\implies |f(x)|\leq \varepsilon.$$

Comme souvent dans ce genre d'exercices, il ne faut pas utiliser ceci pour toutes les valeurs de $\varepsilon$,
mais pour une valeur bien choisie....

F.

Merci bien bonjour. Pour epsilon égal a 1 cette définition montre que, autour de plus l'infini a fonction est bornée. Et autour de moins l'infini que se passe t'il??

Je suis confus.


J'ai pu résoudre le problème pour n=1 mais je n'arrive pas a generaliser

Pied de page des forums