This is the Matlab project post by my professor Actually thi

This is the Matlab project post by my professor. Actually this is my first time use Matlab, I have no idea about that . Please show all detail to teach me how to simulation that by using Matlab. I really appreciate it !

Solution

You may take the Figure 1. mentioned in the document located at https://www.fairchildsemi.com/datasheets/FQ/FQP50N06.pdf

as reference for commercial n-channel mosfet characteristics.

As we know the n_channel mosfet characteristic equation is given by

Id = K(Vgs - Vto)² (1+ lambda * Vdc)

So, we need to draw Id vs Vds for different Vgs values (3V, 3.2V, 4.2V, 6V, 11V)

The matlab program need to be written to draw the simulation waveform.

The following piece of code will give the waveform of Id vs Vgs for given input of Vgs.

Here the Vgs need be entered by the user.

=================================================================

%n-channel enhancement mode MOSFET output characteristics
clear all;
%kn=un*cox = 100 microA/Volts
kn=1e-4;
%Let W/L= 2
%assuming 28nm channel length
W=56*(10^(-9));
L=28*(10^(-9));
beta=kn*W/L;
%Vth is the threshold voltage
Vth=1;
%lambda is the inverse of early voltage in voltage inverse
lambda=1/1000;
%Sweep drain to source voltge from 0 to 10V
vds=0:0.5:10;
%Ask the user to enter gate to source voltage
vgs=input(\'ENTER THE Vgs in volts\');
%Estimate length of the array
m=length(vds);
for i=1:m

if vgs < Vth
current(1,i)=0;
current1(1,i)=0;
elseif vds(i) >= (vgs - Vth)
current(1,i)=0.5* beta * ((vgs - Vth)^2)*(1+lambda*vds(i));
current1(1,i)=0.5* beta * ((vgs - Vth)^2);
elseif vds(i) < (vgs - Vth)
current(1,i)= beta*((vgs-Vth)*vds(i) - 0.5*(vds(i)^2))*(1+lambda*vds(i));
current1(1,i)=beta*((vgs-Vth)*vds(i) - 0.5*(vds(i)^2));
end

end
plot(vds,current(1,:),\'b\',vds,current1(1,:),\'r\')
xlabel(\'Vds, V\')
ylabel(\'Drain Current,A\')
title(\'I-V Characteristics of a MOSFET\')

===============================================

This is an example code and we may tweak the values of Vds, Ids to get better curve as shown in the link.

Similarly, for p_mosfet, the chracteristic equation would become,Id = K(Vgs - Vto)² (1 - lambda * Vdc),

The above matlab code needs to be modified to match the above equation to get the curve as shown in

http://www.physics.csbsju.edu/trace/pMOSFET.CC.html,

If you are a beginner to run matlab code, you may get it here

https://www.youtube.com/watch?v=QOJCvHUkJSs&list=PL2_ZBQcQ1dYjFW9GTAM7iamzSS6bIB2x0