Yhtälön ratkaisu, lineaariset yhtälöryhmät

Ratkaise yhtälörymä
A x-B y-C(x^2+y^2)⩵0
A y-B x⩵D x y
Solve-käskyllä ja mahdollisimman yleisessä muodossa. A ja B ovat vakioita.

Ratkaisu on hyvin monimutkainen:

[Graphics:Images/4_ratkaisut_gr_1.gif]
[Graphics:Images/4_ratkaisut_gr_2.gif]
[Graphics:Images/4_ratkaisut_gr_3.gif]
[Graphics:Images/4_ratkaisut_gr_4.gif]

Ratkaise lineaarinen yhtälöryhmä LinearSolve-komennolla:
3 x1+4x2-x3=5
4x2+x4=2
x1-x4+x2-1=0
x4-x2=2

Matriisit täytyy muodostaa siten, että x1,x2,x3 ja x4 ovat järjestyksessä. Vakiot viedään yhtäsuuruusmerkin vasemmalle puolelle (huomaa, että toiseksi viimeisessä yhtälössä vakion merkki muuttuu).

[Graphics:Images/4_ratkaisut_gr_5.gif]
[Graphics:Images/4_ratkaisut_gr_6.gif]
[Graphics:Images/4_ratkaisut_gr_7.gif]
[Graphics:Images/4_ratkaisut_gr_8.gif]

Taulukoiden muodostaminen ja interpolointi

Muodosta taulukko, jonka arvoina ovat {x,E^(-x)}-lukuparit, missä x saa arvoja 0:sta 5:een ½:n välein (eli 0,½,1,1½...4½,5). Sijoita taulukko muuttujaan tau.

[Graphics:Images/4_ratkaisut_gr_9.gif]
[Graphics:Images/4_ratkaisut_gr_10.gif]

Interpoloi edellä muodostettua taulukkoa

1. interpoloivan polynomin avulla
2. paloittain interpoloivan neliöllisen polynomin avulla

ja plottaa nämä sekä Exp[-x] samaan kuvaajaan väleillä 0-5 ja 0-10. Mitä tapahtuu kuvaajille, kun välinä on 0-10? Miksi? Osaatko piirtää kuvaajat eri värisinä?
Piirrä vielä erotuspolynomi Exp[-x]-InterpolantingPolynomial.

Interpoloiva polynomi muodostetaan käyttämällä InterpolatingPolynomial-funktiota.Ko.funktiota voi käyttää vain,kun interpoloitavia pisteitä on suhteellisen vähän,sillä funktio sovittaa n-1-asteisen polynomin pistejouukkoon,jolloin suurella n:llä suoritusaika nousee huomattavan suureksi.Käytämme N-funktiota tekemään kaavasta puhtaasti numeerisen.

[Graphics:Images/4_ratkaisut_gr_11.gif]
[Graphics:Images/4_ratkaisut_gr_12.gif]

Interpolation sovittaa alemman asteisen polynomin jokaiseen pisteväliin. Polynomin asteen voi asettaa InterpolationOrder-optiolla.

[Graphics:Images/4_ratkaisut_gr_13.gif]
[Graphics:Images/4_ratkaisut_gr_14.gif]

Seuraavassa plotataan kaikki kuvaajat eri väreillä käyttämällä PlotStyle-optiota. RGBColor-funktio puolestaan luo värin käyttämällä punaisen, vihreän ja sinisen sävyjä (jotka asetetaan parametreina).

[Graphics:Images/4_ratkaisut_gr_15.gif]
[Graphics:Images/4_ratkaisut_gr_16.gif]
[Graphics:Images/4_ratkaisut_gr_17.gif]

[Graphics:Images/4_ratkaisut_gr_18.gif]

[Graphics:Images/4_ratkaisut_gr_19.gif]

[Graphics:Images/4_ratkaisut_gr_20.gif]

[Graphics:Images/4_ratkaisut_gr_21.gif]

[Graphics:Images/4_ratkaisut_gr_22.gif]

[Graphics:Images/4_ratkaisut_gr_23.gif]

[Graphics:Images/4_ratkaisut_gr_24.gif]

[Graphics:Images/4_ratkaisut_gr_25.gif]

Ekstrapolointi ei onnistu, koska sovitettavat polynomit karkaavat äärettömyyksiin sovitusalueen ulkopuolella (kuten äärellistä kertalukua olevat polynomit aina tekevät):

[Graphics:Images/4_ratkaisut_gr_26.gif]

[Graphics:Images/4_ratkaisut_gr_27.gif]

[Graphics:Images/4_ratkaisut_gr_28.gif]
[Graphics:Images/4_ratkaisut_gr_29.gif]
[Graphics:Images/4_ratkaisut_gr_30.gif]
[Graphics:Images/4_ratkaisut_gr_31.gif]
[Graphics:Images/4_ratkaisut_gr_32.gif]

[Graphics:Images/4_ratkaisut_gr_33.gif]

[Graphics:Images/4_ratkaisut_gr_34.gif]

[Graphics:Images/4_ratkaisut_gr_35.gif]

[Graphics:Images/4_ratkaisut_gr_36.gif]

[Graphics:Images/4_ratkaisut_gr_37.gif]

[Graphics:Images/4_ratkaisut_gr_38.gif]
[Graphics:Images/4_ratkaisut_gr_39.gif]

[Graphics:Images/4_ratkaisut_gr_40.gif]

[Graphics:Images/4_ratkaisut_gr_41.gif]

Muodosta taulukko tau2, jossa on lukuparit {x,Sin[x]+satunnaislukuja väliltä 0.0-0.5} ja x saa arvoja 0-10 ykkösen välein.

SeedRandom tarvitaan satunnaislukugeneraattorin alustamista varten. Randomille voi antaa parametreja, joista ensimmäinen määrää sen, millaisia satunnaislukuja halutaan (Real=reaalilukuja) ja toinen antaa välin, jolta satunnaislukuja tuotetaan.

[Graphics:Images/4_ratkaisut_gr_42.gif]
[Graphics:Images/4_ratkaisut_gr_43.gif]

Sovita taulukkoon tau2 funktiot x, Sin[x] ja Cos[x]. Plottaa tulos ja taulukon tau pisteet (ListPlot) samaan kuvaajaan.

[Graphics:Images/4_ratkaisut_gr_44.gif]
[Graphics:Images/4_ratkaisut_gr_45.gif]
[Graphics:Images/4_ratkaisut_gr_46.gif]

[Graphics:Images/4_ratkaisut_gr_47.gif]

[Graphics:Images/4_ratkaisut_gr_48.gif]

[Graphics:Images/4_ratkaisut_gr_49.gif]

[Graphics:Images/4_ratkaisut_gr_50.gif]

[Graphics:Images/4_ratkaisut_gr_51.gif]

[Graphics:Images/4_ratkaisut_gr_52.gif]

Tee tilastollinen analyysi yo. sovituksesta.

Luetaan ensin pakkaus, jossa Regress-funktio on:

[Graphics:Images/4_ratkaisut_gr_53.gif]
[Graphics:Images/4_ratkaisut_gr_54.gif]
[Graphics:Images/4_ratkaisut_gr_55.gif]


Converted by Mathematica      February 3, 2005