質問編集履歴

追記

2020/08/18 23:11

投稿

hataki7

スコア4

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -215,3 +215,217 @@
 ValueError: x and y must have same first dimension, but have shapes (10000,) and (1,)
 ```
+書き換えてみたけど別のエラーが出る。
+```python
+from scipy.special import kv
+import matplotlib.pyplot as plt
+from scipy.integrate import quad
+import numpy as np
+import math
+from math import gamma
+from sympy import *
+import os
+xs=np.logspace(-4, 2, 10000)
+m=9.10938356*1e-31
+c=2.99792458*1e+8
+q=1.6021766208*1e-19
+pa=np.pi/2
+sin_pa=math.sin(pa)
+B1=100*1e-6*1e-4
+a = gamma(1/3)
+f = lambda z: kv(5/3,z)
+F = [quad(f,x,np.inf)[0]*x for x in xs]
+G = [(4*np.pi/np.sqrt(3)/a)*(x/2)**(1/3) for x in xs]
+H = [((np.pi/2)**(1/2))*(x**(1/2))*(np.exp(-x))for x in xs]
+def A(x,F,G,H):
+    if x <= 5.0*1e-3:
+        return G(x,y)
+    elif 5.0*1e-3 < x < 30:
+        return F(x,y)
+    elif 30 <= x:
+        return H(x,y)
+vec_A = np.vectorize(A)
+def L(x):
+    return vec_A(xs,F,G,H)
+def P(v,y):
+    return (math.sqrt(3)*q**3*B1*sin_pa/(m*c**2))*L(v/v_c)
+gmin=1
+gmax=10**5
+N=10/(10**5-1)
+p=-2
+def s(x):
+    b=L*x**(-1/2)
+    return b
+def x1(v):
+    return 4*np.pi*m*c*v/(3*q*B1*sin_pa*gmin**2)
+def x2(v):
+    return 4*np.pi*m*c*v/(3*q*B1*sin_pa*gmax**2)
+def b(v):
+    b1=-1/math.sqrt(v)
+    b2=(math.sqrt(3)*q**3*B1*sin_pa/(2*m*c**2))
+    b3=math.sqrt(3*q*B1*sin_pa/(4*np.pi*m*c))
+    return b1*b2*b3*N
+def Pt(v):
+    a=[quad(s,x1(v),x2(v))[0] for x in xs]
+    return b*a[0]
+v=np.logspace(-4,2,10000)
+plt.plot(v,Pt(v))
+plt.show()
+```
+エラー
+```python
+ValueError                                Traceback (most recent call last)
+<ipython-input-14-4d47f52aab3a> in <module>
+     66 v=np.logspace(-4,2,10000)
+     67
+---> 68 plt.plot(v,Pt(v))
+     69 plt.show()
+     70
+<ipython-input-14-4d47f52aab3a> in Pt(v)
+     61
+     62 def Pt(v):
+---> 63     a=[quad(s,x1(v),x2(v))[0] for x in xs]
+     64     return b*a[0]
+     65
+<ipython-input-14-4d47f52aab3a> in <listcomp>(.0)
+     61
+     62 def Pt(v):
+---> 63     a=[quad(s,x1(v),x2(v))[0] for x in xs]
+     64     return b*a[0]
+     65
+C:\ProgramData\Anaconda3\lib\site-packages\scipy\integrate\quadpack.py in quad(func, a, b, args, full_output, epsabs, epsrel, limit, points, weight, wvar, wopts, maxp1, limlst)
+    336
+    337     # check the limits of integration: \int_a^b, expect a < b
+--> 338     flip, a, b = b < a, min(a, b), max(a, b)
+    339
+    340     if weight is None:
+ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()
+```

文字の修正

2020/08/18 23:11

投稿

hataki7

スコア4

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -1,4 +1,4 @@
-x-Ptグラフをプロットしたいがx and y must have same first dimension, but have shapes (10000,) and (1,)のエラーが出る
+v-Ptグラフをプロットしたいがx and y must have same first dimension, but have shapes (10000,) and (1,)のエラーが出る
@@ -28,6 +28,10 @@
 ys=np.logspace(0, 5, 10000)
+vs=np.logspace(0,5,10000)
 m=9.10938356*1e-31
 c=2.99792458*1e+8
@@ -100,7 +104,7 @@
-def R(y):
+def R(v,y):
     return N*y**-p*P
@@ -108,7 +112,7 @@
 def Pt(v):
-    return [quad(R, gmin, gmax)[0] for y in ys]
+    return [quad(lambda y:R, gmin, gmax)[0] for y in ys]
@@ -124,7 +128,7 @@
 ax.set_xscale('log')
-ax.plot(xs,Pt)
+ax.plot(vs,Pt)
@@ -138,17 +142,17 @@
 ValueError                                Traceback (most recent call last)
-<ipython-input-9-c903738de0cc> in <module>
+<ipython-input-5-33ed3b99f09b> in <module>
-     58 ax.set_yscale('log')
+     60 ax.set_yscale('log')
-     59 ax.set_xscale('log')
+     61 ax.set_xscale('log')
----> 60 ax.plot(xs,Pt)
+---> 62 ax.plot(vs,Pt)
-     61
+     63
-     62 plt.show()
+     64 plt.show()

文字の修正

2020/08/18 18:22

投稿

hataki7

スコア4

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -34,7 +34,7 @@
 q=1.6021766208*1e-19
-pa=np.pi/2
+pa=np.pi/2#ピッチ角#
 sin_pa=math.sin(pa)
@@ -56,29 +56,37 @@
 H = [((np.pi/2)**(1/2))*(x**(1/2))*(np.exp(-x))for x in xs]
-def A(x,F,G,H):
+def A(x,y,F,G,H):
     if x <= 5.0*1e-3:
-        g =G
+        return G(x,y)
     elif 5.0*1e-3 < x < 30:
-        g =F
+        return F(x,y)
     elif 30 <= x:
-        g =H
-    return g
+        return H(x,y)
 vec_A = np.vectorize(A)
+def L(x,y):
-L = vec_A(xs,F,G,H)
+    return vec_A(xs,ys,F,G,H)
+def P(v,y):
+    x=v/v_c
-P=(math.sqrt(3)*q**3*B1*sin_pa/(m*c**2))*L
+    return (math.sqrt(3)*q**3*B1*sin_pa/(m*c**2))*L(v/v_c)
@@ -104,6 +112,8 @@
 fig = plt.figure()
 ax = fig.add_subplot(1,1,1)
@@ -128,17 +138,17 @@
 ValueError                                Traceback (most recent call last)
-<ipython-input-1-8a0453e8c3d3> in <module>
+<ipython-input-9-c903738de0cc> in <module>
-     53 ax.set_yscale('log')
+     58 ax.set_yscale('log')
-     54 ax.set_xscale('log')
+     59 ax.set_xscale('log')
----> 55 ax.plot(xs,Pt)
+---> 60 ax.plot(xs,Pt)
-     56
+     61
-     57 plt.show()
+     62 plt.show()