Python魔法函数与两比特量子系统模拟
技术背景
本文主要涵盖两个领域的知识点:python的魔法函数和量子计算模拟,我们可以通过一个实际的案例来先审视一下这两个需求是如何被结合起来的。
量子计算模拟背景
ProjectQ
是一个非常优雅的开源量子计算编程框架,其原作者是来自与瑞士联邦理工的博士Damian和Thomas。该量子计算编程框架是一个从量子计算应用->量子线路编译->哈密顿量模拟->量子计算模拟->量子硬件API对接都有相应实现的、非常全面的量子计算编程框架。其开源地址为:https://github.com/ProjectQ-Framework/ProjectQ,支持使用pip进行安装:python3 -m pip install projectq --upgrade
。
下面来看一个例子,关于如何使用projectq
进行量子计算的模拟:
[dechin@dechin-manjaro simulator]$ ipython
Python 3.8.5 (default, Sep 4 2020, 07:30:14)
Type 'copyright', 'credits' or 'license' for more information
IPython 7.19.0 -- An enhanced Interactive Python. Type '?' for help.In [1]: from projectq import MainEngineIn [2]: from projectq.ops import XIn [3]: eng = MainEngine()In [4]: qubits = eng.allocate_qureg(2)In [5]: X | qubits[0]In [6]: from projectq.ops import CXIn [7]: CX | (qubits[0], qubits[1])In [8]: eng.flush()In [9]: print (eng.backend.cheat()[1])
[0j, 0j, 0j, (1+0j)]
在这个案例中,我们一共分配了2个量子比特,这2个比特的初始状态都是\(\left|0\right>\)态,对应于projectq
输出的amplitude矢量应为[1, 0, 0, 0]
。这个矢量中的4个元素,分别对应00,01,10,11
这四个量子态可能出现的概率幅,如果需要计算某一个态被测量所出现的概率的话,需要对其进行取模平方操作:
注意概率幅是一个复数(Complex Number),因此需要取厄米共轭之后再进行点乘操作。
那么回到上述projectq的使用案例,这个案例在分配了两个比特之后,对其中的第一个比特执行了泡利矩阵\(\sigma^X\)操作,然后又执行了一个纠缠门操作CX
。这里CX(i,j)
量子门操作对应的操作为:如果量子比特i
处于\(\left|0\right>\)态,不进行任何的操作;但是如果量子比特i
出于\(\left|1\right>\)态,则对量子比特j
执行取反操作
,也就是说,如果原来j
是\(\left|0\right>\)就会变成\(\left|1\right>\),如果原来j
是\(\left|1\right>\)就会变成\(\left|0\right>\)。这就是量子纠缠在量子计算中的作用,而多比特的门操作在实际的硬件体系中的高质量实现,目前依旧是一大难题。而量子叠加特性就体现在,一个量子比特可能处于\(\left|0\right>\)态,也可能处于\(\left|1\right>\)态,还有可能处在\(\left|0\right>\)和\(\left|1\right>\)的中间态,这种中间态会以上述提到的概率幅的形式来对\(\left|0\right>\)和\(\left|1\right>\)进行划分:
这些的概率幅就可以用一个矢量的形式组织起来:
最终这个矢量的元素个数会随着比特数的增加而指数增长
,当比特数增长到41
时,所需要存储的内存空间需要32TB
以上!要注意的是,因为计算过程中需要将所有的概率幅加载到内存中,所以这里区别于硬盘存储空间,单指内存就需要到32TB
的大小!因此,使用经典计算机去模拟量子计算,其实是一种非常消耗资源的手段。当然,量子计算模拟器依然有其研究的价值,在现阶段量子芯片规模和质量无法提升的状态下,模拟器就起到了重要的作用。
Python的魔法函数实现
如果读者需要了解详细全面Python的魔法函数的实现方案,可以从本文的参考链接中获取两篇不错的文章。这里我们仅针对上述projectq
的代码用例中所可能使用到的部分功能:__or__
和__str__
,并且可以针对其进行一个简单的复现。
Python的魔法函数可用于定义一个类(class)的特殊运算算符,如:类的加减乘除等,在引入魔法函数之后,就不需要单独对类中的元素进行操作,而可以用魔法函数对操作进行封装。最后的效果,我们可以直接在代码中使用操作符对不同的类进行操作,比如可以自定义class1 + class2
这样的二元操作算符。在本章节我们不详细展开介绍,可以参考下述的具体使用示例或者参考链接中的博文。
量子态定义及实现
根据第一个章节中对量子态矢量的介绍,这里我们可以实现一个简单的量子态的类,我们可以仅考虑两个量子比特的简单系统:
# QubitPair.py
import numpy as npclass QubitPair:def __init__(self):self.state = np.array([1, 0, 0, 0], dtype=complex)def __str__(self):return str(self.state)
这个量子态的类的定义非常简单,就是一个4×1
的矩阵。需要补充说明的是,这里我们定义了一个__str__(self)
的魔法函数,该函数主要用于打印类的字符串表示,如我们这里直接将量子态矢量转化成str
格式之后进行输出。那么我们如果去print
一个自定义的QubitPair
类的话,就会展示当前类所对应的概率幅的字符串表示。
量子门操作定义及实现
关于量子门操作,我们可以将其视作作用在量子态矢量上的矩阵,这里我们可以先展示定义好的门操作的Python类再对其进行展开说明:
# Operator.py
import numpy as npclass QubitOperator:"""Pauli rotations and entanglement on qubit-pair"""def __init__(self, operation=None, theta=0, index=0):self.index = indexself.name = operationpaulix = np.array([[0, 1], [1, 0]], dtype=complex)pauliy = np.array([[0, -1j], [1j, 0]], dtype=complex)pauliz = np.array([[1, 0], [0, -1]], dtype=complex)cnot = np.array([[1, 0, 0, 0],[0, 1, 0, 0],[0, 0, 0, 1],[0, 0, 1, 0]])if operation == 'X' or operation == 'Rx':self.operation = np.cos(theta/2)*np.identity(2)-1j*np.sin(theta/2)*paulixelif operation == 'Y' or operation == 'Ry':self.operation = np.cos(theta/2)*np.identity(2)-1j*np.sin(theta/2)*pauliyelif operation == 'Z' or operation == 'Rz':self.operation = np.cos(theta/2)*np.identity(2)-1j*np.sin(theta/2)*paulizelif operation == 'CX' or operation == 'CNOT':self.operation = cnotdef __or__(self, qubitpair):if self.name == 'CX' or self.name == 'CNOT':qubitpair.state = np.dot(self.operation, qubitpair.state)return Noneelif self.index == 0:operation = np.kron(self.operation, np.identity(2))else:operation = np.kron(np.identity(2), self.operation)qubitpair.state = np.dot(operation, qubitpair.state)
单位矩阵与泡利矩阵的定义
这些是基本的泡利矩阵,这三个两能级体系的泡利矩阵具有非常好的物理性质,如都是酉矩阵且存在特殊的对易关系等:
矩阵指数与旋转门操作
矩阵的指数计算一般采用泰勒级数展开的方法来进行定义:
这里如果我们代入上述介绍的泡利矩阵就会得到这样的结果:
CX门操作的定义
在上述提到的所有的量子门操作中,CX
是唯一的一个两比特量子门操作,也就是同时作用在两个量子比特上面,其矩阵形式的定义如下所示:
使用魔法函数__or__
来实现量子门操作运算
我们首先简单谈一下为什么要用__or__
这个魔法函数而不是其他的二元运算符来实现,这点跟开源库ProjectQ
是同步的,理由是我们在量子力学中的运算,一般写成如下的形式:
将量子态写成狄拉克符号的形式,中文称为"左矢"和"右矢",英文称之为"bra"和"ket"。因此竖线形式的定义,在形式上会更加契合量子力学的思维,当然,就算是换成其他的符号也是无可厚非的。
功能测试验证
在定义了量子态
的类和量子门操作的类
之后,我们可以写如下所示的一个测试脚本来测试程序的执行效果:
# TestQubits.py
from QubitPair import QubitPair
from Operator import QubitOperatorif __name__ == '__main__':qubits = QubitPair()print ('The initial state is: {}'.format(qubits))QubitOperator('X', 3.1415926, 0) | qubitsprint ('Applying X on the 0th qubit...')print ('The new state is: {}'.format(qubits))QubitOperator('CX') | qubitsprint ('Applying entanglement on qubits...')print ('The new state is: {}'.format(qubits))QubitOperator('X', 3.1415926, 0) | qubitsprint ('Applying X on the 0th qubit...')print ('The new state is: {}'.format(qubits))QubitOperator('CX') | qubitsprint ('Applying entanglement on qubits...')print ('The new state is: {}'.format(qubits))
这个程序的测试逻辑为:先定义一个两比特的量子系统,然后对第一个比特执行X
门操作,使得其从\(\left|0\right>\)态变成\(\left|1\right>\)态,再对这两个比特执行纠缠门CX
操作,观察其态的变化情况。之后再将第一个比特的状态变回\(\left|0\right>\)态,再观察作用CX
的态的变化情况,执行结果如下所示:
[dechin@dechin-manjaro simulator]$ python3 TestQubits.py
The initial state is: [1.+0.j 0.+0.j 0.+0.j 0.+0.j]
Applying X on the 0th qubit...
The new state is: [2.67948966e-08+0.j 0.00000000e+00+0.j 0.00000000e+00-1.j0.00000000e+00+0.j]
Applying entanglement on qubits...
The new state is: [2.67948966e-08+0.j 0.00000000e+00+0.j 0.00000000e+00+0.j0.00000000e+00-1.j]
Applying X on the 0th qubit...
The new state is: [ 7.17966483e-16+0.00000000e+00j -1.00000000e+00+0.00000000e+00j0.00000000e+00-2.67948966e-08j 0.00000000e+00-2.67948966e-08j]
Applying entanglement on qubits...
The new state is: [ 7.17966483e-16+0.00000000e+00j -1.00000000e+00+0.00000000e+00j0.00000000e+00-2.67948966e-08j 0.00000000e+00-2.67948966e-08j]
这个结果所展示出来的数字也许比较乱,这是因为在运算过程中的计算精度不足所导致的,这里低于1e-06
的数字其实我们可以认为就是0
。那么我们从这个结果中可以分析总结出量子态的演变历程:
注意:上面的这种写法,其实不太合乎程序语言的逻辑,一般从右到左的方向才是从低位到高位的写法。因此,严格来说写法应该是:\(\left|00\right>\Rightarrow\left|01\right>\Rightarrow\left|11\right>\Rightarrow\left|10\right>\Rightarrow\left|10\right>\)。
这里我们就完成了基于魔法函数的量子计算模拟的过程,感兴趣的读者可以自行尝试更多的玩法,这里就不进行更多的测试了!
总结概要
本文主要尝试了用Python的魔法函数__str__
来定义一个量子态,以及使用__or__
来定义一个量子门操作的运算,我们附带的也简单介绍了一下量子计算模拟的一些背景知识。因为程序有简单而明确的执行逻辑,因此用程序语言的方式来定义和理解科学常识,也能够加深对科学的理解。
版权声明
本文首发链接为:https://www.cnblogs.com/dechinphy/p/magic.html
作者ID:DechinPhy
更多原著文章请参考:https://www.cnblogs.com/dechinphy/
参考链接
- https://zhuanlan.zhihu.com/p/24567545
- https://www.cnblogs.com/jhao/p/6022928.html
Python魔法函数与两比特量子系统模拟相关推荐
- exit函数_全面深入了解 Python 魔法函数
(点击上方公众号,可快速关注一起学Python) 作者:浪子燕青 链接: http://www.langzi.fun/Python魔法函数.html 魔法函数概念 魔法函数是以双下划线开头 ...
- python 魔法函数
Python魔法函数 前言 魔法函数 定义 常用的魔法方法 类构造与初始化 类的表示 控制属性访问 比较.运算等操作 容器类操作 可调用对象 序列化 总结 前言 本篇博客主要介绍Python的魔法函数 ...
- 一文读懂什么是Python魔法函数
一文读懂什么是Python魔法函数 Python的魔法函数是指Python的类中,一系列函数名由双下划线包裹的函数. 笔者最初接触到魔法函数的使用是在Pytorch中,在Pytorch中的Datase ...
- python 魔法函数 运行时_16个python常用魔法函数
==,is的使用 ·is是比较两个引用是否指向了同一个对象(引用比较). ·==是比较两个对象是否相等 1.__ init__(): 所有类的超类object,有一个默认包含pass的__ init ...
- python 魔法函数是什么意思_Python 中的魔法函数
魔法函数是Python中的特性,学习好魔法函数将有助于我们写出优秀的pythonic(优雅的.地道的.整洁的)代码,同时因为Python语言的特性,我们在进行框架设计的时候除了设计模式等高级技能,魔法 ...
- python 魔法函数 __add___python魔法方法
python魔法方法 1.__init__: 初始化函数,在创建实例对象为其赋值时使用,在__new__方法之后.__init__必须至少有一个参数__self__,就是这个__new__返回的实例, ...
- python 魔法函数 运行时_python内置函数和魔法函数
内置方法:Python中声明每一个类系统都会加上一些默认内置方法,提供给系统调用该类的对象时使用.比如需要实例化一个对象时,需要调用该类的init方法:使用print去打印一个类时,其实调用的是str ...
- Python魔法函数
python中定义的以__开头和结尾的的函数.可以随意定制类的特性.魔法函数定义好之后一般不需要我们自己去调用,而是解释器会自动帮我们调用. __getitem__(self, item) 将类编程一 ...
- python 魔法函数 __add___PythonI/O进阶学习笔记_2.魔法函数
前言: 本文一切观点和测试代码是在python3的基础上. Content: 1.什么是魔法函数,魔法函数__getitem__在python中应用. 2.python的数据模型和数据模型这种设计对p ...
最新文章
- 快手推荐系统及 Redis 升级存储
- 如何使用Nodejs爬虫看漫画
- R语言ggplot2可视化:ggplot2中使用element_text函数设置轴标签文本粗体字体(bold text,使x轴和Y轴的标签文本都使用粗体字体)、注意是轴标签而非轴标题
- Apache Kafka之设计
- Coursera公开课笔记: 斯坦福大学机器学习第七课“正则化(Regularization)”
- 专家:电脑看多了掉头发怎么办?
- Calendar类的常用成员方法
- 多线程介绍和多线程模块-lock-互斥锁
- 微信小程序接口测试时appid为空如何解决
- HAOI 2006 受欢迎的牛 (洛谷2341)
- opencv遍历像素输出像素值
- 基于FPGA的图像比例缩放
- 计算机配置内存容量怎么调,如何设置电脑虚拟内存,电脑虚拟内存设置多少最合理?...
- 操作系统:覆盖技术与交换技术
- [转载]网络数据流的java处理
- 自媒体账号十万粉丝如何变现?
- 计算机组成原理是答案,计算机组成原理(上)_答案mooc
- 实验四 201771010101 白玛次仁
- Excel修改证件照图片背景色
- git命令行修改用户名密码
热门文章
- 【新知实验室TRTC】
- 信鸽推送集成采坑之代码混淆报错/XINGE: [Util] please add wup-1.0.0.E-SNAPSHOT.jar in your libs
- 2022暑初二信息竞赛学习成果分享2
- vivo手机怎么添加日程提醒?vivo手机添加日程方法
- 解决电信只能在唯一的电脑上拨号上网的问题(MAC克隆)
- java魔剑天使_大天使之剑大神教你打造牛逼魔剑士
- matplotlib之pyplot模块之坐标轴配置(axis():设置坐标轴外观、设置坐标轴范围)
- 将test1文件下的.java文件复制到copy文件下
- Java 算法 矩阵求和
- 大陆高校毕业IEEE fellow榜单