الضوضاء المضاعفة

في معالجة الإشارات ، يشير مصطلح الضوضاء المضاعفة إلى إشارة عشوائية غير مرغوب فيها يتم ضربها في إشارة ذات صلة أثناء الالتقاط أو الإرسال أو المعالجة الأخرى.

الضوضاء المضاعفة هي نوع من الضوضاء التي تعتمد على الإشارة، حيث يتناسب حجم الضوضاء طرديًا مع شدة الإشارة. وعلى عكس الضوضاء الجمعية، التي لا تعتمد على الإشارة، فإن الضوضاء المضاعفة تُعقّد عملية المعالجة نظرًا لاعتمادها على الإشارة الأساسية.

ومن الأمثلة المهمة على ذلك ضوضاء التبقع الشائعة في صور الرادار . ومن أمثلة الضوضاء المضاعفة التي تؤثر على الصور الرقمية: الظلال الطبيعية الناتجة عن تموجات سطح الأجسام المصورة، والظلال التي تلقيها الأجسام المعقدة مثل أوراق الشجر والستائر المعدنية ، والبقع الداكنة الناتجة عن الغبار في العدسة أو مستشعر الصورة، والاختلافات في كسب العناصر الفردية لمصفوفة مستشعر الصورة . [ 1 ]

الضوضاء المضاعفة في المعادلات التفاضلية العشوائية (SDEs)

في مجال المعادلات التفاضلية العشوائية ، يُستخدم الضجيج المضاعف لنمذجة الأنظمة التي يعتمد فيها سعة التقلبات العشوائية الداخلية للنظام على حالته. ومن أبرز الأمثلة على الضجيج المضاعف في المعادلات التفاضلية العشوائية حركة براون الهندسية . تُستخدم حركة براون الهندسية على نطاق واسع في مجال التمويل لنمذجة أسعار الأسهم، وأسعار صرف العملات، وغيرها من الأصول. ويُستخدم نموذج حركة براون الهندسية على نطاق واسع في الرياضيات المالية لوصف تطور أسعار الأصول، حيث يفترض أن العوائد النسبية للأصل تتبع التوزيع الطبيعي على فترات زمنية متناهية الصغر. وتُعطى معادلة حركة براون الهندسية التفاضلية العشوائية كما يلي: دXت=μXتدت+σXتددبليوت،{\displaystyle dX_{t}=\mu X_{t}\,dt+\sigma X_{t}\,dW_{t},} أين:

  • Xت{\displaystyle X_{t}}هل سعر الأصل في ذلك الوقتت{\displaystyle t}،
  • μ{\displaystyle \mu }هو العائد المتوقع (معدل الانحراف)،
  • σ{\displaystyle \sigma }هو تقلب العوائد،
  • دبليوت{\displaystyle W_{t}}هي حركة براونية قياسية ( عملية وينر ).

نظرية (صيغة إيتو). ليكنXت{\displaystyle X_{t}}يُعطى بواسطة: دXت=ب(ت،ω)دت+σ(ت،ω)ددبليوت.{\displaystyle dX_{t}=b(t,\omega )\,dt+\sigma (t,\omega )\,dW_{t}.} يتركو(ت،x){\displaystyle f(t,x)}كنج1،2{\displaystyle C^{1,2}}دالة (أي،ج1{\displaystyle C^{1}}مع مرور الوقت،ج2{\displaystyle C^{2}}في الفضاء). ثم العمليةYت=و(Xت){\displaystyle Y_{t}=f(X_{t})}يرضي دو(Xت)=(تو(ت،Xت)+xو(ت،Xت)ب(ت،ω)+12x،xو(ت،Xت)σ2(ت،ω))دت+xو(ت،Xت)σ(ت،ω)ددبليوت.{\displaystyle df(X_{t})=\left(\partial _{t}f(t,X_{t})+\partial _{x}f(t,X_{t})b(t,\omega )+{\frac {1}{2}}\partial _{x,x}f(t,X_{t})\sigma ^{2}(t,\omega )\right)dt+\partial _{x}f(t,X_{t})\sigma (t,\omega )dW_{t}.}

تعيين و(ت،x)=سجلx.{\displaystyle f(t,x)=\log x.}

تطبيق صيغة إيتو علىYت=و(ت،x){\displaystyle Y_{t}=f(t,x)}، نقوم بالحساب: د(سجلXت)=د(Yت)=[1Xت(μXت)+12(-1Xت2)(σ2Xت2)]دت+1Xت(σXت)ددبليوت.{\displaystyle d(\log X_{t})=d(Y_{t})=\left[{\frac {1}{X_{t}}}(\mu X_{t})+{\frac {1}{2}}\left(-{\frac {1}{X_{t}^{2}}}\right)(\sigma ^{2}X_{t}^{2})\right]dt+{\frac {1}{X_{t}}}(\sigma X_{t})dW_{t}.}

