σμ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 Academyخرائط الطريقالمسارات
Ctrl K
AO
+1.2k متعلّم
تسجيل الدخولانضم ↗
خرائط الطريقالمسارات
AO
+1.2k متعلّم
تسجيل الدخول
  • 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 Academyخرائط الطريقالمسارات
Ctrl K
AO
+1.2k متعلّم
تسجيل الدخولانضم ↗
خرائط الطريقالمسارات
AO
+1.2k متعلّم
تسجيل الدخول
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
مسار Linear Algebra›الوحدة 2 · حل الأنظمة الخطية›الموضوع 2.2·~50 دمتوسط

تحليل LU

الموضع في المسارالوحدة 2 من 4 · الموضوع 2 من 3
2.2 تحليل LU
—درجة المختبر— / 100إجمالي XP—التتابع—

الحذف محفوظًا

يحل حذف غاوس المعادلة Ax=bA\mathbf{x}=\mathbf{b}Ax=b بطرح مضاعفات صف من الصفوف التي تحته حتى تصبح المصفوفة مثلثية عليا UUU. والمضاعِفات التي استعملتها تستحق الحفظ: ضعها تحت قطر مصفوفة مثلثية سفلى LLL قطرها آحاد فيكون A=LUA=LUA=LU. يسجّل هذا التحليل عملية الحذف كلها، فيُنجَز العمل مرة واحدة ويُعاد استخدامه.

A=LU,Ly=b,Ux=yA=LU,\qquad L\mathbf{y}=\mathbf{b},\qquad U\mathbf{x}=\mathbf{y}A=LU,Ly=b,Ux=y(1)

للحل نحل أولًا Ly=bL\mathbf{y}=\mathbf{b}Ly=b بالتعويض الأمامي (من الأعلى إلى الأسفل) ثم Ux=yU\mathbf{x}=\mathbf{y}Ux=y بالتعويض الخلفي (من الأسفل إلى الأعلى). وكلٌّ منهما يكلّف نحو n2n^2n2 عملية، أما التحليل فيكلّف نحو 23n3\tfrac{2}{3}n^332​n3. ومع كثرة الطرف الأيمن (سيناريوهات كثيرة على المصفوفة نفسها) تحلّل مرة واحدة وتدفع n2n^2n2 فقط لكل سيناريو.

تحليل LU يدويًالتكن A=(2145)A=\begin{pmatrix}2&1\\4&5\end{pmatrix}A=(24​15​). المضاعِف هو 4/2=24/2=24/2=2: اطرح ضعف الصف 1 من الصف 2 فتحصل على U=(2103)U=\begin{pmatrix}2&1\\0&3\end{pmatrix}U=(20​13​) وL=(1021)L=\begin{pmatrix}1&0\\2&1\end{pmatrix}L=(12​01​). تحقّق: LU=(2142+3)=ALU=\begin{pmatrix}2&1\\4&2+3\end{pmatrix}=ALU=(24​12+3​)=A. لنحل Ax=(3,9)A\mathbf{x}=(3,9)Ax=(3,9): Ly=bL\mathbf{y}=\mathbf{b}Ly=b يعطي y1=3y_1=3y1​=3 وy2=9−2⋅3=3y_2=9-2\cdot3=3y2​=9−2⋅3=3؛ ثم Ux=yU\mathbf{x}=\mathbf{y}Ux=y يعطي x2=1x_2=1x2​=1 و2x1+1=32x_1+1=32x1​+1=3 فيكون x=(1,1)\mathbf{x}=(1,1)x=(1,1).
  • LLL: مثلثية سفلى، آحاد على القطر، ومضاعِفات الحذف تحته.
  • UUU: مثلثية عليا ناتج الحذف؛ قطرها يضم المحاور.
  • التمحور: بدّل الصفوف ليكون المحور كبيرًا فيُنتج PA=LUPA=LUPA=LU بمصفوفة تبديل PPP.
المحور الصفري أو الضئيل يحتاج إلى تبديل صفوففي (0111)\begin{pmatrix}0&1\\1&1\end{pmatrix}(01​11​) المحور الأول 0 فيقسم الحذف البسيط على صفر مع أن المصفوفة قابلة للقلب تمامًا. وتبديل الصفين يحل المشكلة. أما المحور الضئيل (لا الصفري) فلا يسبب انهيارًا لكنه ينتج مضاعِفات ضخمة ويفسد الدقة، لذلك تستعمل المكتبات الحقيقية التمحور دائمًا.

