تم النشر اليوم 2024/05/20 | تحليل الأخطاء في التحليل العددي

التقريب (Approximation)

إنّ معظم الأعداد التي نتعامل معها هي أعداد تقريبية؛ لأنها غالبًا ما تمثّل أطوالًا وقياساتٍ أو قيمًا لمقاديرَ فيزيائيّة بنتيجة القياس، وهي بحدّ ذاتها تقريبية. كذلك فإن الكثير من الأعداد الحقيقية لا يمكن التعبير عنها بعدد منتهٍ من الأرقام، فمثلًا العدد π
{displaystyle pi }
يساوي «تقريبًا» 3.14159، كذلك هو الحال بالنسبة لـلعدد eو 7
{displaystyle {sqrt {7}}} مثلًا؛ فإنّنا لا نستطيع كتابتها كأعداد مضبوطة.
إننا ـ كما أوردنا في السابق ـ يتحتّم علينا أن نستخدم التقريب لحل المسائل التي لا تحلّ في الرياضيات التّحليلية، إلا أن هذا التقريب ستنتج عنه أخطاء، أهمها التالي: التدوير (Rounding)
هو واحد من أهمّ مصادر الأخطاء وهو استعمال الأعداد المدوّرة بدلًا من المضبوطة.
قاعدة التدوير: إذا كان لدينا العدد =
0. x 1 x 2
.
.
. x n
−
1 x n
{displaystyle X=0.x_{1}x_{2}…x_{n-1}x_{n}} ، وأردنا الاكتفاء بـ (n-1) عدد على يمين الفاصلة، أي أن العدد المدوّر يكون كما يلي:

0

=
0.

1

2

.
.
.
.
.

n
−
2

¯

n
−
1

{displaystyle X_{0}=0.X_{1}X_{2}…..X_{n-2}{bar {X}}_{n-1}}

فإن

¯

n
−
1

{displaystyle {bar {X}}_{n-1}}

يتم تحديده بالشكل الآتي:
إذا كان

n
{displaystyle X_{n}} (أكبر من 5 أو يساوي)، فإن ¯
n
−
1
=
n
−
1
+
1
{displaystyle {bar {X}}_{n-1}=X_{n-1}+1} إذا كان

n
{displaystyle X_{n}} (أصغر من 5)، فإن
n
−
1
¯ =
n
−
1
{displaystyle {bar {X_{n-1}}}=X_{n-1}} .
على سبيل المثال العددان 0.045651 و 9033 بتقريبهما باستخدام التدوير لثلاث خانات عشرية ينتج العددان 0.0457 و 9033 على التوالي. أما عن منشأ الخطأ العددي هنا فيكون بسبب الاكتفاء بعدد معين من المنازل العشريّة بعد الفاصلة في الحسابات. القطع (chopping)
إن آلات الحاسبة الإلكترونية لا تدور الأعداد غالبًا وإنما تقطعها.
قاعدة القطع: إذا كان لدينا العدد =
0
0
2
.
.
.
.
(
n
−
1
)
n
{displaystyle X=0X_{0}X_{2}….X_{(}n-1)X_{n}} وأردنا الاكتفاء بــ (1-n) عدد على يمين الفاصلة، فإن العدد المتقطع يكون كما يلي:
0
=
0.
1
2
.
.
.
.
(
n
−
2
)
(
n
−
1
)
{displaystyle X_{0}=0.X_{1}X_{2}….X(n-2)X_{(}n-1)} على سبيل المثال العددان 0.045651 و 9033 بتقريبها باستخدام القطع لثلاث خانات عشرية ينتج العددان 0.0456 و9030 على التوالي. ملاحظة: العدد 0 على يمين الفاصلة غير معتبر في عملية التقريب.
ومنشأ الخطأ العددي هنا هو يكون من الاكتفاء بعدد معين من حدود المتسلسلة اللانهائية.

الأخطاء (Error)

تلعب الأخطاء دورًا محوريًا في التحليل العددي وهي تبين مدى دقة وسرعة الطريقة المستخدمة لاحقًا. وفيما يلي أنواع الأخطاء وتأثيرها: الخطأ المطلق (Absolute Error)
وهو الخطأ المطلق المرتكب في العدد المقرب

