学习笔记 -《量子计算与编程入门》- 量子程序

1 量子计算原理

经典计算中，最基本的单元是比特，最基本的控制模式是逻辑门，可以通过逻辑门的组合来达到控制电路的目的。
类似地，处理量子比特的方式就是量子逻辑门，使用量子逻辑门可以有意识地使量子态发生演化，所以量子逻辑门是构成量子算法的基础。

1.1 酉变换

酉变换是一种矩阵，也是一种操作，它作用在量子态上得到的是一个新的量子态。之前博客也提到过：

$U^{†}$ 表示 $U$ 的转置复共轭矩阵，酉矩阵的转置复共轭矩阵也是一个酉矩阵，酉变换是一种可逆变换。

一般酉变换在量子态上的作用是通过左乘以右矢进行的，例如：
$$|\psi ⟩=U|{\phi }_0⟩或⟨\psi |=⟨{\phi }_0|U^{†}$$ 由此可见，两个矢量的内积经过同一个酉变换之后保持不变：
$$⟨\phi |\psi ⟩=⟨\phi |U†U|\psi ⟩$$ 类似地，也可以通过酉变换表示密度矩阵的演化：
$$\rho =U{\rho }_0{U}^{†}$$这样就连混合态的演化也包含在内了。

1.2 矩阵的指数函数

在矩阵乘法的基础上，我们利用函数的幂级数来定义矩阵的指数函数。我们定义矩阵 A 的指数函数形式为：
$$exp(A)=I+A+\frac{A^2}{2\text{！}}+\frac{A^3}{3\text{！}}+...$$

如果 $A$ 是一个对角矩阵，即 $A=diag\left(A_{11},A_{22},A_{33},\dots \right)$，则有：
$$A^n=diag\left(A_{11}^n,A_{22}^n,A_{33}^n,\dots \right)$$ 从而得到：
$$\begin{array}{rl}exp(A)& =I+A+\frac{A^2}{2\text{！}}+\frac{A^3}{3\text{！}}+...\\ & =I+diag\left(A_{11}^{},A_{22}^{},A_{33}^{},\dots \right)+\frac{diag\left(A_{11}^2,A_{22}^2,A_{33}^2,\dots \right)}{2\text{！}}+\frac{diag\left(A_{11}^3,A_{22}^3,A_{33}^3,\dots \right)}{3\text{！}}+...\\ & =diag\left(1+A_{11}^{}+\frac{A_{11}^2}{2!}+\frac{A_{11}^3}{3!}+\dots ,A_{22}^{}+\frac{A_{22}^2}{2!}+\frac{A_{22}^3}{3!}+\dots ,A_{33}^{}+\frac{A_{33}^2}{2!}+\frac{A_{33}^3}{3!}+\dots ,\dots \dots \right)\\ & =\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left(e^{A_{11}},e^{A_{22}},e^{A_{33}},\dots \right)\end{array}$$ 如果 $A$ 不是一个对角矩阵，则利用酉变换可以将它对角化：$A=UD{U}^{†}$。
则：$$\begin{array}{rl}{A}^n& ={\left(UD{U}^{†}\right)}^n\\ & =\left(UD{U}^{†}\right)\left(UD{U}^{†}\right)\left(UD{U}^{†}\right)\dots \dots \\ & =U{D}^n{U}^{†}\end{array}$$ 那么，类似地
$$\begin{array}{rl}\exp \left(A\right)& =\exp \left(UD{U}^{†}\right)\\ & =U{U}^{†}+UD{U}^{†}+\frac{U{D}^2{U}^{†}}{2\text{！}}+\frac{U{D}^3{U}^{†}}{3\text{！}}+...\\ & =U{U}^{†}+UD{U}^{†}+U\frac{D^2}{2\text{！}}{U}^{†}+U\frac{D^3}{3\text{！}}{U}^{†}+\dots \\ & =U\left(I+D+\frac{D^2}{2\text{！}}+\frac{D^3}{3\text{！}}+\dots \right)U^{†}\\ & =U\exp \left(D\right)U^{†}\\ & =U\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left(e^{D_{11}},e^{D_{22}},e^{D_{33}},\dots \right)U^{†}\end{array}$$

注意：
$$exp\left(A+B\right)\ne exp\left(A\right)exp\left(B\right)\ne exp\left(B\right)exp\left(A\right)$$ 我们知道，当 $A$ 是表示数的时候等号是成立的，但是当 A 表示是矩阵时，等式成立需要满足一定条件。

---

通常，下面这种表达形式被称之为以 $A$ 为生成元生成的酉变换。
$$U(\theta )=exp(-i\theta A)$$ 这种矩阵的指数运算可以利用 Matlab 中的 expm 命令进行方便地计算。

