Create bellstates.py

Author: mhjensen

doc/src/week3/bellstates.py

@@ -0,0 +1,48 @@
+from qiskit import QuantumCircuit, transpile, assemble
+from qiskit.visualization import circuit_drawer
+import qiskit_aer
+import matplotlib.pyplot as plt
+
+# Create a quantum circuit with two qubits
+bell_circuit = QuantumCircuit(2, 2)
+
+# Apply Hadamard gate to the first qubit
+bell_circuit.h(0)
+
+# Apply a CNOT gate with the first qubit as the control and the second qubit as the target
+bell_circuit.cx(0, 1)
+
+# Add measurements to the circuit
+bell_circuit.measure([0, 1], [0, 1])
+
+# Visualize the circuit
+print("Quantum Circuit:")
+print(bell_circuit)
+
+# Number of shots
+num_shots = 10000
+
+# Simulate the circuit using the Aer simulator
+simulator = qiskit_aer.Aer.get_backend('qasm_simulator')
+
+# Transpile the circuit for the simulator
+transpiled_circuit = transpile(bell_circuit, simulator)
+
+# Execute the transpiled circuit on the simulator with the specified number of shots
+result = simulator.run(transpiled_circuit, shots=num_shots).result()
+
+# Get measurement counts
+counts = result.get_counts(bell_circuit)
+
+# Get and print the measurement results
+counts = result.get_counts(bell_circuit)
+print("\nMeasurement Results:")
+print(counts)
+
+# Plot the histogram using Matplotlib
+plt.bar(counts.keys(), counts.values())
+plt.title('Bell State Measurement Results')
+plt.ylabel('Counts')
+plt.show()
+
+