質問編集履歴
1
不明点，エラーコードを追記しました
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -1,9 +1,11 @@
|
|
|
1
|
-
###
|
|
1
|
+
### 質問内容
|
|
2
2
|
|
|
3
3
|
下記作成中のコードにて，極座標(r-θ平面)上にベクトル図を作図しようと考えています．
|
|
4
4
|
r,θを変数とし，各座標の値を代入して，その場の速度を算出してベクトル図を描ければと考えましたが上手くいきませんでした．
|
|
5
|
-
改善点等，アドバイスをいただけると助かります．
|
|
6
5
|
|
|
6
|
+
下記エラーコードが出たため，その対策も調べ確認してみましたが，上手くいかずご質問させていただきました．
|
|
7
|
+
アドバイスなど頂けるとありがたいです．
|
|
8
|
+
|
|
7
9
|
### 作成コード
|
|
8
10
|
|
|
9
11
|
```
|
|
@@ -111,4 +113,37 @@
|
|
|
111
113
|
ax.quiver(theta,radius, VelocityX/TotalVelocity, VelocityY/TotalVelocity)
|
|
112
114
|
|
|
113
115
|
plt.show()
|
|
114
|
-
```
|
|
116
|
+
```
|
|
117
|
+
|
|
118
|
+
### エラーコード
|
|
119
|
+
|
|
120
|
+
```
|
|
121
|
+
ValueError Traceback (most recent call last)
|
|
122
|
+
<ipython-input-1-ee1d8f65e98c> in <module>
|
|
123
|
+
86
|
|
124
|
+
87 A_M = M(lambda r: r2 / (r * cmath.sin(beta2)), r2, r1)
|
|
125
|
+
---> 88 A_F_rc = F_rc(lambda r: cmath.exp((-1j* (omega * cmath.pi * (r**2 - r2**2) / Q - n * ganma * cmath.log(r / r2) / Q ))) * (r / radius)**(n + 1), radius, r3)
|
|
126
|
+
89 A_G_rc = G_rc(lambda r: cmath.exp((-1j* (omega * cmath.pi * (r**2 - r2**2) / Q - n * ganma * cmath.log(r / r2) / Q ))) * (radius / r)**(n - 1), r3, radius)
|
|
127
|
+
90 A_F_r2 = F_r2(lambda r: cmath.exp((-1j* (omega * cmath.pi * (r**2 - r2**2) / Q - n * ganma * cmath.log(r / r2) / Q ))) * (r / r2)**(n + 1), r3, r2)
|
|
128
|
+
|
|
129
|
+
<ipython-input-1-ee1d8f65e98c> in F_rc(func, a, b, **kwargs)
|
|
130
|
+
22 def imag_func(r):
|
|
131
|
+
23 return scipy.imag(func(r))
|
|
132
|
+
---> 24 real_integral = quad(real_func, a, b, **kwargs)
|
|
133
|
+
25 imag_integral = quad(imag_func, a, b, **kwargs)
|
|
134
|
+
26 return (real_integral[0] + 1j*imag_integral[0], real_integral[1:], imag_integral[1:])
|
|
135
|
+
|
|
136
|
+
~/opt/anaconda3/lib/python3.7/site-packages/scipy/integrate/quadpack.py in quad(func, a, b, args, full_output, epsabs, epsrel, limit, points, weight, wvar, wopts, maxp1, limlst)
|
|
137
|
+
336
|
|
138
|
+
337 # check the limits of integration: \int_a^b, expect a < b
|
|
139
|
+
--> 338 flip, a, b = b < a, min(a, b), max(a, b)
|
|
140
|
+
339
|
|
141
|
+
340 if weight is None:
|
|
142
|
+
|
|
143
|
+
ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()
|
|
144
|
+
|
|
145
|
+
```
|
|
146
|
+
|
|
147
|
+
### 出力結果
|
|
148
|
+
|
|
149
|
+
|