σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft AcademyFeuilles de routeParcours
Ctrl K
AO
+1,2 k apprenants
ConnexionRejoindre ↗
Feuilles de routeParcours
AO
+1,2 k apprenants
Connexion
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft AcademyFeuilles de routeParcours
Ctrl K
AO
+1,2 k apprenants
ConnexionRejoindre ↗
Feuilles de routeParcours
AO
+1,2 k apprenants
Connexion
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Parcours Linear Algebra›Module 4 · Moindres carrés et factorisations›Sujet 4.3·~50 minIntermédiaire

Cholesky et matrices définies positives

Position dans le parcoursModule 4 sur 4 · Sujet 3 sur 4
4.3 Cholesky et matrices définies positives
—Score du lab— / 100XP total—Série—

Des racines carrées pour les matrices

Une matrice symétrique AAA est définie positive quand xTAx>0\mathbf{x}^{T}A\mathbf{x}>0xTAx>0 pour tout x≠0\mathbf{x}\neq\mathbf{0}x=0 ; de façon équivalente, toutes ses valeurs propres sont positives. Les matrices de covariance sont toujours semi-définies positives, car wTΣw\mathbf{w}^{T}\Sigma\mathbf{w}wTΣw est la variance d'un portefeuille et une variance ne peut être négative. Elles sont définies positives, tout portefeuille ayant une variance strictement positive, exactement quand aucun actif n'est redondant.

A=LLT,Ljj=Ajj−∑k<jLjk2,Lij=Aij−∑k<jLikLjkLjj (i>j)A=LL^{T},\qquad L_{jj}=\sqrt{A_{jj}-\sum_{k<j}L_{jk}^{2}},\qquad L_{ij}=\frac{A_{ij}-\sum_{k<j}L_{ik}L_{jk}}{L_{jj}}\ (i>j)A=LLT,Ljj​=Ajj​−k<j∑​Ljk2​​,Lij​=Ljj​Aij​−∑k<j​Lik​Ljk​​ (i>j)(1)

Toute matrice définie positive admet une factorisation de Cholesky A=LLTA=LL^{T}A=LLT avec LLL triangulaire inférieure à diagonale positive : l'analogue matriciel d'une racine carrée. Elle exploite la symétrie et coûte environ 13n3\tfrac13n^331​n3, la moitié de la LU, sans aucun pivotage. Si l'algorithme rencontre un nombre négatif sous une racine, la matrice n'est pas définie positive : Cholesky est donc un test bon marché de cette propriété.

Un facteur et une simulationPour A=(4225)A=\begin{pmatrix}4&2\\2&5\end{pmatrix}A=(42​25​) : L11=4=2L_{11}=\sqrt4=2L11​=4​=2, L21=2/2=1L_{21}=2/2=1L21​=2/2=1, L22=5−12=2L_{22}=\sqrt{5-1^2}=2L22​=5−12​=2, donc L=(2012)L=\begin{pmatrix}2&0\\1&2\end{pmatrix}L=(21​02​) et LLT=(4225)LL^{T}=\begin{pmatrix}4&2\\2&5\end{pmatrix}LLT=(42​25​). Pour simuler deux rendements de corrélation 0,60{,}60,6 et de variances unitaires, factorisez Σ=(10,60,61)\Sigma=\begin{pmatrix}1&0{,}6\\0{,}6&1\end{pmatrix}Σ=(10,6​0,61​) en L=(100,60,8)L=\begin{pmatrix}1&0\\0{,}6&0{,}8\end{pmatrix}L=(10,6​00,8​) et posez x=Lz\mathbf{x}=L\mathbf{z}x=Lz avec z\mathbf{z}z des normales standard indépendantes : alors Cov⁡(x)=LLT=Σ\operatorname{Cov}(\mathbf{x})=LL^{T}=\SigmaCov(x)=LLT=Σ.
  • Définie positive : xTAx>0\mathbf{x}^{T}A\mathbf{x}>0xTAx>0 pour tout x≠0\mathbf{x}\neq\mathbf{0}x=0 ; toutes les valeurs propres positives.
  • Cholesky : A=LLTA=LL^{T}A=LLT, moitié du coût de la LU, sans pivotage, et test bon marché de la définie-positivité.
  • Tirages corrélés : x=Lz\mathbf{x}=L\mathbf{z}x=Lz transforme des normales indépendantes z\mathbf{z}z en rendements de covariance Σ\SigmaΣ.
Une covariance empirique peut ne pas se factoriserSi vous estimez une covariance avec moins d'observations que d'actifs, ou si deux actifs ont des rendements identiques, la matrice empirique n'est que semi-définie et Cholesky s'arrête sur un pivot nul ou légèrement négatif. Les remèdes : plus de données, supprimer les doublons, un rétrécissement vers une cible bien conditionnée, ou l'ajout à la diagonale d'un petit multiple de l'identité.

Lab : ce qui est noté

  • Factoriser une matrice20 pts
  • Factoriser une matrice de corrélation20 pts
  • Un tirage corrélé15 pts
  • Quelle matrice est définie positive15 pts
  • Quand Cholesky échoue10 pts
  • Simuler des rendements corrélés20 pts
0/ 6 traités
Exercices
1Matrice · vecteur20 ptsSans réponse

Remplissez le facteur de Cholesky LLL de A=(9335)A=\begin{pmatrix}9&3\\3&5\end{pmatrix}A=(93​35​) (triangulaire inférieure, diagonale positive).

L=L =L=
2Matrice · vecteur20 ptsSans réponse

Remplissez le facteur de Cholesky de la matrice de corrélation (10,50,51)\begin{pmatrix}1&0{,}5\\0{,}5&1\end{pmatrix}(10,5​0,51​) (décimaux acceptés, trois chiffres).

L=L =L=
3Réponse numérique15 ptsSans réponse

Avec Σ=(10,60,61)\Sigma=\begin{pmatrix}1&0{,}6\\0{,}6&1\end{pmatrix}Σ=(10,6​0,61​) et L=(100,60,8)L=\begin{pmatrix}1&0\\0{,}6&0{,}8\end{pmatrix}L=(10,6​00,8​), le tirage indépendant z=(1,2)\mathbf{z}=(1,2)z=(1,2) donne x=Lz\mathbf{x}=L\mathbf{z}x=Lz. Que vaut x2x_2x2​ ?

x2=x_2 =x2​=accepté : ± 0.001
4Choix15 ptsSans réponse

Laquelle de ces matrices symétriques est définie positive ?

5Choix10 ptsSans réponse

Cholesky échoue sur une matrice de covariance empirique. Sélectionnez toutes les causes probables.

6Expliquer avec des mots20 ptsSans réponse

En deux phrases, dites comment le facteur de Cholesky d'une matrice de covariance sert à simuler des rendements corrélés.

Mode mathsSaisie au clavier ou en LaTeXExercice 1 sur 6 · Ctrl ⏎ pour vérifier

Connectez-vous pour exécuter votre code.

Tous les cours sont gratuits à lire. Créez un compte ou connectez-vous pour exécuter du code, lancer les tests et soumettre ce lab.

  • IDE dans le navigateur avec terminal bash — rien à installer
  • Correction instantanée, XP et votre place dans le parcours
Créer un compteSe connecter
OU
Continuer avec GitHub