Алгебралық және трансценденттік теңдеулерді шешудің сандық әдістері

Пән: Математика, Геометрия
Жұмыс түрі: Курстық жұмыс
Тегін: Антиплагиат
Көлемі: 34 бет
Таңдаулыға:

Жоспар

Есептеу әдістері - таным үрдісіндегі есептеу математикасының қажетті бөлімі.
Есептеу математикасының пайда болуындағы сандық әдістердің рөлі.
Есептеу математикасының әртүрлі халық шаруашылық аймағында практикалық қолданылуының болашақтылығы.

Әдебиеттер: [1], 6-13 б., [2] -[4], 13-16 б., [6], 8-13 б.

Дәріс №2 Алгебралық және трансценденттік теңдеулерді шешудің сандық әдістері

2. 1 Теңдеулерді шешу есебінің қойылымы

Теңдеуді шешу - оның түбірлері болатынын, егер бар болатын болса нешеу екенін және оларды белгілі дәлдікпен мәндерін анықтау.

\[F(x)=0\]

түріндегі сызықтық емес теңдеулердің түбірлерін табу есебі әртүрлі ғылыми зерттеулерде кездеседі (мұндағы

\[F(x)\]

- анықталған және шектеулі немесе шектеусіз

\[[a,b]\]

аралығында үздіксіз функция) . Сызықтық емес теңдеулерді екі класқа бөлуге болады: алгебралық және трансценденттік. Алгебралық теңдеулер деп тек алгебралық функцияларды ғана (бүтін, рационал, иррационал) қамтитын теңдеулерді айтады. Дербес жағдайда, көпмүше бүтін алгебралық теңдеу болып табылады. Басқадай функцияларды (тригонометриялық, көрсеткіштік, логарифмдік, және т. б. ) қамтитын теңдеулерді трансценденттік деп атайды.

Әрбір

\[\zeta\ \in[a,b]\]

сандар

\[F({\boldsymbol{x}})\]

функциясын нөлге айналдыратын болса, яғни

\[F(x)=0\]

, берілген теңдеудің түбірі деп аталады.

\[\begin{array}{l}{{\stackrel{\sim}{\sim}}}\end{array}\]

саны

\[{\mathcal{N}}\]

еселі түбір деп аталады, егер

\[x=\zeta\]

болғанда

\[F({\boldsymbol{x}})\]

функциясымен бірге оның

\[(k-1)\]

-ші ретті туындылары да нөлге тең болса:

\[{\cal F}(\zeta\ )={\cal F}\phi\varsigma\ J=\cal{K}={\cal F}^{(k-1)}(\zeta\ )\!=\!0\]

Сызықтық емес теңдеулерді шешудің әдістері тура және итерациялық болып бөлінеді. Тура әдістер түбірлерді шекті қатынас (формула) түрінде жазуға мүмкіндік береді. Мектеп курсынан тригонометриялық, логарифмдік, көрсеткіштік, сонымен қатар қарапайым алгебралық теңдеулерді шешу үшін әдістер белгілі. Бірақ та тәжірибеде теңдеулердің мұндай әдістермен шешілмейтіндері де кездеседі. Оларды шешу үшін итерациялық әдістер пайдаланады, яғни тізбектелген жуықтау әдістері (сандық әдістер) .

Теңдеудің түбірлерін сандық әдіспен табу есебі екі кезеңнен тұрады: түбірлерді айыру, яғни түбірдің бір ғана мәнін қамтитын жеткілікті аз (сығылған) аймақтарды табу және түбірлерді анықтау, яғни қандайда бір аймақтағы түбірді белгілі дәлдікпен есептеу.

2. 2 Алгебралық және трансценденттік теңдеулердің түбірлерін айыру

Айталық бізге

\[F(x)=0\]

түріндегі теңдеу берілсін. Мұндағы

\[F({\boldsymbol{x}})\]

