بخشی از مقاله
معادلات دیفرانسیل – روش های تفاضل متناهی
«روشهاي تفاضل متناهي»
روابط واضح يا غيرواضح بين مشتقات و مقادير توابع در نقاط آغازي وجود دارد.
نقاط آغازي بر روي [a,b] مي تواند به وسيله [j= 1,2,…,N] و xj= a+jh به طوريكه ، ، در نظر گرفته شود.
اين عبارت براي مشتقات تحت شرايط مقادير تابعي است.
جواب مسأله مقدار مرزي يك تفاضل متناهي بوسيله جايگذاري معادله ديفرانسيل در هر نقطه آغازين به وسيله يك معادله تفاضلي بدست مي آيد.
با در نظر گرفتن شرايط مرزي در معادلات تفاضلي، سيستم جبري معادلات مورد حصول حل مي شود، اين يك جواب عددي تخميني براي مسأله مقدار مرزي بدست مي دهد.
- Linear Second Order Differential Equations
[معادلات ديفرانسيل خطي مرتبه دوم] [صفحه 5, 4 ]
به معادله ديفرانسيل مرتبه دوم زير توجه مي كنيم:
، (46)
در رابطه با شرايط مرزي نوع اول: ، (47)
مقدار قطعي u(m) از با مشخص شده و مقدار تقريبي آن با ، با استفاده از سريهاي تيلورها مي توانيم مشخص كنيم كه:
( .42)
به طوري كه و
(49)
به طوري كه
ما فرض كرديم كه پيوستگي بدين صورت است:
به طوري كه .
با در نظر گرفتن شرايط در 48 ، 49 و جايگذاري در 46 ، تفاضل تقريبي متناهي معادله ديفرانسيل مذكور در به صورت زير است:
( .50)
شرايط مرزي ( .42) به صورت زير تبديل مي شود:
( .51)
پس از ضرب با ، ( .50) مي تواند به صورت زير نوشته شود:
و ( .52)
به طوري كه:
و و
سيستم ( .52) در نوشتار ماتريسي، پس از لحاظ شرايط مرزي، تبديل ميشود به:
( .53) Au=b
به طوري كه:
حل سيستم معادلات خطي ( .53) جواب تفاضل متناهي معادله ديفرانسيل ( .46) را ارائه مي دهد كه پاسخگوي شرايط مرزي مدنظر است.
اشتباه بريدگي داخلي. (p.565) (خطاي برش)
غلط بريدگي داخلي از معادله ( .52) بوسيله
( .54)
نشان داده مي شود. به طوري كه
بسط هر شرط در طرف اول معادله ( .54) در سري تيلور آن مول ، بدست مي دهد:
( .55)
به طوري كه .
بنابراين روش مذكور، روش حل معادله مرتبه دوم مي باشد.
شرايط مرزي اشتقاقي: (p.596)
هم اكنون توجه خود را به شرايط مرزي نوع سوم معطوف مي كنيم:
( .56)
تفاضل تقريبي معادله ديفرانسيل ( .46) در گرههاي داخلي j=1,2,…,N ، بوسيله معادله ( .52) داده شده كه داراي N+2 مجموع در N معادله ميباشد. هم اكنون ما نياز داريم دو يا چند معادله متناظر براي شرايط مرزي ( .56) بيابيم.
با حذف شرايط در ( .48) ، تفاضل تقريبي متناهي ( .56) به صورت زير مي باشد:
در : يا
( .57)
در يا
( .58)
به طوري كه و ، مقادير تابعي در و مي باشند. گرههاي و خارج از بازه [a,b] قرار دارند و گرههاي غيرواقعي خوانده ميشوند:
ديفرانسيل:
مقادير و مي توانند با اين فرض كه معادله تفاضلي ( .52) براي N+1 و j= 0 در نقاط مرزي و باقي مي ماند و مي تواند ناديده گرفته شود.
جايگذاري مقادير و در ( .57) و ( .58) در معادلات ( .52) به ازاي N+1 و j= 0 ما را مي رساند به:
( .59)
معادلات ، ( .52) ، و يك سيستم سهگانه از معادلات بوجود مي آورند.
تا زماني كه تفاضل تقريبي ( .52) براي معادله ديفرانسيل ( .46) و تفاضلات تقريبي ( .59) براي شرايط مرزي ( .56) ، همگي مرتبه دوم هستند. تمام معادلات براي ، همچنين مرتبه دوم هستند.
به طور متقابل، ما نمي توانيم از نقاط غيرواقعي ، استفاده كنيم. در اين مورد ما مي توانيم از تقريب هاي زير استفاده كنيم:
يا
( .60)
( .61)
يا
تا زماني كه تقريب هاي ( .60) ، ( .61) از نوع اول هستند، تمام معادلات
( .60) ، (7.62) و (7.61) براي j= 0,…,N+1 نمي توانند مرتبه دوم بمانند. اين معادلات همچنين يك دستگاه معادلات تشكيل مي دهند.
يا
( .62)
يا
( .63)
تا زماني كه تقريب هاي ( .62) و ( .63) از مرتبه دوم هستند، تمام معادلات ( .62)، ( .52) و ( .63) براي همچنين از مرتبه دوم هستند. اگر ما را از ( .62) كه از اولين معادله مجموعه ( .52) استفاده مي كند. و را از ( .63) كه از آخرين معادله مجموعه ( .52) حذف كنيم سپس معادلات حاصله يك دستگاه معادلات سهگانه تشكيل ميدهند.
روش مرتبه چهارم در غياب در ( .46) . (p.598)
به معادله ديفرانسيل زير توجه كنيد:
( .64)
كه در ارتباط با شرايط مرزي نوع اول ( .42) است.
براي اين مسئله ما مي توانيم يك روش مرتبه بالاتر يا مرتبه چهارم بسازيم. ما معادله ديفرانسيل را به صورت زير:
( .65)
و يك روش Numeruv براي حل آن مي نويسيم.
( .66)
( .67)
شرايط مرزي اشتقاقي براي ( .56) . (p.598)
بار ديگر توجه خود را به شرايط مرزي نوع سوم معطوف مي كنيم:
( .68)
( .69)
نظر به اينكه روش Numeruv ( .67) براي ( .65) از مرتبه چهارم ميباشد، به تقريبات مرتبه چهارم براي و نياز داريم. با ، و با استفاده از بسط سري تيلور مي نويسيم:
( .70)
با استفاده از قانون سيسمون براي بررسي كران انتگرال طرف راست داريم:
( .71)
به طوري كه و
تخمين خطا از مي باشد.
هم اكنون به يك تخمين براي نياز داريم. با استفاده از سريهاي تيلور مي نويسيم:
( .72)
تخمين خطار از مي باشد. اگر تخمين ( .72) در ( .71) مورد استفاده قرار گيرد، سپس مرتبهاش را با حفظ مي كند. بنابراين با شكل دادن ( .71) ، ( .70) و ( .72) تخمين زير را داريم:
( .73)
( .74)
به طوريكه با حل كردن براي داريم:
كه از مي باشد با جايگذاري در ( .68) ، تقريب تفاضلي صحيح از داريم كه در x =a صحيح است به صورت زير:
( .75)
به طور مشابه مي نويسيم:
( .76)
بار ديگر با استفاده از قانون سيسون براي بررسي طرف راست انتگرال داريم:
( .77)
به طوري كه ؛
تخمين خطار از است.
هم اكنون با استفاده از بسط سري تيلور مي نويسيم:
( .78)
تقريب خطا از مي باشد. بنابراين با شكل دادن ( .76) ، ( .77) و (7.78) تقريب را بدين صورت داريم:
( .79)
به طوري كه ( .80)
با حل كردن براي بدست مي آوريم.
( .81)
كه از است.
با جايگذاري در ( .69) تقريب تفاضل را كه در x =b صحيح است بدست مي آوريم كه بدين صورت است.
( .82)
به جاي تقريبي كه در ( .78) داده شده، همچنين مي توانيم از عبارت زير استفاده كنيم:
پس بجاي تقريب داده شده در (7.80) داريم:
( .83)
مثال 1. حل كنيد مسأله مقدار مرزي زير را با و به كار بردن متد فوق.
؛ ؛
حل. بازه بسته [0,1] را به چهار زير جازه تقسيم مي كنيم، نقاط گرهاي عبارتند از:
، ؛
روش حل معادله ديفرانسيل مرتبه دوم يك دستگاه معادلات زير را بدست مي دهد.
با ضرب در طرفين معادله بالا بدست مي آوريم:
براي داريم.
با به كار بردن شرايط مرزي داريم:
مثال 2: مسأله مقدار مرزي زير را وقتي حل كنيد.
؛ و
حل: با به كار بردن روش حل معادله ديفرانسيل مرتبه دوم
وقتي ؛ ما چهار نقطه گرهاي داريم: و كه عبارتند از 0 ، ، و 1.
سيستم معادلات زير را بدست مي آوريم: