دمج الكونتور

في مجال التحليل المركب ، يُعد التكامل الكفافي طريقةً لتقييم بعض التكاملات على طول مسارات في المستوى المركب . [ 1 ] [ 2 ] [ 3 ] يُستخدم التكامل الكفافي لدراسة الدوال ذات القيم المركبة التي تكون تامة الشكل في منطقة معينة.

يرتبط التكامل الكفافي ارتباطًا وثيقًا بحساب البواقي ، [ 4 ] وهو أسلوب من أساليب التحليل المركب . تكمن قوة التكامل الكفافي في أن تكاملات الدوال الهولمورفية تظل ثابتة عند تشويه الكفاف، شريطة ألا يتقاطع التشوه مع نقطة تفرد أو قطع فرعي . وبالتالي، فإن قيمة التكامل الكفافي بين نقطتي نهاية ثابتتين لا تحددها الشكل الدقيق للكفاف، بل التفافه حول نقاط تفرد الدالة المراد تكاملها.

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

المنحنيات في المستوى المركب

في التحليل المركب ، يُعرف الخط الكفافي بأنه نوع من المنحنيات في المستوى المركب . وفي التكامل الكفافي، توفر الخطوط الكفافية تعريفًا دقيقًا للمنحنيات التي يمكن تعريف التكامل عليها بشكل مناسب. ويُعرَّف المنحنى في المستوى المركب بأنه دالة متصلة من فترة مغلقة على خط الأعداد الحقيقية إلى المستوى المركب.z:[أ،ب]ج{\displaystyle z:[a,b]\to \mathbb {C} }.

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

منحنيات ناعمة موجهة

غالباً ما يتم تعريف الخطوط الكنتورية من حيث المنحنيات الملساء الموجهة. [ 5 ] توفر هذه تعريفاً دقيقاً لـ "جزء" من منحنى أملس، والذي يتم تكوين الخط الكنتوري منه.

المنحنى الأملس هو منحنىz:[أ،ب]ج{\displaystyle z:[a,b]\to \mathbb {C} }بمشتقة متصلة غير معدومة بحيث يتم اجتياز كل نقطة مرة واحدة فقط ( حيث يكون z أحاديًا)، باستثناء المنحنى الذي تتطابق فيه نقاط النهاية (z(أ)=z(ب){\displaystyle z(a)=z(b)}في حالة تطابق نقاط النهاية، يُطلق على المنحنى اسم المنحنى المغلق، ويجب أن تكون الدالة أحادية في كل مكان آخر، ويجب أن تكون المشتقة متصلة عند النقطة المحددة.z(أ)=z(ب){\displaystyle z'(a)=z'(b)}). يُشار غالبًا إلى المنحنى الأملس غير المغلق باسم القوس الأملس. [ 5 ]

توفر عملية تحديد معلمات المنحنى ترتيبًا طبيعيًا للنقاط على المنحنى:z(x){\displaystyle z(x)}يأتي قبلz(y){\displaystyle z(y)}لوx<y{\displaystyle x<y}يؤدي هذا إلى مفهوم المنحنى الأملس الموجه . من المفيد جدًا دراسة المنحنيات بغض النظر عن المعاملات المحددة. يمكن تحقيق ذلك من خلال النظر في فئات التكافؤ للمنحنيات الملساء ذات الاتجاه نفسه. يمكن تعريف المنحنى الأملس الموجه على أنه مجموعة مرتبة من النقاط في المستوى المركب، وهي صورة لمنحنى أملس بترتيبها الطبيعي (وفقًا للمعاملات). تجدر الإشارة إلى أن ليس كل ترتيب للنقاط هو الترتيب الطبيعي للمنحنى الأملس. في الواقع، لا يوجد سوى ترتيبين طبيعيين للمنحنى الأملس. كذلك، يمكن أن يكون للمنحنى المغلق أي نقطة كنقطة نهاية، بينما لا يوجد للقوس الأملس سوى خيارين لنقطتي النهاية.

الخطوط الكنتورية

الخطوط الكنتورية هي فئة المنحنيات التي نُعرّف عليها التكامل الكنتوري. الخط الكنتوري هو منحنى موجه يتكون من سلسلة منتهية من المنحنيات الملساء الموجهة التي تتطابق نقاط نهايتها لتُعطي اتجاهًا واحدًا. وهذا يتطلب أن تكون سلسلة المنحنياتγ1،...،γن{\displaystyle \gamma _{1},\dots ,\gamma _{n}}أن تكون بحيث تكون النقطة النهائية لـγأنا{\displaystyle \gamma _{i}}يتزامن مع نقطة البداية لـγأنا+1{\displaystyle \gamma _{i+1}}للجميعأنا{\displaystyle i}بحيث1أنا<ن{\displaystyle 1\leq i<n}يشمل ذلك جميع المنحنيات الملساء الموجهة. كما تُعتبر نقطة واحدة في المستوى المركب محيطًا. الرمز+{\displaystyle +}يُستخدم هذا المصطلح غالبًا للدلالة على تجميع المنحنيات معًا لتشكيل منحنى جديد. وبالتالي، يمكننا كتابة منحنى محيطي.Γ{\displaystyle \Gamma }ذلك يتكون منن{\displaystyle n}المنحنيات كـ Γ=γ1+γ2++γن.{\displaystyle \Gamma =\gamma _{1}+\gamma _{2}+\cdots +\gamma _{n}.}

التكاملات الكفافية

التكامل الكفافي لدالة مركبةو:جج{\displaystyle f:\mathbb {C} \to \mathbb {C} }يُعدّ التكامل على طول منحنى مُوجَّه تعميمًا للتكامل على الدوال الحقيقية. بالنسبة للدوال المتصلة في المستوى المركب ، يُمكن تعريف التكامل على طول منحنى مُوجَّه أملس بدلالة التكامل على مُعامل حقيقي. يُمكن تقديم تعريف أكثر عمومية بدلالة تجزئة المنحنى على طول المنحنى، قياسًا بتجزئة فترة زمنية وتكامل ريمان . في كلتا الحالتين، يُعرَّف التكامل على طول منحنى ما على أنه مجموع التكاملات على المنحنيات المُوجَّهة الملساء التي تُشكِّل هذا المنحنى.

للدوال المتصلة

لتعريف التكامل الكفافي بهذه الطريقة، يجب أولاً النظر في التكامل، على متغير حقيقي، لدالة ذات قيم مركبة.و:Rج{\displaystyle f:\mathbb {R} \to \mathbb {C} }لتكن دالة ذات قيم مركبة لمتغير حقيقي،ت{\displaystyle t}الأجزاء الحقيقية والخيالية منو{\displaystyle f}يُشار إليها غالبًا باسمu(ت){\displaystyle u(t)}وv(ت){\displaystyle v(t)}، على التوالي، بحيث و(ت)=u(ت)+أناv(ت).{\displaystyle f(t)=u(t)+iv(t).} ثم تكامل الدالة ذات القيم المركبةو{\displaystyle f}خلال الفترة[أ،ب]{\displaystyle [a,b]}يُعطى بواسطة أبو(ت)دت=أب(u(ت)+أناv(ت))دت=أبu(ت)دت+أناأبv(ت)دت.{\displaystyle {\begin{aligned}\int _{a}^{b}f(t)\,dt&=\int _{a}^{b}{\big (}u(t)+iv(t){\big )}\,dt\\&=\int _{a}^{b}u(t)\,dt+i\int _{a}^{b}v(t)\,dt.\end{aligned}}}