- алгебралық немесе трансценденттік функция. Егер біз

\[F(x)\]

функциясының графигін пайдалансақ, онда теңдеудің түбірлері жуықтап алғанда, абсцисса осімен қиылысу нүктелері болмақ. Есепті ықшамдау арқылы, берілген теңдеуді оған мәндес

\[f_{1}(x)=f_{2}(x)\]

теңдеуімен алмастыруға болады. Мұндай жағдайда

\[f_{!}(x)\]

және

\[f_{2}\left(x\right)\]

функцияларының графиктері салынып, Ох осіндегі осы графиктердің қиылысу нүктелерін көрсететін кесінділері белгіленеді.

Мысал 1.

\[\sin2x-\ln x=0\]

теңдеуінің түбірлерін айыру керек.

Түбірлерін графикалық түрде айыру үшін, оны оған мәндес

\[\sin2x=\ln x\]

түрге келтіреміз.

\[y_{1}=\sin2x\]

және

\[y_{2}=\ln x\]

функцияларының графиктерін жеке-жеке саламыз.

Графикке қарап, оның

\[\begin{array}{l}{{\dot{\gamma}}}\\ {{\sim}}\end{array}\]

бір түбірі болатынын көреміз және ол

\[\left[\mathbf{\varepsilon}_{?}\mid\mathbf{\varepsilon}_{3}\otimes\right]\]

кесіндісінде жатады.

Түбірлерді айыру туралы есептерді шешу барысында келесі жайттардың пайдасы бар:

Егеркесіндісінде үздіксізфункциясы, оның шеткі нүктелерінде әртүрлі таңбалы мәндер қабылдаса (яғни) , онда берілген теңдеудің осы кесіндіде кем дегенде бір түбірі бар болады.
Егерфункциясы монотонды (кемімелі немесе өспелі) болса, кесіндісіндегі түбір жалғыз ғана болады.

Тексеру үшін

\[F(x)=\sin2x-\ln x\]

функциясының

\[\left[\pm,\ \mid_{\sim}\mp\right]\]

кесіндісінің шеткі нүктелеріндегі мәндерін есептейік:

\[F(1)=0.909298\]

;

\[F\left(1,5\right)=-0,264344\]

. Байқауымызша,

\[\left[\pm,\ \mid_{\sim}\right]\]

кесіндісінде түбірдің болатынын аламыз.

Қарапайым жағдайда, түбірлерді графикалық айыруды қолмен еептеуге болады, кейде күрделі жағдайларда теңдеудің түбірі берілген кесіндіде болуын (санын) анықтауда компьютердің қолданбалы бағдарламасын пайдалануға немесе программалау тілінде программа құрастыруға болады.

Айталық

\[F(x)=0\]

теңдеуінің барлық түбірлері

\[\left[A;C\right]\]

кесіндісіне тиісті болсын, яғни

\[F(A)\chi F(C)<0\]

. Бізге теңдеудің түбірлерін айыру керек, яғни бір түбірден жататын барлық

\[\left[a;b\right]\subset\left[A;\,C\right]\]

кесінділерді көрсету керек.

\[F(x)\]

-тің мәнін

\[x=A\]

нүктесінен бастап оң жаққа қарай қандайда бір

\[{\mathcal{J}}_{\bar{\boldsymbol{I}}}\]

қадаммен қозғала отырып есептейміз.

\[F({\boldsymbol{x}})\]

-тің көршілес екі әртүрлі таңбалы мәндері пайда болған кезде, алынған кесіндіге түбір тиісті болатынын аламыз.

Теңдеудің шешімін программалау тілі көмегімен қарастырайық. Осыған сәйкес келетін алгоритмнің жалпы схемасын көрсетейік. Қойылған есептің нәтижесі экранда көрсетілген

\[\textstyle{X_{1}}\]

және

\[X_{2}\]

параметрлерінің мәндері (белгіленген кесіндінің шеткі нүктелері) болады.

