python ループ　zip プロット

#うまくループができず、zipの中に落とし込めません。どうやらfor文でエラーが起きているようなので構文が違うのか、式が違うのかよくわからないので、教えてください。

from sympy.geometry import Point, Circle, Triangle, Segment, Line
import numpy as np
import matplotlib.pyplot as plt

fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
ax.set_aspect('equal')
ax.grid()


t=0.63
z=-0.133
L=0.0768

for n in range(15);
     if 15<n:
  p_ｎ= (x_n, y_n) #(n=0,1,2...)
  x_n=t*(n//3)+L*(n%3)  
  if n%3=0,2
  y_n=0
  else:
  y_n=-0.133
  if i == 15:
　　break

ax.plot(*zip(P_n)) #p1, p2, p3を結んでできる折れ線
g=ax.plot
plt.xlim([0,7])
plt.ylim([-0.2,0.1])
plt.show(g)

－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－
ｐの中は下記を想定しています。
p1 = (0, 0)
p2 = (0.0768, -0.133)
p3 = (0.1536, 0)
p4= (0.63, 0)
p5= (0.7068,-0.133)
p6= (0.7836,0)
p7= (1.26,0)
p8= (1.3368,-0.133)
p9= (1.4136,0)
p10=(1.89,0)
p11=(1.9668,-0.133)
p12=(2.0436,0)
p13=(2.52,0)
p14=(2.5968,-0.133)
p15=(2.6736,0)

tiitoi

2019/10/15 05:51 編集

コードは markdown 記法で記載してください。 https://qiita.com/shizuma/items/8616bbe3ebe8ab0b6ca1#2%E3%82%B3%E3%83%BC%E3%83%89%E3%81%AE%E6%8C%BF%E5%85%A5 あと for(int n = 0; n < 15; n++) などコードの中に Python ではなく C 言語？のものが混じっているようですが、いかがでしょうか。

11390

2019/10/15 06:12

参考になります。すいませんC言語は前使っていたので混じっていたようです。 pythonではどのようにすればよいでしょうか？

tiitoi

2019/10/15 06:20

n 番目の点 pn がどのような数式で表されるのか記載していただけますか？

11390

2019/10/15 06:29

p_n=(x_n,y_n) 中身は x_n=0.63*(n//3)+0.0768*(n%3) if n%3=0,2 y_n=0 else: y_n=-0.133

行動規範の内容に同意します

回答1件

ベストアンサー

p_n=(x_n,y_n)
中身は
x_n=0.63*(n//3)+0.0768*(n%3)
if n%3=0,2
y_n=0
else:
y_n=-0.133

各点を表す数式がわかっているのであれば、それをそのままコードにすればよいです。
numpy を使うと、for 文で1項ずつ計算するのではなく、n 項すべてを一度に計算できます。

python
1import matplotlib.pyplot as plt
2import numpy as np
3
4n = np.arange(0, 21)  # 第0項 ～ 第20項
5
6# 各項の x, y の値
7xs = 0.63 * (n // 3) + 0.0768 * (n % 3)
8ys = np.where(n % 3 == 1, -0.133, 0)  # n % 3 == 1 のときは -0.133 それ以外は 0
9
10# 描画する。
11for i, (x, y) in enumerate(zip(xs, ys)):
12    print(f"p[{i}] = ({x:.4f}, {y:.4f})")
13
14fig, ax = plt.subplots()
15ax.plot(xs, ys)
16
17plt.show()