تبسيط كل مصطلح: د(سجلXت)=د(Yت)=(μ-12σ2)دت+σددبليوت.{\displaystyle d(\log X_{t})=d(Y_{t})=\left(\mu -{\frac {1}{2}}\sigma ^{2}\right)dt+\sigma dW_{t}.}

وبدمج البيانات مع مرور الوقت، نحصل على: سجلXت=Yت=Y0+(μ-12σ2)ت+σدبليوت،أينY0=سجلX0.{\displaystyle \log X_{t}=Y_{t}=Y_{0}+\left(\mu -{\frac {1}{2}}\sigma ^{2}\right)t+\sigma W_{t},\quad {\text{حيث}}\quad Y_{0}=\log X_{0}.}

برفع كلا الطرفين إلى الأس نحصل على الحل لـXت{\displaystyle X_{t}}: Xت=خبرة(Yت)=X0خبرة((μ-12σ2)ت+σدبليوت).{\displaystyle X_{t}=\exp(Y_{t})=X_{0}\exp \left(\left(\mu -{\frac {1}{2}}\sigma ^{2}\right)t+\sigma W_{t}\right).}

يمكن كتابة حل هذه المعادلة التفاضلية العشوائية بشكل صريح على النحو التالي: Xت=X0خبرة((μ-12σ2)ت+σدبليوت)،{\displaystyle X_{t}=X_{0}\exp \left(\left(\mu -{\frac {1}{2}}\sigma ^{2}\right)t+\sigma W_{t}\right),} أينX0{\displaystyle X_{0}}هو سعر الأصل الأولي.

تشمل الخصائص الرئيسية لنموذج GBM ما يلي:

  • التوزيع اللوغاريتمي الطبيعي : لأي قيمة ثابتةت>0{\displaystyle t>0}،Xت{\displaystyle X_{t}}يتبع التوزيع اللوغاريتمي الطبيعي.
  • عدم السلبية :Xت>0{\displaystyle X_{t}>0}من شبه المؤكد للجميعت{\displaystyle t}، مما يضمن نمذجة واقعية لأسعار الأصول التي لا يمكن أن تصبح سلبية.
  • الضوضاء المضاعفة : مصطلح التذبذب العشوائيσXتددبليوت{\displaystyle \sigma X_{t}\,dW_{t}}يتناسب معXت{\displaystyle X_{t}}وهذا يعكس الحقيقة التجريبية المتمثلة في أن أسعار الأصول الأكبر تميل إلى إظهار تقلبات مطلقة أكبر.

يشكل نموذج GBM الأساس لنموذج Black–Scholes المستخدم لاستخلاص حلول مغلقة الشكل لتسعير الخيارات الأوروبية.

في الرياضيات المالية، يعكس وجود الضوضاء المضاعفة الملاحظة التجريبية بأن حجم تقلبات أسعار الأصول يميل إلى التناسب مع قيمة الأصل. وتُعد هذه الخاصية أساسية في اشتقاق نماذج مثل نموذج بلاك-شولز لتسعير الخيارات.

يتم وصف نموذج كوكس -إنجرسول-روس (CIR) بواسطة المعادلة التفاضلية العشوائية: درت=أ(ب-رت)دت+σرتددبليوت،{\displaystyle dr_{t}=a(b-r_{t})\,dt+\sigma {\sqrt {r_{t}}}\,dW_{t},} أين:

  • رت{\displaystyle r_{t}}هو سعر الفائدة قصير الأجل،
  • أ>0{\displaystyle a>0}هي سرعة العودة إلى المتوسط،
  • ب>0{\displaystyle b>0}هو المستوى المتوسط ​​على المدى الطويل،
  • σ>0{\displaystyle \sigma >0}هو معامل التقلب،
  • دبليوت{\displaystyle W_{t}}هي حركة براونية قياسية.

لا يوجد حل بسيط مغلق لعملية CIR من حيثرت{\displaystyle r_{t}}ودبليوت{\displaystyle W_{t}}ومع ذلك، فإن توزيعها الشرطي معروف: لقيمة ابتدائية ثابتةر0{\displaystyle r_{0}}المتغيررت{\displaystyle r_{t}}يتبع توزيع كاي تربيع غير مركزي مُقاس .

يمكن تطبيق طريقة أويلر-ماروياما للمحاكاة العددية، وذلك بتقسيم الزمن إلى خطوات زمنية.Δت{\displaystyle \Delta t}: رن+1=رن+أ(ب-رن)Δت+σرنΔدبليون،{\displaystyle r_{n+1}=r_{n}+a(b-r_{n})\Delta t+\sigma {\sqrt {r_{n}}}\,\Delta W_{n},} أينΔدبليون{\displaystyle \Delta W_{n}}هي زيادات طبيعية مستقلة مع Δدبليونشمال(0،Δت).{\displaystyle \Delta W_{n}\sim {\mathcal {N}}(0,\Delta t).}