Мысал 2.

\[\cos x=0,1x\]

теңдеуінің

\[\left[-\mid0,\mid\Theta\right]\]

кесіндідегі түбірлерін 0, 1 қадаммен бөліктерге бөліңдер.

Бұл есепті Excel көмегімен шешуге болады. Компьютер экранында

\[y=\cos x-0,1x\]

функциясының графигінен басқа мәндері анықталған кесте және одан 0, 1 қадаммен жеті бөлінген кесінділердің нәтижесін алуға болады:

\[\cdot9,7\times x_{1}\times\lnot9,6\]

;

\[-\ \ A,3 ;

\[-\ 1_{3}8 ;

\[1_{\cdot}\ 4<\,x_{5}\,<1_{\cdot}5\]

;

\[5,2 ;

\[\mathcal{D}<\,\chi_{7}\,<\,\mathcal{D}_{,1}\]

Түбірлерді айыру алгоритмінің схемасы бойынша программаны оңай құрастыруға болады. Төменде Turbo Pascal тілінде құрылған программаны ұсынамыз:

Program Separat_root;

uses crt;

var a, b, x1, x2, y1, y2, h:real; n, k:integer;

function f(x:real) :real; {уравнение вида F(x) =0}

begin f:=cos(x) -0. 1*x;

end;

begin

clrscr;

writeln (' Введите a, b, h ') ; read (a, b, h) ;

k:=0; x1:=a; x2:=x1+h; y1:=f (x1) ;

while x2<b do

begin y2:=f (x2) ;

if y1*y2<0 then

begin inc(k) ;

writeln (k, '-й корень

[ ', x1:4:1, , ' ; ', x2:4:1, , ' ] ' )

end;

x1:=x2; x2:=x1+h; y1:=y2;

end;

repeat until keypressed

end.

Программаның орындалу нәтижесі:

Ескерту. Теңдеудің түбірлерін айыру алгоритмінің сенімділігі

\[F({\boldsymbol{x}})\]

функциясының сипатына және таңдалған

\[{\mathcal{J}}_{\bar{\boldsymbol{I}}}\]

қадамның шамасына тәуелді.

\[F({\boldsymbol{x}})\]

функциясының

\[\left[{\chi}_{\nu}^{*}\ {\chi}\ +\ {\cal J}_{\zeta}\right]\]

кесіндісінде бір емес бірнеше түбірлері болуы мүмкін. Сондықтан түбірлерді айыру барысында

\[{\mathcal{J}}_{\bar{\boldsymbol{I}}}\]

қадамның шамасын өте аз етіп алған дұрыс.

Тапсырма 1. Берілген теңдеудің түбірлерін графикалық әдісті пайдаланып айырыңдар және осы тапсырманы программалау тілінің көмегімен орындаңдар:

2. 3 Жартылай бөлу әдісімен теңдеудің түбірін анықтау

Сызықтық емес теңдеулердің түбірлерін табудың қарапайым әдістердің бірі болып табылады.

\[F(x)=0\]

теңдеуін шешу барысында, ереже бойынша

\[\begin{array}{l}{{\dot{\gamma}}}\\ {{\sim}}\end{array}\]

түбірі жуық мәнінің алдын ала

\[\begin{array}{c}{{\underline{{{\mathcal{C}}}\Phi}}}\\ {{\underline{{{\langle\Omega\rangle}}}}}\end{array}\]

мүмкін қателігі беріледі. Түбірлерді анықтау үрдісінде олардың

қателігінен үлкен емес жуық мәндерін табуды талап етеді.

Айталық теңдеудің

\[[a,b]\]

кесіндісінде бір түбірі бар болсын, яғни

\[F(x)\]

функциясы осы кесіндіде үздіксіз.

\[[a,b]\]

кесіндісін

\[{\tilde{n}}={\frac{a+b}{2}}\]

нүктесімен жартылай бөлеміз.

\[F({\boldsymbol{x}})\]

функциясының

\[[a,c]\]

және

\[\textstyle\left[{\widetilde{l}},\ J\right]\]

кесінділерінің шеткі нүктелеріндегі мәндерін зерттейміз, яғни

\[a,\,c,\,b\]

нүктелерінде. Егер

\[F(c)\neq0\]

болса, онда екі жағдай болуы мүмкін:

\[F(x)\]

функциясы

\[[a,c]\]

кесіндісінде немесе

\[\textstyle\left[{\widetilde{l}},\ J\right]\]

кесіндісінде таңбасын ауыстырады (2. 6-сурет) .

Функцияның таңбасы ауыспайтын кесіндіні алып тастап, таңбасы ауысатын кесіндіні таңдап алып (өйткені, ол аралықта ізделінді түбір жатыр), оны жаңа кесінді ретінде қарастырамыз. Жартылай бөлу үрдісін жалғастыра отырып, теңдеудің түбірі тиісті болатын ең кіші кесіндіге келеміз.

Айталық анықтық үшін

\[F(a)<0,\,F(b)>0\]

болсын.

Бастапқы жуық түбір ретінде

\[\tilde{n}_{0}=\frac{a+b}{2}\]

аламыз. Қарастырылған жағдайда

\[F({\bar{n}}_{0})<0\]

, онда

\[\tilde{\cal H}_{0}\;<\tilde{n}<\mathcal{I}\]

және

\[\left[{\widetilde{n}}_{0},\,b\right]\]

кесіндісін ғана қарастырамыз. Келесі жуықтау:

\[{\tilde{n}}_{1}={\frac{c_{0}+b}{2}}\]

. Мұнда

\[\textstyle\left[{\tilde{n}}_{1}\,,\,b\right]\]

кесіндісін алып тастаймыз, өйткені

\[F({\bar{n}}_{1})>0\]

және

\[F(b)>0\]

, яғни

\[{\tilde{n}}_{0}\prec{\tilde{n}} . Осыған ұқсас басқа жуықтауларды табамыз:

\[\tilde{n}_{2}=\frac{c_{0}+c_{1}}{2}\]

және т. б.

Жуықтау үрдісін

\[F({\boldsymbol{x}})\]

функциясының модульімен алынған мәні берілген

ең аз санынан аз болғанша жалғастырамыз, яғни

\[|F({\bar{n}}_{n})|<\varepsilon\]

Төмендегі суретте теңдеудің түбірін кесіндіні жартылай бөлу әдісімен табудағы итерациялық үрдістің блок-схемасын ұсынамыз. Мұнда кесіндіні сығу

\[{\mathcal Q}\,\]

және

\[\mathit{\mathcal{I}}\]

шектерін ағымдық

\[{\boldsymbol{C}}\]

түбіріне ауыстырумен жүргізіледі.

Сонымен қатар алынған кесіндінің ұзындығын бағалауға болады: егер ол

\[\begin{array}{c}{{\underline{{{\mathcal{C}}}\Phi}}}\\ {{\underline{{{\langle\Omega\rangle}}}}}\end{array}\]

мүмкін қателігінен кіші болса, есептеу тоқтатылады, яғни

\[b-a<\varepsilon\]

\[{\mathcal Q}\,\]

және

\[\mathit{\mathcal{J}}\]

шектерін ағымдық

\[{\boldsymbol{C}}\]

түбіріне ауыстырумен жүргізіледі.

Мысал 3.

\[\sin2x-\ln x=0\]

теңдеуінің

\[\left[1,3;1,5\right]\]

кесіндісінде бір түбірі бар. Осы теңдеуді

\[1\cdot10^{-4}\]

дәлдігіне дейін жартылай бөлу әдісімен компьютердің программасы көмегімен шешейік.

Алгоритмнің блок-схемасына сәйкес Turbo Pascal тілінде программа мынадай болады:

Программаның орындалу нәтижесі:

Тапсырма 2. Жартылай бөлу әдісімен берілген теңдеудің бір түбірін 10 ^-3 дәлдікпен есепте:

а) Есептеме кестенің көмегімен;

б) Прграммалау тілінің көмегімен.

Есептің берілгендерін 1-тапсырмадан алыңыз.

2. 4 Ньютон әдістері

Жартылай бөлу әдісімен қатар күрделі және тиімді итерациялық әдістер бар. Бұл әдістерге Ньютон есімімен байланысқан әдістердің тобы қатысады. Олардың екеуін қарастырайық: жанама әдісі және хорда (қиюшы) әдісі. Бл әдістердің екеуі де мынадай тәсілге негізделген.

\[F(x)=0\]

теңдеуінің

\[\left[a;\,b\right]\]

кесіндісінде жалғыз түбірі бар болсын. Оны оған мәндес теңдеуге түрлендіреміз:

\[x=x-\,\phi\!\left(x\right)\!\cdot F\!\left(x\right)\]

мұндағы,

\[\phi({\boldsymbol{x}})\]

\[\left[a;\,b\right]\]

кесіндісінде анықталған және осы кесіндіде нөлге айналмайтын кез келген функция.

\[\phi({\boldsymbol{x}})\]

- ті әртүрлі тәсілмен таңдай отырып, көрсетілген әдістерді алуға болады.

Жанама әдісі

а) Бірінші тәсіл

Айталық

\[\phi(x)\!=\!\,{\frac{1}{F^{\prime}(x)}}\]

. Сонымен итерациялық тізбек

\[x_{n+1}=x_{n}\cdot{\frac{F(x)}{F{\dot{Q}}x}}\ (n=0,1,2,\mathrm{K})\]

реккуренттік қатынасының көмегімен құрылады. Бастапқы

\[{\widetilde{O}}_{0}\]

мәнін таңдау мәселесі,

\[F(x)\]

функциясының мынадай шарттарды қанағаттандыруымен шешіледі:

\[\left[a;\,b\right]\]

кесіндісінде екінші рет дифференциалданады;

2) Бірінші және екінші ретті туындылары осы кесіндіде таңбасын сақтайды, яғни

\[F({\boldsymbol{x}})\]

функция монотонды және дөңестік сипатын ауыстырмайды.

Мұндай жағдайда

\[{\widetilde{O}}_{0}\]

мәні ретінде

\[\left[a;\,b\right]\]

кесіндісінің шеткі нүктелерінің бірі алынады және ол нүктеде

\[F({\boldsymbol{x}})\]

функциясы және оның екінші ретті туындысы бірдей таңбалы болуы керек, яғни

\[F(x_{0})\times F\,\phi(x_{0})>0\]

шарты орындалады.

Реккуренттік қатынаспен (

\[n=0\]

) болғанда анықталған

\[\textstyle{X_{1}}\]

нүктесі,

\[\varnothing{}=F{\bigl(}x{\bigr)}\]

функциясының графигіне

\[{\mathcal{X}}_{0}\]

нүктесінде жүргізілген жанамамен абсциссаның қиылысу нүктесі болады.

Итерациялық тізбектің әрбір келесі мүшесіне

\[\varnothing{}=F{\bigl(}x{\bigr)}\]

функциясының графигіне тізбектің алдыңғы мүшесі арқылы жүргізілген жанаманың абсциссамен қиылысу нүктесі сәйкес келетін болады.

Қателікті бағалау мынадай теңсіздіктің көмегімен жүзеге асырылады:

\[F(x_{n})\]

мәндері реккуренттік тізбектің мүшелерін табуда есептелетін болады.

Мысал 4.

\[\mathrm{sin}\,2x-\,\left|\mathrm{n}\,x=0\right.\]

теңдеуінің

кесіндісіндегі түбірін

\[1\cdot10^{-4}\]

-ке дейінгі дәлдікпен жанама әдісі арқылы анықта.

Берілген жағдайда реккуренттік қатынас мынадай түрде болады:

\[x_{n+1}=x_{n}-\,\frac{\sin2x_{n}-\,\ln x_{n}}{2\cos2x_{n}\,-\,1/x_{n}}\quad\left(n=0,1,2,\mathrm{K}\,\right)\]

\[{\mathcal{X}}_{0}\]

нүктесін анықтау үшін

кесіндісінің шеткі нүктелеріндегі

\[F(x)\!=\!\sin2x\!-\!\ln x\]

және

\[F\,\phi(x)=-\,4\,\mathrm{cos}\,2\,x+\frac{1}{x^{2}}\]

таңбаларын табамыз:

\[F(1,3)\sim0,25>0\]

\[F\phi(1,3)\sim4,02>0\]

\[F(1,5)\sim-0,26<0\]

\[F\phi(1.5)\sim4.04>0\]

Сонымен,

\[x_{0}=1,3\]

Итерациялық үрдісті тоқтатуды тексеру үшін (белгілі дәлдік бойынша)

мәндерін анықтаймыз. Кесіндінің шеткі нүктелері үшін

\[F\phi(1,3)\sim-2.48,\;\;F\dot{d}(1.5)\gg-2.65\]

. Кесінді өте аз болғандықтан,

\[m=2,4\]

деп алуға болады.

Итерациялық тізбектің бірнеше мүшелерін есептейік:

\[x_{1}=1.3-\frac{\sin(2\lambda,3),\;\;\ln1,3}{2\cos(2\lambda,3)-\;1/1,3}\sim1,401948\;\;\;\;\;(n=0)\]

Түбірдің жету дәлдігіне тексеру жасайық:

- талап етілген дәлдік жеткілікті емес.

\[x_{2}=1,401948-{\frac{\sin(2\lambda,401948)-\ln1,401948}{2\cos(2\lambda,401948)-1/1,401948}}\sim1,399430\quad\quad{\mathrm{(n=1)}}\]

Тағы тексеру жасаймыз:

- талап етілген дәлдік жеткілікті.

Итерациялық тізбек мүшелерін тізбектеліп есептелуін кесте түрінде көрсетейік:

\[{\cal J}_{\bar{\cal L}}\]

\[X_{\it~l}\]

\[x_{n+1}=x_{n}-\,\frac{\sin2x_{n}-\,\ln x_{n}}{2\cos2x_{n}-1/\,x_{n}}\,\]

\[\left|x_{n}-\xi\right|\]

: 0

: 1, 3

: 1, 401948

: 0, 003

: 1

: 1, 401948

: 1, 399430

: 0, 1

Сонымен,

\[\xi=1.3994430\]

нәтижесінде барлық сандар дұрыс болады.

б) Екінші тәсіл

