effvscj.m.bak

求解导通角的matlab源代码

clear;
echo off;

Omega = 2*pi*9e8;

  S = solve ('Rs=10','Vb=0.2','N=2','Vd=1.5','Rl=5e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');

  ThetaOn = double(S.ThetaOn);
  Rs = double(S.Rs);
  Vb = double(S.Vb);
  N = double(S.N);
  Vd = double(S.Vd);
  Rl = double(S.Rl);
  corner = 360*ThetaOn/2/pi;
  
for i=1:1:100;

  Cj(i) = i*1e-14;
  
  Eta1(i) = (N*Rl/pi/Rs)*((0.5/N+Vb/Vd)^2)*(ThetaOn*(2+1/(cos(ThetaOn)^2))-3*tan(ThetaOn));
  Eta2(i) = (N*Rs*Rl*((Cj(i))^2)*(Omega^2)/pi)*((0.5/N+Vb/Vd)^2)*((pi-ThetaOn)/((cos(ThetaOn))^2)+tan(ThetaOn));
  Eta3(i) = (2*N*Rl/pi/Rs)*(0.5/N+Vb/Vd)*(Vb/Vd)*(tan(ThetaOn)-ThetaOn);
  Eta(i)  = 1/(1+Eta1(i)+Eta2(i)+Eta3(i));

end

plot (Cj, Eta);

grid on;

echo on;