blog.monophile.net

Takaaki Yamamoto

東京工業大学において計算機科学と応用数学を学び、 情報科学芸術大学院大学[IAMAS]においてメディア表現を専攻し修了。 digitiminimi Inc. において、インフラエンジニアとして生計をたててている。

各種環境の構築と管理を承ります。

  • 仮想環境: Openstack, GCP, AWS, Azure, ...
  • アプリケーション: WordPress, GitLab, Redmine, ...

List

2元2次式の可視化

code

import matplotlib.pyplot as plt
import numpy as np

def show2d(func,(sx,ex),(sy,ey),num=20):
    A = np.zeros((num,num) , dtype=np.float)
    for i,x in enumerate(np.linspace(sx,ex,num)):
        for j,y in enumerate(np.linspace(sx,ex,num)):
            A[num-i-1][j] = func(x,y)
    plt.imshow(A, interpolation="nearest", extent=[sx,ex,sy,ey])
    plt.show()

if __name__ == "__main__":
    dim = 2
    mat = np.random.random((dim,dim))
    for i in range(dim):
        for j in range(dim-i):
            mat[i][j] = A[j][i]
    quad = lambda x: x.T.dot(mat).dot(x)
    func = lambda x,y : quad(np.array((x,y)))
    show2d(func,(0,1),(0,1),num=100)