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.
পত্রের ID : 2503.04575শিরোনাম : ভগ্নাংশ ব্রাউনীয় গতির আনুমানিক প্রতিনিধিত্ব সম্পর্কেলেখক : কনস্টান্টিন এ. রাইবাকভশ্রেণীবিভাগ : math.PR (সম্ভাব্যতা তত্ত্ব)প্রকাশনার সময় : ২০২৫ সালের মার্চ (arXiv প্রাক-প্রিন্ট, সংস্করণ v2 ২০২৫ সালের অক্টোবর ১৪ তারিখে)পত্রের লিঙ্ক : https://arxiv.org/abs/2503.04575 এই পত্রটি লিজেন্ড্র বহুপদের সাপেক্ষে ভগ্নাংশ ব্রাউনীয় গতির অর্থোগোনাল সম্প্রসারণ অধ্যয়ন করে। এই সম্প্রসারণটি শুধুমাত্র তাত্ত্বিক গুরুত্বই নয়, বরং ব্যবহারিক মূল্যও রাখে, কারণ এটি ক্রমাগত সময়ে ভগ্নাংশ ব্রাউনীয় গতির আনুমানিক এবং অনুকরণের জন্য ব্যবহার করা যায়। নিবন্ধটি গড় বর্গ আনুমানিক ত্রুটির সম্পর্ক প্রদান করে এবং পূর্বে প্রাপ্ত ফলাফলের সাথে তুলনা করে।
মূল সমস্যা : ক্রমাগত সময়ে ভগ্নাংশ ব্রাউনীয় গতি (FBM) কীভাবে কার্যকরভাবে আনুমানিক এবং অনুকরণ করা যায়গুরুত্ব : ভগ্নাংশ ব্রাউনীয় গতি আর্থিক মডেলিং, পদার্থবিজ্ঞান, প্রকৌশল এবং অন্যান্য অনেক ক্ষেত্রে ব্যাপক প্রয়োগ রয়েছে, এর অনুকরণ তাত্ত্বিক ফলাফল যাচাই এবং ব্যবহারিক সমস্যা সমাধানের জন্য অপরিহার্যবিদ্যমান পদ্ধতির সীমাবদ্ধতা :
ঐতিহ্যবাহী বিচ্ছিন্ন সময় জালিকা পদ্ধতি শুধুমাত্র জালিকা নোডে নির্ভুলতা নিশ্চিত করতে পারে বিদ্যমান সাধারণ গাউসীয় প্রক্রিয়া অনুকরণ পদ্ধতি উচ্চ গণনামূলক জটিলতা রয়েছে আনুমানিক ত্রুটি সঠিকভাবে গণনা করতে পারে এমন ক্রমাগত সময় পদ্ধতির অভাব রয়েছে গবেষণা প্রেরণা : লেখকের পূর্ববর্তী কাজ 2 এর উপর ভিত্তি করে, আরও সহজ সম্প্রসারণ সহগ গণনা পদ্ধতি প্রস্তাব করা, গণনামূলক জটিলতা হ্রাস এবং নির্ভুলতা উন্নত করাতাত্ত্বিক অবদান : ভগ্নাংশ ব্রাউনীয় গতির অবিচ্ছেদ্য প্রতিনিধিত্বের সাথে সম্পর্কিত কার্নেল ফাংশন সম্প্রসারণ সহগের স্পষ্ট সম্পর্ক উদ্ভাবন করা (লিজেন্ড্র বহুপদের সাপেক্ষে)পদ্ধতি উন্নতি : পূর্ববর্তী কাজ 2 এর চেয়ে সহজ সম্প্রসারণ সহগ গণনা পদ্ধতি প্রদান করা, গণনামূলক জটিলতা হ্রাস করাত্রুটি বিশ্লেষণ : ভগ্নাংশ ব্রাউনীয় গতির বহুপদী আনুমানিক গড় বর্গ ত্রুটি সঠিকভাবে গণনা করতে পারে এমন সমীকরণ স্থাপন করাক্রমাগত সময় আনুমানিক : ভগ্নাংশ ব্রাউনীয় গতির ক্রমাগত সময়ে আনুমানিক প্রতিনিধিত্ব বাস্তবায়ন করা, শুধুমাত্র বিচ্ছিন্ন সময় বিন্দুতে নয়প্রদত্ত ভগ্নাংশ ব্রাউনীয় গতি B H ( t ) B_H(t) B H ( t ) (হার্স্ট সূচক 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 (সম্প্রসারণ সহগের স্পষ্ট অভিব্যক্তি):
ম্যাট্রিক্স K H K^H K H এর j j j তম স্তম্ভের জন্য:
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) ব্যবহার করে শক্তি ফাংশনের সম্প্রসারণ সহগ গণনা করাম্যাট্রিক্স নির্মাণ : উপপাদ্য 1 এর মাধ্যমে সম্প্রসারণ সহগ ম্যাট্রিক্স K H K^H K H নির্মাণ করার্যান্ডম ভেরিয়েবল উৎপাদন : স্বাধীন মান স্বাভাবিক র্যান্ডম ভেরিয়েবল 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 ] হার্স্ট সূচক : 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 সংমিশ্রণ হার শূন্যের দিকে প্রবণ, যা সীমানা হার্স্ট সূচকের কম সংমিশ্রণ ব্যাখ্যা করে ক্লাসিক পদ্ধতি পর্যালোচনা :ম্যান্ডেলব্রট-ভ্যান নেস প্রতিনিধিত্ব 25 ডিক্রিউসফন্ড-উস্টুনেল অবিচ্ছেদ্য প্রতিনিধিত্ব 15 বিভিন্ন অনুকরণ পদ্ধতির সমীক্ষা 4-7 এই পত্রের সাথে সম্পর্কিত কাজের সম্পর্ক :লেখকের পূর্ববর্তী বর্ণালী প্রতিনিধিত্ব পদ্ধতির উপর ভিত্তি করে 2 নিয়ন্ত্রণ সিস্টেম বর্ণালী ফর্ম গাণিতিক বর্ণনা থেকে অনুপ্রেরণা 16,17 উইনার প্রক্রিয়ার লিজেন্ড্র বহুপদ প্রতিনিধিত্ব সম্প্রসারিত করা 29,30 প্রয়োগ সম্প্রসারণ :ভগ্নাংশ অর্নস্টাইন-উহলেনবেক প্রক্রিয়ার জন্য ব্যবহারযোগ্য ভগ্নাংশ ব্রাউনীয় সেতু অনুকরণের জন্য উপযুক্ত পুনরাবৃত্তিমূলক র্যান্ডম অবিচ্ছেদ্য প্রতিনিধিত্বে সম্প্রসারণযোগ্য ভগ্নাংশ ব্রাউনীয় গতির কার্নেল ফাংশন সম্প্রসারণ সহগের স্পষ্ট সূত্র সফলভাবে উদ্ভাবন করা নির্ভুল গড় বর্গ আনুমানিক ত্রুটি গণনা পদ্ধতি স্থাপন করা ক্রমাগত সময়ে ভগ্নাংশ ব্রাউনীয় গতির আনুমানিক প্রতিনিধিত্ব বাস্তবায়ন করা পদ্ধতির কার্যকারিতা এবং গণনামূলক সুবিধা যাচাই করা গণনা নির্ভুলতা প্রয়োজনীয়তা : নতুন পদ্ধতি উচ্চ নির্ভুলতা ফ্লোটিং পয়েন্ট গণনা বা প্রতীকী গণনা প্রয়োজনসীমানা প্রভাব : যখন হার্স্ট সূচক 0 বা 1 এর কাছাকাছি থাকে সংমিশ্রণ ধীরগণনা জটিলতা : যদিও সহগ গণনা সরলীকৃত হয়েছে, তবুও ম্যাট্রিক্স অপারেশন প্রয়োজনবহুমাত্রিক ভগ্নাংশ ব্রাউনীয় গতিতে সম্প্রসারণ আরও জটিল র্যান্ডম প্রক্রিয়ায় প্রয়োগ (যেমন ভগ্নাংশ র্যান্ডম ডিফারেনশিয়াল সমীকরণ) আরও দক্ষ সংখ্যাসূচক অ্যালগরিদম উন্নয়ন তাত্ত্বিক কঠোরতা : সম্পূর্ণ গাণিতিক উদ্ভাবন এবং কঠোর ত্রুটি বিশ্লেষণ প্রদান করাপদ্ধতি উদ্ভাবন : সম্প্রসারণ সহগ গণনা সরলীকৃত করা, পদ্ধতির ব্যবহারিকতা উন্নত করাপর্যাপ্ত পরীক্ষা : একাধিক হার্স্ট সূচক এবং ছাঁটাই দৈর্ঘ্যের মাধ্যমে পদ্ধতির কার্যকারিতা যাচাই করাপ্রয়োগ মূল্য : ভগ্নাংশ ব্রাউনীয় গতির ক্রমাগত সময় অনুকরণের জন্য নতুন সরঞ্জাম প্রদান করাগণনা প্রয়োজনীয়তা : সংখ্যাসূচক নির্ভুলতার প্রতি উচ্চ প্রয়োজনীয়তা, ব্যবহারিক প্রয়োগ সীমাবদ্ধ করতে পারেসীমানা কর্মক্ষমতা : চরম হার্স্ট সূচকের অধীনে সংমিশ্রণ কর্মক্ষমতা উন্নত করা প্রয়োজনতুলনা বিশ্লেষণ : অন্যান্য বিদ্যমান পদ্ধতির সাথে বিস্তারিত কর্মক্ষমতা তুলনা অপর্যাপ্তএকাডেমিক অবদান : ভগ্নাংশ ব্রাউনীয় গতি তত্ত্ব এবং গণনা পদ্ধতিতে গুরুত্বপূর্ণ অবদানব্যবহারিক মূল্য : আর্থিক গণিত, পদার্থবিজ্ঞান মডেলিং ইত্যাদি ক্ষেত্রে প্রয়োগ সম্ভাবনা রয়েছেপদ্ধতিগত তাৎপর্য : অন্যান্য র্যান্ডম প্রক্রিয়ার অর্থোগোনাল সম্প্রসারণের জন্য রেফারেন্স ফ্রেমওয়ার্ক প্রদান করাক্রমাগত সময় ভগ্নাংশ ব্রাউনীয় গতি অনুকরণ প্রয়োজনীয় প্রয়োগ আনুমানিক ত্রুটির উপর নির্ভুল নিয়ন্ত্রণ প্রয়োজনীয় পরিস্থিতি ভগ্নাংশ র্যান্ডম ডিফারেনশিয়াল সমীকরণের সংখ্যাসূচক সমাধান আর্থিক ডেরিভেটিভ মূল্য নির্ধারণ এবং ঝুঁকি ব্যবস্থাপনা এই পত্রটি 37টি সম্পর্কিত সংদর্ভ উদ্ধৃত করে, প্রধানত অন্তর্ভুক্ত:
ভগ্নাংশ ব্রাউনীয় গতির মৌলিক তত্ত্ব 22-25 র্যান্ডম প্রক্রিয়া অনুকরণ পদ্ধতি 1,3,4-7 বর্ণালী পদ্ধতি এবং অর্থোগোনাল সম্প্রসারণ 16,17,26 ভগ্নাংশ ক্যালকুলাস তত্ত্ব 32