effvscj.m

求解导通角的matlab源代码

clear;
echo off;

Omega = 2*pi*9e8;

  S = solve ('tan(x)-x=pi*10/(5e4*(0.5/2+0.2/1.5))','x');

  ThetaOn = double(S);
  corner = 360*ThetaOn/2/pi;
  
for i=1:1:100;

  Cj(i) = i*1e-14;
  
  Eta1(i) = (2*5e4/pi/10)*((0.5/2+0.2/1.5)^2)*(ThetaOn*(2+1/(cos(ThetaOn)^2))-3*tan(ThetaOn));
  Eta2(i) = (2*10*5e4*((Cj(i))^2)*(Omega^2)/pi)*((0.5/2+0.2/1.5)^2)*((pi-ThetaOn)/((cos(ThetaOn))^2)+tan(ThetaOn));
  Eta3(i) = (2*2*5e4/pi/10)*(0.5/2+0.2/1.5)*(0.2/1.5)*(tan(ThetaOn)-ThetaOn);
  Eta(i)  = 1/(1+Eta1(i)+Eta2(i)+Eta3(i));

end

plot (Cj, Eta);

grid on;

echo on;