mat3.m

mos transistor, basic simulation

hold on
vgs=0.001:.001:1.8;
vds=1.8;
for vbs=0:.010:.5
    vt=.4+.45*(sqrt(.9-vbs)-sqrt(.9));
    for i=1:1800
        if vgs(i)<vt
        id(i)=vgs(i).*0;
        end
        if vds<vgs(i)-vt && vgs(i)>=vt
        id(i)=.774.*((vgs(i)-vt).*vds-vds.^2./2).*(1+.1.*vds);
        end
        if vds>=vgs(i)-vt && vgs(i)>=vt
        id(i)=.774./2.*((vgs(i)-vt)^2).*(1+.1.*vds);
        end
    end
    plot(vgs,id)
    end