Törmäsin vanhaan tiedostooni vuodelta 2004, siinä oli ratkaistu hankalan vastussysteemin resistanssi käyttäen taulukkolaskentaa. Tiedosto oli alun perin tehty Excelissä ja nytpä se ei enää toiminut oikein. Arvelen, että alun perin olin tehnyt tämän viime vuosituhannella. Päätinpä tehdä yhtälöryhmän ratkaisun uudelleen käyttäen matriisifunktioita.

Esitän menetelmän Google Sheetsissä, menetelmä toimii samalla tavoin Excelissä ja LibreOfficen Calcissa. Lopuksi teen saman GeoGebralla matriisikomennoilla.

resistanssiongelma

Ongelma on oheisessa kuvassa. Pitää ratkaista kyseisen vastussysteemin resistanssi pisteiden A ja B välillä. Muistaakseni ongelma oli alkuperäisessä muodossa sellainen, että kaikkien vastusten resistanssi oli 1 Ω. Sen ratkaiseminen onnistunee lahjakkaalta fyysikolta pelkästään katsomalla tuota kytkentäkaaviota :o)

ratkaisu Kirchoffin laeilla

Itse en ainakaan keksi mitään simppeliä kaavaa resistanssin laskemiseksi. Tehdään Kirchoffin lakien avulla yhtälöt eri osissa kulkevista virroista, ratkaistaan niiden suuruudet ja käytetään sitä tietoa kokonaisresistanssin laskemiseen.

Merkitään 1 Ω:n vastuksen läpi kulkeva virtaa I₁:llä, 2 Ω:n läpi kulkeva virtaa I₂:lla, jne. Virtojen suunta on vasemmalta oikealle. Kuvitellaan, että pisteiden A ja B välillä jännite U = 10 V:n. Näin saadaan Kirchoffin lakien avulla helposti viisi yhtälöä

Kokonaisresistanssi R tulee olemaan 

ratkaisu Google Sheetsillä

Laitan Google Sheetsiin resistanssit soluihin B6:B10 ja jännitteen soluuun B12.

Yhtälöryhmän kertoimet ovat soluissa A17:F21. Laitan ne soluihin viittauksena R:n ja U:n arvoihin, eli solussa A17 on =-B6, solussa C17 on =B8 jne.

Näytä-valikon Näytä kaavat-komennolla kaavat saa näkyviin ja tietysti takaisin näkyville kuten näissä kuvissa.

Soluun A25 lasketaan kerroinmatriisin determinantti komennolla

=MDETERM(A17:E21)

Mikäli determinantti on nollasta poikkeava, niin yhtälöryhmällä on ratkaisu.

Ratkaistaan virtojen arvot soluihin A28:E28. Kirjoitetaan soluun A28 komento

=TRANSPOSE(MMULT(MINVERSE(A17:E21);F17:F21))

Suomenkielisessä Excelissä ja LibreOfficessa tuo rivi pitäisi kirjoittaa muodossa

=TRANSPONOI(MKERRO(MKÄÄNTEINEN(A17:E21);F17:F21))

Kaavassa lasketaan kerroinmatriisin käänteismatriisin ja vakiotermivektorin tulo ja käännetään se vaakasuuntaan.

Kokonaisresistanssi lasketaan soluihin A31 ja A32

Tätä kirjoittaessa löysin vanhemman version alkuperäisestä Excel-tiedostosta, se toimi moitteettomasti. Siinä solun A28 ratkaisu oli tehty LINREGR-funktiolla

=LINREGR(F17:F21;A17:E21;EPÄTOSI;EPÄTOSI)

GeoGebran CAS 6×6 yhtälöryhmä

GeoGebran CAS:issa ongelma olisi ratkennut vaikkapa seuraavasti.

Solun 7 ratkaisu kertoo, että kokonaisresistanssi R on riippumaton U:sta niin kuin pitääkin.

Solussa 8 on sijoitus

Lasketaan samalla kuinka suuri kokonaisresistanssi olisi, jos kaikkien vastusten resistanssi olisi 1.

Nyt kun tietää ratkaisun, niin saattaa olla helpompaa nähdä miksi näin käy tuijottamalla pelkästään kytkentäkaaviota.

GeoGebran matriisikomennot

Edellistä CAS-versiota tehdessä oli aika työlästä kirjoittaa yhtälöitä muuttujineen, mitäpä jos halutaan tehdä  GeoGebralla saman tapainen matriisityylinen ratkaisu kuin Google Sheetsissä.

Kopioidaan Sheetsratkaisusta alkuperäinen matriisi ja sijoitetaan se GeoGebran taulukkolaskentaan alueelle A2:F6. Valitaan alue A2:E2 eli ensimmäisen yhtälön kertoimet ja hiiren oikealla painikkeella valitaan Luo → Lista.

Luodaan samalla tavalla muiden yhtälöiden kertoimista listat l2, l3, l4 ja l5 sekä vakiotermeistä lista l6.

Muodostetaan kerroinlistoista matriisi kirjoittamalla CAS:iin

M:={l1, l2, l3, l4, l5}

Dererminantin arvoksi saadaan 157 komennolla

Determinantti(M)

Annetaan yhtälöryhmän ratkaisun nimeksi Ra ja ratkaistaan se näppäilemällä 

Ra:=M^-1*l6

Kokonaisresistanssi saadaan käyttämällä vakiotermilistan viidettä ja  ratkaisulistan Ra ensimmäistä ja kolmatta jäsentä.

R=l6(5)/(Ra(1)+Ra(3))

linkit

Google Sheets-tiedosto https://docs.google.com/spreadsheets/d/1aW1ZqA-rO2MnbvGXrVf9V8fZVz_58nVNwzjC-jEqWko/edit?usp=sharing

GeoGebra-tiedosto yhtälöryhmänä https://www.geogebra.org/m/zhjndgt7

Vaikuttaa siltä, että alla olevat ei avaudu geogebra.orgissa ja Classic 6:lla. Jos haluat tiedoston omalle koneellesi, niin klikkaa tuota linkkiä ja oikeasta yläkulmasta kolmen pisteetn takaa Tietoa ja Lataa koneellesi

GeoGebra-tiedosto matriisikomentoratkaisuun https://www.geogebra.org/m/jqgmknb7