0

現在
https://www.eidos.ic.i.u-tokyo.ac.jp/~tau/lecture/computational_physics/text/4-matplot-numpy.pdf
こちらのサイトをPythonの自学学習として利用させていただいています。係数行列A、ベクトルfまでは作成できたのですが、この先がどうしてもわかりません。何かアドバイスやヒントを教えていただけると幸いです。

import numpy as np
from sympy import *
from scipy.integrate import odeint
import matplotlib.pyplot as plt


n=10
m=2.0
k=2.0
x0=0.0

X=np.zeros((n,2))
print(X)

StartTime=0.0
EndTime=30.0
dt=0.01

Count=int((EndTime-StartTime)/dt)

A=np.arange(n*n).reshape(n,n)
for i in range(0,n):
    for j in range(0,n):
        if j==i:
            A[i][j]=2
        elif j==i-1 or j==i+1:
            A[i][j]=-1
        else:
            A[i][j]=0
A[9][9]=1
A=A*(-k/m)
print(A)

f=np.zeros((n,1))
f[0][0]=-x0*(-k/m)
print(f)

このようなプログラムを作成し、

[[0. 0.]
 [0. 0.]
 [0. 0.]
 [0. 0.]
 [0. 0.]
 [0. 0.]
 [0. 0.]
 [0. 0.]
 [0. 0.]
 [0. 0.]]
[[-2.  1. -0. -0. -0. -0. -0. -0. -0. -0.]
 [ 1. -2.  1. -0. -0. -0. -0. -0. -0. -0.]
 [-0.  1. -2.  1. -0. -0. -0. -0. -0. -0.]
 [-0. -0.  1. -2.  1. -0. -0. -0. -0. -0.]
 [-0. -0. -0.  1. -2.  1. -0. -0. -0. -0.]
 [-0. -0. -0. -0.  1. -2.  1. -0. -0. -0.]
 [-0. -0. -0. -0. -0.  1. -2.  1. -0. -0.]
 [-0. -0. -0. -0. -0. -0.  1. -2.  1. -0.]
 [-0. -0. -0. -0. -0. -0. -0.  1. -2.  1.]
 [-0. -0. -0. -0. -0. -0. -0. -0.  1. -1.]]
[[0.]
 [0.]
 [0.]
 [0.]
 [0.]
 [0.]
 [0.]
 [0.]
 [0.]
 [0.]]

この結果を得ています。

2
  • Yutoさん、もしかしてアカウントがふたつに分かれてしまっていませんでしょうか? 質問の投稿者であれば自由に内容を編集できます。もし分かれてしまった場合、こちらのヘルプページからアカウントを統合できます :)
    – nekketsuuu
    2019年5月11日 10:47
  • 「この先が分かりません」ではなく、具体的に何をするプログラムが書けないのかを書いて頂けると答えやすいです。
    – nekketsuuu
    2019年5月11日 10:48

0

この質問に回答するには、ログインする必要があります。

のタグが付いた他の質問を参照する。