Wiki

Metodo numerikoen konparatibak

Taula hauek ez dituzte gidak ordezkatzen: metodoa hipotesien, kostuaren eta espero den portaeraren arabera aukeratzeko ikuspegi azkarra dira.

Newton vs sekantea vs bisekzioa

Erro eskalarretarako hiru estrategia: segurtasun globala, abiadura lokala eta iterazio bakoitzeko kostua.

MetodoaHipotesiakOrdenaKostuaNoiz erabili
BisekzioaJarraitutasuna eta zeinu-aldaketa [a,b][a,b] tarteanLineala; tartea 2z zatitzen daff-ren ebaluazio berri bat iterazio bakoitzeanErroa inguratzeko eta hasierako estimazio fidagarria lortzeko
Newton-RaphsonErro sinplea, ff' eskuragarri eta hasierako puntua hurbilKoadratiko lokalaff eta ff' iterazio bakoitzeanHazi ona eta deribatu merkea daudenean
SekanteaBi hasierako estimazio eta ff nahikoa erregularraSuperlineala, p1.618p\approx1.618ff-ren ebaluazio berri bat iterazio bakoitzeanff' ebaluatu nahi edo ezin denean
Irristatu horizontalki taula osoa ikusteko →

Euler vs Heun vs RK4

PVIetarako urrats bakarreko metodoak: zehaztasuna, ff-ren ebaluazioak eta erabilera praktikoa.

MetodoaIdeiaErrore globalaEbaluazioakNoiz erabili
Euler esplizituaTartearen hasierako maldarekin aurrera egiten duO(h)\mathcal{O}(h)1 pauso bakoitzeanEstimazio azkarretarako edo iragarle gisa
HeunEuler-en malda eta zuzendutako malda batezbestekatzen dituO(h2)\mathcal{O}(h^2)2 pauso bakoitzeanEuler baino hobekuntza merkea nahi denean
RK4Lau malda konbinatzen ditu 1,2,2,11,2,2,1 pisuekinO(h4)\mathcal{O}(h^4)4 pauso bakoitzeanEkuazio inpliziturik ebatzi gabe zehaztasun handia behar denean
Irristatu horizontalki taula osoa ikusteko →

Jacobi vs Gauss-Seidel vs SOR

Metodo iteratibo linealak A=MNA=M-N partizio gisa eta beren iterazio-matrizearen bidez.

MetodoaPartizioaEguneratzeaAbantailaArriskua
JacobiM=DM=DAurreko iterazioko balioak bakarrik erabiltzen dituSinplea eta paralelizagarriaAstiro konbergitu dezake
Gauss-SeidelM=DLM=D-LBalio berriak kalkulatu bezain laster berrerabiltzen dituAskotan Jacobi hobetzen du parametro gehigarririk gabeEkuazioen ordenaren menpe dago
SORM=1ωDLM=\frac{1}{\omega}D-LGauss-Seidel erlaxatzen du ω\omega pisu batekinω\omega egokiarekin asko azkar dezakeω\omega txar batek konbergentzia okertu edo hautsi dezake
Irristatu horizontalki taula osoa ikusteko →

Lagrange vs Newton vs Hermite vs splineak

Datu berak interpolatzeko lau modu: zer informazio erabiltzen duten, zenbat kostatzen den eraikitzea eta noiz komeni den bakoitza.

InterpolatzaileaErabiltzen dituen datuakEraikuntzaIndarguneaMuga
Lagrangef(xi)f(x_i) balioak n+1n+1 nodotani(x)\ell_i(x) oinarri esplizitua, ebatzi beharreko sistemarik gabeForma itxia, teoria eta kuadratura-dedukzioetarako aproposaNodo bat gehitzeak oinarri osoa berregitera behartzen du
Newtonf(xi)f(x_i) balioak n+1n+1 nodotanDiferentzia zatituen taula, forma habiaratuaNodo bat gehitzeak diagonal bat gehiago baino ez du kostatzenLagrangeren polinomio bera da: gradu handiko muga berak ditu
Hermitef(xi)f(x_i) balioak eta f(xi)f'(x_i) deribatuakDiferentzia zatituak nodo errepikatuekinNodo bakoitzeko informazioa bikoizten du: 2n+12n+1 gradua eta doitasun handiagoaDeribatuak ezagutzea eskatzen du, eta datuetan ez daude beti
Spline kubikoaf(xi)f(x_i) balioak nodo guztietanSistema tridiagonala momentuentzat, polinomio kubikoa zatikoLeuna (C2C^2) eta Runge fenomenorik gabea nodo askorekinEz du polinomio global bakarra ematen eta sistema bat ebatzi behar da
Irristatu horizontalki taula osoa ikusteko →

Trapezioa vs Simpson vs erdiko puntua vs Gauss

Eguneroko kuadraturak alderatuta: nodoak, errore konposatuaren ordena eta zein egoeratan irabazten duten.

FormulaNodoakErrore konposatuaDoitasun-graduaNoiz erabili
Trapezio konposatuaEkidistanteak, itxia (muturrak erabiltzen ditu)O(h2)\mathcal{O}(h^2)1Datu tabulatuak; sinplea eta sendoa
Simpson konposatuaEkidistanteak, itxia, azpitarte kopuru bikoitiaO(h4)\mathcal{O}(h^4)3Integrakizun leun tabulatua: bi ordena gehiago nodo bakoitzeko kostu berean
Erdiko puntu konposatuaEkidistanteak, irekia (muturrak saihesten ditu)O(h2)\mathcal{O}(h^2)1Integrakizuna muturretan ebaluatu ezin denean
Gauss-LegendreLegendre polinomioen erroak, ez ekidistanteak2n12n-1 gradua nn nodorekin2n12n-1ff nahi duzun tokian ebaluatu badezakezu eta ebaluazioak garestiak badira
Irristatu horizontalki taula osoa ikusteko →

AB2 vs AM2 vs iragarle-zuzentzailea

2. ordenako urrats anitzekoak aurrez aurre: pauso bakoitzeko kostua, errore-konstantea eta zenbat ordaintzen den egonkortasun inplizituagatik.

MetodoaMotaErrore lokalaKostua pauso bakoitzekoOharrak
AB2Esplizitua, bi pausokoa512h3y(ξ)\frac{5}{12}h^3y'''(\xi)ff-ren ebaluazio berri batAbioko balio bat behar du (Heun edo RK4) eta sare uniformea
AM2 (trapezio inplizitua)Inplizitua, pauso berri bat112h3y(ξ)-\frac{1}{12}h^3y'''(\xi)yk+1y_{k+1}-en ekuazio bat ebatziAB2 baino 5 aldiz konstante txikiagoa eta egonkortasun askoz handiagoa
Iragarle-zuzentzailea AB2+AM2Esplizitua praktikanZuzentzailearen ordenakoaff-ren 2 ebaluazioAM2-ren doitasuna ekuazio inplizitua ebatzi gabe; iragarle-zuzentzaile diferentziak errorea estimatzen du
Irristatu horizontalki taula osoa ikusteko →