bjt

function [ z ] = bjt( iceo,vceQ,etavt,beta,va,vce,vbeSat,ib )

%BJT Visualizzazione caratteristiche BJT

 

ic=iceo*(1-exp((vbeSat-vceQ)/etavt))+beta*ib*(1+vceQ/va);

betaF=ic/ib;

 

vcesim=linspace(-0.1,max(vce),100);

ibsim=linspace(0,0.1,10);

ibsim=[ibsim ib];

[vceg,ibg]=meshgrid(vcesim,ibsim);

icsim=iceo*(1-exp((vbeSat-vceg)/etavt))+beta*ibg.*(1+vceg/va);

 

plot(vcesim,icsim,vceQ,ic,'or');

xlabel('Vce [V]');

ylabel('Ic [mA]');

title('Curve caratteristiche del BJT');

 

ax=axis;

axis([ax(1) ax(2) -1 ax(4)]);

xlabel('Vce [V]');

ylabel('Ic [mA]');

riga1=[' IB=' num2str(ib) ' mA'];

riga2=['VCE=' num2str(vceQ) ' V'];

riga3=['BETAF=' num2str(betaF) ' '];

text(ax(1),ax(4)-1,riga1);

text((ax(1)+ax(2))/2,ax(4)-1,riga2,'HorizontalAlignment','center');

text(ax(2),ax(4)-1,riga3,'HorizontalAlignment','right');