\[\mathcal{O}=\cal{F}(x)\]

қисығына

\[X\longrightarrow C_{k}\]

нүктесінде жанама жүргізіліп, жанаманың абсцисса осімен қиылысу нүктесі анықталады. Мұндай жағдайда

\[x=c_{0}\]

бастапқы жуықтау түбірін табу керек.

\[\varnothing{}=F{\bigl(}x{\bigr)}\]

қисығына

\[M_{0}(c_{0}\colon F(c_{0}))\]

нүктесінде жүргізілген жанаманың теңдеуі

\[\mathcal{O}-\ F\big(c_{0}\big)=F\phi c_{0}\big)\big(x-c_{0}\big)\]

Осыдан жанаманың О х осімен (

\[{\hat{\sigma}}=0\]

) қиылысу нүктесі ретіндегі

\[\scriptstyle x\;=c_{1}\]

жуық түбірін табамыз:

\[\tilde{n}_{1}=c_{0}-\frac{F(c_{0})}{F^{\prime}(c_{0})}\]

Осыған ұқсас

\[M_{1}(c_{1}\colon F(c_{1}))\]

\[M_{2}(c_{2}\colon F(c_{2}))\]

және т. б. нүктелерінде жүргізілген жанамалардың абсцисса осімен қиылысу нүктелері ретіндегі жуықтаулар табылады.

