Python NumPy 三角波の生成の仕方が知りたいです

Question

前提

ここに質問の内容を詳しく書いてください。

この式を三角波で生成したいです。

実現したいこと

xn=の後が知りたい

発生している問題・エラーメッセージ

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


t = np.linspace(0,3,400)#0～3秒間の時間を400に分割して配列tを生成



x1 = 8/(1**2*np.pi)*np.sin(2*np.pi*1*t)
x3 = 8/(1**2*np.pi)*np.sin(2*np.pi*1*t)- 8/(3**2*np.pi)*np.sin(2*np.pi*1*t)
x5 = 8/(1**2*np.pi)*np.sin(2*np.pi*1*t)- 8/(3**2*np.pi)*np.sin(2*np.pi*3*t)+ 8/(5**2*np.pi)*np.sin(2*np.pi*5*t)
xn=0 #n個のsin波を重ね合わせる前にリセット
n=100 #重ね合わせるn個の回数を指定
c=1
for i in range(1,n,2): #1～nまで(2つ飛ばし)のsin波の重ね合わせを行う
  xn=　　　　　#ここの部分の一行が知りたい。



plt.title("TRIANGLE WAVE")
plt.plot(t, x1, color=(1.0,0.0,0.0), linewidth=1.0, label="n=1")
plt.plot(t, x3, color=(0.0,1.0,0.0), linewidth=1.0, label="n=3")
plt.plot(t, x5, color=(0.0,0.0,1.0), linewidth=1.0, label="n=5")
plt.plot(t, xn, color=(1.0,0.0,1.0), linewidth=1.0, label="n=100")
plt.legend()
plt.xlabel('t')
plt.ylabel('x(t)')
plt.show()

Answer 1

for内のxnは+=で加算していく必要があります。
加算項は式が書いてある画像の2番めのものを素直に実装すれば大丈夫です。
また、xnを0ではなく0 * tで初期化しておく必要があります。

以下実装例です。

import numpy as np
import matplotlib.pyplot as plt


t = np.linspace(0, 3, 400)  # 0～3秒間の時間を400に分割して配列tを生成


x1 = 8 / (1 ** 2 * np.pi) * np.sin(2 * np.pi * 1 * t)
x3 = 8 / (1 ** 2 * np.pi) * np.sin(2 * np.pi * 1 * t) - 8 / (3 ** 2 * np.pi) * np.sin(
    2 * np.pi * 1 * t
)
x5 = (
    8 / (1 ** 2 * np.pi) * np.sin(2 * np.pi * 1 * t)
    - 8 / (3 ** 2 * np.pi) * np.sin(2 * np.pi * 3 * t)
    + 8 / (5 ** 2 * np.pi) * np.sin(2 * np.pi * 5 * t)
)
xn = 0 * t  # n個のsin波を重ね合わせる前にリセット
n = 100  # 重ね合わせるn個の回数を指定
for i in range(1, n, 2):  # 1～nまで(2つ飛ばし)のsin波の重ね合わせを行う
    xn += 8 / (i ** 2 * np.pi) * np.sin(i * np.pi / 2) * np.sin(2 * np.pi * i * t)


plt.title("TRIANGLE WAVE")
plt.plot(t, x1, color=(1.0, 0.0, 0.0), linewidth=1.0, label="n=1")
plt.plot(t, x3, color=(0.0, 1.0, 0.0), linewidth=1.0, label="n=3")
plt.plot(t, x5, color=(0.0, 0.0, 1.0), linewidth=1.0, label="n=5")
plt.plot(t, xn, color=(1.0, 0.0, 1.0), linewidth=1.0, label="n=100")
plt.legend()
plt.xlabel("t")
plt.ylabel("x(t)")
plt.show()

Comments

[huta_88]

+=で加算させてあげればいいのですね。
ありがとうございました。