以单位矩阵为生成元，可以构建一种特殊的酉变换 $u(\theta )$。
$$u(\theta )=exp(-i\theta I)=\left(\begin{array}{cc}{e}^{-i\theta }& 0\\ 0& e^{-i\theta }\end{array}\right)=exp(-i\theta )I$$ $u(\theta )$ 作用在态矢上面，相当于态矢每个分量同时乘以一个系数 $e^{-i\theta }$。如果将作用后的态矢带入到密度矩阵的表达式中，会发现这一项系数会被消去。

这项系数称为量子态的整体相位。因为任何操作和测量都无法分辨两个相同的密度矩阵，所以量子态的整体相位一般情况下是不会对系统产生任何影响的。

1.3 单量子比特逻辑门

在经典计算机中，单比特逻辑门只有一种——非门（NOT gate），但是在量子计算机中，量子比特情况相对复杂，存在叠加态、相位，所以单量子比特逻辑门会有更加丰富的种类。

量子比特门可以用下面这幅图来表示：

横线表示一个量子比特从左到右按照时序演化的路线，通俗点讲，就是 $|{\psi }_0⟩$ 演变为 $|{\psi }_1⟩$ 的路线；

蓝色方框表示量子逻辑门，这个图标表示一个名为 $U$ 的逻辑门作用在这条路线所代表的量子比特上。

对于一个处于 $|{\psi }_0⟩$ 的量子态，将这个量子逻辑门作用在上面时，相当于将这个代表量子逻辑门的酉矩阵左乘这个量子态的矢量，然后得到下一个时刻的量子态 $|{\psi }_1⟩$：
$$|{\psi }_1⟩=U|{\psi }_0⟩$$ 这个表达式对于所有的单比特门或者多比特门都是适用的。

对于一个有 $n$ 个量子比特的量子系统，它的演化是通过一个 $2^n\times 2^n$ 的酉矩阵来表达。

1.3.1 泡利矩阵

泡利矩阵（Pauli matrices）有时也被称作自旋矩阵（spin matrices）。有以下三种形式，分别是：
$${\sigma }_x=\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right){\sigma }_y=\left(\begin{array}{cc}0& -i\\ i& 0\end{array}\right){\sigma }_z=\left(\begin{array}{cc}1& 0\\ 0& -1\end{array}\right)$$ 三个泡利矩阵所表示的泡利算符代表着对量子态矢量最基本的操作。

如将 ${\sigma }_x$作用到 $|0⟩$ 态上，经过矩阵运算，得到的末态为 $|1⟩$ 态。

泡利矩阵的线性组合是完备的二维酉变换生成元，即所有满足 $U{U}^{†}=I$ 的 $U$ 都能通过下面这种方式得到：
$$U=e^{-i\theta \left(a{\sigma }_x+b{\sigma }_y+c{\sigma }_z\right)}$$

1.3.2 常见逻辑门以及含义

· Hadamard (H) 门

形式：$$H=\frac{1}{\sqrt{2}}\left[\begin{array}{cc}1& 1\\ 1& -1\end{array}\right]$$ 线路表示：

功能：

• 作用在单量子比特上：将基态变为叠加态
  $$\frac{1}{\sqrt{2}}\left[\begin{array}{cc}1& 1\\ 1& -1\end{array}\right]|0⟩=\frac{|0⟩+|1⟩}{\sqrt{2}}$$$$\frac{1}{\sqrt{2}}\left[\begin{array}{cc}1& 1\\ 1& -1\end{array}\right]|1⟩=\frac{|0⟩-|1⟩}{\sqrt{2}}$$
• 作用在任意量子态 $|\psi ⟩=\alpha |0⟩+\beta |1⟩$ 上：
  $$|{\psi }^{\prime }⟩=H|\psi ⟩=\frac{1}{\sqrt{2}}\left[\begin{array}{cc}1& 1\\ 1& -1\end{array}\right]\left[\begin{array}{c}\alpha \\ \beta \end{array}\right]=\frac{1}{\sqrt{2}}\left[\begin{array}{c}\alpha +\beta \\ \alpha -\beta \end{array}\right]=\frac{\alpha +\beta }{\sqrt{2}}|0⟩-\frac{\alpha -\beta }{\sqrt{2}}|1⟩$$

· Pauli-X 门

形式：$$X={\sigma }_x=\left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right]$$ 线路表示：
功能：

• 作用在单量子比特上：是经典计算机 NOT 门的量子等价，即将量子态进行翻转，因此又称 NOT 门。
  $$X|0⟩=\left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right]\left[\begin{array}{c}1\\ 0\end{array}\right]=\left[\begin{array}{c}0\\ 1\end{array}\right]=|1⟩$$ X ∣ 1 > = [