\[n+1\]

- ші ретті жуықтау үшін формула:

\[\tilde{n}_{n+1}=c_{n}-\,\frac{F\big(c_{n}\big)}{F^{\prime}(c_{n})}\]

Мұндай жағдайда,

\[F\phi(c_{n})\neq0\]

болуы керек. Итерациялық үрдісті тоқтату үшін

\[\left|F(c_{n})<\varepsilon\right|\]

шартын немесе тізбектің екі жуықтауының жақындалу

\[\left|{\widetilde n}_{n+1}-\ c_{n}\right|<\varepsilon\]

шартын пайдалануға болады.

Хорда (қиюшы) әдісі

Жанама әдісін жүзеге асыру барысында,

\[F({\boldsymbol{x}})\]

функиясының мәнін ғана емес оның

\[F^{\prime}(x)\]

туындысының мәнінде есептеу қажетті. Бірақ Ньютон әдісінің тек

\[F({\boldsymbol{x}})\]

мәнін есептеумен шектелетін нұсқасы бар.

а) Бірінші тәсіл

Егер

\[\phi\left(x\right)=\frac{x-c}{F(x)-F(c)}\]

деп алып, с мәні ретінде

\[\left[a;\,b\right]\]

кесіндісінің шеткі нүктелерінің бірі алынады және ол нүктеде

