xuchuanxiang.m

晶闸管设计的徐传镶法

function [a,b]=xuchuanxiang(romin,VDSM,tp);
%设计陆晓东
%日期2008.7.18
rox=linspace(romin,400,800);
Xn=0.531*(rox.*VDSM).^(1/2);
if VDSM>=1500;
    Vb=96*rox.^(3/4);
else
    Vb=126*rox.^(0.63);;
end
Vb=100*rox.^(3/4);
Dp=12.40;
Lp=(Dp*tp*10^-6)^(1/2);
alfa1=1-(VDSM./Vb).^4;
Wen1=10^4*Lp*log(1./alfa1+((1./alfa1).^2-1).^(1/2));
Wn1=Xn+Wen1;
a=round(min(Wn1));
for xx=1:length(rox);
   XXn=0.531*(rox(xx)*VDSM)^(1/2);
   Vbb=100*rox(xx)^(3/4);
   alfa1=1-(VDSM./Vbb).^4;
   XWen1=10^4*Lp*log(1/alfa1+((1/alfa1)^2-1)^(1/2));
   WWn1=XXn+XWen1;
   if a==round(WWn1);
       b=rox(xx);
   end
end

if VDSM>1500;
  Vb=96*rox.^(3/4);
else
   Vb=126*rox.^(0.63); 
end
for xx=1:10;
tp=5*xx;
Lp=(Dp*tp*10^-6)^(1/2);
alfa1=1-(VDSM./Vb).^4;
Wen1=10^4*Lp*log(1./alfa1+((1./alfa1).^2-1).^(1/2));
Wn1=Xn+Wen1;;
subplot(2,2,1);
plot(rox,Wn1);
hold on;
axis([romin,400,200,600])
end