3 СЫЗЫҚТЫ ЕМЕС ТЕҢДЕУЛЕРДІ ШЕШУ
3.1 Сызықты емес теңдеулерді шешу кезеңдері
)
( x
f
функцияны қарастыраиқ.
0
)
(
f
шартты қанағаттандыратын кез келген сан ζ функцияның нөлі
деп, немесе (1) теңдеудің шешімі деп аталады:
0
)
(
x
f
.
(1)
Сандық әдістердің бір бөлімі «бір өлшемді сызықты емес теңдеулер»
болып табылады. Физикалық және басқа да құбылыстардың теңдеумен
сипатталатыны белгілі. Сол теңдеуді классикалық математикалық
формуламен шешу мүмкін емес жағдайлар бар. Бұл уақытта практикада
сандық әдістерге жататын әдістермен шешілетінін дәлелдеу керек. Әрине ең
алдымен құрылған теңдеудің қай аралықта анықталғандығын, үзіліссіздігін,
түбірінің барлығын, оның жалғыздығын дәлелдейтін аргументтерді бақылау
керек. Осы этаптан өткеннен кейін ғана есепті осы теңдеуге қолдануға
келетін алгоритм көмегімен шығаруға болады. Сонымен, сызықты емес
теңдеулерді шешудің екі кезеңі бар, ол:
1. Түбір жатқан аралықты анықтау.
2. Түбірін берілген дәлділікпен анықтау.
Сызықты емес теңдеуді сандық шешу екі тәсілден тұрады.
1 Тура тәсіл - есепті математикалық дәлелденген бір формулаға қою
арқылы тікелей шығару;
2 Итерациялық тәсіл – есепті формула көмегімен бастапқы жуықтауды
беру арқылы жуықтап, біртіндеп шығару;
Тура тәсілмен шығарылған есептер дәл мәнді береді. Ал итерациялық
тәсілмен шешілген есептер есептің жуық мәнін береді. Мұның ішінде
итерациялық әдістер сандық әдіске жатады.
Бір өлшемді сызықты емес теңдеуді шешудің келесі әдістері бар.
1 Кесіндіні қақ бөлу - дихотомия әдісі деп аталады;
2 Хорда әдісі;
3 Жанама әдісі немесе Ньютон әдісі;
4 Қарапайым итерациялық әдіс немесе жәй итерация әдісі т.б.;
Түбір жатқан аралықты анықтау әдісі
F(x)=0
(2.1)
Бірөлшемді сызықты емес теңдеу берілген. Мұндағы F(x) функциясы
[a,b] кесіндісінде анықталған және үзіліссіз болсын.
Теорема1.1: [а,в] аралығында анықталған, үзіліссіз F(x) функциясының
екі шеткі нүктелердегі мәндерінің таңбалары әр түрлі болса, яғни мына
шарт орындалса f(a)*f(b)<0, онда осы аралықта (2.1)-теңдеудің түбірі
бар және жалғыз болады.
Практикада кейде теореманың орындалуын функцияның мәндер
кестесін құру арқылы да анықтайды. Функцияның анықталу облысы
бойынша а нүктесін беріп, ол нүктедегі функция мәнін анықтайды, сосын һ
қадаммен келесі нүктеге жылжып, сол нүктедегі функция мәнін анықтайды,
сол сияқты бірнеше нүктедегі функция мәндерін анықтап, таңбасын
салыстырады. Егер көрші нүктелерде функция әр түрлі таңба қабылдаса, сол
аралықта жалғыз түбірі жатыр деп айтады.
1-мысал:
Берілген теңдеудің түбірін анықтау:
)
;
(-
0
10
x
e
x
Теңдеудің түбірі жатқан аралықты аналитикалық тәсілмен табамыз: ол
үшін функция туындысын тауып, оны нөлге теңестіру арқылы
экстремумдарын анықтаймыз:
10
)
(
x
e
x
f
, экстремумы: х
1
=Ln10=2,3;
Экстремум
нүктелеріндегі
функция
таңбасының
1-кестесін
толтырамыз.
1-кесте-
x
e
x
f
x
10
)
(
функциясының таңбасын анықтау
Нүктелер
2,3
sign(f)
+
-
+
Функция таңбасының ауысуы (
; 2,3] және [2,3;
) аралығында
байқалды. Яғни осы аралықта теңдеудің түбірі бар.
Енді графиктік әдісті қарастырайық. Ол үшін теңдеуді мына түрлерге
жіктейміз, себебі функция күрделі, трансцендентті, бірден графигін құруға
болмайды:
x
e
f
x
10
f
;
2
1
. Екі функцияның графигін саламыз, екеуінің
қиылысқан нүктесі теңдеудің түбірі болып табылады (1-сурет). Қиылысу
нүктелерінің аймақтарын анықтаймыз.
1-сурет-
x
e
f
x
10
f
;
2
1
функцияларының графиктері
1
2
3
4
5
6
Ряд1
-1
0
1
2
3
4
Ряд2 0,3679
1
2,7183 7,3891 20,086 54,598
Ряд3
-10
0
10
20
30
40
-20
-10
0
10
20
30
40
50
60
Ряд1
Ряд2
Ряд3
Бірінші түбірі [0,1] аралығында, ал екінші түбірі [2,6] аралығында
жататыны суретте көрініп тұр. Енді осы аралықтағы қай нүкте (2.1)-ші
теңдеуді қанағаттандыратынын анықтаймыз.
3.2 Теңдік түбірін анықтау әдістер
3.2.1 Ньютон әдісі
Алдыңғы әдістерге қарағанда бастапқы жуықтау дұрыс таңдалынып
алынса Ньютон әдісі тез жинақталады. Бұл әдіске қатысты теореманы келтіре
кетейік:
Теорема 1.3.: f(x) функциясы [a,b] аралығында анықталған және екі
ретті туындысы бар, осы аралықта түбір жатыр f(a)*f(b)<0,
туындылардың таңбалары осы аралықта тұрақты болса f(x)*f
'
(x)>0,
онда f(x
0
)*f
''
(x
0
)>0 теңсіздігін қанағаттандыратын
]
,
[
0
b
a
x
бастапқы
жуықтаудан бастап (2.1)-ші теңдеуді қанағаттандыратын [a,b]-да
жататын жалғыз шешімге жинақталатын
0,1,2,...
n
)
(
)
(
'
1
n
n
n
n
x
f
x
f
x
x
итерациялық тізбек құруға болады.
Ньютон әдісінің геометриялық мағынасы: координаталары (x
n
;f(x
n
)) ,
болатын нүктеден қисыққа жанама жүргізсек, оның ох өсімен қиылысу
нүктесі теңдеудің түбіріне х
n+1
– кезекті жуықтау болып табылады.
Түбірге n-ші жуықтаудың қателігін бағалау үшін келесі теңсіздіктің
орындалуын
қадағалау
керек:
2
1
1
2
2
n
n
n
x
x
m
M
x
.
Мұндағы
М
2
–
функцияның екінші ретті туындысының аралықтағы максимумы, m
1
-
минимумы. Егер,
e
x
x
n
n
1
болса, онда
1
2
2
2 m
e
M
x
n
болады, яғни түбірге
дұрыс жуықталынса, әр итерациядан кейін кезекті жуықтаудың ондық таңба
саны екіге артады да процесс тез жинақталады. Егер түбірді берілген е
дәлдікпен табу керек болса, итерациялық процесті
2
1
0
1
2
M
e
m
e
x
x
n
n
шарты
орындалғанша жалғастырамыз.
Бұл біртіндеп жуықтау әдісі деп аталады. Әдістің жинақталу
жылдамдығы х
0
бастапқы нүктені дұрыс таңдауға байланысты. Егер итерация
процесінде функцияның туындысы нөлге тең болса, қисыққа жүргізілген
жанама
x
осіне параллель болса, онда бұл әдісті қолдану қиындайды. Сол
сияқты функцияның екінші ретті туындысының мәні шексіз үлкен болса
және функцияның өзі бірінші ретті туындысы нөлге тең болса, онда шыққан
түбірлер еселі болып, жинақталмауы мүмкін. Бұл әдісті қолдану үшін
)
( x
f
функциясы үзіліссіз және дифференциалданған болуы керек.
0
x
бастапқы жуықтауды таңдалынған уақытта құрылатын тізбек монотонды
кемімелі болуы керек.
)
(
)
(
1
n
n
x
f
x
f
Алгоритмі:
1
Қисықтың бойынан қандайда бір
0
x
нүктесін таңдап, осы нүкте
арқылы қисыққа жанама жүргізіледі.
2
Жанама
x
осін қиған кезде табылған нүктенің мәні мына
формуламен есептелінеді.
)
(
)
(
1
n
n
n
n
x
f
x
f
x
x
,
0
)
(
n
x
f
(2.5)
Табылған нүктедегі функцияның мәні нөлге өте жуық болса, онда сол нүкте
(2.1) теңдеудің түбірі, болмаса процесс жалғасады.
Егер түбірлер еселі болса, ол еселікті
деп белгілейік.
)
(
)
(
1
n
n
n
n
x
f
x
f
x
x
,
0
)
(
n
x
f
(2.6)
3.2.2 Хорда әдісі
Бұл әдіс кесіндіні қаққа бөлу әдісіне қарағанда шешімге тез
жинақталады
.
Алгоритмі:
1
х
n
, x
n+1
аралығында f (x) және f (x
n+1
) функцияларының таңбасы
бір біріне қарама-қарсы және түбірі бар болсын.
2
Осы екі шеткі нүктеден хорда жүргізіп, хорданың х осімен
қиылысқан нүктесін мына формуламен анықтаймыз.
)
(
)
(
)
(
)
(
1
1
*
n
n
n
n
n
n
x
f
x
x
x
x
x
f
x
f
x
(2.4)
Егер f(a)>0 шарты орындалса, а нүктесі тұрақты болады да
)
(
)
(
)
(
*
a
f
b
f
a
b
a
f
a
x
формуласымен есептеледі
Егер f(b)>0 шарты орындалса, b нүктесі тұрақты болады да
)
(
)
(
)
(
*
b
f
a
f
b
a
b
f
b
x
формуласымен есептеледі
3
х
*
нүктесіндегі функция мәнін F(x
*
)-ны есептеу. Оның таңбасын екі
шеткі нүктедегі функцияның таңбасымен салыстырылады. Егер f
(x
n
) және f(x
*
) функциясының таңбасы бірдей болса, онда хорданы
x
n+1
және x
*
нүктесі арқылы жүргізіледі. Оның мәнін (2.4)
формуламен табады. Егер f(x
n+1
) мен f(x
*
) функцияның таңбалары
бірдей болса, онда хорданы x
n
және x
*
нүктесі арқылы жүргізіледі.
Шыққан нүктенің мәні (2.4) формуламен есептелінеді.
4
x
*
нүктедегі мәнін есептеп, мәні нөлге жуық болса
e
x
x
k
k
1
, онда
x
*
нүктесі (2.1) теңдеудің түбірі деп аталады. Егер нөлге жуық
болмаса, онда процесс жалғасады.
Алдындағы мысал үшін программасы келесідей болады:
3.2.3 Жай итерация әдісі
Бұл әдісті қолдану үшін (2.1)-ші теңдеудің сызықты мүшесі
айшықталып мына түрге келтіру керек:
)
( x
x
(2.3)
Сосын теңдеудің түбіріне кез келген Х
0
бастапқы жуықтау беріп
)
(
1
k
k
x
x
k=1,2,… формуласымен х
1
, х
2
,…,х
n
нүктелер тізбегін құрамыз.
Бұл тізбек x=z түбіріне жинақталуы керек. Егер limX
k
=z болса, онда z нүктесі
)
( x
x
теңдеуінің түбірі бола алады. Итерация әдісінің жинақтылық шарты
1
)
(
x
және бастапқы жуықтау кез келген болады. Итерациялық процесс
берілген дәлдікке жетуі үшін
q
q
e
x
x
k
k
1
1
шарты орындалуы керек.
Итерациялық тізбектің жинақтылығы теореманың ([1] қараңыз)
шарттарымен де тексерілуі керек:
Теорема1.2.:
)
( x
x
теңдеуінің [a,b] аралығында жалғыз түбірі бар және келесі
шарттар орындалсын:
1
)
( x
функциясы
[a,b]
аралығында
анықталған
және
дифференциалданады;
2
]
,
[
b
a
x
үшін
]
,
[
)
(
b
a
x
;
3 барлық
]
,
[
b
a
x
үшін
1
)
(
q
x
болатындай q саны табылсын,
онда
)
(
1
k
k
x
x
, (k=1,2,…) итерациялық тізбегі
]
,
[
0
b
a
x
кез келген
бастапқы жуықтауда жинақталады.
Теңдеуді итерациялық түрге келтіру
(2.1)-теңдеуді (2.3)-ші итерациялық түрге келтіру әртүрлі тәсілдермен
орындалады. Қай тәсілмен итерациялық түрге келтірілсе де жоғарыдағы
теореманың шарттары орындалуы керек. Практикада көрінгендей,
теореманың 3)-ші шартының орындалуы күрделірек, сондықтан төмендегі
тәсілдердің бірін қолдануға болады:
1 (2.1)-теңдеуді
)
( x
f
m
x
x
түріне келтірейік, мұндағы m=const нөлден
өзгеше. Бұл уақытта
)
(
)
(
x
f
m
x
x
деуге болады. Оны дифференциалдасақ:
)
(
)
(
x
f
m
x
x
. Теореманың 3)-ші шарты орындалуы үшін
]
,
[
b
a
x
үшін
1
)
(
x
f
m
шарты орындалатындай етіп m-ді таңдап алу керек.
2 (2.1)-теңдеу (2.3)-ші түрге келтірілген, бірақ теореманың 3)-ші шарты
орындалмай тұрса, онда y=
)
( x
функциясының орнына x=g(y) функциясын
қарастыруға болады. Мұндағы g(y) функциясы
)
( x
функциясының кері
функциясы. Енді y=g(y) теңдеуін шешетін боламыз. Немес ескі белгілеулер
бойынша x=g(x) теңдеуін шешеміз деген сөз. Кері функциялардың
туындыларының қасиеттері бойынша
1
)
(
1
)
(
x
x
g
.
2 – мысал
)
;
(-
0
10
x
e
x
теңдеуінің түбірін қарапайым итерация әдісімен
табу керек болсын. Теңдеуді итерациялық түрге келтіреміз:
. Ал және барлық х нүктелері үшін
1
1
,
0
1
,
0
)
(
x
e
x
. Яғни q=0,1 деп алып, бастапқы жуықтауды х
0
=0 десек
e
x
x
k
k
1
шарты орындалғанша итерациялық процесті құрамыз: х
0
=0:
1
,
0
1
,
0
0
1
e
x
;
01
,
0
1
,
0
1
,
0
2
e
x
, т.с.с. Түбір мәні х
6
=0,111833, итерация саны 5-
ке тең.
x
e
x
1
.
0
)
(
0
(x)
,
1.
0
)
(
x
e
x
Достарыңызбен бөлісу: |