回答編集履歴
2
説明の追記
test
CHANGED
|
@@ -156,6 +156,8 @@
|
|
|
156
156
|
|
|
157
157
|
|
|
158
158
|
|
|
159
|
+
# 追加
|
|
160
|
+
|
|
159
161
|
A=np.zeros((f_points),dtype=complex)
|
|
160
162
|
|
|
161
163
|
B=np.zeros((f_points),dtype=complex)
|
1
コード変更
test
CHANGED
|
@@ -58,7 +58,7 @@
|
|
|
58
58
|
|
|
59
59
|
def FF(i):
|
|
60
60
|
|
|
61
|
-
return Q
|
|
61
|
+
return matQ[i]@matW[i]@matE[i]@matW[i]@matQ[i]
|
|
62
62
|
|
|
63
63
|
|
|
64
64
|
|
|
@@ -70,7 +70,7 @@
|
|
|
70
70
|
|
|
71
71
|
#return np.linalg.matrix_power(FF, n)
|
|
72
72
|
|
|
73
|
-
return np.linalg.matrix_power(FF
|
|
73
|
+
return np.linalg.matrix_power(matFF[i], n)
|
|
74
74
|
|
|
75
75
|
|
|
76
76
|
|
|
@@ -156,8 +156,6 @@
|
|
|
156
156
|
|
|
157
157
|
|
|
158
158
|
|
|
159
|
-
# 追加
|
|
160
|
-
|
|
161
159
|
A=np.zeros((f_points),dtype=complex)
|
|
162
160
|
|
|
163
161
|
B=np.zeros((f_points),dtype=complex)
|
|
@@ -206,13 +204,13 @@
|
|
|
206
204
|
|
|
207
205
|
#D[i]=F[i[1,1]]
|
|
208
206
|
|
|
209
|
-
A[i]=F
|
|
207
|
+
A[i]=matF[i][0,0]
|
|
210
|
-
|
|
208
|
+
|
|
211
|
-
B[i]=F
|
|
209
|
+
B[i]=matF[i][0,1]
|
|
212
|
-
|
|
210
|
+
|
|
213
|
-
C[i]=F
|
|
211
|
+
C[i]=matF[i][1,0]
|
|
214
|
-
|
|
212
|
+
|
|
215
|
-
D[i]=F
|
|
213
|
+
D[i]=matF[i][1,1]
|
|
216
214
|
|
|
217
215
|
s11[i]=(A[i]+B[i]/Z0-(Z0*C[i]+D[i]))/(A[i]+B[i]/Z0+Z0*C[i]+D[i])
|
|
218
216
|
|