偶对称离散余弦变换 EDCT
偶对称离散余弦变换 EDCT
元禛慎独 发表于9个月前
偶对称离散余弦变换 EDCT
  • 发表于 9个月前
  • 阅读 6
  • 收藏 0
  • 点赞 0
  • 评论 0

腾讯云 技术升级10大核心产品年终让利>>>   

from __future__ import division
import numpy as np

def fourier_matrix_ec(J):
    def ec_fourier(k,l):
        return np.exp(-1.j*np.pi/J*(k+0.5)*l)

    return np.array([[ec_fourier(k,l) for k in np.arange(-J,J)] for l in np.arange(2*J)])/2/J


U=fourier_matrix_ec(4)

#print np.round(U)

A1=np.array([[1,2,0,1],\
             [1,0,0,0],\
             [0,0,2,2],\
             [1,2,2,0]])

B=np.zeros_like(A1)
for i in range(len(B)):
    B[i][len(B)-1-i]=1

print B

A2=np.dot(A1,B)
A3=np.dot(B,A2)
A4=np.dot(B,A1)
A=np.r_[np.c_[A3,A4],np.c_[A2,A1]]

print A
g=np.dot(U,np.dot(A,U.T))
print np.round(g,4)
print np.round(g.real,4)

print np.round(g.imag,4)

 

共有 人打赏支持
粉丝 3
博文 107
码字总数 54993
×
元禛慎独
如果觉得我的文章对您有用,请随意打赏。您的支持将鼓励我继续创作!
* 金额(元)
¥1 ¥5 ¥10 ¥20 其他金额
打赏人
留言
* 支付类型
微信扫码支付
打赏金额:
已支付成功
打赏金额: