الحساب الفاصل

دالة التسامح ( باللون الفيروزي ) والتقريب ذو القيم الفاصلية ( باللون الأحمر )

الحساب الفتري (المعروف أيضًا باسم رياضيات الفترات ، أو تحليل الفترات ، أو حساب الفترات ) هو أسلوب رياضي يُستخدم للحد من أخطاء التقريب والقياس في العمليات الحسابية ، وذلك بحساب حدود الدوال . تضمن الطرق العددية التي تتضمن الحساب الفتري نتائج موثوقة وصحيحة رياضيًا. فبدلًا من تمثيل القيمة برقم واحد، يُمثل الحساب الفتري كل قيمة بنطاق من الاحتمالات .

رياضيًا، بدلًا من التعامل مع متغير حقيقي غير محدد x ، يتعامل حساب الفترات مع الفترة [ a , b ] التي تحدد نطاق القيم التي يمكن أن يأخذها x . بعبارة أخرى، تقع أي قيمة للمتغير x في الفترة المغلقة بين a و b . عند تطبيق دالة f على x ، ينتج عنها الفترة [ c , d ] التي تشمل جميع القيم الممكنة لـ f ( x ) لكل x ∈ [ a , b ] .

يُعدّ حساب الفترات مناسبًا لأغراضٍ متنوعة؛ وأكثرها شيوعًا في الأعمال العلمية، لا سيما عند إجراء الحسابات بواسطة برامج حاسوبية، حيث يُستخدم لتتبع أخطاء التقريب في الحسابات وعدم اليقين في معرفة القيم الدقيقة للمعايير الفيزيائية والتقنية. غالبًا ما ينشأ هذا الأخير من أخطاء القياس والتفاوتات المسموح بها للمكونات أو بسبب قيود دقة الحساب. كما يُساعد حساب الفترات في إيجاد حلول مضمونة للمعادلات (مثل المعادلات التفاضلية ) ومسائل التحسين .

مقدمة

الهدف الرئيسي من حساب الفترات هو توفير طريقة بسيطة لحساب الحدين الأعلى والأدنى لمدى دالة ما في متغير واحد أو أكثر. لا تُمثل هذه الحدود بالضرورة القيم العليا أو الدنيا الحقيقية للمدى، إذ قد يكون حسابها بدقة أمرًا صعبًا أو مستحيلاً؛ ويكفي أن تحتوي الحدود على مدى الدالة كمجموعة جزئية.

يقتصر هذا العلاج عادةً على فترات زمنية حقيقية، لذا فإن الكميات تأتي على شكل

[أ،ب]={xR|أxب}،{\displaystyle [a,b]=\{x\in \mathbb {R} \mid a\leq x\leq b\},}

حيث يُسمح بقيمتي a = −∞ و b = ∞ . إذا كانت إحدى القيمتين a أو b لانهائية، فإن الفترة ستكون فترة غير محدودة؛ أما إذا كانت كلتاهما لانهائية، فإن الفترة ستكون امتدادًا لخط الأعداد الحقيقية. وبما أن العدد الحقيقي r يُمكن تفسيره على أنه الفترة [ r , r ] ، فإنه يُمكن الجمع بين الفترات والأعداد الحقيقية بحرية.

مثال

مؤشر كتلة الجسم لشخص  طوله 1.80 متر بالنسبة لوزن الجسم ( بالكيلوغرام)

لنفترض حساب مؤشر كتلة الجسم (BMI). يُحسب مؤشر كتلة الجسم بقسمة وزن الجسم بالكيلوغرام على مربع طوله بالمتر. لنفترض أن شخصًا ما يستخدم ميزانًا بدقة كيلوغرام واحد، حيث لا يمكن تمييز القيم المتوسطة، ويتم تقريب الوزن الحقيقي إلى أقرب عدد صحيح. على سبيل المثال، لا يمكن التمييز بين 79.6  كجم و80.3  كجم، لأن الميزان لا يعرض إلا القيم لأقرب كيلوغرام. من غير المرجح أن  يكون وزن الشخص 80.0 كجم بالضبط  عندما يُظهر الميزان 80 كجم. بالتالي، يشير الميزان الذي يُظهر 80  كجم إلى وزن يقع بين 79.5  كجم و80.5  كجم، أو ضمن النطاق [ 79.5، 80.5] .

يبلغ مؤشر كتلة الجسم لرجل يزن 80  كجم وطوله 1.80  متر حوالي 24.7. أما وزن 79.5  كجم ونفس الطول فيعطيان مؤشر كتلة جسم قدره 24.537، بينما وزن 80.5  كجم يعطي مؤشر كتلة جسم قدره 24.846. وبما أن كتلة الجسم متغيرة باستمرار وتزداد دائمًا لجميع القيم ضمن نطاق الوزن المحدد، فإن مؤشر كتلة الجسم الحقيقي يقع ضمن النطاق [ 24.537، 24.846 ] . ولأن النطاق بأكمله أقل من 25، وهو الحد الفاصل بين الوزن الطبيعي والوزن الزائد، يمكن الاستنتاج بيقين أن وزن الرجل طبيعي.

لا يؤثر الخطأ في هذا المثال على النتيجة (الوزن الطبيعي)، ولكن هذا ليس صحيحًا بشكل عام. فلو كان الرجل أثقل وزنًا بقليل، فقد يشمل نطاق مؤشر كتلة الجسم القيمة الحدية 25. في هذه الحالة، لن تكون دقة الميزان كافية للوصول إلى نتيجة قاطعة.

يمكن الإبلاغ عن نطاق أمثلة مؤشر كتلة الجسم على أنه [ 24.5، 24.9 ] لأن هذه الفترة هي مجموعة فرعية من الفترة المحسوبة. لكن لا يمكن الإبلاغ عن النطاق على أنه [ 24.6، 24.8 ] ، لأن هذه الفترة لا تحتوي على قيم محتملة لمؤشر كتلة الجسم.

فترات متعددة

مؤشر كتلة الجسم للأوزان المختلفة بالنسبة للطول L (بالمتر)

يؤثر كل من الطول والوزن على قيمة مؤشر كتلة الجسم (BMI). مع أن المثال السابق اقتصر على دراسة تغير الوزن، إلا أن الطول أيضاً عرضة للخطأ. تُقرّب قياسات الطول بالمتر عادةً إلى أقرب سنتيمتر: فقياس 1.79 متر يُمثل طولاً ضمن النطاق [ 1.785، 1.795] . بما أن مؤشر كتلة الجسم يزداد مع زيادة الوزن وينخفض ​​مع زيادة الطول، يُمكن حساب هامش الخطأ باستبدال أدنى وأعلى قيم كل نطاق، ثم اختيار أدنى وأعلى القيمتين كحدود. بالتالي، يجب أن يقع مؤشر كتلة الجسم ضمن هذا النطاق.

[79.5،80.5)[1.785،1.795)2[24.673،25.266].{\displaystyle {\frac {[79.5,80.5)}{[1.785,1.795)^{2}}}\subseteq [24.673,25.266].}

في هذه الحالة، قد يكون وزن الرجل طبيعياً أو يعاني من زيادة الوزن؛ لم تكن قياسات الوزن والطول دقيقة بما يكفي للتوصل إلى استنتاج نهائي.

عوامل الفترات

تُعرَّف العملية الثنائية على فترتين، مثل الجمع أو الضرب، بواسطة

[x1،x2][y1،y2]={xy|x[x1،x2]y[y1،y2]}.{\displaystyle [x_{1},x_{2}]{\,\star \,}[y_{1},y_{2}]=\{x\star y\,|\,x\in [x_{1},x_{2}]\,\land \,y\in [y_{1},y_{2}]\}.}

بمعنى آخر، هي مجموعة جميع القيم الممكنة لـ xy ، حيث x و y تقعان في فتراتهما المناظرة. إذا كانت ∗ دالة رتيبة لكل معامل في هذه الفترات، كما هو الحال في العمليات الحسابية الأساسية الأربع (باستثناء القسمة عندما يكون المقام صفرًا)، فإن القيم القصوى تقع عند نهايتي فترات المعاملات. إحدى طرق التعبير عن ذلك هي كتابة جميع التوليفات الممكنة.

[x1،x2][y1،y2]=[مين{x1y1،x1y2،x2y1،x2y2}،الأعلى{x1y1،x1y2،x2y1،x2y2}]،{\displaystyle [x_{1},x_{2}]\star [y_{1},y_{2}]=\left[\min\{x_{1}\star y_{1},x_{1}\star y_{2},x_{2}\star y_{1},x_{2}\star y_{2}\},\max\{x_{1}\star y_{1},x_{1}\star y_{2},x_{2}\star y_{1},x_{2}\star y_{2}\}\right],}

بشرط أن يكون xy معرفًا لجميع x[ x 1 , x 2 ] و y[ y 1 , y 2 ] .

أما بالنسبة للتطبيقات العملية، فيمكن تبسيط ذلك أكثر:

  • إضافة :[x1،x2]+[y1،y2]=[x1+y1،x2+y2]{\displaystyle [x_{1},x_{2}]+[y_{1},y_{2}]=[x_{1}+y_{1},x_{2}+y_{2}]}
  • الطرح :[x1،x2]-[y1،y2]=[x1-y2،x2-y1]{\displaystyle [x_{1},x_{2}]-[y_{1},y_{2}]=[x_{1}-y_{2},x_{2}-y_{1}]}
  • الضرب :[x1،x2][y1،y2]=[مين{x1y1،x1y2،x2y1،x2y2}،الأعلى{x1y1،x1y2،x2y1،x2y2}]{\displaystyle [x_{1},x_{2}]\cdot [y_{1},y_{2}]=[\min\{x_{1}y_{1},x_{1}y_{2},x_{2}y_{1},x_{2}y_{2}\},\max\{x_{1}y_{1},x_{1}y_{2},x_{2}y_{1},x_{2}y_{2}\}]}
  • قسم :[x1،x2][y1،y2]=[x1،x2]1[y1،y2]،{\displaystyle {\frac {[x_{1},x_{2}]}{[y_{1},y_{2}]}}=[x_{1},x_{2}]\cdot {\frac {1}{[y_{1},y_{2}]}},}أين1[y1،y2]=[1y2،1y1]لو0[y1،y2]1[y1،0]=[-،1y1]1[0،y2]=[1y2،]1[y1،y2]=[-،1y1][1y2،][-،]لو0(y1،y2){\displaystyle {\begin{aligned}{\frac {1}{[y_{1},y_{2}]}}&=\left[{\frac {1}{y_{2}}},{\frac {1}{y_{1}}}\right]&&{\textrm {if}}\;0\notin [y_{1},y_{2}]\\{\frac {1}{[y_{1},0]}}&=\left[-\infty ,{\frac {1}{y_{1}}}\right]\\{\frac {1}{[0,y_{2}]}}&=\left[{\frac {1}{y_{2}}},\infty \right]\\{\frac {1}{[y_{1},y_{2}]}}&=\left[-\infty ,{\frac {1}{y_{1}}}\right]\cup \left[{\frac {1}{y_{2}}},\infty \right]\subseteq [-\infty ,\infty ]&&{\textrm {if}}\;0\in (y_{1},y_{2})\end{aligned}}}

في الحالة الأخيرة ، تُفقد معلومات مفيدة حول استبعاد الفترة ( 1 / y1 , 1 / y2 ) . لذا، من الشائع التعامل مع الفترتين [-∞, 1/y1] و[1 /y2, ∞ ] كفترات منفصلة . وبشكل أعم ، عند التعامل مع الدوال غير المتصلة ، يكون من المفيد أحيانًا إجراء الحساب باستخدام ما يُسمى بالفترات المتعددة من الشكل ⋃i [ ai , bi ] . تحافظ العمليات الحسابية للفترات المتعددة على مجموعة من الفترات (عادةً ما تكون منفصلة) وتتيح أيضًا دمج الفترات المتداخلة. [ 1 ]

ضرب الفترات الموجبة

غالبًا ما يتطلب ضرب الفترات عمليتي ضرب فقط. إذا كانت x1 و y1 غير سالبتين ،

[x1،x2][y1،y2]=[x1y1،x2y2]، لو x1،y10.{\displaystyle [x_{1},x_{2}]\cdot [y_{1},y_{2}]=[x_{1}\cdot y_{1},x_{2}\cdot y_{2}],\qquad {\text{ if }}x_{1},y_{1}\geq 0.}

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

بمساعدة هذه التعريفات، أصبح من الممكن بالفعل حساب مدى الدوال البسيطة، مثل f ( a , b , x ) = ax + b . على سبيل المثال، إذا كانت a = [ 1, 2 ] و b = [ 5, 7 ] و x = [ 2, 3 ] :

و(أ،ب،x)=[1،2][2،3]+[5،7]=[12،23]+[5،7]=[7،13].{\displaystyle f(a,b,x)=[1,2]\cdot [2,3]+[5,7]=[1\cdot 2,2\cdot 3]+[5,7]=[7,13].}

الترميز

يمكن استخدام الأقواس لتقصير تدوين الفترات.

يمكن استخدام [ x ][ x1 , x2 ] لتمثيل فترة. لاحظ أنه في هذا الترميز المختصر، يجب عدم الخلط بين [ x ] وفترة النقطة الواحدة [ x1 , x2 ] والفترة العامة. بالنسبة لمجموعة جميع الفترات، يمكننا استخدام

[R]:={[x1،x2]|x1x2 و x1،x2R{-،}}{\displaystyle [\mathbb {R} ]:=\left\{\,[x_{1},x_{2}]\,|\,x_{1}\leq x_{2}{\text{ and }}x_{1},x_{2}\in \mathbb {R} \cup \{-\infty ,\infty \}\right\}}

كاختصار. بالنسبة لمتجه الفترات، يمكننا استخدام خط غامق:

[x]=([x]1،...،[x]ن)[R]ن.{\displaystyle [\mathbf {x} ]=\left([x]_{1},\ldots ,[x]_{n}\right)\in [\mathbb {R} ]^{n}.}

الدوال الأولية

قيم دالة رتيبة

يمكن أيضاً تعريف دوال الفترات التي تتجاوز العوامل الأساسية الأربعة.

بالنسبة للدوال الرتيبة في متغير واحد، يكون حساب نطاق القيم بسيطًا. إذا كانت الدالة f  : ℝ → ℝ متزايدة بشكل رتيب في الفترة [ x 1 , x 2 ] ، فإنه لكل y 1 , y 2[ x 1 , x 2 ] بحيث y 1 < y 2 ، فإن f ( y 1 ) ≤ f ( y 2 ) (وإذا كانت f متناقصة، فإن f ( y 1 ) ≥ f ( y 2 ) ).

وبالتالي، يمكن حساب المدى المقابل للفترة [ y 1 , y 2 ][ x 1 , x 2 ] عن طريق تطبيق الدالة على نقاط نهايتها:

و([y1،y2])=[مين{و(y1)،و(y2)}،الأعلى{و(y1)،و(y2)}].{\displaystyle f([y_{1},y_{2}])=\left[\min \left\{f(y_{1}),f(y_{2})\right\},\max \left\{f(y_{1}),f(y_{2})\right\}\right].}

ومن هذا، يمكن بسهولة تحديد الخصائص الأساسية التالية لدوال الفترات:

  • الدالة الأسية :أ[x1،x2]=[أx1،أx2]ل أ>1،{\displaystyle a^{[x_{1},x_{2}]}=[a^{x_{1}},a^{x_{2}}]\qquad {\text{for }}a>1,}
  • اللوغاريتم :سجلأ[x1،x2]=[سجلأx1،سجلأx2]للفترات الموجبة [x1،x2] و أ>1،{\displaystyle \log _{a}[x_{1},x_{2}]=[\log _{a}{x_{1}},\log _{a}{x_{2}}]\quad {\text{for positive intervals }}[x_{1},x_{2}]{\text{ and }}a>1,}
  • قوى غريبة:[x1،x2]ن=[x1ن،x2ن]للفردي نشمال.{\displaystyle [x_{1},x_{2}]^{n}=[x_{1}^{n},x_{2}^{n}]\qquad {\text{for odd }}n\in \mathbb {N} .}

