MNNE:Algebra

Z Skrypty dla studentów Ekonofizyki UPGOW

(Różnice między wersjami)
m
m (Wektory)
Linia 14: Linia 14:
Zmiana bazy  
Zmiana bazy  
 +
 +
 +
==  Zmiana bazy i iloczyn skalarny==
 +
Mamy <math>(x,y)=\sum_{i=1}^{N}x_i y_j</math>
 +
<source lang='matlab'>
 +
b1=[1;1;1]
 +
b2=[1;0;1]
 +
b3=[1;0;-1]
 +
C=[b1,b2,b3]
 +
C'*C
 +
M=[b1'*b1,b1'*b2,b1'*b3;b2'*b1,b2'*b2,b2'*b3;b3'*b1,b3'*b2,b3'*b3]
 +
</source>
 +
 +
Ortogonalizacja Grama-Schmidta
 +
 +
<source lang='matlab'>
 +
e1=b1
 +
e2=b2-(e1'*b2)/(e1'*e1)*e1
 +
e3=b3-( (e1'*b3)/(e1'*e1)*e1 + (e2'*b3)/(e2'*e2)*e2 )
 +
</source>
 +
 +
<source lang='matlab'>
 +
E=[e1,e2,e3]
 +
quiver3([0,0,0],[0,0,0],[0,0,0], C(1,:), C(2,:), C(3,:), 0,'r')
 +
hold on
 +
quiver3([0,0,0],[0,0,0],[0,0,0], E(1,:), E(2,:), E(3,:), 0,'b')
 +
xlim([-1,1])
 +
ylim([-1,1])
 +
zlim([-1,1])
 +
hold off
 +
 +
</source>
 +
  Ćwiczenia:  
  Ćwiczenia:  

Wersja z 21:34, 23 lut 2011

Wektory

Wektor w przestrzeni euklidesowej N wymiarowej jest reprezentowany przez N liczb. Iloczyn skalarny dwóch wektorów \(\mathbf a\) oraz \(\mathbf b\) o oznaczany symbolem \(\mathbf a \cdot \mathbf b\) określony jest jako:

\(\mathbf a \cdot \mathbf b = \sum_{i=1}^{N}a_i b_i\),

gdzie \(a_i\) to i-ty element wektora a. Można pokazać, że iloczyn ten jest też dany przez

\(\mathbf a \cdot \mathbf b = \|\mathbf a\| \|\mathbf b\| \cos \theta\),

gdzie \(\theta\) jest kątem między \(\mathbf a\) a \(\mathbf b\).

Jeśli jeden wektorów jest wektorem o długości jeden to mnoże go przez dowolny inny wektor może być interpretowane jako rzutowanie na kierunek wyznaczony przez pierwszy wektor.

Baza

Zmiana bazy


Zmiana bazy i iloczyn skalarny

Mamy \((x,y)=\sum_{i=1}^{N}x_i y_j\)

b1=[1;1;1]
b2=[1;0;1]
b3=[1;0;-1]
C=[b1,b2,b3]
C'*C
M=[b1'*b1,b1'*b2,b1'*b3;b2'*b1,b2'*b2,b2'*b3;b3'*b1,b3'*b2,b3'*b3]

Ortogonalizacja Grama-Schmidta

e1=b1
e2=b2-(e1'*b2)/(e1'*e1)*e1
e3=b3-( (e1'*b3)/(e1'*e1)*e1 + (e2'*b3)/(e2'*e2)*e2 )
E=[e1,e2,e3]
quiver3([0,0,0],[0,0,0],[0,0,0], C(1,:), C(2,:), C(3,:), 0,'r')
hold on
quiver3([0,0,0],[0,0,0],[0,0,0], E(1,:), E(2,:), E(3,:), 0,'b')
xlim([-1,1])
ylim([-1,1])
zlim([-1,1])
hold off


Ćwiczenia: 
  • obliczyć współrzędne wektora \(a=(1,2,3)\) w bazie \((e_1+e_2...\)
  • obliczyc rząd macierzy