MathCad_Complex_Matrices_Tutorial_Krowiak.pdf

(215 KB) Pobierz
167500076 UNPDF
6.05.2008
T WORZENIE ORAZ OPERACJE NA MACIERZACH
O ELEMENTACH BĘDĄCYMI LICZBAMI ZESPOLONYMI W
PROGRAMIE M ATH C AD ®
Autor:
Radosław Krowiak
Spis treści
1.Macierze liczb zespolonych w programie MathCad........................................................................2
1.1.Tworzenie macierzy liczb zespolonych krok po kroku [1,2]....................................................2
1.2.Zapis moduł-faza liczb zespolonych w programie MathCad...................................................3
2.Operacje na macierzach....................................................................................................................4
3.Przykłady zastosowań oraz interakcji z systemem operacyjnym [1]...............................................5
Bibliografia:................................................................................................................................7
1
1. Macierze liczb zespolonych w programie MathCad.
1.1. Tworzenie macierzy liczb zespolonych krok po kroku [1,2].
1. Wcisnąć kombinację klawiszy [CTR M] lub kliknąć lewym
klawiszem myszki na ikonkę Matrix or Vector na pasku
Macierzy, jak pokazuje zrzut ekranu obok:
2. Wpisać pożądane wymiary macierzy w okienku,
które ukaże się po przejściu kroku pierwszego, jak
na zrzucie obok:
3. Wypełnić macierz pożądanymi elementami macierz używając klawisza [TAB] do
przesuwania się między pozycjami, jak na załączonym zrzucie:
Uwaga!
Należy pamiętać o symbolicznym zapisie liczb zespolonych w programie MathCad, otóż
jednostkę urojoną i zapisujemy:
DOBRZE BŁĄD
1 i 1* i
Jest tak ponieważ w drugim przypadku MathCad odczyta i jako zupełnie nową zmienną.
2
167500076.048.png 167500076.051.png 167500076.052.png
1.2. Zapis moduł-faza liczb zespolonych w programie MathCad.
Dla przypomnienia podaję równanie Euler'a opisujące przejście między dwoma zapisami:
z = r e i = r ⋅cos isin 
r = a 2 b 2
= arctg b
a
W MathCadzie nie istnieje domyślny operator zapisu liczb zespolonych moduł-faza, jednak
bardzo łatwo takowy utworzyć:
1. Zgodnie z przyjętymi powyżej oznaczeniami tworzymy nową funkcję operującą
znakiem (sam znak kopiujemy z Extra Math Symbols w dziale quicksheets):
∢ r,  : = r ⋅cos i ⋅sin
2. Wywołanie funkcji odbywać będzie się poprzez pasek Evaluation:
Kliknięcie na zaznaczony na czerwono przycisk spowoduje
pojawienie się wypełniacza, który należy wypełnić zgodnie z przyjętą
notacją funkcji:
Na przykład:
3 Ð
è
3
ø
=
1.5 2.598i
,
lub: 1 Ð
è
4
ø
=
0.707 0.707i
Adnotacja:
Niekoniecznie musimy używać operatora - równie dobrze możemy jawnie nazwać
definiowaną funkcję, np. „phasor” czy też inną pasująca nam nazwę a funkcjonalność jej
się nie zmieni, co ukazuje poniższy przykład:
( ) phasor p
3
è
ø
=
0.5 0.866i
3
{
{ z = a ib
+
+
1
+
167500076.053.png 167500076.001.png 167500076.002.png 167500076.003.png 167500076.004.png
2. Operacje na macierzach [1].
W niniejszym rozdziale przedstawię sposoby operacji na macierzach w programie
MathCad 13. Operacje przedstawiam w przejrzystej formie tabelki:
# Wywołanie Wygląd Komentarz odnośnie wyniku lub działania
1
|
|M|
Zwraca wartość wyznacznika macierzy (dla wektora
zwraca normę, wartość bezwzględną dla skalara)
2 ^−1
M^−1
Zwraca macierz będącą odwrotnością M, jeśli macierz
taka nie istnieje wysłany zostaje odpowiedni komunikat o
błędzie
3 [Ctrl 1]
M^T
Zwraca macierz transponowaną do M NxM (zamienia
kolumny z wierszami) o wymiarze MxN
4
+
x + y
Wynik jest:
– skalarem – sumą x oraz y, jeśli x i y to skalary
– macierzą – jeśli X to macierz a y – skalar, każdy element
wynikowej macierzy jest sumą y oraz odpowiedniego el.
X
– macierzą – sumą X oraz Y, jeśli x i y to macierze o tych
samych wymiarach
5
x – y Patrz komentarz wyżej
6 [Ctrl 8]
x × y Wynikiem jest iloczyn wektorowy macierzy X oraz Y
7
*
x * y Wynikiem jest iloczyn skalarny macierzy X oraz Y
8 [Ctrl /]
x / y Wynikiem jest iloczyn skalarny macierzy: x * y^-1
9
^
y ^ x
Dla macierzy wynik zależy od x:
– x = 0, zwraca macierz jednostkowa o wymiarach takich
samych co Y
– x = –1, zwraca odwrotność macierzy Y
– x>0, wynikiem jest macierz Y pomnożona przez siebie x
razy
– x<0 (bez –1), wynikiem jest odwrotność macierzy Y
pomnożona przez siebie x razy
10 eigenvals eigenvals(M)
Zwraca wektor, którego elementami są posortowane
wartości własne macierzy M
11 eigenvec
eigenvec(M,
z)
Zwraca pojedynczy, znormalizowany wektor, skojarzony z
wartością własną z macierzy M
12
[
M[m,n
Zwraca element macierzy będący przecięciem m-tego
wiersza oraz n-tej kolumny
Na następnej stronie zamieszczam przykłady operowania na macierzach o
współczynnikach symbolicznych. Dla uproszczenia ilustracji zakładam iż wymiary
macierzy wynoszą 2x2. Operacje nie obsługiwane przez MathCada symbolicznie są
zilustrowane przykładami numerycznymi.
4
167500076.005.png
# Wygląd
Przykład operacji:
ç
ç
è
2 1i
+
100
54
1i
2
7
2
3
1
÷
÷
ø
ç
ç
è
16.732 3.589i
+
÷
÷
ø
0 Założenia:
M
:=
eigenvals M
( )
=
-
8.038
-
7.458i
-
3.694
+
4.869i
1
|M|
ç
è
a 1
a 3
a 2
a 4
÷
ø
®
a 1 a 4
×
-
a 2 a 3
×
ç
è
a 1
a 3
a 2
a 4
÷
ø
-
1
ç
è
a 4
a 3
-
a 2
÷
ø
1
2
M^−1
®
×
a 1 a 4
×
-
a 2 a 3
×
-
a 1
ç
è
a 1
a 3
a 2
a 4
÷
ø
T
ç
è
a 1
a 2
a 3
a 4
÷
ø
3
M^T
®
ç
è
a 1
a 3
a 2
a 4
÷
ø
ç
è
b 1
b 3
b 2
b 4
÷
ø
ç
è
a 1
+
b 1
a 2
+
b 2
÷
ø
4
x + y
+
®
a 3
+
b 3
a 4
+
b 4
ç
è
a 1
a 3
a 2
a 4
÷
ø
ç
è
b 1
b 3
b 2
b 4
÷
ø
ç
è
a 1
-
b 1
a 2
-
b 2
÷
ø
5
x – y
-
®
a 3
-
b 3
a 4
-
b 4
ç
ç
è
13
3
÷
÷
ø
ç
ç
è
33
6
÷
÷
ø
ç
ç
è
54 100 1i
-
×
÷
÷
ø
6
x × y
-
50
´
-
+
82
2i
®
584
21 26 1i
+
×
ç
è
a 1
a 3
a 2
a 4
÷
ø
ç
è
b 1
b 3
b 2
b 4
÷
ø
ç
è
a 1 b 1
×
+
a 2 b 3
×
a 1 b 2
×
+
a 2 b 4
×
÷
ø
7
x * y
×
®
a 3 b 1
×
+
a 4 b 3
×
a 3 b 2
×
+
a 4 b 4
×
ç
è
a 1
a 3
a 2
a 4
÷
ø
ê
ê
ê
ê
ë
( )
-
a 1
×
b 4
+
a 2
×
b 3
a 1
×
b 2
-
a 2
×
b 1
ú
ú
ú
ú
û
( )
( )
( )
( )
-
b 1
×
b 4
+
b 2 b 3
×
-
b 1
×
b 4
+
b 2 b 3
×
-
b 1
×
b 4
+
b 2 b 3
×
-
b 1
b 4
+
b 2 b 3
×
8
x / y
®
ç
è
b 1
b 3
b 2
b 4
÷
ø
b 4
b 3
b 2
b 1
( )
-
a 3
×
+
a 4
×
a 3
×
-
a 4
×
( )
( )
( )
( )
-
b 1
×
b 4
+
b 2 b 3
×
-
b 1
×
b 4
+
b 2 b 3
×
-
b 1
×
b 4
+
b 2 b 3
×
-
b 1
b 4
+
b 2 b 3
×
9
y ^ x
ê
ê
ê
ê
ê
ê
ë
1
2 b 1
1
2 b 4
1
2
ë
( )
2
( )
2
û
1
2
ú
ú
ú
ú
ú
ú
û
ç
è
ç
è
b 1
b 3
b 2
b 4
÷
ø
÷
ø
×
+
×
+
×
b 1
-
2 b 1 b 4
×
×
+
b 4
+
4 b 2 b 3
×
×
10 eigenvals(M) eigenvals
®
1
2
1
2 b 1
1
2 b 4
1
2
ë
( )
2
( )
2
û
×
+
×
-
×
b 1
-
2 b 1 b 4
×
×
+
b 4
+
4 b 2 b 3
×
×
ç
ç
è
0.044 0.056i
-
÷
÷
ø
11 eigenvec(M, z)
eigenvec M 3.694
(
,
-
+
4.869i
)
=
-
0.771
-
0.018
4.794i 10
´
-
4
-
+
0.632i
12
M[m,n
M 1 2
,
=
3
5
×
×
167500076.006.png 167500076.007.png 167500076.008.png 167500076.009.png 167500076.010.png 167500076.011.png 167500076.012.png 167500076.013.png
 
167500076.014.png 167500076.015.png 167500076.016.png
 
167500076.017.png
 
167500076.018.png 167500076.019.png 167500076.020.png 167500076.021.png 167500076.022.png 167500076.023.png 167500076.024.png 167500076.025.png 167500076.026.png 167500076.027.png 167500076.028.png 167500076.029.png 167500076.030.png 167500076.031.png 167500076.032.png 167500076.033.png 167500076.034.png 167500076.035.png 167500076.036.png 167500076.037.png 167500076.038.png 167500076.039.png 167500076.040.png 167500076.041.png 167500076.042.png 167500076.043.png 167500076.044.png 167500076.045.png 167500076.046.png 167500076.047.png 167500076.049.png 167500076.050.png
 
Zgłoś jeśli naruszono regulamin