Numeeriset (laskennalliset) menetelmät - menetelmät matemaattisten ongelmien ratkaisemiseksi numeerisessa muodossa [1] .
Sekä alkutietojen esitys tehtävässä että sen ratkaisu - luvun tai lukujoukon muodossa .
Monet numeeriset menetelmät ovat osa matemaattisia ohjelmakirjastoja [2] . Koulutusjärjestelmässä teknisten erikoisalojen insinöörit ovat tärkeä osa.
Laskennallisten menetelmien perusteet ovat:
Kaikki laskennallisen matematiikan tehtävät ratkaistaan seuraavassa järjestyksessä [3] :
Symbolisesti tuntemattoman suuren etsimisen ongelma kirjoitetaan muodossa . Laskennallisen matematiikan hakuun käytetään yhtä tai useampaa avaruuden substituutiota, jossa määritetään suuret , , tai funktiot laskennan helpottamiseksi. Tuloksena olevalla uudella ongelmalla tulisi olla ratkaisu, joka on lähellä alkuperäisen ongelman ratkaisua. Esimerkiksi integraalia laskettaessa janan jatkuva funktio voidaan aina korvata polynomilla , jonka integraali on helppo määrittää; tai korvaa integraali äärellisellä summalla ja ratkaise tuloksena oleva ongelma. Tällaisen korvaamisen suorittamiseksi on löydettävä rajallinen joukko elementtejä, jotka lähentyvät hyvin päätilaa. Viimeinen ehto asettaa rajoituksia metriavaruuteen . Suurin rajoitus on -verkon olemassaolo, josta tila on itsessään kompakti ja erotettavissa . Tämä rajoitus ei kuitenkaan ole pakollinen. Nykyaikaiset funktionaalisen analyysin menetelmät mahdollistavat ongelman olosuhteisiin sopivimpien metristen avaruuksien valitsemisen [7] .
Numeerisia menetelmiä käytettäessä syntyy monenlaisia virheitä. Kun yhtä lukua lähestyy toinen, tapahtuu pyöristysvirhe, virheelliseen alkutietoon liittyvää virhettä kutsutaan kohtalokkaaksi, lisäksi alkuperäisen ongelman korvaamisen vuoksi likimääräisellä menetelmässä on virhe. Kokonaisvirhe on tässä tapauksessa menetelmän virheen ja laskelmien virheen summa, toisin sanoen yhtälön sijasta ratkaistaan yhtälö, jonka ratkaisun tarkkuus määräytyy kaavalla [8]
Virheen suuruuden määrittämiseen käytetään absoluuttisen ja suhteellisen virheen käsitteitä sekä rajoittavaa absoluuttista ja suhteellista virhettä, kun taas virheteoria määrittää virheiden suuruuden muutoksen eri aritmeettisten operaatioiden aikana [9] . Virheiden tarkan arvioinnin menetelmien lisäksi, joiden seurauksena virheiden raja-arvot määritetään, käytetään tilastollisia menetelmiä yksittäisten virheiden saavuttamisen mahdollisuuden määrittämiseen [10] ja myös satunnaisten virheiden matemaattiset ominaisuudet. liittyy poikkeamiseen määritellyistä koeolosuhteista, kun useat mittaustulokset fyysinen määrä määritetään sen likimääräisellä arvolla [11] .
Arvotaulukon antaman funktion arvon saamiseksi argumentin väliarvoille rakennetaan likimääräinen funktio , joka tietyissä pisteissä , joita kutsutaan interpolaatiosolmuiksi, ottaa arvot ja muissa pisteissä. kuuluu funktion alueeseen. Useimmiten likimääräinen funktio konstruoidaan algebrallisena polynomina, joka sisältää lineaarisesti riippumattoman järjestelmän ensimmäiset alkiot. Käytännössä lineaarisesti riippumattoman järjestelmän elementteinä potenssien sarja : , trigonometriset funktiot : , eksponentiaaliset funktiot : [12] .
Interpoloivan funktion muodostamiseksi tässä tapauksessa on tarpeen ratkaista yhtälöjärjestelmä tuntemattomilla. Järjestelmän tuloksena olevalle matriisille asetetaan tietyt ehdot: matriisin asteen tulee olla yhtä suuri kuin , ja — lineaarisen riippumattomuuden ehdon takaamiseksi — jotta ongelman ratkaisu on yksiselitteinen, matriisin determinantti — niin että ratkaisu on olemassa ja lisäksi ainutlaatuinen [13] . Lagrangen interpolaatiopolynomin rakentaminen on perusmenetelmä tällaisten ongelmien ratkaisemiseksi, erittäin resurssiintensiivinen ja vaikeasti laajennettavissa [14] .
Seuraava askel on ottaa käyttöön - :nnen kertaluvun jaetun erotuksen käsite, joka perustuu naapurisolmujen funktion arvon eron suhteeseen solmujen väliseen etäisyyteen, jolla on määritelmänsä perusteella luku. hyödyllisistä ominaisuuksista erityisesti jaetuilla kertalukueroilla astepolynomista on aste , eli järjestyserot ovat vakioita, kun taas korkeamman kertaluvun erot ovat [15] . Jaetut erot mahdollistavat Lagrangen interpolaatiopolynomin kirjoittamisen uudelleen muotoon, joka on helpompi laskea. Uutta kaavaa kutsutaan Newtonin interpolaatiopolynomiksi [16] , yhtäläisten intervallien tapauksessa kaava yksinkertaistuu huomattavasti [17] . Jaettujen erojen avulla muodostetaan Gaussin , Stirlingin , Besselin , Everettin interpolaatiokaavat [18] . Yleisessä tapauksessa jaetut erot ensin pienenevät järjestyksen kasvaessa ja alkavat sitten taas kasvaa, eli laskelmissa ei ole järkevää käyttää korkean kertaluvun eroja [19] . Tämä herättää kysymyksen interpolointiprosessin konvergenssista, jonka ratkaisemiseen käytetään erilaisia matemaattisen analyysin menetelmiä [20] .
Käytännön ongelmia ratkaistaessa on tarpeen laskea toistuvasti tietyn funktion arvot, mikä yleensä on resurssivaltaista toimintaa. On tarpeen löytää parhaan yhtenäisen approksimoinnin funktio [21] . Approksimaatiota varten funktiot lineaarisessa normiavaruudessa muodostavat aliavaruuden kaikkien mahdollisten lineaariyhdistelmien dimensioista, joille normi on määritelty ja sen infimumi on olemassa . Elementtiä, jossa tämä reuna saavutetaan, kutsutaan parhaan approksimaation elementiksi eli projektioksi [22] . Voidaan todistaa, että aliavaruudessa on aina olemassa parhaan approksimation omaava elementti [23] , ja tilan tiukan normalisoinnin ehdolla tällainen elementti on ainutlaatuinen [24] . Jatkuvien funktioiden avaruudessa normin kanssa
on myös parhaan approksimaatioelementti [25] , mutta sen ainutlaatuisuuden ehtona on korkeintaan erotettavien yleisen polynomin nollien läsnäolo välissä ( Tšebyševin polynomit ) [26] .
Funktioteoria soveltuu tehofunktioiden järjestelmään, koska se on Tšebyševin järjestelmä millä tahansa aikavälillä [27] . Weierstrassin lauseen mukaan aliavaruuden ( ) mittasuhteen kasvaessa projektion ja annetun funktion välinen ero pyrkii nollaan [28] . Tämän approksimoinnin järjestys riippuu funktion rakenteellisista ominaisuuksista, se voidaan määrittää Bernsteinin polynomeilla [29] . Trigonometristen funktioiden järjestelmällä on myös Tšebyševin järjestelmän ominaisuuksia välissä , sillä projektion ja annetun funktion välinen ero pyrkii myös nollaan [30] .
Huolimatta esitetystä parhaan approksimaatiopolynomin olemassaolosta, ei ole tapoja rakentaa sitä tarkasti. Sen sijaan käytetään useita menetelmiä parhaan yhtenäisen approksimation omaavien polynomien konstruktion approksimoimiseksi [31] .
Monissa tapauksissa yhtenäisen approksimoinnin vaatimus on tarpeeton ja funktioiden "integraalinen" läheisyys riittää, lisäksi kokeesta saadut likimääräisten funktioiden arvot sisältävät satunnaisia virheitä, eikä ole suositeltavaa vaatia funktioiden yhteensopivuutta. approksimoivat ja approksimoivat funktiot, jos jälkimmäinen sisältää epätarkkuuksia. Neliön keskiarvon likiarvomenetelmä ottaa seuraavan arvon läheisyyden mittana
mikä mahdollistaa integrandin interpoloinnin ja jatkuvuuden vaatimuksen luopumisen säilyttäen vain neliöintegroitavuuden vaatimukset [32] .
Funktionavaruuteen määritelty yhtälö, jonka muoto on , voi sisältää differentiaatio- ja integrointioperaattoreita , joille on mahdotonta löytää tarkkaa ratkaisua. Numeerisen differentioinnin ja integroinnin menetelmät perustuvat interpolointiin [33] .
Pääfunktion derivaatan katsotaan olevan suurin piirtein yhtä suuri kuin interpolointifunktion derivaatta, kun taas interpolointikaavan lopputermin derivaatta voi olla suuri, erityisesti korkeamman kertaluvun derivaatan kohdalla [34] . Numeeriset differentiaatiokaavat perustuvat suurelta osin Newtonin [35] , Gaussin, Stirlingin ja Besselin [36] interpolaatiokaavojen suoraan differentiaatioon, joka on rakennettu hajautettujen erojen varaan, mutta on myös erottomia kaavoja. Erityisesti kun numeeriselle differentiaalille käytetään suoraan Lagrangen kaavaa yhtäläisille intervalleille [37] , määrittelemättömien kertoimien menetelmää ja muita [38] .
Integraalin tapauksessa jo integraalin määritelmä osoittaa mahdollisuuden korvata se integraalisummalla , mutta tämä tekniikka konvergoituu hitaasti ja siitä on vähän hyötyä. Pääfunktion integraalin katsotaan olevan suunnilleen yhtä suuri kuin interpoloivan funktion integraali, ja jatkossa käytetään interpolointikaavoja, joissa on useita solmuja [39] . Lagrangen interpolaatiopolynomin käyttö yhtäläisille intervalleille integrandina johtaa Newton-Cotesin kaavoihin [40] ja sen erityistapauksiin, puolisuunnikkaan kaavaan , kun integrandikäyrä korvataan sointeella ja integraali on yhtä suuri kuin . puolisuunnikas ja Simpsonin kaava, kun integrandikäyrä korvataan kolmen pisteen läpi kulkevalla paraabelilla [ 41] . Luopumalla yhtäläisten välien vaatimuksesta Lagrangen interpolaatiopolynomin avulla voidaan saada tarkempia kaavoja numeeriseen integrointiin, erityisesti Gaussin kaavat [42] , Hermite-kaavat [43] , Markovin kaavat [44] , Chebyshev-kaavat [45 ] ] . Gaussin interpolaatiokaavojen pohjalta rakennetut kvadratuuriprosessit konvergoivat aina, kun taas Newton-Cotesin kaavoilla ei yleensä ole näitä ominaisuuksia [46] .
Numeeriseen integrointiin on muitakin tapoja, joista tärkein on Eulerin kaavojen käyttö, joissa muuttujien muutos ja sitä seuraava integrointi osien mukaan johtaa kaavaan numeerista integrointia varten puolisuunnikkaan ja korjaustermiin, johon muuttujien muutos ja osien integrointi otetaan uudelleen käyttöön. Yleisessä tapauksessa Eulerin kaava käyttää kertoimina lukuja ja Bernoullin polynomeja [47] . Kysymys yhden tai toisen numeerisen integroinnin menetelmän soveltamisesta riippuu sellaisista tekijöistä kuin laskentatyökalut, vaadittava tarkkuus ja integrandin määrittelytapa. Manuaalisissa laskelmissa on suositeltavaa käyttää eroja sisältäviä kaavoja, kun taas automaattisissa laskelmissa - ei-erotuskaavoja, erityisesti Gaussin kaavoja [48] .
Useiden integraalien likimääräiseen laskemiseen käytetään toistuvasti yksittäisten integraalien numeerisen integroinnin kaavoja, kun taas funktion ominaisuuksista riippuen eri integraaleille voidaan käyttää erilaisia kaavoja. Tätä menetelmää käytettäessä integrandi on laskettava suuressa määrässä pisteitä, joten on suositeltavaa käyttää Gaussin ja Chebyshev-kaavoja, jotka ovat tarkempia [49] . Toinen tapa on korvata integrandi interpolaatiopolynomilla kahdessa tai useammassa muuttujassa [50] . Lyusternik ja Ditkin ehdottivat Maclaurin-kaavojen käyttöä moniintegraalin likimääräiseen laskemiseen [51] . Samaan aikaan integraalin moninkertaisuuden kasvaessa niiden pisteiden määrä, joille on tarpeen tietää integrandin arvot, jotta voidaan käyttää interpolointiin perustuvia menetelmiä, kasvaa jyrkästi. Useampien integraalien laskemiseen käytetään useammin todennäköisyyspohjaisia Monte Carlo -menetelmiä , kun taas tarve saada yhtä mahdollisia sekvenssejä aiheuttaa ylimääräisiä virheitä, joita on vaikea arvioida [52] .
Lineaaristen algebrallisten yhtälöiden järjestelmien ratkaisemiseen on olemassa kaksi menetelmäryhmää: tarkat menetelmät mahdollistavat äärellisen määrän operaatioiden avulla saada tarkat tuntemattomien arvot ja sisältävät järjestelmän muuntamisen yksinkertaiseen muotoon ja yksinkertaistettu järjestelmä; Alkuperäisiin approksimaatioihin perustuvat peräkkäiset approksimaatiomenetelmät mahdollistavat "parannettujen" likimääräisten arvojen saamisen, joille "parannus"-toiminto tulee toistaa peräkkäin; Monte Carlon menetelmät mahdollistavat satunnaismuuttujien matemaattisten odotusten perusteella ratkaisun saamisen järjestelmään [53] .
Algebran koulukurssilta tunnettu eliminointimenetelmä mahdollistaa järjestelmän matriisin pelkistämisen diagonaali- tai kolmiomuotoon [54] . Gaussin eliminointimenetelmä pääelementin valinnalla, joka on välttämätön laskennallisen virheen vähentämiseksi, sisältää liikkeen eteenpäin (itse eliminointiprosessi) ja taaksepäin liikkeen (kolmiomatriisin järjestelmän ratkaisu) [55] . Sen kompaktia versiota käytetään käänteismatriisin määrittämiseen, josta voi olla hyötyä, jos vain oikea puoli muuttuu lineaariyhtälöjärjestelmässä [56] ja determinanttien laskemiseen [57] . Jordan-kaavio mahdollistaa käänteisen liikkeen helpotuksen [58] , ja järjestelmässä ilman käänteistä liikettä, joka perustuu solumatriisin muuntamiseen , jälkimmäistä ei vaadita [59] . Matriisisymmetriaehdon avulla voimme tehdä useita yksinkertaistuksia ja käyttää neliöjuurimenetelmää, jossa järjestelmämatriisi esitetään alemman kolmiomatriisin tulona sen suhteen transponoidulla matriisilla, jossa kolmiomatriisit määritetään kaavoilla alkuperäisen matriisin alkioiden tulojen kautta (jos positiivisten määrällisten matriisien ehtoa ei ole, jotkut kaavat voivat sisältää imaginaarielementtejä), ja sitten järjestelmä ratkaistaan kahdessa vaiheessa apuratkaisun kautta. kolmiomatriiseille rakennetut järjestelmät [60] . On olemassa myös skalaaritulon ominaisuuksiin perustuva ortogonalisointimenetelmä [61] , konjugaattigradienttimenetelmä, jossa konstruoidaan apufunktio, joka muodostaa ellipsoidien perheen, joilla on yhteinen keskus ja jolle on tarpeen löytää vektori. jolle se ottaa vähimmäisarvon [62] . Korkean asteen matriiseille käytetään solun osiointimenetelmää, jolloin ongelma rajoittuu alemman asteen matriisien ongelmien ratkaisemiseen [63] .
Peräkkäisten approksimaatioiden tapauksessa käytetään toistuvaa kaavaa
missä on funktio, joka riippuu systeemimatriisista, oikeasta, approksimaatioluvusta ja aiemmista approksimaatioista , missä on alkuvektori. Tässä tapauksessa menetelmää pidetään ensim- mäisenä, jos funktio riippuu vain viimeisestä aikaisemmista approksimaatioista. Tässä tapauksessa kaava voidaan kirjoittaa muodossa , missä . Laskelmien mukavuuden vuoksi on toivottavaa käyttää diagonaalista tai kolmiomatriisia , joka on kätevä kääntää. Tämän matriisin valinnasta riippuen menetelmiä kutsutaan täysvaiheiseksi ja yksivaiheiseksi [64] . Lineaariset täysivaiheiset menetelmät sisältävät yksinkertaisen iteroinnin [65] , Richardsonin menetelmän [66] ; lineaarisiin yksivaiheisiin menetelmiin - Seidelin menetelmä [67] , relaksaatiomenetelmä [68] ; epälineaarisiin menetelmiin - jyrkimmän laskun menetelmään [69] .
Algebrallisen yhtälön ratkaisu , jossa todellisen tai kompleksisen argumentin funktio on vasemmalla puolella, on kompleksitasossa [70] . Sen määrittämiseksi on ensinnäkin tarpeen sulkea jokainen juuri riittävän pienelle alueelle, eli erottaa se, johon usein käytetään graafisia menetelmiä [71] . Reaalijuurille käytetään myös yleistettyä Descartes-sääntöä, Sturmin lausetta [72] ja Fourier-menetelmää [73] . Neliöjuurimenetelmä tai Lobachevsky-menetelmä [74] on löytänyt laajan sovelluksen . Perusmuodossaan se koskee todellisia juuria [75] , jotka ovat kaukana toisistaan, mutta yleistyksiä on sekä kompleksisiin [76] että todellisiin yhtäläisiin tai läheisiin juuriin [77] .
Iteratiiviset menetelmät algebrallisten yhtälöiden ratkaisemiseksi jaetaan stationäärisiin, kun funktio liittyy toiseen funktioon, jolla on samat juuret, iteraatioluvusta riippumatta [78] , ja ei-stationaarisiin, jolloin funktio voi riippua iteraatioluvusta. Yksinkertaisimpia stationaarisia iteratiivisia menetelmiä ovat sekanttimenetelmä (tai lineaarinen interpolointimenetelmä) ja tangenttimenetelmä (tai Newtonin menetelmä), jotka ovat vastaavasti ensimmäisen ja toisen asteen menetelmiä. Näiden menetelmien yhdistelmä, jossa peräkkäiset approksimaatiot sijaitsevat juuren vastakkaisilla puolilla, mahdollistaa nopeamman konvergenssin saavuttamisen [79] . Tšebyshevin menetelmä, joka perustuu käänteisfunktion laajentamiseen Taylorin kaavalla, mahdollistaa korkeamman asteen menetelmien rakentamisen erittäin nopealla konvergenssilla [80] . On olemassa myös Koenigin lauseeseen [81] ja Aitkenin menetelmään [82] perustuva menetelmä . Iteratiivisten menetelmien konvergenssin osoittamiseksi käytetään pakatun kuvauksen periaatetta [83] .
Sanakirjat ja tietosanakirjat | |
---|---|
Bibliografisissa luetteloissa |
Matematiikan alat | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
Portaali "Tiede" | ||||||||||
Matematiikan perusteet joukko teoria matemaattinen logiikka logiikan algebra | ||||||||||
Lukuteoria ( aritmetiikka ) | ||||||||||
| ||||||||||
| ||||||||||
| ||||||||||
| ||||||||||
| ||||||||||
|