effvsrl.m

求解导通角的matlab源代码

clear;
echo off;

Cj = 1e-13;
Omega = 2*pi*9e8;

  S(1) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.2','Rl=3e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');
  S(2) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.2','Rl=6e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');
  S(3) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.2','Rl=9e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');
  S(4) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.2','Rl=12e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');
  S(5) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.2','Rl=15e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');
  S(6) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.2','Rl=18e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');
  S(7) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.2','Rl=21e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');
  S(8) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.2','Rl=24e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');
  S(9) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.2','Rl=27e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');
  S(10) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.2','Rl=30e4','tan(ThetaOn)-ThetaOn=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','ThetaOn');

for i=1:1:10;

  ThetaOn(i) = S(i).ThetaOn;
  Rs(i) = S(i).Rs;
  Vb(i) = S(i).Vb;
  N(i) = S(i).N;
  Vd(i) = S(i).Vd;
  Rl(i) = S(i).Rl;
  
  corner(i) = double (ThetaOn(i)/2/pi);
  x(i) = double (Rl(i));
  
  Eta1(i) = (N(i)*Rl(i)/pi/Rs(i))*((1/2/N(i)+Vb(i)/Vd(i))^2)*(ThetaOn(i)*(2+1/(cos(ThetaOn(i))^2))-3*tan(ThetaOn(i)));
  Eta2(i) = (N(i)*Rs(i)*Rl(i)*(Cj^2)*(Omega^2)/pi)*((1/2/N(i)+Vb(i)/Vd(i))^2)*((pi-ThetaOn(i))/((cos(ThetaOn(i)))^2)+tan(ThetaOn(i)));
  Eta3(i) = (2*N(i)*Rl(i)/pi/Rs(i))*(1/2/N(i)+Vb(i)/Vd(i))*(Vb(i)/Vd(i))*(tan(ThetaOn(i))-ThetaOn(i));
  Eta(i)  = 1/(1+Eta1(i)+Eta2(i)+Eta3(i));
  
  y(i) = double (Eta(i));

end

plot (x, y);

grid on;

echo on;