MATLAB Answers

Plotting Schottkly diode I-V chracteristics

13 views (last 30 days)
I am trying to plot current density vs voltage for the schottky diode.
This is my code so far:
clear all
n = 1.02; %Ideality factor
A = 33.65; %Richardson constant
q = 1.602*10^-19; % electron charge
K = 1.38*10^-23; %Boltzmann constant
T = 300; % Absolute temperature
phi_b = 0.28; %Barrier Height
V = 0:0.2:2;
J_0 = A* (T* T)* exp((-q* phi_b)./(K* T));
J = J_0 * exp((q* V)./(n* K* T))* (1 - (exp(-q* V)./(K* T)));
plot(V, log(J ./(1 - exp(-q* V)/(k *T))))

Accepted Answer

Srijith Kasaragod
Srijith Kasaragod on 2 Sep 2021
Edited: Srijith Kasaragod on 7 Sep 2021
Assuming correctness of parameters and equations, error exists only in final 2 lines of code. Multiplication to compute J must be elementwise multiplication rather than matrix multiplication. Also Botzmann constant is saved as variable 'K' and not 'k'. Making these two changes generates the plot.
J = J_0* exp((q* V)./(n* K* T)).* (1 - (exp(-q* V)./(K* T)));
plot(V, log(J ./(1 - exp(-q* V)/(K *T))));
Srijith Kasaragod
Srijith Kasaragod on 3 Sep 2021
You have 11 voltage values in array 'V'. Your requirement is to find corresponding 11 current values and plot the I-V characteristics. The two terms being multiplied in line 11 are both of size 1-by-11. To get the required current values, you have to perform element-wise multiplication between these two arrays. Matrix multiplication would give either 1-by-1 or 11-by-11 matrix according to the order of matrix multiplication, which is not the desired output.

Sign in to comment.

More Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!