文档章节

偶对称离散余弦变换 EDCT

元禛慎独
 元禛慎独
发布于 2017/03/26 22:24
字数 123
阅读 8
收藏 0
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
博文 209
码字总数 60366
作品 0
朝阳
程序员
私信 提问
【工具使用系列】关于 MATLAB 图像变换,你需要知道的事

如何进行图像变换(频域处理) 正交变换 傅里叶变换 二维连续傅里叶变换 二维离散傅里叶变换(DFT) 快速傅里叶变换(FFT) 离散余弦变换 一维离散余弦变换 二维离散余弦变换 离散正弦变换 ...

AllenMoore
01/26
8
0
离散余弦变换(DCT)

DCT变换、DCT反变换、分块DCT变换 一、引言 DCT变换的全称是离散余弦变换(Discrete Cosine Transform),主要用于将数据或图像的压缩,能够将空域的信号转换到频域上,具有良好的去相关性的性...

li_wen01
2017/06/05
0
0
离散余弦变换(DCT)的来龙去脉

1. 图像的二维离散变换   与一维的有限长离散非周期信号存在傅里叶变换(DFT)一样,图像作为一个二维离散信号同样存在着二维离散变换(注意这里是介绍一个通用的概念,二维离散变换,包括...

dugudaibo
2017/11/01
0
0
【DCT】OPENCV python 离散余弦变换问题

因为项目需要,目前在研究《基于离散余弦变换和区域生长的白粉虱图像分割算法》,使用的OPENCV python。 论文地址:%e5%9f%ba%e4%ba%8e%e7%a6%bb%e6%95%a3%e4%bd%99%e5%bc%a6%e5%8f%98%e6%8d...

majisong
2016/06/02
1K
0
用Python学《微积分B》(Fourier级数)

  除了“幂级数”,另一类非常重要的函数项级数级数就是“三角(函数)级数”(trigonometric series)。《微积分B》第8-6节主要介绍:什么是三角级数?三角函数系的正交性?三角级数与Fou...

sagittarius_warrior
2017/09/30
0
0

没有更多内容

加载失败,请刷新页面

加载更多

初识flask

文档 0.10.1版本 http://www.pythondoc.com/flask/index.html 1.0.2版本 https://dormousehole.readthedocs.io/en/latest/ 安装flask $ pip3 install flaskCollecting flask Downloading......

yimingkeji
昨天
1
0
Akka系统《sixteen》译

Actor是一个封装状态(state)和行为(behavior)的对象,它们只通过交换消息通信(放入收件人邮箱的邮件)。从某种意义上说,Actor是最严格的面向对象编程形式,但它更适合将他们视为人:在与Act...

woshixin
昨天
1
0
技术工坊|如何开发一款以太坊钱包(深圳)

【好消息!】HiBlock区块链技术工坊已经成功举办了26期,其中北京1期,西安1期,成都2期,上海22期。经常有社区的小伙伴问定期举办技术工坊的除了上海以外,其他城市有没有?现在区块链技术工...

HiBlock
昨天
1
0
Redis 梳理笔记

安装 安装gccyum install gcc-c++下载传输到服务器上解压tar -xzvf *.tar.gzcd redis-3.2.9编译make安装 make PREFIX=/usr/local/redis install将配置文件拷贝出来cp redis...

晨猫
昨天
0
0
聊聊storm TridentWindowManager的pendingTriggers

序 本文主要研究一下storm TridentWindowManager的pendingTriggers TridentBoltExecutor.finishBatch storm-core-1.2.2-sources.jar!/org/apache/storm/trident/topology/TridentBoltExecut......

go4it
昨天
1
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部