بسبب مصطلح الانتشار ذي الجذر التربيعي، يجب توخي الحذر لضمانرن0{\displaystyle r_{n}\geq 0}أثناء المحاكاة. تُستخدم عدة طرق لمعالجة هذه المشكلة:

  • مخطط الاقتطاع الكامل : تعيين القيم السالبة إلى الصفر.
  • مخطط الانعكاس: عكس القيم السالبة إلى قيم موجبة.
  • مخطط شبه صريح:

رن+1=(رن+أ(ب-رن)Δت+σ2Δدبليون)2،{\displaystyle r_{n+1}=\left({\sqrt {r_{n}+a(b-r_{n})\Delta t}}+{\frac {\sigma }{2}}\Delta W_{n}\right)^{2},} مما يحافظ بشكل أفضل على الإيجابية ويحسن الاستقرار العددي .

بدلاً عن ذلك،رت{\displaystyle r_{t}}يمكن أخذ عينات منها بدقة عن طريق توليد متغير عشوائي من توزيع كاي تربيع غير المركزي المناسب .

يُعدّ نموذج هيستون نموذجًا للتقلبات العشوائية يُستخدم في التمويل الرياضي لوصف تطور أسعار الأصول وتقلباتها. وهو يُوسّع إطار عمل بلاك-شولز من خلال السماح للتقلبات بالتغير عشوائيًا بمرور الوقت.

يُعرَّف نموذج هيستون بواسطة نظام المعادلات التفاضلية العشوائية التالي: دSت=μSتدت+vتSتددبليوت(1)،دvت=κ(θ-vت)دت+ξvتددبليوت(2)،{\displaystyle {\begin{aligned}dS_{t}&=\mu S_{t}\,dt+{\sqrt {v_{t}}}S_{t}\,dW_{t}^{(1)},\\dv_{t}&=\kappa (\theta -v_{t})\,dt+\xi {\sqrt {v_{t}}}\,dW_{t}^{(2)},\end{aligned}}} أين:

  • Sت{\displaystyle S_{t}}هل سعر الأصل في ذلك الوقتت{\displaystyle t}،
  • vت{\displaystyle v_{t}}يمثل التباين اللحظي (أي مربع التقلب)،
  • μ{\displaystyle \mu }هو معدل انحراف الأصل،
  • κ{\displaystyle \kappa }هو المعدل الذيvت{\displaystyle v_{t}}يعود إلى متوسطه طويل الأجلθ{\displaystyle \theta }،
  • ξ{\displaystyle \xi }إنها تقلبات التقلبات،
  • دبليوت(1){\displaystyle W_{t}^{(1)}}ودبليوت(2){\displaystyle W_{t}^{(2)}}هي حركات براونية قياسية ذات ارتباطρ{\displaystyle \rho }

السمة الرئيسية لنموذج هيستون هي أن التقلبvت{\displaystyle {\sqrt {v_{t}}}}هي في حد ذاتها عملية عشوائية مدفوعة بانتشار الجذر التربيعي (على غرار عملية CIR). وهذا يسمح للنموذج برصد السمات التجريبية المهمة للأسواق المالية، مثل:

  • تجميع التقلبات،
  • تأثير الرافعة المالية (عبرρ<0{\displaystyle \rho <0})
  • تبتسم التقلبات الضمنية وتتشوه.

The Heston model admits a closed-form solution for European option prices using characteristic functions and Fourier transform methods, which makes it both tractable and flexible for calibration to market data.

General Mathematical Form of Multiplicative Noise

In general, a stochastic differential equation with multiplicative noise can be written as: dXt=f(Xt,t)dt+g(Xt,t)dWt,{\displaystyle dX_{t}=f(X_{t},t)\,dt+g(X_{t},t)\,dW_{t},} where:

  • f(Xt,t){\displaystyle f(X_{t},t)} is the drift term,
  • g(Xt,t){\displaystyle g(X_{t},t)} is the diffusion coefficient,
  • Wt{\displaystyle W_{t}} is a standard Brownian motion.

When the diffusion coefficient g{\displaystyle g} depends explicitly on the state variable Xt{\displaystyle X_{t}}, the noise is said to be multiplicative. This contrasts with additive noise, where g{\displaystyle g} is independent of Xt{\displaystyle X_{t}}. Multiplicative noise introduces complexities in both analytical and numerical treatments of SDEs, including the need to carefully choose between interpretations such as the Itô calculus and the Stratonovich calculus.

In particular, under the Itô interpretation, the presence of state-dependent noise can induce additional drift terms when transforming variables, a phenomenon known as the Itô correction.

References

  1. Maria Petrou, Costas Petrou (2010) Image Processing: The Fundamentals. John Wiley & Sons. 818 pages. ISBN 9780470745861