MatheAss 10.0 − Linjär algebra

Linjära ekvationssystem

Programmet bestämmer lösningsvektorn för ett system av linjära ekvationer med n ekvationer och n obekanta.

Exempel: Söker man en parabel genom punkterna P(1/3), Q(2/1) och R(4/9), leder detta till ekvationssystemet

  1·x1 + 1·x2 + 1·x3 = 3
  4·x1 + 2·x2 + 1·x3 = 1
 16·x1 + 4·x2 + 1·x3 = 9   
            
 L = ( 2; -8; 9 )

Parabeln har alltså ekvationen y = 2x² - 8x + 9.

Exempel med tvådimensionellt lösningsrum:

  0·x1 + 0·x2 + 2·x3 - 1·x4  =  1
  1·x1 + 1·x2 + 1·x3 + 1·x4  =  4
  2·x1 + 2·x2 - 4·x3 + 5·x4  =  5
  1·x1 + 1·x2 - 7·x3 + 5·x4  =  0

  L = { ( 3,5-s-1,5t; s; 0,5+0,5t; t ) | s,t ∈ R }

Linjär optimering (från februari 2022)

Programmet bestämmer den optimala lösningen för en målfunktion med två variabler med linjära olikheter som randvillkor.

Målfunktion:   
  ƒ(x,y) = 140·x + 80·y → Maximum

Randvillkor:
  x ≥ 0
  y ≥ 0
  x ≤ 600
  y ≤ 700
  x + y ≤ 750
  3·x + y ≤ 1200

Maximum
  x = 225   y = 525
  ƒ(x,y) = 73500

Programmet bestämmer linjärkombinationen av en vektor från tre givna vektorer. Rutinen kan också användas för att pröva den linjära oberoendet av tre vektorer i rummet, det vill säga om de ligger i ett plan.

    ⎧ 1 ⎫     ⎧ 2 ⎫     ⎧ 0 ⎫   ⎧ 2 ⎫
  a·⎪ 2 ⎪ + b·⎪ 1 ⎪ + c·⎪ 1 ⎪ = ⎪ 3 ⎪   
    ⎩ 0 ⎭     ⎩ 1 ⎭     ⎩ 0 ⎭   ⎩ 7 ⎭

Lösning :
  a = -12  b = 7  c = 20

Skalärprodukt

Programmet beräknar för två vektorer deras skalärprodukt, längden av de båda vektorerna och den inneslutna vinkeln.

   ->  ⎧ 1 ⎫   ->  ⎧ 5 ⎫
   a = ⎪ 3 ⎪   b = ⎪ 0 ⎪
       ⎩ 1 ⎭       ⎩ 3 ⎭

Skalärprodukt av vektorerna = 8

Längd av första vektorn    = √11 = 3,3166248
Längd av andra vektorn     = √34 = 5,8309519

Innesluten vinkel  α = 65,564402°

Vektorprodukt

Programmet beräknar för två vektorer deras vektorprodukt samt dess belopp. Vektorprodukten står vinkelrät mot det parallellogram som de spänner upp, och dess belopp är lika med parallellogrammets ytinnehåll.

     ->  ⎧ 1 ⎫     ->  ⎧ 7 ⎫
     a = ⎪ 2 ⎪     b = ⎪ 1 ⎪
         ⎩ 3 ⎭         ⎩ 4 ⎭

 ->  ->  ⎧  5 ⎫    ->  ->  
 a x b = ⎪ 17 ⎪   |a x b|= √483 = 21,977261
         ⎩-13 ⎭

Rumsprodukt

Programmet beräknar för tre vektorer rumsprodukten. Dess belopp anger volymen av den förskjutna rätblock (spat), som spänns upp av de tre vektorerna.

  ->  ⎧ 2 ⎫     ->  ⎧ 2 ⎫     ->  ⎧ 3 ⎫   
  a = ⎪ 3 ⎪     b = ⎪-1 ⎪     c = ⎪ 9 ⎪
      ⎩ 5 ⎭         ⎩ 7 ⎭         ⎩ 2 ⎭

    ->  ->    ->  
  ( a x b ) · c = 26

Matrisinversion

Programmet beräknar för en kvadratisk matris av ordning n determinanten, rangen och den inversa matrisen.

Matris: 

  ⎧  0  1  1  ⎫
  ⎪  0  1  3  ⎪
  ⎩  2  0  1  ⎭

Invers matris

  ⎧  0,25 -0,25  0,5  ⎫
  ⎪   1,5  -0,5    0  ⎪
  ⎩  -0,5   0,5    0  ⎭

Ordning = 3,  Rang = 3,  Determinant = 4

Pseudoinvers matris

Om kolumnerna i en matris A är linjärt oberoende, är A^T· A inverterbar och man erhåller med följande formel pseudoinversen:

A⁺ = ( A^T· A )^-1· A^T

Där är A⁺ en vänsterinvers till A , det vill säga det gäller: A⁺ · A = E .

Matris A
¯¯¯¯¯¯¯¯
  ⎧  1  1  1  1  ⎫
  ⎩  5  7  7  9  ⎭

A^T· A
¯¯¯¯
  ⎧  26  36  36  46  ⎫
  ⎪  36  50  50  64  ⎪
  ⎪  36  50  50  64  ⎪
  ⎩  46  64  64  82  ⎭

A^T· A ej inverterbar

A · A^T
¯¯¯¯¯¯
  ⎧   4   28  ⎫
  ⎩  28  204  ⎭

( A · A^T )^-1
¯¯¯¯¯¯¯¯¯¯¯¯¯
  ⎧  6,375 -0,875  ⎫
  ⎩ -0,875  0,125  ⎭

Högerinvers:  A^T · ( A · A^T )^-1
¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯
  ⎧     2  -0,25  ⎫
  ⎪  0,25      0  ⎪
  ⎪  0,25      0  ⎪
  ⎩  -1,5   0,25  ⎭

Matrismultiplikation

Programmet beräknar produktmatrisen för två matriser.

 1. Matris           2. Matris

  ⎧  1  2  3  0  ⎫   ⎧  1   2   3   4  ⎫
  ⎪  4  5  6  0  ⎪   ⎪  5   6   7   8  ⎪
  ⎪  0  0  0  0  ⎪   ⎪  9  10  11  12  ⎪
  ⎩  0  0  0  0  ⎭   ⎩  0   0   0   0  ⎭

 Produktmatris
  
  ⎧  38  44   50   56  ⎫
  ⎪  83  98  113  128  ⎪
  ⎪   0   0   0    0   ⎪
  ⎩   0   0   0    0   ⎭

MatheAss 10.0 − Linjär algebra

A+ = ( AT· A )-1· AT

A⁺ = ( A^T· A )^-1· A^T