والآن، لتعريف التكامل الكفافي، لنفترضو:جج{\displaystyle f:\mathbb {C} \to \mathbb {C} }لتكن دالة متصلة على المنحنى الأملس الموجهγ{\displaystyle \gamma }. يتركz:[أ،ب]ج{\displaystyle z:[a,b]\to \mathbb {C} }أي معلمات لـγ{\displaystyle \gamma }وهذا يتوافق مع ترتيبه (اتجاهه). ثم التكامل على طولγ{\displaystyle \gamma }يُشار إليه بـ γو(z)دz{\displaystyle \int _{\gamma }f(z)\,dz\,} ويتم تحديده بواسطة [ 5 ]γو(z)دz:=أبو(z(ت))z(ت)دت.{\displaystyle \int _{\gamma }f(z)\,dz:=\int _{a}^{b}f{\big (}z(t){\big )}z'(t)\,dt.}

هذا التعريف دقيق. أي أن النتيجة مستقلة عن المعاملات المختارة. [ 5 ] في حالة عدم وجود التكامل الحقيقي على الجانب الأيمن، يكون التكامل على طولγ{\displaystyle \gamma }يقال إنه غير موجود.

كتعميم لتكامل ريمان

يتم تعميم تكامل ريمان ليشمل الدوال ذات المتغير المركب بنفس الطريقة التي تم بها تعريفه للدوال من مجموعة الأعداد الحقيقية. تجزئة منحنى أملس موجهγ{\displaystyle \gamma }يُعرَّف بأنه مجموعة منتهية ومرتبة من النقاط علىγ{\displaystyle \gamma }التكامل على المنحنى هو نهاية المجاميع المحدودة لقيم الدالة، المأخوذة عند النقاط الموجودة على التقسيم، في الحد الذي تصبح فيه أقصى مسافة بين أي نقطتين متتاليتين على التقسيم (في المستوى المركب ثنائي الأبعاد)، والمعروف أيضًا باسم الشبكة، تؤول إلى الصفر.

الأساليب المباشرة

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

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

مثال

من النتائج الأساسية في التحليل المركب أن التكامل الكفافي لـ 1 / z يساوي 2πi ، حيث يُفترض أن مسار الكفاف هو دائرة الوحدة التي تم اجتيازها عكس اتجاه عقارب الساعة (أو أي منحنى جوردان موجب الاتجاه حول الصفر ). في حالة دائرة الوحدة ، توجد طريقة مباشرة لحساب التكامل. ج1zدz.{\displaystyle \oint _{C}{\frac {1}{z}}\,dz.}

عند حساب هذا التكامل، استخدم دائرة الوحدة | z | = 1 كمسار، معطاة بالمعادلة z ( t ) = e<sup> it</sup> ، حيث t ∈ [0, 2π] ، عندئذٍ يكون dz / dt = ie <sup> it </sup> و ج1zدz=02π1هـأناتأناهـأناتدت=أنا02π1دت=أنات|02π=(2π-0)أنا=2πأنا،{\displaystyle \oint _{C}{\frac {1}{z}}\,dz=\int _{0}^{2\pi }{\frac {1}{e^{it}}}ie^{it}\,dt=i\int _{0}^{2\pi }1\,dt=i\,t{\Big |}_{0}^{2\pi }=\left(2\pi -0\right)i=2\pi i,} وهي قيمة التكامل. تنطبق هذه النتيجة فقط في حالة رفع z إلى القوة -1 . إذا لم تكن القوة تساوي -1 ، فإن النتيجة ستكون دائمًا صفرًا.

تطبيقات النظريات التكاملية

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

تُستخدم نظريات التكامل مثل صيغة كوشي التكاملية أو نظرية البواقي بشكل عام في الطريقة التالية:

  • يتم اختيار شكل محيطي محدد:
    يتم اختيار المسار بحيث يتبع جزء المستوى المركب الذي يصف التكامل ذي القيم الحقيقية، ويحيط أيضًا بنقاط التفرد للدالة المراد تكاملها، مما يسمح بتطبيق صيغة كوشي التكاملية أو نظرية البواقي.
  • تطبيق نظرية كوشي التكاملية
    يتم اختزال التكامل إلى مجرد تكامل حول دائرة صغيرة حول كل قطب.
  • تطبيق صيغة كوشي التكاملية أو نظرية الباقي
    إن تطبيق هذه الصيغ التكاملية يعطينا قيمة للتكامل حول كامل المحيط.
  • تقسيم المحيط إلى محيط على طول الجزء الحقيقي وجزء تخيلي
    يمكن تقسيم المسار بأكمله إلى مسار يتبع الجزء من المستوى المركب الذي يصف التكامل الحقيقي المحدد مسبقًا (لنسميه R )، وتكامل يقطع المستوى المركب (لنسميه I ). والتكامل على كامل المسار هو مجموع التكاملات على كل من هذين المسارين.
  • إثبات أن التكامل الذي يعبر المستوى المركب لا يلعب أي دور في المجموع
    إذا أمكن إثبات أن التكامل I يساوي صفرًا، أو إذا كان التكامل الحقيقي المطلوب غير مناسب، فإذا أثبتنا أن التكامل I كما هو موضح أعلاه يميل إلى 0، فإن التكامل على طول R سيميل إلى التكامل حول المحيط R + I.
  • خاتمة
    إذا استطعنا إظهار الخطوة المذكورة أعلاه، فسنتمكن من حساب R مباشرة ، وهو التكامل ذو القيم الحقيقية.

المثال 1

ضع في اعتبارك التكامل -1(x2+1)2دx،{\displaystyle \int _{-\infty }^{\infty }{\frac {1}{\left(x^{2}+1\right)^{2}}}\,dx,}

لتقييم هذا التكامل، ننظر إلى الدالة ذات القيم المركبة و(z)=1(z2+1)2{\displaystyle f(z)={\frac {1}{\left(z^{2}+1\right)^{2}}}}

والتي لها نقاط تفرد عند i و -i . نختار مسارًا محيطًا يحيط بالتكامل ذي القيم الحقيقية، وهنا يكون نصف دائرة قطره على خط الأعداد الحقيقية (يمتد، على سبيل المثال، من -a إلى a ) مناسبًا . لنسمي هذا المسار المحيط C.

هناك طريقتان للمضي قدماً، إما باستخدام صيغة كوشي التكاملية أو بطريقة البواقي:

باستخدام صيغة كوشي التكاملية

لاحظ أن: جو(z)دz=-أأو(z)دz+قوسو(z)دz{\displaystyle \oint _{C}f(z)\,dz=\int _{-a}^{a}f(z)\,dz+\int _{\text{Arc}}f(z)\,dz} هكذا -أأو(z)دz=جو(z)دz-قوسو(z)دz{\displaystyle \int _{-a}^{a}f(z)\,dz=\oint _{C}f(z)\,dz-\int _{\text{Arc}}f(z)\,dz}

علاوة على ذلك، لاحظ أن و(z)=1(z2+1)2=1(z+أنا)2(z-أنا)2.{\displaystyle f(z)={\frac {1}{\left(z^{2}+1\right)^{2}}}={\frac {1}{(z+i)^{2}(z-i)^{2}}}.}

بما أن النقطة الشاذة الوحيدة في المسار هي تلك الموجودة عند النقطة i ، فيمكننا كتابة  و(z)=1(z+أنا)2(z-أنا)2،{\displaystyle f(z)={\frac {\frac {1}{(z+i)^{2}}}{(z-i)^{2}}},}

مما يضع الدالة في صيغة تسمح بتطبيق الصيغة مباشرةً. ثم، باستخدام صيغة كوشي التكاملية، جو(z)دz=ج1(z+أنا)2(z-أنا)2دz=2πأناددz1(z+أنا)2|z=أنا=2πأنا[-2(z+أنا)3]z=أنا=π2{\displaystyle \oint _{C}f(z)\,dz=\oint _{C}{\frac {\frac {1}{(z+i)^{2}}}{(z-i)^{2}}}\,dz=2\pi i\,\left.{\frac {d}{dz}}{\frac {1}{(z+i)^{2}}}\right|_{z=i}=2\pi i\left[{\frac {-2}{(z+i)^{3}}}\right]_{z=i}={\frac {\pi }{2}}}

نأخذ المشتقة الأولى، في الخطوات السابقة، لأن القطب من الرتبة الثانية. أي أن ( z - i ) مرفوع للأس الثاني، لذا نستخدم المشتقة الأولى لـ f ( z ) . لو كان ( z - i ) مرفوعًا للأس الثالث، لاستعملنا المشتقة الثانية وقسمنا على 2!، وهكذا. أما حالة ( z - i ) المرفوع للأس الأول فتُقابل مشتقة من الرتبة الصفرية، وهي f ( z ) نفسها.

نحتاج إلى إثبات أن التكامل على قوس نصف الدائرة يؤول إلى الصفر عندما a → ∞ ، باستخدام مبرهنة التقدير|قوسو(z)دz|مL{\displaystyle \left|\int _{\text{Arc}}f(z)\,dz\right|\leq ML}

حيث M هو حد أعلى لـ | f ( z ) | على طول القوس و L هو طول القوس. الآن، |قوسو(z)دz|أπ(أ2-1)20 مثل أ.{\displaystyle \left|\int _{\text{Arc}}f(z)\,dz\right|\leq {\frac {a\pi }{\left(a^{2}-1\right)^{2}}}\to 0{\text{ as }}a\to \infty .} لذا -1(x2+1)2دx=-و(z)دz=ليمأ+-أأو(z)دz=π2.{\displaystyle \int _{-\infty }^{\infty }{\frac {1}{\left(x^{2}+1\right)^{2}}}\,dx=\int _{-\infty }^{\infty }f(z)\,dz=\lim _{a\to +\infty }\int _{-a}^{a}f(z)\,dz={\frac {\pi }{2}}.\quad \square }

باستخدام طريقة البقايا

لنعتبر متسلسلة لوران للدالة f ( z ) حول i ، وهي النقطة الشاذة الوحيدة التي نحتاج إلى أخذها في الاعتبار. عندئذٍ لدينا و(z)=-14(z-أنا)2+-أنا4(z-أنا)+316+أنا8(z-أنا)+-564(z-أنا)2+{\displaystyle f(z)={\frac {-1}{4(z-i)^{2}}}+{\frac {-i}{4(z-i)}}+{\frac {3}{16}}+{\frac {i}{8}}(z-i)+{\frac {-5}{64}}(z-i)^{2}+\cdots }

(انظر إلى مثال حساب لوران من متسلسلة لوران لاشتقاق هذه المتسلسلة.)

يتضح من الفحص أن الباقي هو i / 4 ، لذا، وفقًا لنظرية الباقي ، لدينا جو(z)دz=ج1(z2+1)2دz=2πأناريسz=أناو(z)=2πأنا(-أنا4)=π2{\displaystyle \oint _{C}f(z)\,dz=\oint _{C}{\frac {1}{\left(z^{2}+1\right)^{2}}}\,dz=2\pi i\,\operatorname {Res} _{z=i}f(z)=2\pi i\left(-{\frac {i}{4}}\right)={\frac {\pi }{2}}\quad \square }

وهكذا نحصل على نفس النتيجة كما في السابق.

ملاحظة كونتورية

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

هذا لا يؤثر على استخدام طريقة البقايا عن طريق السلسلة.

المثال 2 – توزيع كوشي

التكامل -هـأناتxx2+1دx{\displaystyle \int _{-\infty }^{\infty }{\frac {e^{itx}}{x^{2}+1}}\,dx}

محيط
محيط

(الذي يظهر في نظرية الاحتمالات كمضاعف قياسي للدالة المميزة لتوزيع كوشي ) يقاوم تقنيات حساب التفاضل والتكامل الأساسي . سنقوم بتقييمه من خلال التعبير عنه كحد لتكاملات الكفاف على طول الكفاف C الذي يمتد على طول خط الأعداد الحقيقية من −a إلى ثم عكس اتجاه عقارب الساعة على طول نصف دائرة مركزها 0 من a إلى −a . لنفترض أن a أكبر من 1 ، بحيث تكون الوحدة التخيلية i محصورة داخل المنحنى. تكامل الكفاف هو جهـأناتzz2+1دz.{\displaystyle \int _{C}{\frac {e^{itz}}{z^{2}+1}}\,dz.}

بما أن e itz دالة كاملة (لا تحتوي على نقاط شاذة في أي نقطة في المستوى المركب)، فإن هذه الدالة تحتوي على نقاط شاذة فقط عندما يكون المقام + 1 يساوي صفرًا. وبما أن + 1 = ( z + i )( z - i ) ، فإن ذلك يحدث فقط عندما z = i أو z = -i . تقع إحدى هاتين النقطتين فقط في المنطقة المحصورة بهذا المسار. باقي f ( z ) عند z = i هو ليمzأنا(z-أنا)و(z)=ليمzأنا(z-أنا)هـأناتzz2+1=ليمzأنا(z-أنا)هـأناتz(z-أنا)(z+أنا)=ليمzأناهـأناتzz+أنا=هـ-ت2أنا.{\displaystyle \lim _{z\to i}(z-i)f(z)=\lim _{z\to i}(z-i){\frac {e^{itz}}{z^{2}+1}}=\lim _{z\to i}(z-i){\frac {e^{itz}}{(z-i)(z+i)}}=\lim _{z\to i}{\frac {e^{itz}}{z+i}}={\frac {e^{-t}}{2i}}.}

وفقًا لنظرية البواقي ، إذن، لدينا جو(z)دz=2πأناريسz=أناو(z)=2πأناهـ-ت2أنا=πهـ-ت.{\displaystyle \int _{C}f(z)\,dz=2\pi i\operatorname {Res} _{z=i}f(z)=2\pi i{\frac {e^{-t}}{2i}}=\pi e^{-t}.}

يمكن تقسيم المحيط C إلى جزء "مستقيم" وقوس منحني، بحيث مستقيم+قوس=πهـ-ت،{\displaystyle \int _{\text{straight}}+\int _{\text{arc}}=\pi e^{-t},} وبالتالي -أأ=πهـ-ت-قوس.{\displaystyle \int _{-a}^{a}=\pi e^{-t}-\int _{\text{arc}}.}

وفقًا لفرضية جوردان ، إذا كان t > 0 فإنقوسهـأناتzz2+1دz0 مثل أ.{\displaystyle \int _{\text{arc}}{\frac {e^{itz}}{z^{2}+1}}\,dz\rightarrow 0{\mbox{ as }}a\rightarrow \infty .}

لذلك، إذا كانت t > 0 فإن-هـأناتxx2+1دx=πهـ-ت.{\displaystyle \int _{-\infty }^{\infty }{\frac {e^{itx}}{x^{2}+1}}\,dx=\pi e^{-t}.}

تُظهر حجة مماثلة باستخدام قوس يلتف حول i بدلاً من i أنه إذا كان t < 0 فإن-هـأناتxx2+1دx=πهـت،{\displaystyle \int _{-\infty }^{\infty }{\frac {e^{itx}}{x^{2}+1}}\,dx=\pi e^{t},} وأخيراً لدينا هذا: -هـأناتxx2+1دx=πهـ-|ت|.{\displaystyle \int _{-\infty }^{\infty }{\frac {e^{itx}}{x^{2}+1}}\,dx=\pi e^{-|t|}.}

(إذا كانت t = 0 فإن التكامل يخضع مباشرة لطرق حساب التفاضل والتكامل ذات القيم الحقيقية وقيمته هي π .)

مثال 3 – التكاملات المثلثية

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

على سبيل المثال، انظر -ππ11+3(كوست)2دت.{\displaystyle \int _{-\pi }^{\pi }{\frac {1}{1+3(\cos t)^{2}}}\,dt.}

نسعى إلى استبدال z = e it . الآن، تذكر كوست=12(هـأنات+هـ-أنات)=12(z+1z){\displaystyle \cos t={\frac {1}{2}}\left(e^{it}+e^{-it}\right)={\frac {1}{2}}\left(z+{\frac {1}{z}}\right)} و دzدت=أناz، دت=دzأناz.{\displaystyle {\frac {dz}{dt}}=iz,\ dt={\frac {dz}{iz}}.}

باعتبار C دائرة الوحدة، نستبدل لنحصل على: ج11+3(12(z+1z))2دzأناz=ج11+34(z+1z)21أناzدz=ج-أناz+34z(z+1z)2دz=-أناجدzz+34z(z2+2+1z2)=-أناجدzz+34(z3+2z+1z)=-أناجدz34z3+52z+34z=-أناج43z3+10z+3zدz=-4أناجدz3z3+10z+3z=-4أناجz3z4+10z2+3دz=-4أناجz3(z+3أنا)(z-3أنا)(z+أنا3)(z-أنا3)دz=-4أنا3جz(z+3أنا)(z-3أنا)(z+أنا3)(z-أنا3)دz.{\displaystyle {\begin{aligned}\oint _{C}{\frac {1}{1+3\left({\frac {1}{2}}\left(z+{\frac {1}{z}}\right)\right)^{2}}}\,{\frac {dz}{iz}}&=\oint _{C}{\frac {1}{1+{\frac {3}{4}}\left(z+{\frac {1}{z}}\right)^{2}}}{\frac {1}{iz}}\,dz\\&=\oint _{C}{\frac {-i}{z+{\frac {3}{4}}z\left(z+{\frac {1}{z}}\right)^{2}}}\,dz\\&=-i\oint _{C}{\frac {dz}{z+{\frac {3}{4}}z\left(z^{2}+2+{\frac {1}{z^{2}}}\right)}}\\&=-i\oint _{C}{\frac {dz}{z+{\frac {3}{4}}\left(z^{3}+2z+{\frac {1}{z}}\right)}}\\&=-i\oint _{C}{\frac {dz}{{\frac {3}{4}}z^{3}+{\frac {5}{2}}z+{\frac {3}{4z}}}}\\&=-i\oint _{C}{\frac {4}{3z^{3}+10z+{\frac {3}{z}}}}\,dz\\&=-4i\oint _{C}{\frac {dz}{3z^{3}+10z+{\frac {3}{z}}}}\\&=-4i\oint _{C}{\frac {z}{3z^{4}+10z^{2}+3}}\,dz\\&=-4i\oint _{C}{\frac {z}{3\left(z+{\sqrt {3}}i\right)\left(z-{\sqrt {3}}i\right)\left(z+{\frac {i}{\sqrt {3}}}\right)\left(z-{\frac {i}{\sqrt {3}}}\right)}}\,dz\\&=-{\frac {4i}{3}}\oint _{C}{\frac {z}{\left(z+{\sqrt {3}}i\right)\left(z-{\sqrt {3}}i\right)\left(z+{\frac {i}{\sqrt {3}}}\right)\left(z-{\frac {i}{\sqrt {3}}}\right)}}\,dz.\end{aligned}}}

توجد النقاط الشاذة التي يجب أخذها في الاعتبار عند±أنا3.{\displaystyle {\tfrac {\pm i}{\sqrt {3}}}.}لتكن C1 دائرة صغيرة حولأنا3،{\displaystyle {\tfrac {i}{\sqrt {3}}},}و C 2 تكون دائرة صغيرة حول-أنا3.{\displaystyle {\tfrac {-i}{\sqrt {3}}}.}ثم نصل إلى ما يلي: -4أنا3[ج1z(z+3أنا)(z-3أنا)(z+أنا3)z-أنا3دz+ج2z(z+3أنا)(z-3أنا)(z-أنا3)z+أنا3دz]=-4أنا3[2πأنا[z(z+3أنا)(z-3أنا)(z+أنا3)]z=أنا3+2πأنا[z(z+3أنا)(z-3أنا)(z-أنا3)]z=-أنا3]=8π3[أنا3(أنا3+3أنا)(أنا3-3أنا)(أنا3+أنا3)+-أنا3(-أنا3+3أنا)(-أنا3-3أنا)(-أنا3-أنا3)]=8π3[أنا3(43أنا)(-2أنا3)(23أنا)+-أنا3(23أنا)(-43أنا)(-23أنا)]=8π3[أنا3أنا(43)(23)(23)+-أنا3-أنا(23)(43)(23)]=8π3[13(43)(23)(23)+13(23)(43)(23)]=8π3[131633+131633]=8π3[316+316]=π.{\displaystyle {\begin{aligned}&-{\frac {4i}{3}}\left[\oint _{C_{1}}{\frac {\frac {z}{\left(z+{\sqrt {3}}i\right)\left(z-{\sqrt {3}}i\right)\left(z+{\frac {i}{\sqrt {3}}}\right)}}{z-{\frac {i}{\sqrt {3}}}}}\,dz+\oint _{C_{2}}{\frac {\frac {z}{\left(z+{\sqrt {3}}i\right)\left(z-{\sqrt {3}}i\right)\left(z-{\frac {i}{\sqrt {3}}}\right)}}{z+{\frac {i}{\sqrt {3}}}}}\,dz\right]\\={}&-{\frac {4i}{3}}\left[2\pi i\left[{\frac {z}{\left(z+{\sqrt {3}}i\right)\left(z-{\sqrt {3}}i\right)\left(z+{\frac {i}{\sqrt {3}}}\right)}}\right]_{z={\frac {i}{\sqrt {3}}}}+2\pi i\left[{\frac {z}{\left(z+{\sqrt {3}}i\right)\left(z-{\sqrt {3}}i\right)\left(z-{\frac {i}{\sqrt {3}}}\right)}}\right]_{z=-{\frac {i}{\sqrt {3}}}}\right]\\={}&{\frac {8\pi }{3}}\left[{\frac {\frac {i}{\sqrt {3}}}{\left({\frac {i}{\sqrt {3}}}+{\sqrt {3}}i\right)\left({\frac {i}{\sqrt {3}}}-{\sqrt {3}}i\right)\left({\frac {i}{\sqrt {3}}}+{\frac {i}{\sqrt {3}}}\right)}}+{\frac {-{\frac {i}{\sqrt {3}}}}{\left(-{\frac {i}{\sqrt {3}}}+{\sqrt {3}}i\right)\left(-{\frac {i}{\sqrt {3}}}-{\sqrt {3}}i\right)\left(-{\frac {i}{\sqrt {3}}}-{\frac {i}{\sqrt {3}}}\right)}}\right]\\={}&{\frac {8\pi }{3}}\left[{\frac {\frac {i}{\sqrt {3}}}{\left({\frac {4}{\sqrt {3}}}i\right)\left(-{\frac {2}{i{\sqrt {3}}}}\right)\left({\frac {2}{{\sqrt {3}}i}}\right)}}+{\frac {-{\frac {i}{\sqrt {3}}}}{\left({\frac {2}{\sqrt {3}}}i\right)\left(-{\frac {4}{\sqrt {3}}}i\right)\left(-{\frac {2}{\sqrt {3}}}i\right)}}\right]\\={}&{\frac {8\pi }{3}}\left[{\frac {\frac {i}{\sqrt {3}}}{i\left({\frac {4}{\sqrt {3}}}\right)\left({\frac {2}{\sqrt {3}}}\right)\left({\frac {2}{\sqrt {3}}}\right)}}+{\frac {-{\frac {i}{\sqrt {3}}}}{-i\left({\frac {2}{\sqrt {3}}}\right)\left({\frac {4}{\sqrt {3}}}\right)\left({\frac {2}{\sqrt {3}}}\right)}}\right]\\={}&{\frac {8\pi }{3}}\left[{\frac {\frac {1}{\sqrt {3}}}{\left({\frac {4}{\sqrt {3}}}\right)\left({\frac {2}{\sqrt {3}}}\right)\left({\frac {2}{\sqrt {3}}}\right)}}+{\frac {\frac {1}{\sqrt {3}}}{\left({\frac {2}{\sqrt {3}}}\right)\left({\frac {4}{\sqrt {3}}}\right)\left({\frac {2}{\sqrt {3}}}\right)}}\right]\\={}&{\frac {8\pi }{3}}\left[{\frac {\frac {1}{\sqrt {3}}}{\frac {16}{3{\sqrt {3}}}}}+{\frac {\frac {1}{\sqrt {3}}}{\frac {16}{3{\sqrt {3}}}}}\right]\\={}&{\frac {8\pi }{3}}\left[{\frac {3}{16}}+{\frac {3}{16}}\right]\\={}&\pi .\end{aligned}}}

مثال 3أ - التكاملات المثلثية، الإجراء العام

يمكن تطبيق الطريقة المذكورة أعلاه على جميع التكاملات من النوع 02πP(الخطيئة(ت)،الخطيئة(2ت)،...،كوس(ت)،كوس(2ت)،...)سؤال(الخطيئة(ت)،الخطيئة(2ت)،...،كوس(ت)،كوس(2ت)،...)دت{\displaystyle \int _{0}^{2\pi }{\frac {P{\big (}\sin(t),\sin(2t),\ldots ,\cos(t),\cos(2t),\ldots {\big )}}{Q{\big (}\sin(t),\sin(2t),\ldots ,\cos(t),\cos(2t),\ldots {\big )}}}\,dt} حيث P و Q كثيرتا حدود، أي يتم تكامل دالة كسرية بدلالة الدوال المثلثية. لاحظ أن حدود التكامل قد تكون π و ، كما في المثال السابق، أو أي زوج آخر من نقاط النهاية يفصل بينهما .

يكمن السر في استخدام التعويض z = e it حيث dz = ie it dt وبالتالي 1أناzدz=دت.{\displaystyle {\frac {1}{iz}}\,dz=dt.}

هذا الاستبدال يحوّل الفترة [ 0، 2π ] إلى دائرة الوحدة. علاوة على ذلك، الخطيئة(كت)=هـأناكت-هـ-أناكت2أنا=zك-z-ك2أنا{\displaystyle \sin(kt)={\frac {e^{ikt}-e^{-ikt}}{2i}}={\frac {z^{k}-z^{-k}}{2i}}} و كوس(كت)=هـأناكت+هـ-أناكت2=zك+z-ك2{\displaystyle \cos(kt)={\frac {e^{ikt}+e^{-ikt}}{2}}={\frac {z^{k}+z^{-k}}{2}}} وبذلك ينتج عن التعويض دالة كسرية f ( z ) في z ، ويصبح التكامل |z|=1و(z)1أناzدz{\displaystyle \oint _{|z|=1}f(z){\frac {1}{iz}}\,dz} والتي يتم حسابها بدورها عن طريق جمع بقايا f ( z ) 1 / iz داخل دائرة الوحدة.

توضح الصورة على اليمين هذا الأمر بالنسبة لـ أنا=0π211+(الخطيئةت)2دت،{\displaystyle I=\int _{0}^{\frac {\pi }{2}}{\frac {1}{1+(\sin t)^{2}}}\,dt,} والتي سنقوم بحسابها الآن. الخطوة الأولى هي إدراك ذلك. أنا=1402π11+(الخطيئةت)2دت.{\displaystyle I={\frac {1}{4}}\int _{0}^{2\pi }{\frac {1}{1+(\sin t)^{2}}}\,dt.}

ينتج عن الاستبدال 14|z|=14أناzz4-6z2+1دz=|z|=1أناzz4-6z2+1دz.{\displaystyle {\frac {1}{4}}\oint _{|z|=1}{\frac {4iz}{z^{4}-6z^{2}+1}}\,dz=\oint _{|z|=1}{\frac {iz}{z^{4}-6z^{2}+1}}\,dz.}

تقع أقطاب هذه الدالة عند 1 ± √2 و −1 ± √2 . من بين هذه الأقطاب، يقع 1 + √2 و −1√2 خارج دائرة الوحدة (الموضحة باللون الأحمر، وليست وفقًا للمقياس)، بينما يقع 1√2 و −1 + √2 داخل دائرة الوحدة ( الموضحة باللون الأزرق). والباقي المقابل يساوي i √2 / 16 ، وبالتالي فإن قيمة التكامل هي أنا=2πأنا2(-216أنا)=π24.{\displaystyle I=2\pi i\;2\left(-{\frac {\sqrt {2}}{16}}i\right)=\pi {\frac {\sqrt {2}}{4}}.}

المثال 4 - تقليم الفروع

لننظر إلى التكامل الحقيقي 0xx2+6x+8دx.{\displaystyle \int _{0}^{\infty }{\frac {\sqrt {x}}{x^{2}+6x+8}}\,dx.}

يمكننا البدء بصياغة التكامل المركب جzz2+6z+8دz=أنا.{\displaystyle \int _{C}{\frac {\sqrt {z}}{z^{2}+6z+8}}\,dz=I.}

يمكننا استخدام صيغة كوشي التكاملية أو نظرية البواقي مرة أخرى للحصول على البواقي المطلوبة. مع ذلك، من المهم ملاحظة أن z 1/2 = e (Log z )/2 ، لذا فإن z 1/2 له قطع فرعي . يؤثر هذا على اختيارنا للمسار C. عادةً ما يُعرَّف القطع الفرعي اللوغاريتمي بأنه المحور الحقيقي السالب، ولكن هذا يجعل حساب التكامل أكثر تعقيدًا بعض الشيء، لذا نُعرِّفه بأنه المحور الحقيقي الموجب.

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

لاحظ أن z = −2 و z = −4 تقعان داخل الدائرة الكبيرة. وهما القطبان المتبقيان، واللذان يمكن استنتاجهما بتحليل مقام الدالة التكاملية. وقد تم تجنب نقطة التفرع عند z = 0 بالالتفاف حول نقطة الأصل.

لتكن γ الدائرة الصغرى بنصف قطر ε ، و Γ الدائرة الكبرى بنصف قطر R ، عندئذٍ ج=εR+Γ+Rε+γ.{\displaystyle \int _{C}=\int _{\varepsilon }^{R}+\int _{\Gamma }+\int _{R}^{\varepsilon }+\int _{\gamma }.}

يمكن إثبات أن التكاملات على Γ و γ تؤول إلى الصفر عندما ε → 0 و R → ∞ ، وذلك باستخدام حجة تقديرية سابقة، مما يترك حدين. الآن، بما أن z 1/2 = e (Log z )/2 ، على المسار خارج قطع الفرع، فقد كسبنا 2π في الموضع على طول γ . (بحسب متطابقة أويلر ، يمثل e i π متجه الوحدة ، والذي بالتالي يكون لوغاريتمه π . هذا π هو المقصود بموضع z . معامل 1/2 يجبرنا على استخدام . )Rεzz2+6z+8دz=Rεهـ12سجلzz2+6z+8دz=Rεهـ12(سجل|z|+أناargz)z2+6z+8دz=Rεهـ12سجل|z|هـ12(2πأنا)z2+6z+8دz=Rεهـ12سجل|z|هـπأناz2+6z+8دz=Rε-zz2+6z+8دz=εRzz2+6z+8دz.{\displaystyle {\begin{aligned}\int _{R}^{\varepsilon }{\frac {\sqrt {z}}{z^{2}+6z+8}}\,dz&=\int _{R}^{\varepsilon }{\frac {e^{{\frac {1}{2}}\operatorname {Log} z}}{z^{2}+6z+8}}\,dz\\[6pt]&=\int _{R}^{\varepsilon }{\frac {e^{{\frac {1}{2}}(\log |z|+i\arg {z})}}{z^{2}+6z+8}}\,dz\\[6pt]&=\int _{R}^{\varepsilon }{\frac {e^{{\frac {1}{2}}\log |z|}e^{{\frac {1}{2}}(2\pi i)}}{z^{2}+6z+8}}\,dz\\[6pt]&=\int _{R}^{\varepsilon }{\frac {e^{{\frac {1}{2}}\log |z|}e^{\pi i}}{z^{2}+6z+8}}\,dz\\[6pt]&=\int _{R}^{\varepsilon }{\frac {-{\sqrt {z}}}{z^{2}+6z+8}}\,dz\\[6pt]&=\int _{\varepsilon }^{R}{\frac {\sqrt {z}}{z^{2}+6z+8}}\,dz.\end{aligned}}}

لذلك: جzz2+6z+8دz=20xx2+6x+8دx.{\displaystyle \int _{C}{\frac {\sqrt {z}}{z^{2}+6z+8}}\,dz=2\int _{0}^{\infty }{\frac {\sqrt {x}}{x^{2}+6x+8}}\,dx.}

باستخدام نظرية البواقي أو صيغة كوشي التكاملية (بعد استخدام طريقة الكسور الجزئية لاستنتاج مجموع تكاملين بسيطين على مستوى الكفاف) نحصل على πأنا(أنا2-أنا)=0xx2+6x+8دx=π(1-12).{\displaystyle \pi i\left({\frac {i}{\sqrt {2}}}-i\right)=\int _{0}^{\infty }{\frac {\sqrt {x}}{x^{2}+6x+8}}\,dx=\pi \left(1-{\frac {1}{\sqrt {2}}}\right).\quad \square }

المثال 5 – مربع اللوغاريتم

يتناول هذا القسم نوعًا من التكاملات التي 0سجلx(1+x2)2دx{\displaystyle \int _{0}^{\infty }{\frac {\log x}{\left(1+x^{2}\right)^{2}}}\,dx} هذا مثال.

لحساب هذا التكامل، يتم استخدام الدالة و(z)=(سجلz1+z2)2{\displaystyle f(z)=\left({\frac {\log z}{1+z^{2}}}\right)^{2}} وفرع اللوغاريتم المقابل لـ −π < arg z ≤ π .

سنحسب تكامل الدالة f ( z ) على طول محيط ثقب المفتاح الموضح على اليمين. وكما اتضح، فإن هذا التكامل هو مضاعف للتكامل الأولي الذي نرغب في حسابه، وبحسب نظرية البواقي لكوشي، لدينا (R+م+شمال+ر)و(z)دz= 2πأنا(ريسz=أناو(z)+ريسz=-أناو(z))= 2πأنا(-π4+116أناπ2-π4-116أناπ2)= -أناπ2.{\displaystyle {\begin{aligned}\left(\int _{R}+\int _{M}+\int _{N}+\int _{r}\right)f(z)\,dz=&\ 2\pi i{\big (}\operatorname {Res} _{z=i}f(z)+\operatorname {Res} _{z=-i}f(z){\big )}\\=&\ 2\pi i\left(-{\frac {\pi }{4}}+{\frac {1}{16}}i\pi ^{2}-{\frac {\pi }{4}}-{\frac {1}{16}}i\pi ^{2}\right)\\=&\ -i\pi ^{2}.\end{aligned}}}

ليكن R نصف قطر الدائرة الكبيرة، و r نصف قطر الدائرة الصغيرة. سنرمز للخط العلوي بـ M ، وللخط السفلي بـ N. وكما في السابق، نأخذ النهاية عندما R → ∞ و r → 0. تتلاشى مساهمات الدائرتين. على سبيل المثال، لدينا الحد الأعلى التالي باستخدام مبرهنة ML : |Rو(z)دz|2πR(سجلR)2+π2(R2-1)20.{\displaystyle \left|\int _{R}f(z)\,dz\right|\leq 2\pi R{\frac {(\log R)^{2}+\pi ^{2}}{\left(R^{2}-1\right)^{2}}}\to 0.}

لحساب مساهمات M و N ، نضع z = − x + على M و z = − x على N ، مع 0 < x < ∞ : -أناπ2=(R+م+شمال+ر)و(z)دz=(م+شمال)و(z)دzR،ر يختفي=-0(سجل(-x+أناε)1+(-x+أناε)2)2دx-0(سجل(-x-أناε)1+(-x-أناε)2)2دx=0(سجل(-x+أناε)1+(-x+أناε)2)2دx-0(سجل(-x-أناε)1+(-x-أناε)2)2دx=0(سجلx+أناπ1+x2)2دx-0(سجلx-أناπ1+x2)2دxε0=0(سجلx+أناπ)2-(سجلx-أناπ)2(1+x2)2دx=04πأناسجلx(1+x2)2دx=4πأنا0سجلx(1+x2)2دx{\displaystyle {\begin{aligned}-i\pi ^{2}&=\left(\int _{R}+\int _{M}+\int _{N}+\int _{r}\right)f(z)\,dz\\[6pt]&=\left(\int _{M}+\int _{N}\right)f(z)\,dz&&\int _{R},\int _{r}{\mbox{ vanish}}\\[6pt]&=-\int _{\infty }^{0}\left({\frac {\log(-x+i\varepsilon )}{1+(-x+i\varepsilon )^{2}}}\right)^{2}\,dx-\int _{0}^{\infty }\left({\frac {\log(-x-i\varepsilon )}{1+(-x-i\varepsilon )^{2}}}\right)^{2}\,dx\\[6pt]&=\int _{0}^{\infty }\left({\frac {\log(-x+i\varepsilon )}{1+(-x+i\varepsilon )^{2}}}\right)^{2}\,dx-\int _{0}^{\infty }\left({\frac {\log(-x-i\varepsilon )}{1+(-x-i\varepsilon )^{2}}}\right)^{2}\,dx\\[6pt]&=\int _{0}^{\infty }\left({\frac {\log x+i\pi }{1+x^{2}}}\right)^{2}\,dx-\int _{0}^{\infty }\left({\frac {\log x-i\pi }{1+x^{2}}}\right)^{2}\,dx&&\varepsilon \to 0\\&=\int _{0}^{\infty }{\frac {(\log x+i\pi )^{2}-(\log x-i\pi )^{2}}{\left(1+x^{2}\right)^{2}}}\,dx\\[6pt]&=\int _{0}^{\infty }{\frac {4\pi i\log x}{\left(1+x^{2}\right)^{2}}}\,dx\\[6pt]&=4\pi i\int _{0}^{\infty }{\frac {\log x}{\left(1+x^{2}\right)^{2}}}\,dx\end{aligned}}} مما يعطي 0سجلx(1+x2)2دx=-π4.{\displaystyle \int _{0}^{\infty }{\frac {\log x}{\left(1+x^{2}\right)^{2}}}\,dx=-{\frac {\pi }{4}}.}

المثال 6 - اللوغاريتمات والباقي عند اللانهاية

نسعى إلى التقييم أنا=03x34(3-x)145-xدx.{\displaystyle I=\int _{0}^{3}{\frac {x^{\frac {3}{4}}(3-x)^{\frac {1}{4}}}{5-x}}\,dx.}

يتطلب هذا دراسة متأنية لـ و(z)=z34(3-z)14.{\displaystyle f(z)=z^{\frac {3}{4}}(3-z)^{\frac {1}{4}}.}

سنقوم بإنشاء الدالة f ( z ) بحيث يكون لها فرع قطع على الفترة [ 0, 3 ] ، كما هو موضح باللون الأحمر في الرسم التخطيطي. وللقيام بذلك، نختار فرعين من اللوغاريتم، ونضع z34=خبرة(34سجلz)أين -πargz<π{\displaystyle z^{\frac {3}{4}}=\exp \left({\frac {3}{4}}\log z\right)\quad {\mbox{where }}-\pi \leq \arg z<\pi } و (3-z)14=خبرة(14سجل(3-z))أين 0arg(3-z)<2π.{\displaystyle (3-z)^{\frac {1}{4}}=\exp \left({\frac {1}{4}}\log(3-z)\right)\quad {\mbox{where }}0\leq \arg(3-z)<2\pi .}

وبالتالي ، فإن قطع z ¾ هو [−∞, 0 ] ، وقطع (3 − z ) ¼ هو [−∞, 3 ] . من السهل ملاحظة أن قطع حاصل ضربهما، أي f ( z ) ، هو [0, 3] ، لأن f ( z ) متصلة فعليًا عبر (−∞, 0) . هذا لأنه عندما z = −r < 0 ونقترب من القطع من الأعلى، فإن f ( z ) تأخذ القيمة ر34هـ34πأنا(3+ر)14هـ24πأنا=ر34(3+ر)14هـ54πأنا.{\displaystyle r^{\frac {3}{4}}e^{{\frac {3}{4}}\pi i}(3+r)^{\frac {1}{4}}e^{{\frac {2}{4}}\pi i}=r^{\frac {3}{4}}(3+r)^{\frac {1}{4}}e^{{\frac {5}{4}}\pi i}.}

عندما نقترب من الأسفل، فإن f ( z ) تأخذ القيمة ر34هـ-34πأنا(3+ر)14هـ04πأنا=ر34(3+ر)14هـ-34πأنا.{\displaystyle r^{\frac {3}{4}}e^{-{\frac {3}{4}}\pi i}(3+r)^{\frac {1}{4}}e^{{\frac {0}{4}}\pi i}=r^{\frac {3}{4}}(3+r)^{\frac {1}{4}}e^{-{\frac {3}{4}}\pi i}.}

لكن هـ-34πأنا=هـ54πأنا،{\displaystyle e^{-{\frac {3}{4}}\pi i}=e^{{\frac {5}{4}}\pi i},}

بحيث يكون لدينا استمرارية عبر القطع. ويتضح هذا في الرسم التوضيحي، حيث تم تسمية الدائرتين السوداوين الموجهتين بالقيمة المقابلة لوسيط اللوغاريتم المستخدم في z 3 4 و (3 − z ) 1/4 .

سنستخدم الخط المحيطي الموضح باللون الأخضر في الرسم التخطيطي. وللقيام بذلك، يجب علينا حساب قيمة f ( z ) على طول القطع المستقيمة أعلى وأسفل القطع مباشرةً.

لنفترض أن z = r (في النهاية، أي عندما تتقلص الدائرتان الخضراوان إلى نصف قطر صفر)، حيث 0 ≤ r ≤ 3. على طول القطعة العلوية، نجد أن f ( z ) لها القيمة التالية: ر34هـ04πأنا(3-ر)14هـ24πأنا=أنار34(3-ر)14{\displaystyle r^{\frac {3}{4}}e^{{\frac {0}{4}}\pi i}(3-r)^{\frac {1}{4}}e^{{\frac {2}{4}}\pi i}=ir^{\frac {3}{4}}(3-r)^{\frac {1}{4}}} وعلى طول الجزء السفلي، ر34هـ04πأنا(3-ر)14هـ04πأنا=ر34(3-ر)14.{\displaystyle r^{\frac {3}{4}}e^{{\frac {0}{4}}\pi i}(3-r)^{\frac {1}{4}}e^{{\frac {0}{4}}\pi i}=r^{\frac {3}{4}}(3-r)^{\frac {1}{4}}.}

ويترتب على ذلك أن تكامل f ( z ) / 5 − z على طول الجزء العلوي هو iI في النهاية، وعلى طول الجزء السفلي، I .

إذا استطعنا إثبات أن التكاملات على طول الدائرتين الخضراوين تتلاشى في النهاية، فسنحصل أيضًا على قيمة I ، وفقًا لنظرية البواقي لكوشي . لنفترض أن نصف قطر الدائرتين الخضراوين هو ρ ، حيث ρ < 0.001 و ρ → 0 ، ونطبق متباينة ML . بالنسبة للدائرة CL على اليسار ، نجد |جLو(z)5-zدz|2πρρ343.001144.999يا(ρ74)0.{\displaystyle \left|\int _{C_{\mathrm {L} }}{\frac {f(z)}{5-z}}dz\right|\leq 2\pi \rho {\frac {\rho ^{\frac {3}{4}}3.001^{\frac {1}{4}}}{4.999}}\in {\mathcal {O}}\left(\rho ^{\frac {7}{4}}\right)\to 0.}

وبالمثل، بالنسبة للدائرة C R على اليمين، لدينا |جRو(z)5-zدz|2πρ3.00134ρ141.999يا(ρ54)0.{\displaystyle \left|\int _{C_{\mathrm {R} }}{\frac {f(z)}{5-z}}dz\right|\leq 2\pi \rho {\frac {3.001^{\frac {3}{4}}\rho ^{\frac {1}{4}}}{1.999}}\in {\mathcal {O}}\left(\rho ^{\frac {5}{4}}\right)\to 0.}

باستخدام نظرية كوشي للبواقي ، لدينا (-أنا+1)أنا=-2πأنا(ريسz=5و(z)5-z+ريسz=و(z)5-z).{\displaystyle (-i+1)I=-2\pi i\left(\operatorname {Res} _{z=5}{\frac {f(z)}{5-z}}+\operatorname {Res} _{z=\infty }{\frac {f(z)}{5-z}}\right).} حيث تُعزى الإشارة السالبة إلى اتجاه عقارب الساعة حول البواقي. باستخدام فرع اللوغاريتم من السابق، من الواضح ريسz=5و(z)5-z=-534هـ14سجل(-2).{\displaystyle \operatorname {Res} _{z=5}{\frac {f(z)}{5-z}}=-5^{\frac {3}{4}}e^{{\frac {1}{4}}\log(-2)}.}

يُظهر الرسم التخطيطي القطب باللون الأزرق. وتُبسط القيمة إلى -534هـ14(سجل2+πأنا)=-هـ14πأنا534214.{\displaystyle -5^{\frac {3}{4}}e^{{\frac {1}{4}}(\log 2+\pi i)}=-e^{{\frac {1}{4}}\pi i}5^{\frac {3}{4}}2^{\frac {1}{4}}.}

نستخدم الصيغة التالية لحساب الباقي عند اللانهاية: ريسz=ح(z)=ريسz=0(-1z2ح(1z)).{\displaystyle \operatorname {Res} _{z=\infty }h(z)=\operatorname {Res} _{z=0}\left(-{\frac {1}{z^{2}}}h\left({\frac {1}{z}}\right)\right).}

بالتعويض، نجد 15-1z=-z(1+5z+52z2+53z3+){\displaystyle {\frac {1}{5-{\frac {1}{z}}}}=-z\left(1+5z+5^{2}z^{2}+5^{3}z^{3}+\cdots \right)} و (1z3(3-1z))14=1z(3z-1)14=1zهـ14πأنا(1-3z)14،{\displaystyle \left({\frac {1}{z^{3}}}\left(3-{\frac {1}{z}}\right)\right)^{\frac {1}{4}}={\frac {1}{z}}(3z-1)^{\frac {1}{4}}={\frac {1}{z}}e^{{\frac {1}{4}}\pi i}(1-3z)^{\frac {1}{4}},} حيث استخدمنا حقيقة أن −1 = e π i للفرع الثاني من اللوغاريتم. ثم نطبق مفكوك ذات الحدين، فنحصل على 1zهـ14πأنا(1-(1/41)3z+(1/42)32z2-(1/43)33z3+).{\displaystyle {\frac {1}{z}}e^{{\frac {1}{4}}\pi i}\left(1-{1/4 \choose 1}3z+{1/4 \choose 2}3^{2}z^{2}-{1/4 \choose 3}3^{3}z^{3}+\cdots \right).}

والخلاصة هي أن ريسz=و(z)5-z=هـ14πأنا(5-34)=هـ14πأنا174.{\displaystyle \operatorname {Res} _{z=\infty }{\frac {f(z)}{5-z}}=e^{{\frac {1}{4}}\pi i}\left(5-{\frac {3}{4}}\right)=e^{{\frac {1}{4}}\pi i}{\frac {17}{4}}.}

وأخيرًا، يترتب على ذلك أن قيمة I هي أنا=2πأناهـ14πأنا-1+أنا(174-534214)=2π2-12(174-534214){\displaystyle I=2\pi i{\frac {e^{{\frac {1}{4}}\pi i}}{-1+i}}\left({\frac {17}{4}}-5^{\frac {3}{4}}2^{\frac {1}{4}}\right)=2\pi 2^{-{\frac {1}{2}}}\left({\frac {17}{4}}-5^{\frac {3}{4}}2^{\frac {1}{4}}\right)} مما ينتج عنه أنا=π22(17-534294)=π22(17-4034).{\displaystyle I={\frac {\pi }{2{\sqrt {2}}}}\left(17-5^{\frac {3}{4}}2^{\frac {9}{4}}\right)={\frac {\pi }{2{\sqrt {2}}}}\left(17-40^{\frac {3}{4}}\right).}

التقييم باستخدام نظرية البواقي

باستخدام نظرية البواقي ، يمكننا حساب التكاملات الكفافية المغلقة. فيما يلي أمثلة على حساب التكاملات الكفافية باستخدام نظرية البواقي.

باستخدام نظرية البقايا، دعونا نحسب هذا التكامل الكفافي. جهـzz3دz{\displaystyle \oint _{C}{\frac {e^{z}}{z^{3}}}\,dz}

تذكر أن نظرية الباقي تنص على جو(z)دz=2πأناريس(و،أك)،{\displaystyle \oint _{C}f(z)dz=2\pi i\cdot \sum \operatorname {Res} (f,a_{k}),} أينريس{\displaystyle \operatorname {Res} }هو بقاياو(z){\displaystyle f(z)}وأك{\displaystyle a_{k}}هي نقاط التفرد فيو(z){\displaystyle f(z)}يقع داخل المحيطج{\displaystyle C}(مع عدم وجود أي منها يقع مباشرة علىج{\displaystyle C}).

و(z){\displaystyle f(z)}له قطب واحد فقط،0{\displaystyle 0}ومن ذلك، نحدد أن بقاياو(z){\displaystyle f(z)}يكون12{\displaystyle {\tfrac {1}{2}}}جو(z)دz=جهـzz3دz=2πأناريسz=0و(z)=2πأناريسz=0هـzz3=2πأنا12=πأنا{\displaystyle {\begin{aligned}\oint _{C}f(z)dz&=\oint _{C}{\frac {e^{z}}{z^{3}}}dz\\&=2\pi i\cdot \operatorname {Res} _{z=0}f(z)\\&=2\pi i\operatorname {Res} _{z=0}{\frac {e^{z}}{z^{3}}}\\&=2\pi i\cdot {\frac {1}{2}}\\&=\pi i\end{aligned}}}

وبالتالي، باستخدام نظرية البواقي ، يمكننا تحديد ما يلي: جهـzz3دz=πأنا.{\displaystyle \oint _{C}{\frac {e^{z}}{z^{3}}}dz=\pi i.}

تكاملات الكفاف متعددة المتغيرات

لحل تكاملات الكفاف متعددة المتغيرات (أي تكاملات السطح ، وتكاملات الحجم المركب ، والتكاملات ذات الرتب العليا )، يجب علينا استخدام نظرية التباعد . في الوقت الحالي، لنفترض أن{\displaystyle \nabla \cdot }يمكن استبدالها بـdiv{\displaystyle \operatorname {div} }سيمثل هذان العنصران تباعد حقل المتجهات المشار إليه بـF{\displaystyle \mathbf {F} }تنص هذه النظرية على ما يلي: يونdiv(F)دV=يون-1FندS{\displaystyle \underbrace {\int \cdots \int _{U}} _{n}\operatorname {div} (\mathbf {F} )\,dV=\underbrace {\oint \cdots \oint _{\partial U}} _{n-1}\mathbf {F} \cdot \mathbf {n} \,dS}

بالإضافة إلى ذلك، نحتاج أيضًا إلى التقييمF{\displaystyle \nabla \cdot \mathbf {F} }أينF{\displaystyle \nabla \cdot \mathbf {F} }هي صيغة بديلة لـdiv(F){\displaystyle \operatorname {div} (\mathbf {F} )}يمكن وصف تباعد أي بُعد على النحو التالي :div(F)=F=(u،x،y،z،...)(Fu،Fx،Fy،Fz،...)=(Fuu+Fxx+Fyy+Fzz+){\displaystyle {\begin{aligned}\operatorname {div} (\mathbf {F} )&=\nabla \cdot \mathbf {F} \\&=\left({\frac {\partial }{\partial u}},{\frac {\partial }{\partial x}},{\frac {\partial }{\partial y}},{\frac {\partial }{\partial z}},\dots \right)\cdot (F_{u},F_{x},F_{y},F_{z},\dots )\\&=\left({\frac {\partial F_{u}}{\partial u}}+{\frac {\partial F_{x}}{\partial x}}+{\frac {\partial F_{y}}{\partial y}}+{\frac {\partial F_{z}}{\partial z}}+\cdots \right)\end{aligned}}}

المثال 1

ليكن حقل المتجهاتF=الخطيئة(2x)هـx+الخطيئة(2y)هـy+الخطيئة(2z)هـz{\displaystyle \mathbf {F} =\sin(2x)\mathbf {e} _{x}+\sin(2y)\mathbf {e} _{y}+\sin(2z)\mathbf {e} _{z}}وأن تكون محصورة بما يلي 0x10y3-1z4{\displaystyle {0\leq x\leq 1}\quad {0\leq y\leq 3}\quad {-1\leq z\leq 4}}

سيتم إعداد التكامل المزدوج المقابل على النحو التالي:

{\displaystyle }\oiintS{\displaystyle {\scriptstyle S}}FندS{\displaystyle \mathbf {F} \cdot \mathbf {n} \,dS}

نقوم الآن بالتقييمF{\displaystyle \nabla \cdot \mathbf {F} }في هذه الأثناء، قم بإعداد التكامل الثلاثي المقابل: =V(Fxx+Fyy+Fzz)دV=V(الخطيئة(2x)x+الخطيئة(2y)y+الخطيئة(2z)z)دV=V2(كوس(2x)+كوس(2y)+كوس(2z))دV=0103-142(كوس(2x)+كوس(2y)+كوس(2z))دxدyدz=0103(10كوس(2y)+الخطيئة(8)+الخطيئة(2)+10كوس(z))دyدz=01(30كوس(2z)+3الخطيئة(2)+3الخطيئة(8)+5الخطيئة(6))دz=18الخطيئة(2)+3الخطيئة(8)+5الخطيئة(6){\displaystyle {\begin{aligned}&=\iiint _{V}\left({\frac {\partial F_{x}}{\partial x}}+{\frac {\partial F_{y}}{\partial y}}+{\frac {\partial F_{z}}{\partial z}}\right)dV\\[6pt]&=\iiint _{V}\left({\frac {\partial \sin(2x)}{\partial x}}+{\frac {\partial \sin(2y)}{\partial y}}+{\frac {\partial \sin(2z)}{\partial z}}\right)dV\\[6pt]&=\iiint _{V}2\left(\cos(2x)+\cos(2y)+\cos(2z)\right)dV\\[6pt]&=\int _{0}^{1}\int _{0}^{3}\int _{-1}^{4}2(\cos(2x)+\cos(2y)+\cos(2z))\,dx\,dy\,dz\\[6pt]&=\int _{0}^{1}\int _{0}^{3}(10\cos(2y)+\sin(8)+\sin(2)+10\cos(z))\,dy\,dz\\[6pt]&=\int _{0}^{1}(30\cos(2z)+3\sin(2)+3\sin(8)+5\sin(6))\,dz\\[6pt]&=18\sin(2)+3\sin(8)+5\sin(6)\end{aligned}}}

المثال 2

ليكن حقل المتجهاتF=u4هـu+x5هـx+y6هـy+z-3هـz{\displaystyle \mathbf {F} =u^{4}\mathbf {e} _{u}+x^{5}\mathbf {e} _{x}+y^{6}\mathbf {e} _{y}+z^{-3}\mathbf {e} _{z}}لاحظ أن هناك 4 معلمات في هذه الحالة. ليكن حقل المتجهات هذا محدودًا بما يلي: 0x1-10y2π4z5-1u3{\displaystyle {0\leq x\leq 1}\quad {-10\leq y\leq 2\pi }\quad {4\leq z\leq 5}\quad {-1\leq u\leq 3}}

لتقييم ذلك، يجب علينا استخدام نظرية التباعد كما ذكرنا سابقًا، ويجب علينا تقييمF{\displaystyle \nabla \cdot \mathbf {F} }. يتركدV=دxدyدzدu{\displaystyle dV=dx\,dy\,dz\,du}

\oiiintS{\displaystyle {\scriptstyle S}}FندS{\displaystyle \mathbf {F} \cdot \mathbf {n} \,dS}

=V(Fuu+Fxx+Fyy+Fzz)دV=V(u4u+x5x+y6y+z-3z)دV=V4u3z4+5x4z4+5y4z4-3z4دV=V4u3z4+5x4z4+5y4z4-3z4دV=01-102π45-134u3z4+5x4z4+5y4z4-3z4دV=01-102π45(4(3u4z3+3y6+91z3+3)3z3)دyدzدu=01-102π(4u4+74344021+4z3)دzدu=01(-12π2+1486880π21+8πu4+40u4+3717200211050)دu=3717284211050+14869136π3-105210π2576468.77{\displaystyle {\begin{aligned}&=\iiiint _{V}\left({\frac {\partial F_{u}}{\partial u}}+{\frac {\partial F_{x}}{\partial x}}+{\frac {\partial F_{y}}{\partial y}}+{\frac {\partial F_{z}}{\partial z}}\right)\,dV\\[6pt]&=\iiiint _{V}\left({\frac {\partial u^{4}}{\partial u}}+{\frac {\partial x^{5}}{\partial x}}+{\frac {\partial y^{6}}{\partial y}}+{\frac {\partial z^{-3}}{\partial z}}\right)\,dV\\[6pt]&=\iiiint _{V}{\frac {4u^{3}z^{4}+5x^{4}z^{4}+5y^{4}z^{4}-3}{z^{4}}}\,dV\\[6pt]&=\iiiint _{V}{\frac {4u^{3}z^{4}+5x^{4}z^{4}+5y^{4}z^{4}-3}{z^{4}}}\,dV\\[6pt]&=\int _{0}^{1}\int _{-10}^{2\pi }\int _{4}^{5}\int _{-1}^{3}{\frac {4u^{3}z^{4}+5x^{4}z^{4}+5y^{4}z^{4}-3}{z^{4}}}\,dV\\[6pt]&=\int _{0}^{1}\int _{-10}^{2\pi }\int _{4}^{5}\left({\frac {4(3u^{4}z^{3}+3y^{6}+91z^{3}+3)}{3z^{3}}}\right)\,dy\,dz\,du\\[6pt]&=\int _{0}^{1}\int _{-10}^{2\pi }\left(4u^{4}+{\frac {743440}{21}}+{\frac {4}{z^{3}}}\right)\,dz\,du\\[6pt]&=\int _{0}^{1}\left(-{\frac {1}{2\pi ^{2}}}+{\frac {1486880\pi }{21}}+8\pi u^{4}+40u^{4}+{\frac {371720021}{1050}}\right)\,du\\[6pt]&={\frac {371728421}{1050}}+{\frac {14869136\pi ^{3}-105}{210\pi ^{2}}}\\[6pt]&\approx {576468.77}\end{aligned}}}

وبالتالي، يمكننا حساب التكامل الكفافي باستخدامن=4{\displaystyle n=4}يمكننا استخدام نفس الطريقة لتقييم تكاملات الكفاف لأي حقل متجهي معن>4{\displaystyle n>4}أيضًا.

التمثيل التكاملي

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

ومن أهم الأمثلة على ذلك صيغة كوشي التكاملية ، التي توفر طريقة لإعادة بناء دالة تحليلية من قيمها على محيط محيط بها:

و(z)=12πأناγو(ζ)ζ-zدζ{\displaystyle f(z)={\frac {1}{2\pi i}}\oint _{\gamma }{\frac {f(\zeta )}{\zeta -z}}\,d\zeta }

أينو{\displaystyle f}هي دالة هولومورفية على وداخل المحيط المغلق البسيطγ{\displaystyle \gamma }،z{\displaystyle z}هي نقطة في الداخلγ{\displaystyle \gamma }، وζ{\displaystyle \zeta }يمثل متغير التكامل. توضح هذه الصيغة أن قيمو{\displaystyle f}يتم تحديد القيم داخل المحيط من خلال قيمه على طول المحيط.

أمثلة

تحويل لابلاس العكسي

يُعرَّف تحويل لابلاس العكسي بواسطة تكامل محيطي معقد يُعرف باسم تكامل برومويتش : و(ت)=12πأناγ-أناγ+أناهـsتF(s)دs{\displaystyle f(t)={\frac {1}{2\pi i}}\int _{\gamma -i\infty }^{\gamma +i\infty }e^{st}F(s)\,ds}

هذا التكامل يعبر عن دالةو(ت){\displaystyle f(t)}من حيث تحويل لابلاس الخاص بهF(s){\displaystyle F(s)}.

تمثيل دالة الجيب

التكامل التالي يعطي تمثيلاً لدالة sinc :

الخطيئةxx=1π0كوس(xت)دت{\displaystyle {\frac {\sin x}{x}}={\frac {1}{\pi }}\int _{0}^{\infty }\cos(xt)\,dt}

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

دالة غاما

دالة غاما لها التمثيل التكاملي التالي:

Γ(z)=0تz-1هـ-تدتل يكرر(z)>0{\displaystyle \Gamma (z)=\int _{0}^{\infty }t^{z-1}e^{-t}\,dt\quad {\text{for }}\operatorname {Re} (z)>0}

تتضمن امتدادات هذا التعريف التكاملات الكفافية في المستوى المركب .

دالة زيتا لريمان

التعريف الأصلي لدالة زيتا لريمانζ(s){\displaystyle \zeta (s)}عبر سلسلة ديريشليه ، ζ(s)=ك=11كs،{\displaystyle \zeta (s)=\sum _{k=1}^{\infty }{\frac {1}{k^{s}}},} صالحة فقط لـيكرر(s)>1{\displaystyle \operatorname {Re} (s)>1}، لكن ζ(s)=Γ(1-s)2πأناح(-ت)s-1هـت-1دت،{\displaystyle \zeta (s)={\dfrac {\Gamma (1-s)}{2\pi i}}\int _{H}{\dfrac {(-t)^{s-1}}{e^{t}-1}}dt,} حيث يتم إجراء التكامل على محيط هانكلح{\displaystyle H}، صالحة لجميع المركبات المعقدةs{\displaystyle s}لا يساوي1{\displaystyle 1}.

التطبيقات

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

انظر أيضاً

مراجع

  1. ستالكر، جون (1998). التحليل المركب: أساسيات النظرية الكلاسيكية للدوال . سبرينغر. ص  77. ISBN 0-8176-4038-X.
  2. باك، جوزيف؛ نيومان، دونالد ج. (1997). "الفصلان 11 و12" . التحليل المركب . سبرينغر. ص 130-156 . ISBN  0-387-94756-6.
  3. كرانز، ستيفن جورج (1999). "الفصل 2" . دليل المتغيرات المركبة . سبرينغر. ISBN 0-8176-4011-8.
  4. ميترينوفيتش، دراغوسلاف س.؛ كيتشكيتش، يوفان د. (1984). "الفصل 2". طريقة كوشي للبواقي: النظرية والتطبيقات . سبرينغر. ISBN 90-277-1623-4.
  5. 1 2 3 4 5 ساف، إدوارد ب.؛ سنايدر، آرثر ديفيد (2003). "الفصل 4". أساسيات التحليل المركب مع تطبيقات في الهندسة والعلوم والرياضيات ( الطبعة الثالثة). برنتيس هول. ISBN  0-1390-7874-6.

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