In [2]:

    

import numpy as np
import matplotlib.pyplot as plt
import os
import scipy.optimize as opt
%matplotlib inline


    

In [13]:

    

def fermi(x,a,b):
    return 1/(1+np.exp(b*(x-a)))
def fermi1(x,a,b):
    return 1/(1+np.exp((a-x)/100))
data=np.arange(0,40,0.1)

y1=fermi(data,20,0.5)
# y2=fermi(data,20,1)
# y3=fermi(data,20,2)
# y4=fermi(data,20,5)
# y5=fermi(data,20,10)
y6=fermi(data,20,100)

plt.figure(figsize=(8,8))
ax=plt.subplot(111)
ax.plot(data,y1,'#FF0000',data,y6,'blue',linewidth=3)#,data,y3,'yellow',data,y4,'green',data,y5,'blue',data,y6,'purple',linewidth=3)
ax.axis([0,max(data),-0.1,max(y1)+0.1])
ax.set_xlabel("Vg",fontsize=20)
ax.set_ylabel("Ids",fontsize=20)


    

    




/home/leo/anaconda3/lib/python3.5/site-packages/ipykernel/__main__.py:2: RuntimeWarning: overflow encountered in exp
  from ipykernel import kernelapp as app






    
Out[13]:





<matplotlib.text.Text at 0x7fbdf8940908>






    

In [16]:

    

data=np.loadtxt('235RC/235RC_pin03-13_8_RvsT_300s_-15Vg_2016_06_28',skiprows=1)
def f2(x,a,b,c):
    return a*np.exp((b*x))+c
# fit=opt.curve_fit(f2,np.log(data[1:,0]),data[1:,1]/1e7,p0=(1,-0.05,4))
# print(fit,data[:10,1])

plt.plot(data[1:,0],data[1:,1]/1e7,"r-")#,data[1:,0], f2( data[1:,0],1,-0.05,4), 'b-')
plt.xlabel("Time (s)",fontsize=15)
plt.ylabel("Resistance $\Omega$",fontsize=15)


    

    
Out[16]:





<matplotlib.text.Text at 0x7fd1060cbf98>






    

In [ ]:

    




    

In [41]:

    

data2=np.loadtxt('235RC/235RC_pin03-13_15_RvsT_600s_updown_15Vg_2016_06_28',skiprows=1)
data1=np.loadtxt('235RC/235RC_pin03-13_14_RvsT_600s_updown_15Vg_2016_06_28',skiprows=1)
plt.figure(figsize=(16,7))
ax1=plt.subplot(121)
ax1.plot(data1[:,0],data1[:,1])

ax1.set_xlabel("Time (s)",fontsize=15)
ax1.set_ylabel("Resistance ($\Omega$)",fontsize=15)
ax2=plt.subplot(122)
ax2.plot(data2[:,0],data2[:,1])

ax2.set_xlabel("Time (s)",fontsize=15)
ax2.set_ylabel("Resistance ($\Omega$)",fontsize=15)
fig.suptitle("235RC Device 1 (03-13) R(t) with $\pm$15V pulses",fontsize=20)
offset_text = ax1.yaxis.get_offset_text()
offset_text.set_size(15)
offset_text = ax2.yaxis.get_offset_text()
offset_text.set_size(15)


    

In [ ]: