bonjour à tous !!!
tous d'abord je m'excuse pour se long silence.
en fait, je suis confronté à un problème qui me torture depuis plus de 4 mois.
le problème c'est: je n'arrive pas à tracer la dérivé de la courbe comme le montre le code matlab ci-joint.
Code : Sélectionner tout - Visualiser dans une fenêtre à part
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)
à chaque fois ça me donne une droite verticale. je voudrais que vous y jetiez un cout d'oeil please