% 2010 EE 482 Homework 7 Question 4 plots using Matlab – By Leo close all, clear all, clc %Constants q = 1.6e-19; esc = 11.7; eox = 3.9; e0 = 8.85e-14; Na = 4e18; xi = 0.5e-7; xox = 1e-7; VT = 25.9e-3; ni = 1e10; phims = 1.07; k = 8.62e-5; T = 300; %Set Vsc range Vsc = 1:.01:4.5; %Depletion charge vs Vsc @ Vcb=0 Vcb = 0; Qdep1 = -sqrt(2*esc*e0*q*Na*abs(Vsc)); %Inversion charge vs Vsc @ Vcb=0 ns = (ni^2/Na)*exp((Vsc-Vcb)/(k*T)); Qinv = -q*ns*xi; %Same two things but with Vcb=1 Vcb=1; Qdep2 = -sqrt(2*esc*e0*q*Na*abs(Vsc)); ns2 = (ni^2/Na)*exp((Vsc-Vcb)/(k*T)); Qinv2 = -q*ns2*xi; %Plot result figure(1); plot(Vsc(1:350),abs(Qdep1(1:350)), 'c'), hold plot(Vsc(1:350),abs(Qdep2(1:350)), 'r') title('Depletion charge vs. Vs') xlabel('Vs (V)') ylabel('|Q| (C)') legend('Cyan: Vcb=0', 'Red: Vcb=1') figure(2); plot(Vsc(1:250),abs(Qinv(1:250)), 'g'), hold plot(Vsc(1:350),abs(Qinv2(1:350)), 'b') title('Inversion charge vs. Vs') xlabel('Vs (V)') ylabel('|Q|Inv (C)')

legend('Green: Vcb=0', 'Blue: Vcb=1') %Vox vs Vs Cox = eox*e0/xox; Vox = -(Qdep1+Qinv)./Cox; Vgc = phims + Vox + Vsc; Vox2 = -(Qdep2+Qinv2)./Cox; Vgc2 = phims + Vox2 + Vsc; %Range of data plotted was modified to fit physical criteria figure(3); plot(Vsc(1:20),Vox(1:20), 'r'), hold plot(Vsc(1:120),Vox2(1:120)), axis tight title('Vox vs Vs') xlabel('Vs (V)') ylabel('Vox (V)') legend('Vcb=0', 'Vcb=1') %Vgc vs Vs figure(4) plot(Vsc(1:20),Vgc(1:20)), hold plot(Vsc(1:120),Vgc2(1:120), 'r') title('Vgc vs. Vs') xlabel('Vs (V)') ylabel('Vgc (V)') legend('Vcb=0', 'Vcb=1',2)

Plots:

1 1.5 2 2.5 3 3.5 4 4.51

1.5

2

2.5x 10

-6 Depletion charge vs. Vs

Vs (V)

|Q| (

C)

Cyan: Vcb=0

Red: Vcb=1

1 1.5 2 2.5 3 3.5 4 4.50

1

2

3

4

5

6

7

8

9x 10

33 Inversion charge vs. Vs

Vs (V)

|Q|In

v (C

)

Green: Vcb=0

Blue: Vcb=1

1 1.2 1.4 1.6 1.8 2

1

2

3

4

5

6

Vox vs Vs

Vs (V)

Vox

(V

)

Vcb=0

Vcb=1

1 1.2 1.4 1.6 1.8 2 2.2 2.42

3

4

5

6

7

8

9

10Vgc vs. Vs

Vs (V)

Vgc

(V

)

Vcb=0

Vcb=1