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 33 34 35 36 37 38
| gp=0
B1=5e-6;
B2=7e-6;
C0n=110e-12;
%Cs=95e-12;
Cso=0;
B=B1.*B2./(B1+B2);
j=1;
for kp=0:0.001:pi/2
Q1=-(-2.*Cso^2+(2.*Cso+C0n)^2)^2;
Q2=-8.*cos(2.*kp).^2.*Cso^3./B-(2.*(-2.*C0n./B-2.*Cso./B-gp^2));
Q3=4.*cos(2.*kp).^2.*(2.*Cso^2./(B1.*B2)+Cso^2./B^2)-(2.*(1./B^2-1./B1^2-1./B2^2)).*(-2.*Cso^2+(2.*Cso+C0n)^2);
Q4=-8.*cos(2.*kp).^2.*Cso./(B1.*B2.*B);
Q5=-(1./B^2-1./B1^2-1./B2^2)^2+4.*cos(2.*kp).^2./(B1^2.*B2^2);
pol=[Q1 Q2 Q3 Q4 Q5];
eigen=roots(pol)
for n=1:length(eigen)
if real(eigen(4))>=0
OME(j,1)=kp ;
OME1(j,2) =real(sqrt(eigen(4)));
j=j+1
end
end
end
%\\\\\\\\\\\\\\\\\\\\\courbe proprement dite\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
plot(OME(:,1),OME1(:,2),'.')
%/////////////////////tracé de la dérivé//////////////////////////////////
uk=diff(OME1(:,2))./diff(OME(:,1))% définition de la dérivé
plot(OME1(length(OME1)-1),uk) |
Partager