質問編集履歴
4
マークダウンで累乗**が太字になっていたので^で書き直しました。
test
CHANGED
File without changes
|
test
CHANGED
@@ -34,6 +34,6 @@
|
|
34
34
|
|
35
35
|
欲しい行列
|
36
36
|
|
37
|
-
m.T * m = [a, b].T * [a, b] = [[a
|
37
|
+
m.T * m = [a, b].T * [a, b] = [[a^2, ab],
|
38
38
|
|
39
|
-
[ab, b
|
39
|
+
[ab, b^2]]
|
3
壊滅的に質問が悪いようなので、質問を全面的に刷新いたしました。回答いただいた方々、申し訳ございません。
test
CHANGED
@@ -1 +1 @@
|
|
1
|
-
python:列ベクトルx行ベクトルの
|
1
|
+
python:列ベクトルx行ベクトルの積を行列として出力する方法
|
test
CHANGED
@@ -1,43 +1,39 @@
|
|
1
1
|
### 前提・実現したいこと
|
2
|
+
|
3
|
+
最終目標
|
4
|
+
|
5
|
+
鏡の面方向ベクトルmと光の進行方向ベクトルnが与えられたときに反射後の光の進行方向を算出すること
|
2
6
|
|
3
7
|
|
4
8
|
|
5
|
-
|
9
|
+
反射後の光の進行方向を算出するには、以下の手順が必要になります。
|
6
10
|
|
7
|
-
|
11
|
+
1. 鏡の面方向ベクトルmから反射用の回転行列を作成する。
|
8
12
|
|
9
|
-
よろしくお願いいたします。
|
10
|
-
|
11
|
-
|
13
|
+
2. 回転行列と光の進行方向ベクトルnの積によって反射後の進行方向ベクトルを算出する。
|
12
14
|
|
13
15
|
|
14
16
|
|
15
|
-
|
17
|
+
課題
|
18
|
+
|
19
|
+
手順1で回転行列が作成できずに困っています。
|
20
|
+
|
21
|
+
ある行ベクトル(m行1列)とその転置ベクトル(1行m列)の積で行列(m行m列)ができれば問題は解決します。
|
22
|
+
|
23
|
+
*後学のために、その方法が、列ベクトル(m行1列)と行ベクトル(1行n列)をそのまま行列(m行n列)にできる方法かどうか教えていただけると幸いです。
|
16
24
|
|
17
25
|
|
18
26
|
|
19
|
-
|
27
|
+
具体例:
|
20
28
|
|
21
|
-
|
29
|
+
ベクトル
|
22
30
|
|
23
|
-
|
31
|
+
m = [a, b]
|
24
|
-
|
25
|
-
[[1/√2],と[1/√2, -1/√2]の積が[[1/2, -1/2],
|
26
|
-
|
27
|
-
[-1/√2]] [-1/2, 1/2]]
|
28
32
|
|
29
33
|
|
30
34
|
|
31
|
-
|
35
|
+
欲しい行列
|
32
36
|
|
33
|
-
|
37
|
+
m.T * m = [a, b].T * [a, b] = [[a ** 2, ab],
|
34
38
|
|
35
|
-
n = np.array([np.sqrt(0.5), - np.sqrt(0.5)])
|
36
|
-
|
37
|
-
|
39
|
+
[ab, b ** 2]]
|
38
|
-
|
39
|
-
n_matrix = np.dot(n, n.T)
|
40
|
-
|
41
|
-
print(m_matrix, n_matrix)
|
42
|
-
|
43
|
-
```
|
2
与えられるベクトルと欲しい行列の一例を追加しました。
test
CHANGED
File without changes
|
test
CHANGED
@@ -16,7 +16,15 @@
|
|
16
16
|
|
17
17
|
|
18
18
|
|
19
|
-
以下、
|
19
|
+
以下、ベクトルと欲しい行列です
|
20
|
+
|
21
|
+
[[1/√2],と[1/√2, 1/√2]の積が[[1/2, 1/2],
|
22
|
+
|
23
|
+
[1/√2]] [1/2, 1/2]]
|
24
|
+
|
25
|
+
[[1/√2],と[1/√2, -1/√2]の積が[[1/2, -1/2],
|
26
|
+
|
27
|
+
[-1/√2]] [-1/2, 1/2]]
|
20
28
|
|
21
29
|
|
22
30
|
|
1
具体的なコード追加
test
CHANGED
File without changes
|
test
CHANGED
@@ -9,3 +9,27 @@
|
|
9
9
|
よろしくお願いいたします。
|
10
10
|
|
11
11
|
調べてもあんまり見つけられなくて...
|
12
|
+
|
13
|
+
|
14
|
+
|
15
|
+
追記
|
16
|
+
|
17
|
+
|
18
|
+
|
19
|
+
以下、具体的なコードです
|
20
|
+
|
21
|
+
|
22
|
+
|
23
|
+
```ptyhon3
|
24
|
+
|
25
|
+
m = np.array([np.sqrt(0.5), np.sqrt(0.5)])
|
26
|
+
|
27
|
+
n = np.array([np.sqrt(0.5), - np.sqrt(0.5)])
|
28
|
+
|
29
|
+
m_matrix = np.dot(m, m.T)
|
30
|
+
|
31
|
+
n_matrix = np.dot(n, n.T)
|
32
|
+
|
33
|
+
print(m_matrix, n_matrix)
|
34
|
+
|
35
|
+
```
|