بالنسبة للقوى الزوجية، يُعد نطاق القيم المُراد دراستها مهمًا ويجب تحديده قبل إجراء أي عملية ضرب. على سبيل المثال، xⁿ حيث x [ −1, 1 ] يجب أن يُنتج الفترة [ 0, 1 ] لأي عدد زوجي n ∈ ℕ . ولكن إذا تم الحصول على [ −1, 1 ] بتكرار ضرب الفترات بالشكل [ −1, 1 ] · [ −1, 1 ] · ⋯ · [ −1, 1 فإن النتيجة ستكون [ −1, 1 ] ، وهي أوسع من اللازم.

بشكل عام ، يمكن القول إنه بالنسبة للدوال الرتيبة جزئيًا، يكفي النظر إلى نقطتي النهاية x₁ و x₂ لفترة ما، بالإضافة إلى ما يُسمى بالنقاط الحرجة داخل هذه الفترة، وهي النقاط التي يتغير عندها اتجاه رتابة الدالة. بالنسبة لدالتي الجيب وجيب التمام ، تقع النقاط الحرجة عند ( n + 1/2 ) π و nπ على التوالي ، حيث n ∈ ℤ . وبالتالي، يكفي النظر إلى خمس نقاط فقط داخل الفترة، لأن الفترة الناتجة هي [ -1 , 1 ] إذا كانت تحتوي على قيمتين قصوى على الأقل. بالنسبة لدالتي الجيب وجيب التمام، لا يلزم سوى حساب نقاط النهاية بالكامل، لأن النقاط الحرجة تؤدي إلى قيم يمكن حسابها مسبقًا بسهولة، وهي -1 و 0 و 1.

امتدادات الفترات للدوال العامة

بشكل عام، قد لا يكون من السهل إيجاد وصف بسيط لفترة الإخراج للعديد من الدوال. ولكن قد يظل من الممكن توسيع الدوال لتشمل حساب الفترات. إذا كانت f  : ℝ n → ℝ دالة من متجه حقيقي إلى عدد حقيقي، فإن [ f ]  : [ℝ] n → [ℝ] تُسمى امتدادًا لفترة f إذا

[و]([x]){و(y)|y[x]}.{\displaystyle [f]([\mathbf {x} ])\supseteq \{f(\mathbf {y} )\mid \mathbf {y} \in [\mathbf {x} ]\}.}

لا يُعطي هذا التعريف لتمديد الفترة نتيجة دقيقة. على سبيل المثال، يُعد كل من [f]([x1, x2 ] ) = [ e ^ x1 , e ^ x2 ] و [ g ] ( [ x1 , x2 ] ) = [ −∞ , ] تمديدين مقبولين للدالة الأسية. يُفضّل استخدام تمديدات أدق، مع مراعاة التكاليف النسبية للحساب وعدم الدقة؛ في هذه الحالة، يُنصح باختيار [ f ] لأنه يُعطي أدق نتيجة ممكنة.

بالنظر إلى تعبير حقيقي، يتم تحقيق امتداده الطبيعي للفترة باستخدام امتدادات الفترة لكل من تعبيراته الفرعية ووظائفه ومعاملاته.

امتداد فترة تايلور ( من الدرجة k ) هو دالة قابلة للتفاضل k + 1 مرة f معرفة بـ

[و]([x]):=و(y)+أنا=1ك1أنا!دأناو(y)([x]-y)أنا+[ر]([x]،[x]،y)،{\displaystyle [f]([\mathbf {x} ]):=f(\mathbf {y} )+\sum _{i=1}^{k}{\frac {1}{i!}}\mathrm {D} ^{i}f(\mathbf {y} )\cdot ([\mathbf {x} ]-\mathbf {y} )^{i}+[r]([\mathbf {x} ],[\mathbf {x} ],\mathbf {y} ),}

لبعض y ∈ [ x ] ، حيث D i f ( y ) هو التفاضل من الرتبة i لـ f عند النقطة y و [ r ] هو امتداد فترة لباقي تايلور .

ر(x،ξ،y)=1(ك+1)!دك+1و(ξ)(x-y)ك+1.{\displaystyle r(\mathbf {x} ,{\boldsymbol {\xi }},\mathbf {y} )={\frac {1}{(k+1)!}}\mathrm {D} ^{k+1}f({\boldsymbol {\xi }})\cdot (\mathbf {x} -\mathbf {y} )^{k+1}.}
صيغة القيمة المتوسطة

يقع المتجه ξ بين x و y حيث x و y ∈ [ x ] ؛ ξ محمي بواسطة [ x ] . عادةً ما يتم اختيار y ليكون نقطة منتصف الفترة، ويتم استخدام امتداد الفترة الطبيعي لتقييم الباقي.

تُعرف الحالة الخاصة لتمديد فترة تايلور من الدرجة k = 0 أيضًا باسم شكل القيمة المتوسطة .

الحساب الفترات المركبة

يمكن تعريف الفترة بأنها مجموعة من النقاط ضمن مسافة محددة من المركز، ويمكن توسيع هذا التعريف من الأعداد الحقيقية إلى الأعداد المركبة . [ 2 ] كما يُعرّف امتداد آخر الفترات على أنها مستطيلات في المستوى المركب. وكما هو الحال في الحساب بالأعداد الحقيقية، فإن الحساب بالأعداد المركبة ينطوي على بيانات غير مؤكدة. لذا، ونظرًا لأن عدد الفترة هو فترة حقيقية مغلقة، والعدد المركب هو زوج مرتب من الأعداد الحقيقية ، فلا يوجد سبب لحصر تطبيق حساب الفترات على قياس حالات عدم اليقين في العمليات الحسابية بالأعداد الحقيقية. [ 3 ] وبالتالي، يمكن توسيع حساب الفترات، عبر أعداد الفترات المركبة، لتحديد مناطق عدم اليقين في الحساب بالأعداد المركبة. ويمكن تعريف حساب الفترات المركبة إما باستخدام المستطيلات أو الأقراص، ولكل منهما مزاياه وعيوبه. [ 3 ]

يمكن تعميم العمليات الجبرية الأساسية للأعداد الحقيقية (الفترات المغلقة الحقيقية) على الأعداد المركبة. لذا، ليس من المستغرب أن يكون حساب الفترات المركبة مشابهًا لحساب الأعداد المركبة العادي، ولكنه ليس مطابقًا له تمامًا. [ 3 ] يمكن إثبات أنه، كما هو الحال في حساب الفترات الحقيقية، لا توجد خاصية التوزيع بين جمع وضرب أعداد الفترات المركبة إلا في حالات خاصة معينة، كما أن العناصر المعكوسة لا توجد دائمًا لأعداد الفترات المركبة. [ 3 ] هناك خاصيتان مفيدتان أخريان لحساب الأعداد المركبة العادي لا تنطبقان على حساب الفترات المركبة: خصائص الجمع والضرب، الخاصة بالمرافقات المركبة العادية، لا تنطبق على مرافقات الفترات المركبة. [ 3 ]

يمكن توسيع نطاق الحساب الفتري، بطريقة مماثلة، ليشمل أنظمة عددية متعددة الأبعاد أخرى مثل الكواترنيونات والأوكتونيونات ، ولكن على حساب التضحية بخصائص أخرى مفيدة للحساب العادي. [ 3 ]

طرق الفترات

لا يمكن نقل أساليب التحليل العددي الكلاسيكي بشكل مباشر إلى خوارزميات ذات قيم فترية، حيث لا يتم عادةً أخذ التبعيات بين القيم العددية في الاعتبار.

حساب الفترات المقربة

الحدود الخارجية عند مستويات تقريب مختلفة

لضمان فعالية التطبيق العملي، يجب أن تكون الفترات الزمنية متوافقة مع حسابات الفاصلة العائمة. كانت العمليات السابقة تعتمد على الحساب الدقيق، ولكن بشكل عام، قد لا تتوفر طرق حل عددية سريعة لها. في هذا المثال، تتراوح قيم الدالة f ( x , y ) = x + y عندما x[ 0.1, 0.8 ] و y[ 0.06, 0.08 ] بين [ 0.16, 0.88 ] . أما عند إجراء الحساب نفسه بدقة رقم واحد، فستكون النتيجة عادةً [ 0.2, 0.9 ] . لكن [ 0.2, 0.9 ][ 0.16, 0.88 ] ، لذا فإن هذا الأسلوب يتعارض مع المبادئ الأساسية لحساب الفترات، حيث سيُفقد جزء من مجال الدالة f ( [ 0.1, 0.8 ] , [ 0.06, 0.08 ] ) . لذا، يُستخدم الحل المُقَرَّب للخارج [ 0.1, 0.9 ] .

يحدد معيار IEEE 754 للحسابات الثنائية ذات الفاصلة العائمة إجراءات لتنفيذ التقريب. يسمح النظام المتوافق مع معيار IEEE 754 للمبرمجين بالتقريب إلى أقرب عدد ذي فاصلة عائمة؛ وتشمل البدائل التقريب نحو الصفر (التقريب إلى أقرب عدد صحيح)، أو التقريب نحو اللانهاية الموجبة (التقريب إلى الأعلى)، أو التقريب نحو اللانهاية السالبة (التقريب إلى الأسفل).

يمكن تحقيق التقريب الخارجي المطلوب لحسابات الفترات بتغيير إعدادات التقريب في المعالج عند حساب الحد الأعلى (up) والحد الأدنى (down). أو بدلاً من ذلك، يمكن إضافة فترة صغيرة مناسبة [ ε1 , ε2 ] .

تم أيضًا تطبيق الحساب الفاصل باستخدام حسابات الفاصلة العائمة ذات الدقة العشوائية ، حيث تكون الدقة أعلى بكثير (لكنها لا تزال محدودة، وبالتالي لا تزال تتطلب التقريب). تستغرق الأعداد ذات الدقة العشوائية الكبيرة وقتًا طويلاً في الحساب وتستهلك مساحة تخزين كبيرة. يُطلق على أحد أنواع الحساب الفاصل اسميخزن حساب الكرة نقطة منتصف عالية الدقة ونصف قطر منخفض الدقة بدلاً من حدين أدنى وأعلى عاليي الدقة كما هو الحال في حساب الفترات التقليدي (الحد الأدنى - الحد الأعلى)، [ 4 ] مما يقلل من تكلفة الحساب. [ 5 ]

مشكلة التبعية

تقدير تقريبي لنطاق القيمة

تُعدّ مشكلة التبعية عائقًا رئيسيًا أمام تطبيق الحساب الفتري. فعلى الرغم من قدرة طرق الفترات على تحديد نطاق العمليات والدوال الحسابية الأساسية بدقة عالية، إلا أن هذا لا ينطبق دائمًا على الدوال الأكثر تعقيدًا. فإذا تكررت فترة ما عدة مرات في عملية حسابية باستخدام معاملات، وتمّ أخذ كل تكرار على حدة، فقد يؤدي ذلك إلى توسع غير مرغوب فيه في الفترات الناتجة.

التعامل مع كل ظهور لمتغير بشكل مستقل

كمثال توضيحي، لنأخذ الدالة f المعرفة بالعلاقة f ( x ) = + x . قيم هذه الدالة على الفترة [ -1 , 1 ] هي [ -1/4 , 2 ] . وباعتبارها امتدادًا طبيعيًا للفترة ، تُحسب على النحو التالي :

[-1،1]2+[-1،1]=[0،1]+[-1،1]=[-1،2]،{\displaystyle [-1,1]^{2}+[-1,1]=[0,1]+[-1,1]=[-1,2],}

وهو أكبر قليلاً؛ بدلاً من ذلك ، قمنا بحساب الحد الأدنى والحد الأعلى للدالة h ( x , y ) = + y على الفترة [ -1, 1 ] . يوجد تعبير أفضل للدالة f حيث يظهر المتغير x مرة واحدة فقط، وذلك بإعادة كتابة المعادلة التربيعية f ( x ) = + x بإكمال المربع .

و(x)=(x+12)2-14.{\displaystyle f(x)=\left(x+{\tfrac {1}{2}}\right)^{2}-{\tfrac {1}{4}}.}

لذا فإن حساب الفترة المناسب هو

([-1،1]+12)2-14=[-12،32]2-14=[0،94]-14=[-14،2]{\displaystyle \left([-1,1]+{\tfrac {1}{2}}\right)^{2}-{\tfrac {1}{4}}=\left[-{\tfrac {1}{2}},{\tfrac {3}{2}}\right]^{2}-{\tfrac {1}{4}}=\left[0,{\tfrac {9}{4}}\right]-{\tfrac {1}{4}}=\left[-{\tfrac {1}{4}},2\right]}

ويعطي القيم الصحيحة.

بشكل عام، يمكن إثبات إمكانية الوصول إلى النطاق الدقيق للقيم، إذا ظهر كل متغير مرة واحدة فقط وكانت الدالة f متصلة داخل المربع. مع ذلك، لا يمكن إعادة كتابة كل دالة بهذه الطريقة.

تأثير التغليف

قد يصل اعتماد المشكلة التي تسبب المبالغة في تقدير نطاق القيمة إلى حد تغطية نطاق واسع، مما يمنع التوصل إلى استنتاجات أكثر جدوى.

ينتج عن حل المناطق التي لا تأخذ شكل متجه فاصل زيادة إضافية في المدى. مجموعة حلول النظام الخطي

{x=صy=صص[-1،1]{\displaystyle {\begin{cases}x=p\\y=p\end{cases}}\qquad p\in [-1,1]}

هو بالضبط الخط الواصل بين النقطتين (−1, −1) و (1, 1) . باستخدام طرق الفترات، نحصل على المربع ذي الوحدة، [ −1, 1 ] × [ −1, 1 ] . يُعرف هذا بتأثير الالتفاف .

أنظمة الفترات الخطية

يتكون نظام الفترات الخطية من امتداد فترة المصفوفة [ A ] ∈ [ℝ] n × m ومتجه فترة [ b ] ∈ [ℝ] n . نريد أصغر متوازي مستطيلات [ x ] ∈ [ℝ] m ، لكل متجه x ∈ ℝ m يوجد زوج ( A , b ) حيث A ∈ [ A ] و b ∈ [ b ] يحقق الشرط التالي:

أx=ب{\displaystyle \mathbf {A} \cdot \mathbf {x} =\mathbf {b} }.

بالنسبة للأنظمة التربيعية (أي عندما n = m )، يمكن إيجاد متجه فترات [ x ] يغطي جميع الحلول الممكنة، وذلك ببساطة باستخدام طريقة جاوس للفترات. تحل هذه الطريقة محل العمليات العددية، حيث تصبح طريقة الجبر الخطي المعروفة باسم حذف جاوس هي النسخة الخاصة بها من طريقة الفترات. مع ذلك، ولأن هذه الطريقة تستخدم عناصر الفترات [ A ] و [ b ] بشكل متكرر في الحساب، فقد تُنتج نتائج غير دقيقة لبعض المسائل. لذا، فإن استخدام نتيجة جاوس ذات القيم الفتراتية لا يوفر سوى تقديرات أولية تقريبية، لأنه على الرغم من احتوائها على مجموعة الحلول كاملة، إلا أنها تحتوي أيضًا على مساحة كبيرة خارجها.

يمكن تحسين الحل التقريبي [ x ] غالبًا باستخدام نسخة فاصلية من طريقة جاوس-سيدل . والدافع وراء ذلك هو أن الصف i من امتداد الفاصل للمعادلة الخطية

([أ11][أ1ن][أن1][أنن])(x1xن)=([ب1][بن]){\displaystyle {\begin{pmatrix}{[a_{11}]}&\cdots &{[a_{1n}]}\\\vdots &\ddots &\vdots \\{[a_{n1}]}&\cdots &{[a_{nn}]}\end{pmatrix}}\cdot {\begin{pmatrix}{x_{1}}\\\vdots \\{x_{n}}\end{pmatrix}}={\begin{pmatrix}{[b_{1}]}\\\vdots \\{[b_{n}]}\end{pmatrix}}}

يمكن تحديدها بواسطة المتغير xᵢ إذا سُمح بالقسمة 1 / [ aᵢᵢ ] . وبالتالي ، فهي في آن واحد

xج[xج]وxج[بأنا]-كج[أأناك][xك][أأناأنا].{\displaystyle x_{j}\in [x_{j}]\quad {\text{and}}\quad x_{j}\in {\frac {[b_{i}]-\sum \limits _{k\neq j}[a_{ik}]\cdot [x_{k}]}{[a_{ii}]}}.}

لذا يمكننا الآن استبدال [ x j ] بـ

[xج][بأنا]-كج[أأناك][xك][أأناأنا]،{\displaystyle [x_{j}]\cap {\frac {[b_{i}]-\sum \limits _{k\neq j}[a_{ik}]\cdot [x_{k}]}{[a_{ii}]}},}

وبالتالي فإن المتجه [ x ] لكل عنصر.

بما أن الإجراء أكثر كفاءة للمصفوفة ذات المصفوفة القطرية المهيمنة ، فبدلاً من النظام [ A ] · x = [ b ] ، يمكن للمرء غالبًا محاولة ضربه بمصفوفة كسرية مناسبة M مع معادلة المصفوفة الناتجة

(م[أ])x=م[ب]{\displaystyle (\mathbf {M} \cdot [\mathbf {A} ])\cdot \mathbf {x} =\mathbf {M} \cdot [\mathbf {b} ]}

بقي حلها. إذا اخترنا، على سبيل المثال، M = A −1 للمصفوفة المركزية A ∈ [ A ] ، فإن M · [ A ] هي امتداد خارجي لمصفوفة الوحدة.

لا تُجدي هذه الطرق نفعًا إلا إذا كانت أطوال الفترات صغيرة بما يكفي. أما في حالة الفترات الأوسع، فقد يكون من المفيد استخدام نظام خطي فتري على أنظمة خطية مكافئة للأعداد الحقيقية المحدودة (وإن كانت كبيرة). إذا كانت جميع المصفوفات A ∈ [ A ] قابلة للعكس، يكفي النظر في جميع التوليفات الممكنة (العليا والسفلى) لنقاط النهاية الواقعة في الفترات. ويمكن حل المسائل الناتجة باستخدام الطرق العددية التقليدية. ولا يزال حساب الفترات يُستخدم لتحديد أخطاء التقريب.

هذا الأسلوب مناسب فقط للأنظمة ذات الأبعاد الصغيرة، لأنه مع مصفوفة n × n ممتلئة تمامًا ، يلزم عكس 2n² مصفوفة حقيقية، مع 2n متجهًا للطرف الأيمن . وقد طوّر جيري رون هذا الأسلوب، ولا يزال قيد التطوير . [ 6 ]

طريقة نيوتن الفاصلة

تقليل مساحة البحث في خطوة نيوتن الفاصلة في الدوال "السميكة".

يمكن اشتقاق صيغة فترة من طريقة نيوتن لإيجاد أصفار متجه فترة [ x ] من امتداد القيمة المتوسطة. [ 7 ] بالنسبة لمتجه مجهول z ∈ [ x ] مطبق على y ∈ [ x ] ، نحصل على

و(z)و(y)+[جو]([x])(z-y).{\displaystyle f(\mathbf {z} )\in f(\mathbf {y} )+[J_{f}](\mathbf {[x]} )\cdot (\mathbf {z} -\mathbf {y} ).}

بالنسبة لقيمة z تساوي صفرًا ، أي f ( z ) = 0 ، يجب أن تحقق الشرط التالي:

و(y)+[جو]([x])(z-y)=0.{\displaystyle f(\mathbf {y} )+[J_{f}](\mathbf {[x]} )\cdot (\mathbf {z} -\mathbf {y} )=0.}

هذا يعادل

zy-[جو]([x])-1و(y).{\displaystyle \mathbf {z} \in \mathbf {y} -[J_{f}](\mathbf {[x]} )^{-1}\cdot f(\mathbf {y} ).}

يمكن تحديد تقدير خارجي لـ [ J f ]([ x ]) −1 · f ( y ) باستخدام الطرق الخطية.

في كل خطوة من خطوات طريقة نيوتن الفاصلية، يتم استبدال قيمة ابتدائية تقريبية [ x ] ∈ [ℝ] n بـ

[x](y-[جو]([x])-1و(y)){\displaystyle [\mathbf {x} ]\cap \left(\mathbf {y} -[J_{f}](\mathbf {[x]} )^{-1}\cdot f(\mathbf {y} )\right)}

وبالتالي يمكن تحسين النتيجة. على عكس الطرق التقليدية، تقترب طريقة الفترات من النتيجة باحتوائها على الأصفار. وهذا يضمن أن النتيجة تُنتج جميع الأصفار في النطاق الأولي. وعلى العكس، تُثبت هذه الطريقة أنه لم تكن هناك أي أصفار للدالة f في النطاق الأولي [ x ] إذا أنتجت خطوة نيوتن المجموعة الفارغة.

تتقارب هذه الطريقة على جميع الأصفار في المنطقة الابتدائية. يمكن أن تؤدي القسمة على صفر إلى فصل الأصفار المتميزة، على الرغم من أن الفصل قد لا يكون كاملاً؛ ويمكن استكماله بطريقة التنصيف .

كمثال، لنفترض الدالة f ( x ) = - 2 ، والمدى الابتدائي [ x ] = [ -2, 2 ] ، والنقطة y = 0. عندئذٍ ، يكون لدينا Jf ( x ) = 2x ، وتعطينا خطوة نيوتن الأولى

[-2،2](0-12[-2،2](0-2))=[-2،2]([-،-12][12،])=[-2،-12][12،2].{\displaystyle [-2,2]\cap \left(0-{\frac {1}{2\cdot [-2,2]}}(0-2)\right)=[-2,2]\cap \left(\left[-\infty ,-{\tfrac {1}{2}}\right]\cup \left[{\tfrac {1}{2}},\infty \right]\right)=\left[-2,-{\tfrac {1}{2}}\right]\cup \left[{\tfrac {1}{2}},2\right].}

تُستخدم خطوات نيوتن الإضافية بشكل منفصل على x[ −2, − 1 / 2 ] و x[ 1 / 2 , 2 ] . وتتقارب هذه الخطوات إلى فترات صغيرة بشكل تعسفي حول 2 و + 2 على التوالي.

يمكن استخدام طريقة نيوتن للفترات مع الدوال السميكة مثل g ( x ) = - [ 2, 3 ] ، والتي ستكون نتائجها فترات في جميع الأحوال . ينتج عن ذلك فترات تحتوي على [ −√3 , −√2 ][ √2 , √3 ] .

التنصيف والتغطية

تقدير تقريبي ( باللون الفيروزي ) وتقديرات محسّنة من خلال "التقطيع" ( باللون الأحمر )

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

تغطية متجه الفترة [ x ] بمربعات أصغر [ x1 ] ، ...، [ xk ] ، بحيث

[x]=أنا=1ك[xأنا]،{\displaystyle [\mathbf {x} ]=\bigcup _{i=1}^{k}[\mathbf {x} _{i}],}

ويكون صالحاً ضمن نطاق القيم.

و([x])=أنا=1كو([xأنا]).{\displaystyle f([\mathbf {x} ])=\bigcup _{i=1}^{k}f([\mathbf {x} _{i}]).}

لذا، بالنسبة لتمديدات الفترات الموصوفة أعلاه، ينطبق ما يلي:

[و]([x])أنا=1ك[و]([xأنا]).{\displaystyle [f]([\mathbf {x} ])\supseteq \bigcup _{i=1}^{k}[f]([\mathbf {x} _{i}]).}

بما أن [ f ]([ x ]) غالباً ما تكون مجموعة فائقة حقيقية للجانب الأيمن، فإن هذا عادة ما يؤدي إلى تقدير محسّن.

يمكن توليد مثل هذا الغطاء باستخدام طريقة التنصيف، مثل العناصر السميكة [ xi1 , xi2 ] لمتجه الفاصل الزمني

[x]=([x11،x12]،...،[xن1،xن2]){\displaystyle [\mathbf {x} ]=([x_{11},x_{12}],\ldots ,[x_{n1},x_{n2}])}

بتقسيمها في المنتصف إلى فترتين

[xأنا1،xأنا1+xأنا22]و[xأنا1+xأنا22،xأنا2].{\displaystyle \left[x_{i1},{\frac {x_{i1}+x_{i2}}{2}}\right]\quad {\text{and}}\quad \left[{\frac {x_{i1}+x_{i2}}{2}},x_{i2}\right].}

إذا لم تكن النتيجة مناسبة، فمن الممكن إجراء تقسيم تدريجي إضافي. ينتج عن تقسيم عناصر المتجه إلى r فترات تغطية تساوي 2r ، مما يزيد بشكل كبير من تكاليف الحساب.

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

طلب

يمكن استخدام حساب الفترات في مجالات متنوعة (مثل عكس المجموعة ، وتخطيط الحركة ، وتقدير المجموعة ، أو تحليل الاستقرار) لمعالجة التقديرات التي لا تحتوي على قيمة عددية دقيقة. [ 8 ]

تحليل أخطاء التقريب

تُستخدم حسابات الفترات مع تحليل الأخطاء للتحكم في أخطاء التقريب الناتجة عن كل عملية حسابية. وتكمن ميزة حسابات الفترات في وجود فترة بعد كل عملية تتضمن النتيجة الحقيقية بدقة. وتعطي المسافة بين حدود هذه الفترة حساب أخطاء التقريب الحالي بشكل مباشر.

الخطأ = القيمة المطلقة ( أ - ب ) لفترة معينة [ أ ، ب ] .

يُضيف تحليل الفترات إلى الطرق التقليدية لتقليل الأخطاء، مثل التمحور ، بدلاً من أن يحل محلها .

تحليل التسامح

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

إذا كان سلوك مثل هذا النظام المتأثر بالتفاوتات يحقق، على سبيل المثال، f ( x , p ) = 0 ، لـ p ∈ [ p ] و x غير معروف ، فإن مجموعة الحلول الممكنة.

{x|ص[ص]،و(x،ص)=0}{\displaystyle \{\mathbf {x} \,|\,\exists \mathbf {p} \in [\mathbf {p} ],f(\mathbf {x} ,\mathbf {p} )=0\}}،

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

الحساب الفاصل الضبابي

تقريب التوزيع الطبيعي بواسطة سلسلة من الفترات

يمكن استخدام حساب الفترات مع دوال الانتماء للكميات الضبابية كما هو الحال في المنطق الضبابي . إلى جانب العبارات الصارمة x ∈ [ x ] و x ∉ [ x ] ، توجد قيم وسيطة تُسند إليها أعداد حقيقية μ[ 0, 1 ] . μ = 1 تعني انتماءً مؤكدًا، بينما μ = 0 تعني عدم الانتماء. تُحدد دالة التوزيع عدم اليقين، والذي يمكن فهمه على أنه فترة إضافية.

في الحساب الضبابي [ 9 ] ، يُؤخذ في الاعتبار عدد محدود فقط من مراحل الانتماء المنفصلة μ i[ 0, 1 ] . ويمكن تمثيل شكل هذا التوزيع لقيمة غير محددة بسلسلة من الفترات.

[x(1)][x(2)][x(ك)].{\displaystyle \left[x^{(1)}\right]\supset \left[x^{(2)}\right]\supset \cdots \supset \left[x^{(k)}\right].}

تتوافق الفترة [ x ( i ) ] تمامًا مع نطاق التذبذب للمرحلة μ i .

التوزيع المناسب للدالة f ( x1 , …, xn ) فيما يتعلق بالقيم غير المميزة x1 , …, xn والمتتابعات المقابلة لها

[x1(1)][x1(ك)]،...،[xن(1)][xن(ك)]{\displaystyle \left[x_{1}^{(1)}\right]\supset \cdots \supset \left[x_{1}^{(k)}\right],\ldots ,\left[x_{n}^{(1)}\right]\supset \cdots \supset \left[x_{n}^{(k)}\right]}

يمكن تقريبها بواسطة المتتالية

[y(1)][y(ك)]،{\displaystyle \left[y^{(1)}\right]\supset \cdots \supset \left[y^{(k)}\right],}

أين

[y(أنا)]=و([x1(أنا)]،...[xن(أنا)]){\displaystyle \left[y^{(i)}\right]=f\left(\left[x_{1}^{(i)}\right],\ldots \left[x_{n}^{(i)}\right]\right)}

ويمكن حسابها باستخدام طرق الفترات. القيمة [ y (1) ] تتوافق مع نتيجة حساب الفترة.

التدقيق بمساعدة الحاسوب

استخدم وارويك تاكر حساب الفترات لحل المسألة الرابعة عشرة من مسائل سميل ، أي لإثبات أن جاذب لورنز هو جاذب غريب . [ 10 ] استخدم توماس هيلز حساب الفترات لحل حدسية كبلر .

يسرد موقع Fast Library for Number Theory الإلكتروني العديد من استخدامات مكتبة FLINT/Arb للحساب الكروي في الأوراق المنشورة والمطبوعات الأولية. [ 11 ]

تاريخ

لا يُعدّ الحساب الفتري ظاهرةً جديدةً تمامًا في الرياضيات؛ فقد ظهر عدة مرات بأسماء مختلفة عبر التاريخ. على سبيل المثال، قام أرخميدس بحساب الحدين الأدنى والأعلى 223 / 71 < π < 22 / 7 في القرن الثالث قبل الميلاد. لم يحظَ الحساب الفعلي باستخدام الفترات بشعبيةٍ واسعةٍ كغيره من التقنيات العددية، ولم يُنسَ تمامًا.

نُشرت قواعد الحساب باستخدام الفترات ومجموعات فرعية أخرى من الأعداد الحقيقية في كتابٍ لروزاليند سيسلي يونغ عام ١٩٣١. [ ١٢ ] ثم نُشرت دراسة حسابية حول أعداد النطاقات لتحسين موثوقية الأنظمة الرقمية في كتابٍ مدرسيٍّ عن الجبر الخطي لبول س. دوير عام ١٩٥١ ؛ [ ١٣ ] استُخدمت الفترات لقياس أخطاء التقريب المرتبطة بالأعداد العشرية. ونُشرت ورقة بحثية شاملة حول جبر الفترات في التحليل العددي لتيروو سوناغا (١٩٥٨). [ ١٤ ]

شهد ظهور كتاب " تحليل الفترات " لرامون إي. مور عام 1966 ميلاد الحساب الفتري الحديث. [ 15 ] [ 16 ] وقد راودته الفكرة في ربيع عام 1958، وبعد عام نشر مقالًا عن الحساب الفتري الحاسوبي. [ 17 ] وتكمن ميزته في أنه انطلاقًا من مبدأ بسيط، قدّم طريقة عامة لتحليل الأخطاء آليًا، وليس فقط الأخطاء الناتجة عن التقريب.

بشكل مستقل في عام 1956، اقترح ميتشيسلاف وارموس صيغًا للحسابات مع الفترات، [ 18 ] على الرغم من أن مور وجد التطبيقات غير التافهة الأولى.

في العشرين عامًا التالية، قامت مجموعات من الباحثين الألمان بأعمال رائدة بقيادة أولريش دبليو كوليش [ 19 ] [ 20 ] وغوتز أليفيلد [ 21 ] في جامعة كارلسروه، ولاحقًا في جامعة بيرغيش في فوبرتال . على سبيل المثال، استكشف كارل نيكل تطبيقات أكثر فعالية، بينما تُعزى إجراءات الاحتواء المحسّنة لمجموعة حلول أنظمة المعادلات إلى أرنولد نوماير، من بين آخرين. في ستينيات القرن الماضي، تناول إلدون آر هانسن امتدادات الفترات للمعادلات الخطية، ثم قدم مساهمات حاسمة في التحسين الشامل، بما في ذلك ما يُعرف الآن بطريقة هانسن، والتي ربما تكون خوارزمية الفترات الأكثر استخدامًا. [ 7 ] غالبًا ما تواجه الطرق الكلاسيكية في هذا المجال مشكلة تحديد أكبر (أو أصغر) قيمة شاملة، ولكنها لا تستطيع سوى إيجاد قيمة مثلى محلية، ولا تستطيع إيجاد قيم أفضل. قام هيلموت راتشيك وجون جورج روكني بتطوير أساليب التفرع والتقييد ، والتي كانت حتى ذلك الحين تنطبق فقط على القيم الصحيحة، وذلك باستخدام الفترات الزمنية لتوفير تطبيقات للقيم المستمرة.

في عام 1988، قام رودولف لونر بتطوير برنامج قائم على لغة فورتران لإيجاد حلول موثوقة لمسائل القيمة الأولية باستخدام المعادلات التفاضلية العادية . [ 22 ]

تُصدر مجلة "الحوسبة الموثوقة " (التي كانت تُعرف سابقًا باسم "حسابات الفترات ") منذ تسعينيات القرن الماضي، وهي مُخصصة لموثوقية الحسابات المُساعدة بالحاسوب. وبصفته رئيس التحرير، ساهم ر. بيكر كيرفوت، بالإضافة إلى عمله في مجال التحسين الشامل، بشكلٍ كبير في توحيد الرموز والمصطلحات المُستخدمة في حساب الفترات. [ 23 ]

في السنوات الأخيرة، تركز العمل بشكل خاص على تقدير الصور العكسية للدوال ذات المعاملات وعلى نظرية التحكم القوي من قبل مجموعة عمل COPRIN التابعة لـ INRIA في صوفيا أنتيبوليس بفرنسا. [ 24 ]

التطبيقات

تتوفر العديد من حزم البرامج التي تتيح تطوير تطبيقات حسابية باستخدام حساب الفترات. [ 25 ] وعادةً ما تُقدم هذه الحزم على شكل مكتبات برمجية. كما توجد مُترجمات للغتين C++ وFortran تدعم أنواع بيانات الفترات والعمليات المناسبة كامتداد للغة، مما يجعل حساب الفترات مدعومًا بشكل مباشر.

منذ عام 1967، طُوِّرت امتدادات الحوسبة العلمية (XSC) في جامعة كارلسروه للغات برمجة متنوعة ، مثل C++ وFortran و Pascal . [ 26 ] كانت المنصة الأولى هي Zuse Z23 ، والتي أُتيح لها نوع بيانات فاصل زمني جديد مع عوامل تشغيل أولية مناسبة. تبع ذلك في عام 1976، Pascal-SC ، وهو إصدار من Pascal مبني على Zilog Z80 ، مما أتاح إنشاء إجراءات سريعة ومعقدة للتحقق الآلي من النتائج. ثم جاء ACRITH-XSC المبني على Fortran 77 لبنية System/370 (FORTRAN-SC)، والذي قدمته لاحقًا شركة IBM. بدءًا من عام 1991، أصبح بالإمكان إنتاج شفرة لمترجمات C باستخدام Pascal-XSC ؛ وبعد عام، دعمت مكتبة فئات C++ لغة C-XSC على العديد من أنظمة الحاسوب المختلفة. في عام 1997، أُتيحت جميع إصدارات XSC بموجب رخصة جنو العمومية العامة . في بداية عام 2000، تم إصدار C-XSC 2.0 تحت قيادة مجموعة العمل الخاصة بالحساب العلمي في جامعة بيرغيش في فوبرتال ليتوافق مع معيار C++ المحسن.

تم إنشاء مكتبة أخرى من فئة C++ في عام 1993 في جامعة هامبورغ للتكنولوجيا تُسمى Profil/BIAS (مكتبة الفترات السريعة المُحسّنة لوقت تشغيل المبرمج، والحسابات الأساسية للفترات)، والتي جعلت عمليات الفترات المعتادة أكثر سهولة في الاستخدام. وقد ركزت على الاستخدام الفعال للأجهزة، وقابلية النقل، واستقلالية عرض الفترات.

تحتوي مجموعة مكتبات Boost للغة C++ على فئة نموذجية للفترات. ويهدف مطوروها إلى تضمين العمليات الحسابية للفترات في لغة C++ القياسية. [ 27 ]

تتميز لغة البرمجة Frink بتطبيق لحساب الفترات يدعم التعامل مع الأعداد ذات الدقة العشوائية . ويمكن للبرامج المكتوبة بلغة Frink استخدام الفترات دون الحاجة إلى إعادة كتابتها أو إعادة تجميعها.

GAOL [ 28 ] هي مكتبة حسابية أخرى للفترات في لغة C++ وهي فريدة من نوعها لأنها توفر عوامل التشغيل العلائقية للفترات المستخدمة في برمجة قيود الفترات .

تُعد مكتبة مور [ 29 ] تطبيقًا فعالًا للحسابات الفاصلية في لغة C++. وهي توفر فترات ذات نقاط نهاية بدقة عشوائية وتعتمد على ميزة المفاهيم في لغة C++ .

تتضمن لغة البرمجة جوليا [ 30 ] تطبيقًا للحسابات الفاصلية إلى جانب ميزات عالية المستوى، مثل إيجاد الجذور (للدوال الحقيقية والمركبة) وبرمجة القيود الفاصلية ، عبر الحزمة [ 31 ]ValidatedNumerics.jl .

بالإضافة إلى ذلك، تستطيع أنظمة الجبر الحاسوبي، مثل Euler Mathematical Toolbox و FriCAS و Maple و Mathematica و Maxima [ 32 ] و MuPAD ، التعامل مع الفترات. ويعتمد امتداد Matlab المسمى Intlab [ 33 ] على إجراءات BLAS ، ويُوفر هذا الصندوق b4mواجهة Profil/BIAS. [ 33 ] [ 34 ]

تمت كتابة مكتبة للغة الوظيفية OCaml بلغة التجميع ولغة C. [ 35 ]

MPFI هي مكتبة للحسابات الفاصلية ذات الدقة العشوائية؛ وهي مكتوبة بلغة C وتعتمد على MPFR . [ 36 ]

معيار IEEE 1788

تمت الموافقة على معيار حساب الفترات، IEEE Std 1788-2015، في يونيو 2015. [ 37 ] يتوفر تطبيقان مرجعيان مجانًا. [ 38 ] وقد طُوِّر هذان التطبيقان من قِبَل أعضاء فريق عمل المعيار: مكتبة libieeep1788 [ 39 ] للغة C++، وحزمة interval [ 40 ] لبرنامج GNU Octave .

تمت الموافقة على مجموعة فرعية مصغرة من المعيار، IEEE Std 1788.1-2017، في ديسمبر 2017 ونُشرت في فبراير 2018. ومن المتوقع أن يكون تطبيقها أسهل وقد يُسرّع من إنتاج التطبيقات. [ 41 ]

المؤتمرات وورش العمل

تُعقد العديد من المؤتمرات وورش العمل الدولية سنوياً في جميع أنحاء العالم. ولعلّ أبرزها مؤتمر SCAN (الندوة الدولية حول الحوسبة العلمية، والحساب الحاسوبي، والحساب العددي المُدقّق)، بالإضافة إلى مؤتمرات أخرى مثل SWIM (ورشة العمل المصغّرة حول أساليب الفترات)، وPPAM (المؤتمر الدولي حول المعالجة المتوازية والرياضيات التطبيقية)، وREC (ورشة العمل الدولية حول الحوسبة الهندسية الموثوقة).

انظر أيضاً

مراجع

  1. 1 2 دراير، ألكسندر (2003). تحليل الفترات للدوائر التناظرية ذات التفاوتات في المكونات . آخن، ألمانيا: شاكر فيرلاغ . ص  15. ISBN 3-8322-4555-3.
  2. الحساب الفتراتي المركب وتطبيقاته ، ميودراغ س. بيتكوفيتش ، ليليانا د. بيتكوفيتش، وايلي-في سي إتش ، 1998، رقم ISBN 978-3-527-40134-5
  3. 1 2 3 4 5 6 هند داوود (2011). نظريات الحساب الفتري: الأسس الرياضية والتطبيقات . ساربروكن: دار لامبرت للنشر الأكاديمي. ISBN 978-3-8465-0154-2.
  4. فان دير هوفن، جوريس. "حساب الكرة" . www.texmacs.org . نركز بشكل أساسي على الحساب التلقائي والفعال لحدود الخطأ عالية الجودة، استنادًا إلى نوع من حساب الفترات الذي نود أن نسميه "حساب الكرة".
  5. "نظرة عامة على الميزات - وثائق FLINT 3.5.0-dev" . flintlib.org . في الحساب التقليدي للفترات (inf-sup)، تكون كلتا نقطتي نهاية الفترة [ a , b ] أعدادًا كاملة الدقة، مما يجعل حساب الفترات أكثر تكلفة بمرتين من حساب الفاصلة العائمة. في حساب الكرة، تكون نقطة المنتصف m فقط من الفترة [ m ± r ] عددًا كامل الدقة، ويكفي عدد قليل من البتات لنصف القطر r . عند الدقة العالية، لا يكون حساب الكرة أكثر تكلفة من حساب الفاصلة العائمة العادي.
  6. "جيري رون، قائمة المنشورات" . مؤرشف من الأصل بتاريخ 23-11-2008 . تم الاطلاع عليه بتاريخ 26-05-2008 .
  7. 1 2 والستر، ج. ويليام؛ هانسن، إلدون روبرت (2004). التحسين العالمي باستخدام تحليل الفاصل الزمني (الطبعة الثانية ). نيويورك، الولايات المتحدة الأمريكية: مارسيل ديكر. رقم ISBN  0-8247-4059-9.
  8. ^ جولين، لوك؛ كيفير، ميشيل. ديدريت، أوليفييه؛ والتر، اريك (2001). تحليل الفاصل الزمني التطبيقي . برلين: سبرينغر. رقم ISBN 1-85233-219-0.
  9. تطبيق الحساب الضبابي لتحديد آثار معلمات النموذج غير المؤكدة، مايكل هانس ، جامعة شتوتغارت
  10. ^ تاكر ، وارويك (1999). جاذب لورينز موجود. Comptes Rendus de l'Académie des Sciences-Series I-Mathematics, 328(12), 1197-1202.
  11. "التطبيقات والمعايير: الحساب العددي" . FLINT: مكتبة سريعة لنظرية الأعداد (flintlib.org) .
  12. يونغ، روزاليند سيسلي (1931). جبر الكميات متعددة القيم. حوليات الرياضيات، 104(1)، 260-290. (ملاحظة: مرشحة لنيل درجة الدكتوراه في جامعة كامبريدج .)
  13. دوير، بول سومنر (1951). الحسابات الخطية. أكسفورد، إنجلترا: وايلي. ( جامعة ميشيغان )
  14. ^ سوناجا ، تيرو (1958). “نظرية الجبر الفاصل وتطبيقها على التحليل العددي”. مذكرات راج (2): 29- 46.
  15. مور، رامون إدغار (1966). تحليل الفترات . إنجلوود كليف، نيو جيرسي، الولايات المتحدة الأمريكية: برنتيس هول . ISBN 0-13-476853-1.
  16. كلاود، مايكل جيه؛ مور، رامون إدغار ؛ كيرفوت، آر. بيكر (2009). مقدمة في تحليل الفترات . فيلادلفيا: جمعية الرياضيات الصناعية والتطبيقية (SIAM). ISBN 978-0-89871-669-6.
  17. هانسن، إلدون روبرت (13 أغسطس 2001). "منشورات متعلقة بأعمال آر إي مور المبكرة في فترة ما بين الفصول" . مطبعة جامعة لويزيانا في لافاييت . تاريخ الاسترجاع: 29 يونيو 2015 .
  18. أوراق تمهيدية حول تحليل الفترات بقلم ميتشيسلاف فارموس، مؤرشفة في 18 أبريل 2008 على موقع Wayback Machine
  19. ^ كوليش، أولريش دبليو (1989). Wissenschaftliches Rechnen mit Ergebnis Verifikation. عين اينفورونج (في المانيا). فيسبادن: Vieweg-Verlag . رقم ISBN 3-528-08943-1.
  20. ^ كوليش، أولريش دبليو (1969). "Grundzüge der Intervallrechnung". في لاوغويتز، ديتليف (محرر). Jahrbuch Überblicke Mathematik (باللغة الألمانية). المجلد. 2. مانهايم، ألمانيا: معهد الببليوغرافيات . ص 51 – 98.  
  21. ^ ألفيلد، جوتز [بالألمانية] ؛ هيرزبيرجر، يورغن (1974). Einführung in die Intervallrechnung . ريهي إنفورماتيك (باللغة الألمانية). المجلد. 12. مانهايم، فيينا، زيورخ: BI-Wissenschaftsverlag . رقم ISBN  3-411-01466-0.
  22. حدود المعادلات التفاضلية العادية لرودولف لونهر، مؤرشفة في 11 مايو 2018 في أرشيف الإنترنت (باللغة الألمانية)
  23. قائمة مراجع ر. بيكر كيرفوت ، جامعة لويزيانا في لافاييت
  24. فيلم تمهيدي (mpeg) لفرق COPRIN في INRIA ، صوفيا أنتيبوليس
  25. برنامج لحسابات الفترات مؤرشف في 2006-03-02 في Wayback Machine جمعه فلاديك كرينوفيتش ، جامعة تكساس في إل باسو .
  26. تاريخ لغات XSC، مؤرشف بتاريخ 29-09-2007 في أرشيف الإنترنت (Wayback Machine)
  27. اقتراح لإضافة العمليات الحسابية الفاصلية إلى مكتبة C++ القياسية
  28. Gaol ليست مجرد مكتبة حسابية أخرى للفترات
  29. مور: الحساب الفاصل في لغة C++ الحديثة
  30. لغة البرمجة جوليا
  31. ValidatedNumerics.jl
  32. حساب الفترات لإيجاد القيم العظمى: ملخص موجز بقلم ريتشارد ج. فيتمان.
  33. 1 2 "مختبر الفترات الزمنية Intlab" . مؤرشف من الأصل بتاريخ 30 يناير 2020. تم الاطلاع عليه بتاريخ 7 نوفمبر 2012 .
  34. b4m
  35. أليوت، جان مارك؛ غوتيلاند، جان باتيست؛ فاناريه، تشارلي؛ دوراند، نيكولا؛ جيانازا، ديفيد (2012). تنفيذ مكتبة حساب الفترات للغة OCaml على معمارية x86/amd64 . المؤتمر الدولي السابع عشر لجمعية ACM SIGPLAN حول البرمجة الوظيفية.
  36. MPFI، مكتبة حسابية متعددة الدقة تعتمد على MPFR
  37. معيار IEEE للحساب الفاصل
  38. ناتالي ريفول (2015). معيار IEEE 1788 (القريب) المستقبلي للحسابات الفاصلية، شرائح العرض التقديمي // ورشة العمل المصغرة الثامنة في أساليب الفترات SWIM 2015. براغ، 9-11 يونيو 2015
  39. تطبيق C++ لمعيار IEEE P1788 الأولي للحساب الفاصل
  40. حزمة فواصل GNU Octave
  41. "معيار IEEE 1788.1-2017 - معيار IEEE للحساب الفاصل (مبسط)" . معيار IEEE . جمعية معايير IEEE. 2017. مؤرشف من الأصل بتاريخ 7 فبراير 2018. تم الاطلاع عليه بتاريخ 6 فبراير 2018 .

للمزيد من القراءة