¡ï ¡ï ¡ï ¡ï ¡ï jjdg: ½ð±Ò+5, ¸Ðл·ÖÏí¾Ñé 2017-04-12 00:20:02
Â¥ÉÏÔõô³öÂÒÂëÁË£¡£¿£¡£¿CODE: ################################################################
#### A 3D Drawing Program
#### Version of April 11, 2017
#### Lu Niu
################################################################
import matplotlib.pyplot as plt
import numpy as np
from matplotlib import pyplot as plt
from matplotlib import rc
from mpl_toolkits.mplot3d import Axes3D
fig = plt.figure(figsize=(8, 6))
ax = Axes3D(fig)
plt.rc('font', family='serif')
ax.text2D(0.20, 0.72, 'Left Current', fontsize = 16, transform=ax.transAxes)
ax.text2D(0.22, 0.68, '$\\tau_s=0.1ps$', fontsize = 16, transform=ax.transAxes)
#### J10
with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC00.dat', 'r') as fig:
x = 0.0
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC02.dat', 'r') as fig:
x = 0.2
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC04.dat', 'r') as fig:
x = 0.4
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC06.dat', 'r') as fig:
x = 0.6
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC08.dat', 'r') as fig:
x = 0.8
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC10.dat', 'r') as fig:
x = 1.0
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True)
#### J01
with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC00.dat', 'r') as fig:
x = 0.0
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC02.dat', 'r') as fig:
x = 0.2
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC04.dat', 'r') as fig:
x = 0.4
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC06.dat', 'r') as fig:
x = 0.6
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC08.dat', 'r') as fig:
x = 0.8
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC10.dat', 'r') as fig:
x = 1.0
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True)
#### J00
with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC00.dat', 'r') as fig:
x = 0.0
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC02.dat', 'r') as fig:
x = 0.2
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC04.dat', 'r') as fig:
x = 0.4
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC06.dat', 'r') as fig:
x = 0.6
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC08.dat', 'r') as fig:
x = 0.8
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True)
with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC10.dat', 'r') as fig:
x = 1.0
y = []
z = []
for line in fig:
data = line.split()
y.append((float(data[1])))
z.append((float(data[2])))
ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True)
ax.view_init(30, 30)
ax.set_xlim3d(0, 1)
ax.set_ylim3d(0, 4)
ax.set_zlim3d(-0.10, 0.25)
ax.set_xlabel(r'Voltage $(V)$', family='serif', fontsize=14)
ax.set_ylabel(r'Time $(ps)$', family='serif', fontsize=14)
ax.set_zlabel(r'Current $(nA)$', family='serif', fontsize=14)
plt.xticks([0, 0.2, 0.4, 0.6, 0.8, 1.0], fontsize=12)
plt.yticks([0.0, 0.5, 1.0, 1,5, 2.0, 2.5, 3.0, 3.5, 4.0], fontsize=12)
ax.set_zticks([-0.10,-0.05,0,0.05,0.10,0.15,0.20,0.25])
#ax.set_zticks(np.arange(0,6,1.0))
plt.gca().invert_xaxis()
plt.savefig('GV-Ts001-L.png', dpi=512)
plt.savefig('GV-Ts001-L.eps', dpi=512)
plt.show()