Python RotMeasTools.gen_RMtau_kerの例

コード例 #1

ファイル: test_deconv.py プロジェクト: nkern/capo

import numpy as np
import RotMeasTools as RMT
from pylab import *

N = 15

fq = np.linspace(0.12, 0.18, N)

W = RMT.gen_RMtau_ker(fq)
Wt = RMT.gen_RMtau_ker(fq, inv=True).T

print W.shape, Wt.shape

KER = np.dot(Wt, W)
TEST = np.dot(W, Wt)

figure()
subplot(311)
imshow(np.abs(KER))
subplot(312)
INV = np.linalg.inv(KER)
imshow(np.abs(INV))
subplot(313)
one = np.dot(INV, KER)
imshow(np.abs(one))
draw()

figure()
TEST2 = np.zeros((N, N))
TEST2[N / 2, N / 2] = 1.
TEST2 = TEST2.flatten()

コード例 #2

ファイル: test_deconv.py プロジェクト: SaulAryehKohn/capo

import numpy as np
import RotMeasTools as RMT
from pylab import *

N = 15

fq = np.linspace(0.12,0.18,N)

W = RMT.gen_RMtau_ker(fq)
Wt= RMT.gen_RMtau_ker(fq,inv=True).T

print W.shape,Wt.shape

KER = np.dot(Wt,W)
TEST = np.dot(W,Wt)

figure()
subplot(311)
imshow(np.abs(KER))
subplot(312)
INV = np.linalg.inv(KER)
imshow(np.abs(INV))
subplot(313)
one = np.dot(INV,KER)
imshow(np.abs(one))
draw()

figure()
TEST2 = np.zeros((N,N))
TEST2[N/2,N/2] = 1.
TEST2 = TEST2.flatten()