0
{displaystyle X_{0}} هو القيمة المطلقة للفرق بينه وبين القيمة الفعلية Δ
= | Δ
.
x | = |
x 0
−
x | = | x
− x 0 | {displaystyle Delta =|Delta .x|=|x_{0}-x|=|x-x_{0}|}
غالبًا مايكون العدد الفعلي x غير معلوم، عندئذ لا يمكن تعيين الخطأ المطلق، من أجل ذلك نلجأ إلى إيجاد حد أعلى لهذا الخطأ مثل
ϵ x
{displaystyle epsilon _{x}} ويحقق المتباينة: | Δ
.
x | = | x
− x 0 | ≤ ε x
{displaystyle |Delta .x|=|x-x_{0}|leq varepsilon _{x}} .
وهكذا نحصل على:
( x 0
− ϵ x
)
≤
x
≤
( x 0
+ ϵ x
)
{displaystyle (x_{0}-epsilon _{x})leq xleq (x_{0}+epsilon _{x})}
أي أن العدد الفعلي يقع بين العددين
x 0
− ϵ x
{displaystyle x_{0}-epsilon _{x}} ,
x 0
+ ϵ x
{displaystyle x_{0}+epsilon _{x}} , حيث القيمة
x 0
+ ϵ x
{displaystyle x_{0}+epsilon _{x}} تمثل تقريب العدد x بالزيادة والقيمة
x 0
− ϵ x
{displaystyle x_{0}-epsilon _{x}} تمثل تقريبه بالنقصان. مثال
لتكن x=3.257 ولتكن القيمة التقريبية لـx هي 3.26 بتعيين الخطأ المطلق:
| Δ
.
x | = |
x 0
−
x | {displaystyle |Delta .x|=|x_{0}-x|}
| 3.26
−
3.257 | =
0.003
{displaystyle |3.26-3.257|=0.003} الخطأ النسبي (Relative)
الخطأ R النسبي المرتكب في عدد ما تقريبي
x 0
{displaystyle x_{0}} هو نسبة الخطأ المطلق لمرتكب بهذا العدد إلى القيمة المطلقة للعدد الفعلي x≠0) X) أي أن:
R
=
Δ
| x |
{displaystyle R={Delta over |x|}}
وفي حالة عدم معرفة القيمة الفعلية X نلجأ إلى الحد الأعلى لهذا الخطأ أي أن:

Δ

|

x

|

≤

S

x

{displaystyle {Delta over |x|}leq {S_{x}}}

،

R

≤

S

x

{displaystyle {R}leq {S_{x}}}

ومنه: Δ
≤ | x |
S x
{displaystyle Delta leq |x|S_{x}}
لذلك يمكن اعتبار:
ϵ x
= | x |
δ x
{displaystyle epsilon _{x}=|x|delta _{x}}
فيكون:
S x
= ϵ x
| x |
{displaystyle S_{x}={epsilon _{x} over |x|}}
الحد الأعلى للخطأ النسبي المرتكب في العدد المقرب
x 0
{displaystyle x_{0}} وبما أن
x ≅
x 0 {displaystyle {x}cong {x_{0}}} فيمكن أن يكتب ϵ x ≅
|
x 0 |
δ x {displaystyle {epsilon _{x}}cong {|x_{0}|delta _{x}}} مثال
لتكن x=3.257 ولتكن
x 0
=
3.026
{displaystyle x_{0}=3.026} قيمة تقريبية لــ x فإن الخطأ النسبي:

R
=

Δ

|

x

|

=

0.003
3.257

=
0.000921

{displaystyle R={Delta over |x|}={0.003 over 3.257}=0.000921}

ويستفاد من الخطأ النسبي في حساب الدقة المعنوية (significant digits) حيث العلاقة التالية:
Δ
| x |
≤ 5
× 10 −
t {displaystyle {Delta over |x|}leq {5times 10^{-}t}}
حيث t أكبر عدد صحيح غير سالب نستطيع أن نوجد من خلاله أكبر فترة تقريب للخطأ الفعلي.
ملاحظات عامّة: من المهمّ أن يُعلم أنّ الخطأ النسبي لا يتأثر بالخطأ كثيرًا بينما الخطأ المطلق يتأثّر بشكل ملحوظ؛ لذا نعتمد في الأغلب على الخطأ النسبي.
أيضًا كلّما كانت قيمة الخطأ الناتجة صغيرة جدًا كلما كان هذا أفضل، وعلى العكس من ذلك فكلما كبرت قيمة الخطأ احتجنا لأن نبحث عن إجراء رياضيّ يعمل على تصغير (تحسين) قيمة الخطأ الناتج.

المقدمة

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

أكبر خطأ نسبي في الدالة في عدة متغيرات (Max relative error in a function of several variables)

نفرض أن دالة f في متغيرn:

x

1

,

x

2

,
.
.
.
,

x

n

{displaystyle x_{1},x_{2},…,x_{n}}