\[F\left({\bar{n}}\right)\times F\phi({\bar{n}})>0\]

шарты орындалады. Осыдан итерациялық әдіс

\[x_{n+1}={\frac{{\tilde{n}}F(x_{n})_{-}\ x_{n}F(c)}{F(x_{n})_{-}\ F(c)}}\qquad\left(n=0,1,2,\mathrm{K}\right)\]

реккуренттік қатынаспен анықталатын хорда әдісіне (қиюшы әдісіне) келеміз.

\[{\widetilde{O}}_{0}\]

мәні ретінде

\[\left[a;\,b\right]\]

кесіндісінен с мәні таңдағаннан қалған екінші шеткі нүктесі алынады (яғни, егер

\[c=a\]

болса, онда

\[x_{0}=b\]

немесе керісінше) .

Тізбек реккуренттік қатынастың формуласы бойынша құрылады. Жуықтау түбірінің бағалауы

теңсіздігінің көмегімен анықталады.

Әдістің геометриялық мағынасы төмендегі суретте көрсетілген. Берілген жағдайда

\[c=b,\;x_{0}=a\]

\[{\widetilde{O}}_{1}\]

мәніне қисықтың шеттерін қосатын хорданың абсцисса осімен қиылысу нүктесіне сәйкес келеді. Кейін қисықтың бойынан абсцисасы