المختبر: ما الذي يُقيَّم

  • المصفوفة L15 نقطة
  • المصفوفة U15 نقطة
  • محدّد من U15 نقطة
  • الحل بـ L وU15 نقطة
  • لماذا التمحور15 نقطة
  • خطوات حل LU10 نقطة
  • لماذا التحليل مرة واحدة15 نقطة
0/ 7 تمت الإجابة
المهام
1مصفوفة · متجه15 نقطةبلا إجابة

لتكن A=(2145)A=\begin{pmatrix}2&1\\4&5\end{pmatrix}A=(24​15​). املأ العامل المثلثي السفلي LLL (آحاد على القطر) في A=LUA=LUA=LU دون تبديل صفوف.

L=L =L=
2مصفوفة · متجه15 نقطةبلا إجابة

لتكن A=(2145)A=\begin{pmatrix}2&1\\4&5\end{pmatrix}A=(24​15​). املأ العامل المثلثي العلوي UUU.

U=U =U=
3إجابة رقمية15 نقطةبلا إجابة

تتحلل المصفوفة A=(211433879)A=\begin{pmatrix}2&1&1\\4&3&3\\8&7&9\end{pmatrix}A=​248​137​139​​ دون تبديل إلى L=(100210431)L=\begin{pmatrix}1&0&0\\2&1&0\\4&3&1\end{pmatrix}L=​124​013​001​​ وU=(211011002)U=\begin{pmatrix}2&1&1\\0&1&1\\0&0&2\end{pmatrix}U=​200​110​112​​. أوجد det⁡A\det AdetA.

det⁡A=\det A =detA=المقبول: ± 0.001
4إجابة رقمية15 نقطةبلا إجابة

لتكن A=(2145)A=\begin{pmatrix}2&1\\4&5\end{pmatrix}A=(24​15​). استخدم L=(1021)L=\begin{pmatrix}1&0\\2&1\end{pmatrix}L=(12​01​) وU=(2103)U=\begin{pmatrix}2&1\\0&3\end{pmatrix}U=(20​13​) لحل Ax=(5,13)A\mathbf{x}=(5,13)Ax=(5,13). ما قيمة x1x_1x1​؟

x1=x_1 =x1​=المقبول: ± 0.001
5اختيار15 نقطةبلا إجابة

ما الذي يحدث إذا طبّقت الحذف البسيط دون تبديل صفوف على (0111)\begin{pmatrix}0&1\\1&1\end{pmatrix}(01​11​)؟

6رتّب الخطوات10 نقطةبلا إجابة

رتّب خطوات حل Ax=bA\mathbf{x}=\mathbf{b}Ax=b بتحليل LU مع التمحور.

اسحب لإعادة الترتيب، أو ركّز على عنصر واضغط مسافة لرفعه والأسهم لتحريكه ومسافة لإفلاته.

  1. 1التعويض الأمامي: حل Ly=PbL\mathbf{y}=P\mathbf{b}Ly=Pb
  2. 2تحليل PA=LUPA=LUPA=LU
  3. 3التعويض الخلفي: حل Ux=yU\mathbf{x}=\mathbf{y}Ux=y
  4. 4تطبيق تبديلات الصفوف على الطرف الأيمن: PbP\mathbf{b}Pb
7اشرح بكلماتك15 نقطةبلا إجابة

عليك حل Ax=bA\mathbf{x}=\mathbf{b}Ax=b لألف متجه b\mathbf{b}b مختلف بالمصفوفة AAA نفسها. في جملتين بيّن كيف يفيد تحليل LU.

وضع الرياضياتإدخال بلوحة المفاتيح أو LaTeXالمهمة 1 من 7 · Ctrl ⏎ للتحقق

سجّل الدخول لتشغيل الكود الخاص بك.

قراءة كل المحاضرات مجانية. أنشئ حسابًا أو سجّل الدخول لتشغيل الكود وتنفيذ الاختبارات وتسليم هذا المختبر.

  • بيئة تطوير داخل المتصفح مع طرفية bash — لا حاجة لتثبيت أي شيء
  • تصحيح فوري ونقاط خبرة ومكانك في المسار
إنشاء حسابتسجيل الدخول
أو
المتابعة عبر GitHub