こんにちは。
仕事の自動化にやりがいと達成感を感じるガッくんです。

この記事の目次

背景・目的

積分をやったんで、微分も記事にまとめます。

２波形について試してみました。

動作環境

・Windows 10
・winpython 64bit 3.4.4

プログラム

ソースコード


###############################################################################
# 離散データの微分を計算するプログラム
###############################################################################
#==============================================================================
# ライブラリインポート
#==============================================================================
import numpy as np
import matplotlib.pyplot as plt
 
#==============================================================================
# テスト用離散データの作成
#==============================================================================
flag = 1
 
if flag == 0:
    x = np.linspace(0, 2.2, num=2**5+1)
    y = 2.1*x**4 - 3.3*x**3 - 2*x + 1
elif flag == 1:
    x = np.arange(-7.0, 8.0, 1)
    y = x**2
 
#==============================================================================
# 微分計算
#==============================================================================
dy = np.gradient(y)
 
#==============================================================================
# グラフ
#==============================================================================
# 微分対象波形
plt.plot(x, y, 'r-', label='wave')
plt.plot(x, y, 'r.')
 
# 微分波形
plt.plot(x, dy, 'b-', label='differential')
plt.plot(x, dy, 'b.')
 
# グラフのタイトル・目盛設定
plt.title("Differential")
 
# 凡例の位置設定
plt.legend(loc='upper left')
 
# 罫線
plt.grid(which='major', color='black', linestyle='-')
plt.grid(which='minor', color='black', linestyle='-')
 
# 軸ラベル
plt.xlabel('x')
plt.ylabel('y')
 
# y 軸目盛
if flag == 1:
    plt.ylim([-20, 70])
 
file_name = 'differential.jpg'
plt.savefig(file_name)
plt.show()