\[{\widetilde{O}}_{1}\]

болатын нүкте табылып, хорда жүргізіледі және т. б.

Мысал 5.

\[\sin2x-\ln x=0\]

теңдеуінің

кесіндісіндегі түбірін

\[1\cdot10^{-4}\]

-ке дейінгі дәлдікпен хорда әдісі арқылы анықта.

\[{\widetilde{\mathcal{I}}}\]

нүктесін анықтау үшін

кесіндісінің шеткі нүктелеріндегі

\[F(x)=\sin2x-\ln x\]

және

\[F\,\phi(x)=-\,4\,\mathrm{cos}\,2\,x+\frac{1}{x^{2}}\]

таңбаларын табамыз:

\[F(1,3)\sim0,25>0\]

\[F\phi(1,3)\sim4,02>0\]

\[F(1.5)\backslash-0.26<0\]

\[F\phi(1.5)\sim4.04>0\]

Сонымен,

\[{\tilde{n}}=1,3\]

. Осыдан

\[x_{0}=1,5\]

\[{\widetilde{O}}_{1}\]

мәнін табамыз:

Берілген дәлдіктің жеткіліктілігін тексерейік:

\[\frac{\left|F{\bigl(}x_{1}\right)}{m}=\frac{\left|\sin(2\times397834)-\ \ln1.397834\right|}{2.4}\ \approx0.002\]

