pyhton matplotlib を用いて円を描きたいのですが

Question

コードはエラーなく通るのですがうまく円がかけずに写真のようになってしまいます．
コード上のミスなどがあれば教えていただきたいです．よろしくお願いします．

import matplotlib.pyplot as plt
import matplotlib.ticker as tick

x = np.arange(-190,190,10)
f1 = np.sqrt(36100-x**2)
f2 = -np.sqrt(36100-x**2)

plt.rc("font",family="serif")
plt.figure()
plt.xlim([-200,200])
plt.ylim([-200,200])
plt.plot(f1,color="black")
plt.plot(f2,color="black")

plt.gca().xaxis.set_minor_locator(tick.MultipleLocator(1))
plt.gca().yaxis.set_minor_locator(tick.MultipleLocator(1))
plt.grid(which="minor")

plt.show()
コード

Comments

[nmoa]

https://teratail.com/help/question-tips#questionTips3-5
上記リンクのアドバイスに書かれているように、ソースをシンタックスハイライトしてください。
「**」の部分がmarkdownとして解釈されてコードが正しく表示されていません。

[naoto8888]

申し訳ありません．ご指摘ありがとうございます．

Answer 1

質問者様と書き方は違いますが、例えばこういった書き方はダメなのでしょうか？
極座標で書きました。

import numpy as np
import matplotlib.pyplot as plt
import math as m


theta = np.arange(0,2*m.pi+1,0.1)
xlist=[];ylist=[]
r = 1 #radius
for i in theta:
    xlist.append(r*m.cos(i))
    ylist.append(r*m.sin(i))
plt.figure(figsize=(3,3))
plt.plot(xlist,ylist)
plt.grid(True)

r =1 で単位ですので、半径１００にしたいなら、r =100です

Comments

[nmoa]

listをループで回すのはpythonでは効率が悪いです。numpy.cos(),numpy.sin()を使えばnumpyのndarrayに対して直接演算ができ、処理が速くなります。一部書き換えるなら以下のような感じでしょうか。

r = 1
theta = np.arange(0,2*m.pi+1,0.1)
plt.plot(r*np.cos(theta), r*np.sin(theta))

[Fallout_18]

アドバイスありがとうございます

[naoto8888]

Fallout_18さん，nmoaさん，回答，それに関する追記ありがとうございます．こちらの方法で行いたいと思います．ありがとうございました

Answer 2

単純に円を描きたいだけであれば、matplotlib.patches.Circle を使ったほうが簡単ですが、いかがでしょうか。

matplotlib - 円、長方形などの描画方法と設定項目

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.patches as mpatches

fig, ax = plt.subplots(figsize=(8, 8), facecolor='w')
ax.set_aspect('equal')

circle = mpatches.Circle(xy=(100, 100), radius=90, fill=False, color='b')
ax.add_patch(circle)
ax.set_xlim(0, 200), ax.set_ylim(0, 200)
ax.set_xticks(np.arange(0, 200), minor=True)
ax.set_yticks(np.arange(0, 200), minor=True)
ax.grid(which='minor', alpha=0.5)

plt.show()

Comments

[naoto8888]

背景にある1mm×1mmのグリッドを残したまま円を描きたいです．もしその方法を知っていれば回答よろしくお願いします．

[tiitoi]

1目盛りごとのグリッドということでしょうか？
円を描画するのと、グリッドを作成するのは別々のタスクであるため、単にその処理を書いてあげればよいです。
ちなみに単位は mm ではないと思いますよ。

[naoto8888]

なるほど．ありがとうございます．そのような方向で行いたいと思います！