This paper considers the orthogonal expansion of the fractional Brownian motion relative to the Legendre polynomials. Such an expansion has not only theoretical but also practical interest, since it can be applied to approximate and simulate the fractional Brownian motion in continuous time. The relations for the mean square approximation error are presented, and a comparison with the previously obtained result is carried out.
معرّف الورقة : 2503.04575العنوان : On Approximate Representation of Fractional Brownian Motionالمؤلف : Konstantin A. Rybakovالتصنيف : math.PR (نظرية الاحتمالات)تاريخ النشر : مارس 2025 (نسخة أولية على arXiv، الإصدار v2 في 14 أكتوبر 2025)رابط الورقة : https://arxiv.org/abs/2503.04575 تدرس هذه الورقة التوسع المتعامد لحركة براونيان الكسرية بالنسبة إلى متعددات حدود لجندر. يتمتع هذا التوسع بأهمية نظرية وعملية، حيث يمكن استخدامه للتقريب والمحاكاة في الوقت المستمر لحركة براونيان الكسرية. تقدم الورقة صيغاً لخطأ التقريب بمتوسط المربعات وتقارنها مع النتائج المحققة سابقاً.
المشكلة الأساسية : كيفية تقريب ومحاكاة حركة براونيان الكسرية (FBM) بفعالية في الوقت المستمرالأهمية : لحركة براونيان الكسرية تطبيقات واسعة في النمذجة المالية والفيزياء والهندسة وغيرها، وتعتبر محاكاتها حاسمة للتحقق من النتائج النظرية وحل المشاكل العمليةقيود الطرق الموجودة :
تستطيع طرق الشبكة الزمنية المنفصلة التقليدية ضمان الدقة فقط عند نقاط الشبكة تتمتع طرق محاكاة العمليات الغاوسية العامة الموجودة بتعقيد حسابي عالي تفتقر إلى طرق الوقت المستمر التي تحسب خطأ التقريب بدقة الدافع البحثي : بناءً على الأعمال السابقة للمؤلف 2 ، اقتراح طريقة أبسط لحساب معاملات التوسع، مما يقلل التعقيد الحسابي ويحسن الدقةالمساهمة النظرية : اشتقاق علاقات صريحة لمعاملات توسع الدالة النواة المرتبطة بالتمثيل التكاملي لحركة براونيان الكسرية (بالنسبة إلى متعددات حدود لجندر)تحسين الطريقة : توفير طريقة أبسط لحساب معاملات التوسع مقارنة بالعمل السابق 2 ، مما يقلل التعقيد الحسابيتحليل الخطأ : إنشاء معادلات تمكن من حساب خطأ التقريب بمتوسط المربعات لحركة براونيان الكسرية بدقةالتقريب في الوقت المستمر : تحقيق تمثيل تقريبي لحركة براونيان الكسرية في الوقت المستمر، وليس فقط عند نقاط زمنية منفصلةبالنظر إلى حركة براونيان الكسرية B H ( t ) B_H(t) B H ( t ) (مع مؤشر Hurst H ∈ ( 0 , 1 ) H \in (0,1) H ∈ ( 0 , 1 ) )، البحث عن معاملات توسعها بالنسبة إلى الأساس المتعامد لمتعددات حدود لجندر، لتحقيق التقريب في الوقت المستمر:
B H ( t ) = ∑ i = 0 ∞ B H i P ^ ( i , t ) B_H(t) = \sum_{i=0}^{\infty} B_H^i \hat{P}(i,t) B H ( t ) = ∑ i = 0 ∞ B H i P ^ ( i , t )
حيث P ^ ( i , t ) \hat{P}(i,t) P ^ ( i , t ) هي متعددات حدود لجندر المعيارية.
تعريف حركة براونيان الكسرية :
دالة التغاير: R H ( t , τ ) = t 2 H + τ 2 H − ∣ t − τ ∣ 2 H 2 R_H(t,\tau) = \frac{t^{2H} + \tau^{2H} - |t-\tau|^{2H}}{2} R H ( t , τ ) = 2 t 2 H + τ 2 H − ∣ t − τ ∣ 2 H التمثيل التكاملي: B H ( t ) = ∫ 0 t k H ( t , τ ) d B ( τ ) B_H(t) = \int_0^t k_H(t,\tau) dB(\tau) B H ( t ) = ∫ 0 t k H ( t , τ ) d B ( τ ) الدالة النواة :
k H ( t , τ ) = a H ( t − τ ) H − 1 / 2 2 F 1 ( 1 2 − H , H − 1 2 , H + 1 2 , 1 − t τ ) 1 ( t − τ ) k_H(t,\tau) = a_H(t-\tau)^{H-1/2} \,_2F_1\left(\frac{1}{2}-H, H-\frac{1}{2}, H+\frac{1}{2}, 1-\frac{t}{\tau}\right)\mathbf{1}(t-\tau) k H ( t , τ ) = a H ( t − τ ) H − 1/2 2 F 1 ( 2 1 − H , H − 2 1 , H + 2 1 , 1 − τ t ) 1 ( t − τ )
حيث a H = π H ( 1 − 2 H ) Γ ( 2 − 2 H ) cos π H a_H = \sqrt{\frac{\pi H(1-2H)}{\Gamma(2-2H)\cos\pi H}} a H = Γ ( 2 − 2 H ) c o s π H π H ( 1 − 2 H ) (عندما H ≠ 1 / 2 H \neq 1/2 H = 1/2 ).
النظرية 1 (التعبير الصريح لمعاملات التوسع):
للعمود j j j من المصفوفة K H K^H K H :
K ∗ j H = a H Γ ( 3 / 2 − H ) 2 j + 1 T ∑ k = 0 j l j k ( 3 / 2 − H ) k T k ( H + 1 / 2 + k ) k ! F H + 1 / 2 + k K^H_{*j} = a_H\Gamma(3/2-H)\sqrt{\frac{2j+1}{T}} \sum_{k=0}^j \frac{l_{jk}(3/2-H)_k}{T^k(H+1/2+k)k!} F^{H+1/2+k} K ∗ j H = a H Γ ( 3/2 − H ) T 2 j + 1 ∑ k = 0 j T k ( H + 1/2 + k ) k ! l jk ( 3/2 − H ) k F H + 1/2 + k
بالنسبة إلى H ≠ 1 / 2 H \neq 1/2 H = 1/2 ، معاملات التوسع تحقق:
K i j H = a H Γ ( 1 / 2 − H ) 2 j + 1 T F i H + 1 / 2 ∑ k = 0 j ( − 1 ) j − k 1 / 2 − H + k H + 1 / 2 + k Π k ( i j ) K^H_{ij} = a_H\Gamma(1/2-H)\sqrt{\frac{2j+1}{T}} F^{H+1/2}_i \sum_{k=0}^j (-1)^{j-k} \frac{1/2-H+k}{H+1/2+k} \Pi^{(ij)}_k K ij H = a H Γ ( 1/2 − H ) T 2 j + 1 F i H + 1/2 ∑ k = 0 j ( − 1 ) j − k H + 1/2 + k 1/2 − H + k Π k ( ij )
حساب المعاملات : استخدام العلاقات التكرارية (12)-(13) لحساب معاملات توسع الدوال الأسيةبناء المصفوفة : بناء مصفوفة معاملات التوسع K H K^H K H من خلال النظرية 1توليد المتغيرات العشوائية : توليد متغيرات عشوائية طبيعية معيارية مستقلة V i V_i V i بناء التقريب : B ~ H ( t ) = ∑ i = 0 L − 1 B ~ i H P ^ ( i , t ) \tilde{B}_H(t) = \sum_{i=0}^{L-1} \tilde{B}^H_i \hat{P}(i,t) B ~ H ( t ) = ∑ i = 0 L − 1 B ~ i H P ^ ( i , t ) الفترة الزمنية : T = [ 0 , 1 ] T = [0,1] T = [ 0 , 1 ] مؤشر Hurst : H = 0.1 , 0.2 , … , 0.9 H = 0.1, 0.2, \ldots, 0.9 H = 0.1 , 0.2 , … , 0.9 طول القطع : L = 4 , 8 , 16 , 32 , 64 , 128 L = 4, 8, 16, 32, 64, 128 L = 4 , 8 , 16 , 32 , 64 , 128 خطأ التقريب بمتوسط المربعات :
ε = E ∫ T ( B H ( t ) − B ~ H ( t ) ) 2 d t \varepsilon = E\int_T (B_H(t) - \tilde{B}_H(t))^2 dt ε = E ∫ T ( B H ( t ) − B ~ H ( t ) ) 2 d t
الصيغة النظرية :
ε = ∥ k H ( ⋅ ) ∥ L 2 ( T 2 ) 2 − ∥ K ˉ H ∥ 2 \varepsilon = \|k_H(\cdot)\|^2_{L^2(T^2)} - \|\bar{K}^H\|^2 ε = ∥ k H ( ⋅ ) ∥ L 2 ( T 2 ) 2 − ∥ K ˉ H ∥ 2
حيث ∥ k H ( ⋅ ) ∥ L 2 ( T 2 ) 2 = T 2 H + 1 2 H + 1 \|k_H(\cdot)\|^2_{L^2(T^2)} = \frac{T^{2H+1}}{2H+1} ∥ k H ( ⋅ ) ∥ L 2 ( T 2 ) 2 = 2 H + 1 T 2 H + 1
الجدول 1: خطأ التقريب بمتوسط المربعات ε \varepsilon ε
H L=4 L=8 L=16 L=32 L=64 L=128 0.1 0.3842 0.3229 0.2720 0.2299 0.1950 0.1659 0.3 0.1035 0.0655 0.0423 0.0275 0.0180 0.0118 0.5 0.0357 0.0167 0.0081 0.0040 0.0020 0.0010 0.7 0.0132 0.0049 0.0019 0.0008 0.0003 0.0001 0.9 0.0812 0.0617 0.0469 0.0356 0.0270 0.0205
تحليل التقارب :بالنسبة إلى H = 0.5 H = 0.5 H = 0.5 ، يكون التقارب أسرع عندما يقترب H H H من القيم الحدية (0 أو 1)، تنخفض سرعة التقارب بشكل ملحوظ الخطأ النسبي هو ε / ( 2 H + 1 ) − 1 \varepsilon/(2H+1)^{-1} ε / ( 2 H + 1 ) − 1 مقارنة الطرق :النتائج من الطريقة الجديدة قريبة نوعياً من نتائج المرجع 2 حتى بالنسبة إلى قيم صغيرة من L L L ، يكون الخطأ من الطريقتين قريباً جداً يشير إلى أن تمثيل حاصل الضرب المصفوفي له تأثير ضئيل على الدقة، لكنه يوفر استقراراً حسابياً تحليل معدل التقارب :معدل تقارب الدالة الأسية t α t^\alpha t α هو O ( L − ( 2 α + 1 ) ) O(L^{-(2\alpha+1)}) O ( L − ( 2 α + 1 ) ) عندما α → − 1 / 2 \alpha \to -1/2 α → − 1/2 يميل معدل التقارب إلى الصفر، مما يفسر انخفاض التقارب عند مؤشرات Hurst الحدية مراجعة الطرق الكلاسيكية :تمثيل Mandelbrot-Van Ness 25 التمثيل التكاملي Decreusefond-Üstünel 15 مسح طرق المحاكاة المختلفة 4-7 علاقة هذه الورقة بالأعمال ذات الصلة :بناءً على طريقة التمثيل الطيفي السابقة للمؤلف 2 الاستفادة من الوصف الرياضي للشكل الطيفي لأنظمة التحكم 16,17 توسيع تمثيل متعددات حدود لجندر لعملية Wiener 29,30 التوسعات التطبيقية :يمكن تطبيقها على عملية Ornstein-Uhlenbeck الكسرية مناسبة لمحاكاة جسر براونيان الكسري قابلة للتوسع إلى التمثيل التكاملي العشوائي المتكرر اشتقاق ناجح للصيغ الصريحة لمعاملات توسع دالة النواة لحركة براونيان الكسرية إنشاء طريقة دقيقة لحساب خطأ التقريب بمتوسط المربعات تحقيق تمثيل تقريبي لحركة براونيان الكسرية في الوقت المستمر التحقق من فعالية الطريقة والمزايا الحسابية متطلبات الدقة الحسابية : تتطلب الطريقة الجديدة حسابات بفاصلة عائمة عالية الدقة أو حسابات رمزيةالتأثيرات الحدية : عندما يقترب مؤشر Hurst من 0 أو 1، يكون التقارب أبطأالتعقيد الحسابي : على الرغم من تبسيط حساب المعاملات، لا تزال تتطلب عمليات مصفوفيةالتوسع إلى حركة براونيان الكسرية متعددة الأبعاد التطبيق على عمليات عشوائية أكثر تعقيداً (مثل معادلات براونيان الكسرية العشوائية التفاضلية) تطوير خوارزميات عددية أكثر كفاءة الصرامة النظرية : توفير اشتقاق رياضي كامل وتحليل خطأ صارمابتكار الطريقة : تبسيط حساب معاملات التوسع، مما يحسن الجدوى العملية للطريقةالتجارب الشاملة : التحقق من فعالية الطريقة من خلال مؤشرات Hurst وأطوال قطع متعددةالقيمة التطبيقية : توفير أداة جديدة لمحاكاة حركة براونيان الكسرية في الوقت المستمرالمتطلبات الحسابية : متطلبات دقة عددية عالية قد تحد من التطبيق العمليالأداء الحدي : أداء التقارب عند مؤشرات Hurst القصوى تحتاج إلى تحسينتحليل المقارنة : المقارنة التفصيلية للأداء مع الطرق الموجودة الأخرى غير كافيةالمساهمة الأكاديمية : تقديم مساهمة مهمة لنظرية حركة براونيان الكسرية وطرق الحسابالقيمة العملية : لها إمكانيات تطبيقية في الرياضيات المالية والنمذجة الفيزيائية وغيرهاالأهمية المنهجية : توفير إطار مرجعي لتوسع متعامد لعمليات عشوائية أخرىالتطبيقات التي تتطلب محاكاة حركة براونيان الكسرية في الوقت المستمر السيناريوهات التي تتطلب تحكماً دقيقاً في خطأ التقريب الحل العددي لمعادلات براونيان الكسرية العشوائية التفاضلية تسعير المشتقات المالية وإدارة المخاطر تستشهد الورقة بـ 37 مرجعاً ذا صلة، تشمل بشكل أساسي:
النظرية الأساسية لحركة براونيان الكسرية 22-25 طرق محاكاة العمليات العشوائية 1,3,4-7 الطرق الطيفية والتوسعات المتعامدة 16,17,26 نظرية التفاضل والتكامل الكسري 32