私は機能を持っています:
def g(R, r):
return (np.sqrt(2.0 * (R + r) / (r * R)) - (1 + np.sqrt(R)) / np.sqrt(R) -
np.sqrt(2.0 / (r * (1 + r))) * (1 - r) -
(1.0 / np.sqrt(R) - np.sqrt(2.0) * (1 - R) / np.sqrt(R * (1 + R))
- 1))
関数は、delta v_B = delta v_H
場所delta v_B
を設定することによって定義されます
np.sqrt(2.0 * (R + r) / (r * R)) - (1 + np.sqrt(R)) / np.sqrt(R) -
np.sqrt(2.0 / (r * (1 + r))) * (1 - r)
そしてdelta v_H
_
1.0 / np.sqrt(R) - np.sqrt(2.0) * (1 - R) / np.sqrt(R * (1 + R)) - 1
ということで、 と書きg
ましdelta v_b - delta v_H
た。
これが私の関数であり、以下で使用しているコードです。
import pylab
import numby as np
def g(R, r):
return (np.sqrt(2.0 * (R + r) / (r * R)) - (1 + np.sqrt(R)) / np.sqrt(R) -
np.sqrt(2.0 / (r * (1 + r))) * (1 - r) -
(1.0 / np.sqrt(R) - np.sqrt(2.0) * (1 - R) / np.sqrt(R * (1 + R))
- 1))
r = np.linspace(11.9, 16, 500000)
R = np.linspace(1, 20, 500000)
fig2 = pylab.figure()
ax2 = fig2.add_subplot(111)
ax2.plot(R, g(R, r), 'r')
pylab.xlabel('$R_1 = \\frac{r_C}{r_A}$')
pylab.ylabel('$R_2 = \\frac{r_B}{r_A}$')
pylab.xlim((0, 25))
pylab.ylim((0, 100))
pylab.show()
関数は、約で無限大に漸近し、約で11.94
線と交差する必要がありますy = x
15.58
どうすればそのようなプロットを作成できますか? 私はこれを行う方法に慣れておらず、そのような関数をプロットする方法も知りません。
私の定義はg
asには不適切g(R, r)
ですか? そうでない場合、どのように定義する必要がありますか?