- талап етілген дәлдік жеткілікті емес.

Келесі жуықтауды табамыз:

- талап етілген дәлдік жеткілікті.

Итерациялық тізбек мүшелерін тізбектеліп есептелуін кесте түрінде көрсетейік:

\[{\mathcal{N}}\]

\[\textstyle{X_{n}}\]

\[x_{n+1}={\frac{{\tilde{n}}F(x_{n})_{\circ}\ x_{n}F(c)}{F(x_{n})-F(c)}}\]

\[\left|x_{n}-\xi\right|\]

: 0

: 1, 3

: 1, 397834

: 0, 002

: 1

: 1, 397834

: 1, 399410

: 0, 2

Сонымен,

\[\xi=1,3994\]

б) Екінші тәсіл

Айталық

\[\mathcal{G}=\cal{F}(x)\]

функциясы таңбасын ауыстыратын

\[\left[a;\,b\right]\]

кесіндісін бар болсын. Анықтық үшін

\[F(a)>0,\ F(b)<0\]

деп алайық. Берілген әдісте итерация үрдісі бойынша, теңдеудің түбірлеріне жуықтау ретінде хорданың абсцисса осімен қиылысқан нүктелердің

\[C_{0},\;C_{1},\mathbb{K}\]

мәндері алынады.

Алдымен АВ хордасының теңдеуін табамыз:

\[{\frac{y-F(a)}{F(b).\ F(a)}}={\frac{x-a}{b-a}}\]

Оның абсцисса осімен қиылысу нүктелері үшін (

\[x=c_{0},\ \delta=0\]

)

\[\tilde{n}_{0}=a-\frac{b-a}{F(b)-F(a)}F(a)\]

теңдеуін аламыз.

Қарастырылған жағдай үшін

\[{\cal F}\left(\cal G\right)\]

және

\[F({\tilde{n}}_{0})\]

шамаларының таңбаларын салыстыра отырып, түбірдің

\[\left(Q\colon\widetilde{J}_{0}\right)\]

аралығында болатынын аламыз, өйткені

\[F(a)F(\bar{n}_{v})<0\]

\[\left[{\widetilde{n}}_{0}:D\right]\]

кесіндісін қарастырмаймыз. Келесі итерация АВ ₁ хордасымен абсцисса осінің қиылысу нүктесі болатын жаңа

\[{\mathcal{C}}_{1}\]

жуықтауын анықтаудан тұрады және т. с. с. Итерациялық үрдіс

шарты орындалғанша жалғастырылады.

Хорда әдісінің блок-схемасы

Тапсырма 3. Берілген теңдеудің бір түбірін 10 ^-6 дәлдікпен Ньютон әдістерінің бірін пайдаланып, программалау тілінің көмегімен есепте. Есептің берілгендерін 1-тапсырмадан алыңыз.

Дәріс №3 Сызықтық теңдеулер жүйесін шешудің сандық әдістері

3. 1 Сызықтық жүйелер

Көптеген қолданбалы және таза математикалық есептер жиыны сызықтық алгебралық жүйелерді шешудің қажеттілігіне әкеледі. Сызықтық теңдеулер жүйесін шешу есептеу математикасының қажетті есептерінің бірі болып табылады.

... жалғасы

Сіз бұл жұмысты біздің қосымшамыз арқылы толығымен тегін көре аласыз.