![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?8acaa2e)
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
です。
![「Pythonの3Dグラフ表示に関する質問」の回答画像1](http://oshiete.xgoo.jp/_/bucket/oshietegoo/images/media/8/991608_63914bd958edf/M.png)
お探しの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ランキング
-
パリティチェッカ47180の回路図...
-
クレーン
-
低圧絶縁抵抗測定について教え...
-
【工学博士に質問です】SNSで「...
-
電気磁気の問題です。 半径 a[m...
-
先端がマグネットになってるド...
-
計測アンプについて
-
製図です、ハッチング狭すぎで...
-
「最もローラン展開公式は、コ...
-
乾電池3本で動くファンがあり...
-
制御工学の問題です!G(S)=8...
-
オペアンプを使用した電圧ー電...
-
製図の断面図についてです。図1...
-
ヒトの可聴周波数は20Hz~20kHz...
-
一般的な充電器で開放型バッテ...
-
「やさしく学べる材料力学」の...
-
電気回路でこの問題でIAとI1が...
-
アブソモータのエンコーダリセ...
-
この突起はなんと呼ばれ何の為...
-
はしご回路 この赤で囲んだ抵抗...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
VBA ListViewの選択Itemを削除する
-
エクセルで検索結果をユーザー...
-
エクセルマクロ(VBA)でテ...
-
xcode 画像をランダムに並べ順...
-
This is the man who they say ...
-
レインボーブリッジ遊歩道近く...
-
平均滞留時間
-
座椅子
-
大阪であんかけスパゲティを食...
-
PCの周辺機器でマニアックなお...
-
この花の名前を教えて下さい
-
アルマダとはなんのこと?
-
「★」を境に文字を分割し隣のセ...
-
sumif関数で列を数字で指定でき...
-
読み方教えてください
-
アルミ電解コンデンサで、20...
-
この問題の解き方を教えてください
-
ひにちって漢字で書くと
-
blowki11er4.6 というソフト
-
横浜市内の格安の理容室を教え...
おすすめ情報