Q1)下記のコードをPython3.11.0で実行しますが、表示画面の座標等が
表示出来ません。
正しいコードをお教え頂けますと大変有難いです。
from mpl_toolkits.mplot3d import Axes3D
from matplotlib import cm
import pylab
import numpy as np
import mpmath
mpmath.dps = 5
# Use instead of arg for a continuous phase
#figsize=(4.0, 4.0)
def arg2(x):
return mpmath.sin(mpmath.arg(x))
#f = lambda z: abs(mpmath.loggamma(z))
#f = lambda z: arg2(mpmath.exp(z))
#f = lambda z: abs(mpmath.besselj(3,z))
f = lambda z: arg2(mpmath.cos(z))
#fig = figure(figsize=(6.0, 6.0))
fig = pylab.figure()
ax = Axes3D(fig)
X = np.arange(-5, 5, 0.125)
Y = np.arange(-5, 5, 0.125)
X, Y = np.meshgrid(X, Y)
xn, yn = X.shape
W = X*0
for xk in range(xn):
for yk in range(yn):
try:
z = complex(X[xk,yk],Y[xk,yk])
w = float(f(z))
if w != w:
raise ValueError
W[xk,yk] = w
except (ValueError, TypeError, ZeroDivisionError):
# can handle special values here
pass
#print( xk, xn)
# can comment out one of these
ax.plot_surface(X, Y, W, rstride=1, cstride=1, cmap=cm.jet)
ax.plot_wireframe(X, Y, W, rstride=5, cstride=5)
pylab.show()
以上、宜しくお願いします。
No.1ベストアンサー
- 回答日時:
私の環境では全然動きませんでした。
真っ白な絵が表示されるだけ。私は pyplot派なんで、pyplotでちょい書き直したらが動きました。
因みに最近は pylab は非推奨と聞いてます。
メンテされてるのかな~
matplotlib は変化が早いので、pylab が追い付いていないのかも。
nativeなメソッドを使うことを推奨します。
from matplotlib import cm
import matplotlib.pyplot as plt
import numpy as np
import mpmath
mpmath.dps = 5
def arg2(x):
return mpmath.sin(mpmath.arg(x))
f = lambda z: arg2(mpmath.cos(z))
fig, ax = plt.subplots(subplot_kw=dict(projection="3d"))
X = np.arange(-5, 5, 0.125)
Y = np.arange(-5, 5, 0.125)
X, Y = np.meshgrid(X, Y)
xn, yn = X.shape
W = X*0
for xk in range(xn):
for yk in range(yn):
try:
z = complex(X[xk,yk],Y[xk,yk])
w = float(f(z))
if w != w:
raise ValueError
W[xk,yk] = w
except (ValueError, TypeError, ZeroDivisionError):
pass
ax.plot_surface(X, Y, W, rstride=1, cstride=1, cmap=cm.jet)
ax.plot_wireframe(X, Y, W, rstride=5, cstride=5)
plt.show()
私の実行環境はこれ
python 3.9.6
matplotlib 3.6.2
mpmath 1.2.1
numpy 1.23.5
です。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 工学 Pythonの3次元描画に関する質問です 3 2022/12/07 20:07
- その他(プログラミング・Web制作) python 2 2022/12/23 09:06
- その他(プログラミング・Web制作) pythonのグローバル変数 2 2022/11/25 18:02
- 英語 L-PRF can be obtained by manual or automated metho 1 2022/04/08 09:39
- その他(プログラミング・Web制作) listへのappendが出来ない件 1 2022/12/06 21:44
- YouTube youtubeからの著作権メールについて 2 2023/02/14 16:12
- 英語 できるだけ直訳で英語の翻訳をお願いします。(英語→日本語) 1 2022/10/15 20:59
- Visual Basic(VBA) ファイル全てを .xlsm に変更したところ、プログラムが途中で落ちてしまっています 17 2022/12/07 12:03
- FX・外国為替取引 1.The U.S. FI sells $100 million for pounds at the 1 2022/05/03 00:08
- Ruby vscode 文字化け 1 2022/05/21 19:17
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
照明器具の水銀
-
電気電子材料についてです。 c...
-
オペアンプを使用したアナログ...
-
ガスタービンとガスエンジン、...
-
仮に自動車がAI化されても、私...
-
低圧絶縁抵抗測定について教え...
-
放射線取扱主任者令和5年度 物...
-
よろしくお願いします! 下の図...
-
電検王さんの問題について
-
電気工学の学び方
-
ラップに表裏はあるのでしょうか?
-
配線の色の意味
-
プラスチックの棒 (太さ1cm) を...
-
水と空気でガソリンを作れるっ...
-
九州大学の電気情報工学科より...
-
双方向通話システムの自作(電...
-
オープンコレクタの使い方
-
電気回路 写真の定常状態の回路...
-
一線地絡時のEVTについて 添付...
-
写真のように電池とコイルと鉄...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
三相交流回路の意味がチンプン...
-
燃料電池を三次電池と言わない...
-
なんで
-
低圧絶縁抵抗測定について教え...
-
「Fパラメータ」が出力端「開放」「...
-
よろしくお願いします。 乾電池...
-
電子回路の記号
-
直流回路の電力測定について 電...
-
【工学】プランジャーとシリン...
-
IPCPネゴシエーションについて
-
皆さん、おはようございます♪ ...
-
電動糸紡ぎ機とドラムカーダー...
-
電気工学と電子工学の学科選択...
-
板の表面の平面度が2005.7μm、...
-
ラップに表裏はあるのでしょうか?
-
絶縁抵抗試験と耐電圧試験
-
九州大学工学部gpaについて
-
プリウスにphevやソーラーpanel...
-
5段階スピードエレキモーターに...
-
しかしながら使用可能なGPIOは...
おすすめ情報