أي أن:
f
=
f
( x 1
, x 2
,
.
.
.
, x n
)
{displaystyle f=f(x_{1},x_{2},…,x_{n})}
ونفرض أن الأخطاء المفردة في هذه المتغيرات هي على الترتيب:
Δ
. x 1
,
Δ
. x 2
,
.
.
.
Δ
. x n
{displaystyle Delta .x_{1},Delta .x_{2},…Delta .x_{n}}
فيكون الخطأ الكلي في الدالة مع اهمال حدود الرتبة الثانية وما فوقها هو:
Δ
.
f
= d
f
d x 1 Δ
. x 1
+ d
f
d x 2 Δ
. x 2
+
.
.
.
.
+ d
f
d x n Δ
. x n
{displaystyle Delta .f={df over {dx_{1}}}Delta .x_{1}+{df over {dx_{2}}}Delta .x_{2}+….+{df over {dx_{n}}}Delta .x_{n}} إذا
| Δ
.
f | =
| d
f |
|
d x 1
| | Δ
. x 1 | +
| d
f |
|
d x 2
| | Δ
. x 2 | +
.
.
.
.
+
| d
f |
|
d x n
| | Δ
. x n | {displaystyle |Delta .f|={|df| over |{dx_{1}}|}|Delta .x_{1}|+{|df| over |{dx_{2}}|}|Delta .x_{2}|+….+{|df| over |{dx_{n}}|}|Delta .x_{n}|} وبالتالي يكون أقصى خطأ نسبي هو:
| Δ
.
f
|
max | f |
=
1
| f |
| d
f |
|
d x 1
| | Δ
. x 1 | +
| d
f |
|
d x 2
| | Δ
. x 2 | +
.
.
.
.
+
| d
f |
|
d x n
| | Δ
. x n | {displaystyle {|Delta .f|_{max } over {|f|}}={1 over |f|}{|df| over |{dx_{1}}|}|Delta .x_{1}|+{|df| over |{dx_{2}}|}|Delta .x_{2}|+….+{|df| over |{dx_{n}}|}|Delta .x_{n}|} مثال
أكبير خطأ نسبي في كل من الدالتين التاليتيين بدلالة المتغيرات
x 1
, x 2
,
.
.
.
, x n
{displaystyle x_{1},x_{2},…,x_{n}}
وأخطأها Δ
. x 1
,
Δ
. x 2
,
.
.
.
, Δ n
{displaystyle Delta .x_{1},Delta .x_{2},…,Delta _{n}} f
= x 1
, x 2
,
.
.
.
, x n
{displaystyle f=x_{1},x_{2},…,x_{n}} g
= x 1 m 1
, x 2 m 2
,
.
.
.
, x n m n
{displaystyle g={x_{1}}^{m_{1}},{x_{2}}^{m_{2}},…,{x_{n}}^{m_{n}}} الحل:
| Δ
.
f
|
max
=
1. | Δ
. x 1 | +
1. | Δ
. x 2 | +
.
.
.
+
1. | Δ
. x n | {displaystyle |Delta .f|_{max }=1.|Delta .x_{1}|+1.|Delta .x_{2}|+…+1.|Delta .x_{n}|}
| Δ
.
f |
| f | max
=
| Δ
. x 1 | + | Δ
. x 2 | +
.
.
.
+ | Δ
. x n |
|
x 1
+ x 2
+
.
.
.
+ x n |
{displaystyle {|Delta .f| over {|f|}}_{max }={|Delta .x_{1}|+|Delta .x_{2}|+…+|Delta .x_{n}| over {|x_{1}+x_{2}+…+x_{n}|}}}
| Δ
.
g |
| g | max
= |
m 1 | | Δ
. x 1 |
|
x 1 |
+ |
m 2 | | Δ
. x 2 |
|
x 2 |
+
.
.
.
+ |
m n | | Δ
. x n |
|
x n |
{displaystyle {|Delta .g| over {|g|}}_{max }=|m_{1}|{|Delta .x_{1}| over {|x_{1}|}}+|m_{2}|{|Delta .x_{2}| over {|x_{2}|}}+…+|m_{n}|{|Delta .x_{n}| over {|x_{n}|}}} يُلاحظ عمومًا أن الخطأ الفعلي في قيمة يكون أقل بكثير من الحد الأقصى للخطأ؛ وذلك لأن الأخطاء من الناحية العلمية تتجه لأن يلغي بعضها بعضًا (حيث بعضها موجب والبعض الآخر سالب) فمثلًا إذا أضيف 200,000 عدد، وقرب كل عدد إلى 4 أرقام عشرية، فإن الحد الأقصى للخطأ يساوي
10
×
.5
×
20
,
000
{displaystyle 10times .5times 20,000} أو يساوي 1 وهذه قيمة كبيرة، بينما الخطأ الكلي المتوقع يكون عادة في حدود 0,005.

شرح مبسط

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