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