The Heston stochastic-local volatility model, consisting of a asset price process and a Cox--Ingersoll--Ross-type variance process, offers a wide range of applications in the financial industry. The pursuit for efficient model evaluation has been assiduously ongoing and central to which is the numerical simulation of CIR process. Different from the weakly convergent noncentral chi-squared approximation used in 25, this paper considers two strongly convergent and positivity-preserving methods for CIR process under Lamperti transformation, namely, the truncated Euler method and the backward Euler method. It should be noted that these two methods are completely different. The explicit truncated Euler method is computationally effective and remains robust under high volatility, while the implicit backward Euler method provides high computational accuracy and stable performance. Numerical experiments on European call options are presented to show the superiority of different methods.
academic- معرّف الورقة البحثية: 2509.24449
- العنوان: محاكاة فعالة لأسعار خيارات الاستدعاء الأوروبية تحت نموذج Heston للتقلب العشوائي-المحلي: مقارنة الطرق
- المؤلفون: Meng Cai, Tianze Li
- المؤسسات: كلية الإحصاء والرياضيات والاقتصاد، جامعة المالية والاقتصاد المركزية
- التصنيف: q-fin.CP (التمويل الحسابي)
- تاريخ النشر: 16 أكتوبر 2025
- رابط الورقة: https://arxiv.org/abs/2509.24449
تدرس هذه الورقة نموذج Heston للتقلب العشوائي-المحلي المكون من عملية سعر الأصل وعملية التباين من نوع Cox-Ingersoll-Ross، وهو نموذج يتمتع بتطبيقات واسعة في الصناعة المالية. بخلاف التقريب الضعيف للتقارب باستخدام توزيع كاي-تربيع غير المركزي المستخدم في المرجع 25، تدرس هذه الورقة طريقتين يتمتعان بتقارب قوي وحفظ الإيجابية تحت تحويل Lamperti لعملية CIR: طريقة Euler المقطوعة وطريقة Euler العكسية. تتمتع طريقة Euler المقطوعة الصريحة بكفاءة حسابية عالية وتحافظ على الاستقرار عند التقلبات العالية، بينما توفر طريقة Euler العكسية الضمنية دقة حسابية عالية واستقراراً متفوقاً. تُظهر التجارب الرقمية على خيارات الاستدعاء الأوروبية تفوق الطرق المختلفة.
- متطلبات النموذج: يتطلب الهندسة المالية تحقيق توازن بين المعايرة المثالية والديناميكيات الصحيحة، حيث يوفر نموذج التقلب العشوائي-المحلي (SLV) حلاً لهذه المشكلة
- قيود النماذج الموجودة:
- تتمتع نماذج التقلب العشوائي البحت بمعاملات محدودة، مما يجعل من الصعب معايرة التقلب الضمني بشكل مثالي لجميع آجال الاستحقاق وأسعار التنفيذ
- على الرغم من أن نماذج التقلب المحلي يمكنها معايرة جميع الخيارات الحالية في السوق، إلا أنها تواجه مشاكل في التنبؤ بسلوك التقلب المستقبلي، مما يؤدي إلى تسطح ابتسامة التقلب الأمامية
محاكاة عملية CIR رقمياً تمثل نقطة الاختناق الحرجة في تقييم النموذج، لأن:
- الحد الانتشاري غير محدود والمعاملات ليست Lipschitz عالمياً
- قد تنتج طريقة Euler-Maruyama التقليدية تباينات سالبة، وهو ما لا معنى له رياضياً وماليّاً
- يمكن لتقريبات كاي-تربيع غير المركزية الموجودة فقط تحقيق تقارب ضعيف
تطوير طرق تتمتع بتقارب قوي وحفظ الإيجابية، مما يضمن أن تكون المسارات بين الحل الدقيق والحل التقريبي قريبة من بعضها، وبالتالي تحسين تسعير الخيارات الإجمالي.
- اقتراح طريقتين جديدتين للحساب الرقمي: طريقة Euler المقطوعة وطريقة Euler العكسية، وكلاهما يتمتع بالتقارب القوي وحفظ الإيجابية
- التحليل النظري: إثبات التقارب القوي لطريقة Euler المقطوعة (رتبة التقارب 1/2)
- مقارنة رقمية شاملة: من خلال تجارب تسعير خيارات الاستدعاء الأوروبية، مقارنة منهجية لأداء الطرق المختلفة
- إرشادات عملية: توفير توصيات اختيار الطريقة لسيناريوهات التطبيق المختلفة
حل مشكلة تسعير خيارات الاستدعاء الأوروبية تحت نموذج Heston للتقلب العشوائي-المحلي، حيث يكون المحور الأساسي هو محاكاة فعالة وسقيمة لعملية التباين CIR.
dS_t = rS_t dt + \sqrt{V_t} S_t \sigma(t, S_t)\left(\rho dW_t + \sqrt{1-\rho^2} d\tilde{W}_t\right) \\
dV_t = \kappa(\theta - V_t)dt + \gamma\sqrt{V_t} dW_t \\
dW_t d\tilde{W}_t = \rho dt
\end{cases}$$
حيث:
- $S_t$: عملية سعر الأصل
- $V_t$: عملية التباين (ديناميكيات CIR)
- $\sigma(t,s)$: دالة التقلب المحلي
- $\kappa, \theta, \gamma$: معاملات عملية CIR
#### تحويل Lamperti
للتغلب على الصعوبات الرقمية لعملية CIR، يتم إدخال التحويل $L_t = \sqrt{V_t}$:
$$dL_t = \frac{1}{2}\kappa\left(\frac{\theta}{L_t} - L_t\right)dt + \frac{1}{2}\gamma dW_t$$
### الطرق الرقمية الأساسية
#### 1. طريقة Euler المقطوعة
تعريف دالة القطع: $\pi_\tau(x) = (b\tau^{1/4}) \vee x$
الصيغة الرقمية:
$$L^N_{t_{n+1}} = L^N_{t_n} + \frac{\kappa}{2}\tau\left(\frac{\theta}{\pi_\tau(L^N_{t_n})} - \pi_\tau(L^N_{t_n})\right) + \frac{1}{2}\gamma\Delta W_n$$
**الخصائص**:
- صيغة صريحة، حسابياً فعالة
- تضمن الإيجابية من خلال القطع
- تحافظ على الاستقرار عند التقلبات العالية
#### 2. طريقة Euler العكسية
$$L^N_{t_{n+1}} = L^N_{t_n} + \frac{\kappa}{2}\tau\left(\frac{\theta}{L^N_{t_{n+1}}} - L^N_{t_{n+1}}\right) + \frac{1}{2}\gamma\Delta W_n$$
**الخصائص**:
- صيغة ضمنية، تتطلب حلاً تكرارياً
- توفر دقة حسابية عالية
- أداء مستقرة لكن التكلفة الحسابية عالية
### نقاط الابتكار التقني
1. **التقارب القوي**: بخلاف طريقة كاي-تربيع غير المركزية ذات التقارب الضعيف، تضمن تقارب على مستوى المسار
2. **تصميم حفظ الإيجابية**: ضمان أن يكون التباين موجباً دائماً من خلال التحويلات الرياضية وتقنيات القطع
3. **التكامل بين الطرق**: تتمتع الطريقتان بمزايا متبادلة في الكفاءة والدقة، مما يناسب سيناريوهات مختلفة
## إعداد التجارب
### معاملات النموذج
- **معاملات Heston الأساسية**: $\gamma = 0.95$, $\kappa = 1.05$, $\rho = -0.315$, $\theta = 0.0855$, $S_0 = 1$, $r = 0$, $v_0 = 0.0945$, $T = 5$
- **معاملات تصحيح HSLV**: إدخال معامل التصحيح $p = 0.25$، لتعديل المعاملات المختلفة
### مؤشرات التقييم
- نسبة الخطأ النسبي المئوية
- وقت الحساب
- الاستقرار في ظروف السوق المختلفة
### طرق المقارنة
- طريقة Euler القياسية
- طريقة AES (محاكاة Andersen الدقيقة)
- طريقة Euler المقطوعة
- طريقة Euler العكسية
### تصميم التجارب
1. **حساب التوقع الشرطي**: استخدام طريقة 2D-COS كمعيار، 20 صندوقاً، 10,000 مسار، حجم الخطوة 0.001
2. **تحليل حساسية المعاملات**: اختبار الأداء تحت قيم $\theta$ و $p$ مختلفة
3. **مقارنة أسعار التنفيذ**: $K = 70\%, 100\%, 150\%$ في ثلاث حالات
## نتائج التجارب
### النتائج الرئيسية
#### مقارنة الأخطاء تحت عدد خطوات مختلف
| عدد الخطوات N | K=70% | K=100% | K=150% |
|-------|--------|---------|---------|
| **طريقة Euler المقطوعة** |
| 5 | 13.49%(0.24%) | 14.08%(0.24%) | 13.96%(0.85%) |
| 40 | 4.66%(0.53%) | 3.86%(0.74%) | 4.80%(1.35%) |
| **طريقة Euler العكسية** |
| 5 | 15.49%(0.28%) | 15.97%(0.10%) | 17.32%(1.58%) |
| 40 | 0.06%(0.66%) | 0.46%(0.68%) | 1.44%(1.15%) |
#### مقارنة وقت الحساب (40 خطوة)
- طريقة AES: 0.1114 ثانية
- طريقة Euler: 0.0939 ثانية
- طريقة Euler المقطوعة: 0.0733 ثانية
- طريقة Euler العكسية: 137.9476 ثانية
### الاكتشافات الرئيسية
1. **عند عدد خطوات قليل**: عادة ما تحقق طريقة Euler المقطوعة أقل خطأ
2. **عند عدد خطوات كبير**: تحقق طريقة Euler العكسية أعلى دقة
3. **استقرار المعاملات**: تحافظ طريقة Euler المقطوعة على استقرار ممتاز مع نمو بطيء للخطأ عند زيادة $\theta$
4. **الكفاءة الحسابية**: تتمتع طريقة Euler المقطوعة بأسرع حساب، بينما تستغرق طريقة Euler العكسية وقتاً طويلاً بسبب الحل التكراري
### التجارب الاستئصالية
#### تأثير معامل التقلب $\theta$
تُظهر التجارب أنه مع زيادة $\theta$ (محاكاة ظروف السوق القاسية):
- يزداد خطأ طرق AES و Euler القياسية بشكل ملحوظ
- تحافظ طريقة Euler المقطوعة على الاستقرار، مع معدل نمو خطأ منخفض
- يشير هذا إلى أن طريقة القطع أكثر ملاءمة لاختبارات الضغط ونمذجة المخاطر الطرفية
#### تأثير معامل التصحيح $p$
مع زيادة $p$، تحقق طريقة Euler العكسية أفضل أداء دائماً، لكنها تتمتع بحساسية أعلى لحجم الخطوة.
## الأعمال ذات الصلة
### الاتجاهات البحثية الرئيسية
1. **معايرة نموذج SLV**: العمل الرائد لـ Guyon و Henry-Labordère (2011)
2. **طرق رقمية لعملية CIR**: تشمل الفروقات المحدودة وطرق جزيئات Monte Carlo
3. **تطبيقات التعلم العميق**: استخدام الشبكات العدائية التوليدية للمعايرة بواسطة Cuchiero وآخرين (2020)
### مزايا هذه الورقة
بالمقارنة مع الأعمال الموجودة، توفر هذه الورقة طرقاً ذات تقارب قوي وحفظ إيجابية، مما يسد الفجوة التي تركتها الطرق ذات التقارب الضعيف، ويوفر أدوات رقمية أكثر موثوقية للتطبيقات العملية.
## الخلاصة والنقاش
### الاستنتاجات الرئيسية
1. **التكامل بين الطرق**: تناسب طريقة Euler المقطوعة المحاكاة واسعة النطاق وإدارة المخاطر في الوقت الفعلي؛ تناسب طريقة Euler العكسية اختبارات الضغط والسيناريوهات التي تتطلب دقة عالية
2. **مزايا الاستقرار**: تُظهر طريقة Euler المقطوعة أداء ممتازة في ظروف السوق القاسية
3. **المقايضة بين الدقة والكفاءة**: يتطلب اختيار حجم الخطوة المناسب والطريقة المناسبة وفقاً لمتطلبات التطبيق المحددة
### القيود
1. التكلفة الحسابية العالية لطريقة Euler العكسية
2. يتطلب اختيار حجم الخطوة المقايضة بين الدقة والكفاءة
3. لا تزال استقرار الطريقة بحاجة إلى التحقق الإضافي لبعض مجموعات المعاملات القاسية
### الاتجاهات المستقبلية
1. تطوير استراتيجيات حجم خطوة تكيفية
2. استكشاف التنفيذ المتوازي
3. التوسع إلى متغيرات نموذج SLV الأكثر تعقيداً
## التقييم المتعمق
### المزايا
1. **الصرامة النظرية**: توفير إثبات تقارب كامل
2. **القيمة العملية العالية**: توفر الطريقتان مزايا متبادلة، مما يوفر خيارات مرنة للتطبيقات العملية
3. **التجارب الشاملة**: مقارنة منهجية لأداء الطرق المختلفة في ظروف متنوعة
4. **الكتابة الواضحة**: التعبير الرياضي الدقيق والهيكل المنطقي الواضح
### أوجه القصور
1. **عمق التحليل النظري**: التحليل النظري لطريقة Euler العكسية نسبياً أقل
2. **التحقق من السوق الفعلية**: نقص التحقق باستخدام بيانات السوق الحقيقية
3. **تحليل التعقيد الحسابي**: عدم توفير تحليل نظري تفصيلي للتعقيد
### التأثير
1. **المساهمة الأكاديمية**: توفير طرق تقارب قوية جديدة لمحاكاة عملية CIR
2. **القيمة العملية**: قابلة للتطبيق المباشر على تسعير الخيارات وإدارة المخاطر في المؤسسات المالية
3. **قابلية إعادة الإنتاج**: وصف الطريقة مفصل، سهل التنفيذ
### السيناريوهات المناسبة
1. **التداول الكمي**: تسعير الخيارات عالي التردد واستراتيجيات التحوط
2. **إدارة المخاطر**: اختبارات الضغط وتقييم مخاطر الذيل
3. **معايرة النموذج**: تقدير المعاملات والتحقق من نموذج SLV
## المراجع
تستشهد الورقة بـ 25 مرجعاً ذا صلة، تغطي النظرية الأساسية لنموذج SLV وطرق رقمية لعملية CIR وتسعير الخيارات وغيرها من المجالات الرئيسية، مما يوفر أساساً نظرياً متيناً للبحث.