1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
|
%function A=optim(x);
Vgs=-0.9:0.1:0;
Vds=0:1:6;
Vpk=-0.2;
p=[1.561 1.409 5.379];
alp=1.156;
L=0.00739;
Ipk=0.761;
%Vds=0:1:6;
a0=(-p(1)*Vpk)+(p(2)*Vpk.^2)-(p(3)*Vpk.^3);
a1=p(1)-(2*p(2)*Vpk)+(3*p(3)*Vpk.^2);
a2=p(2)-(3*p(3)*Vpk);
a3=p(3);
hold on;
for i=1:10,
psi(i)=a0+(a1*Vgs(i))+(a2*Vgs(i)^2)+(a3*Vgs(i)^3);
for j=1:7 ,
Ids(i,j)=Ipk*(1+tanh(psi(i)))*(1+(Vds(j)*L))*tanh(Vds(j)*alp);
end;
figure(1);
plot(Vds,Ids(i,:));
pause;
end;
hold off;
figure(2);
plot(Vgs,psi);
%%%%%%%%%%%
function G=optimal(L);
G=Ipk*(1+tanh(psi(i)))*(1+(Vds(j)*L))*tanh(Vds(j)*alp); |
Partager