>>> import qutip >>> import numpy as np >>> psi = qutip.Qobj ([[ 0.6 ] , [ 0.8 ]]) # สร้างสถานะควอนตัมจากลิสต์>>> psi วัตถุควอนตัม: dims = [[2], [1]], shape = (2, 1), type = ket ข้อมูลQobj = [[0.6] [ 0.8 ]] >>> phi = qutip.Qobj ( np.array ([ 0.8 , -0.6 ])) # สร้างสถานะควอนตัมจากอาร์เรย์ numpy >>> phi วัตถุควอนตั ม : dims = [[2], [1]], shape = ( 2, 1), type = ket ข้อมูล Qobj = [[0.8] [-0.6]] >>> e0 = qutip . basis ( 2 , 0 ) # สร้างเวกเตอร์ฐาน>>> e0 วัตถุควอนตัม: dims = [[2], [1]], shape = (2, 1), type = ket ข้อมูล Qobj = [[1.] [0.]] >>> A = qutip . Qobj ( np . array ([[ 1 , 2 j ], [ - 2 j , 1 ]])) # สร้างตัวดำเนินการควอนตัมจากอาร์เรย์ numpy >>> A วัตถุควอนตัม: dims = [[2], [2]], shape = (2, 2), type = oper, isherm = True ข้อมูล Qobj = [[1.+0.j 0.+2.j] [0.-2.j 1.+0.j]] >>> qutip . sigmay () # วัตถุควอนตัมทั่วไปบางอย่าง เช่น เมทริกซ์ Pauli ถูกกำหนดไว้ล่วงหน้าในแพ็คเกจ qutip วัตถุควอนตัม: dims = [[2], [2]], shape = (2, 2), type = oper, isherm = True ข้อมูล Qobj = [[0.+0.j 0.-1.j] [0.+1.j 0.+0.j]]

>>> A * qutip . sigmax () + qutip . sigmay () # เราสามารถบวกและคูณวัตถุควอนตัมที่มีรูปร่างและมิติที่เข้ากันได้วัตถุควอนตัม: dims = [[2], [2]], shape = (2, 2), type = oper, isherm = False ข้อมูล Qobj = [[0.+2.j 1.-1.j] [1.+1.j 0.-2.j]] >>> psi . dag () # คอนจูเกตเฮอร์มิเชียนวัตถุควอนตัม: dims = [[1], [2]], shape = (1, 2), type = bra ข้อมูล Qobj = [[0.6 0.8]] >>> psi . proj () # โปรเจกเตอร์ไปยังสถานะควอนตัมวัตถุควอนตัม: dims = [[2], [2]], shape = (2, 2), type = oper, isherm = True ข้อมูล Qobj = [[0.36 0.48] [0.48 0.64]] >>> A . tr () # ร่องรอยของตัวดำเนินการ2.0 >>> A . eigenstates () # แปลงตัวดำเนินการให้เป็นแนวทแยง(array([-1., 3.]), array([Quantum object: dims = [[2], [1]], shape = (2, 1), type = ket  Qobj data =  [[-0.70710678+0.j ]  [ 0. -0.70710678j]] ,  Quantum object: dims = [[2], [1]], shape = (2, 1), type = ket  Qobj data =  [[-0.70710678+0.j ]  [ 0. +0.70710678j]] ],  dtype=object)) >>> ( 1 j * A ) . expm () # เมทริกซ์เอกซ์โพเนนเชียลของตัวดำเนินการวัตถุควอนตัม: dims = [[2], [2]], shape = (2, 2), type = oper, isherm = False ข้อมูล Qobj = [[-0.2248451-0.35017549j -0.4912955-0.7651474j ] [ 0.4912955+0.7651474j -0.2248451-0.35017549j]] >>> qutip . เทนเซอร์( qutip . sigmaz (), qutip . sigmay ()) # ผลคูณเทนเซอร์ วัตถุควอนตัม: dims = [[2, 2], [2, 2]], shape = (4, 4), type = oper, isherm = True ข้อมูล Qobj = [[0.+0.j 0.-1.j 0.+0.j 0.+0.j] [0.+1.j 0.+0.j 0.+0.j 0.+0.j] [0.+0.j 0.+0.j 0.+0.j 0.+1.j] [0.+0.j 0.+0.j 0.-1.j 0.+0.j]]

>>> Hamiltonian = qutip.sigmay ( ) >>> times = np.linspace ( 0 , 2 , 10 ) >>> result = qutip.sesolve ( Hamiltonian , psi , times , [ psi.proj ( ) , phi.proj ( ) ] ) # วิวัฒนาการเวลา แบบเอกภาพของระบบตามสมการชโรดิงเจอร์>>> expectpsi , expectphi = result.expect # ค่าคาดหวัง ของตัว ฉายภาพบน psi และ phi >>> plt.figure ( dpi = 200 ) >>> plt.plot ( times , expectpsi ) >>> plt.plot ( times , expectphi ) >>> plt.legend ( [ r " $ \ psi$" , r " $\phi $ " ] ) >>> plt.show ( )