tracé de la dérivée d'une courbe ou d'un ensemble de points

Version imprimable

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:

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 :calim2:

29/01/2014, 18h31
Invité
Bonsoir,

La partie OME1(length(OME1)-1) dans la ligne

Code:

plot(OME1(length(OME1)-1),uk)

ne représente qu'une seule valeur, tu traces donc toutes les valeurs dans uk sur la même abscisse.

Ne voulais-tu pas quelque chose comme :
Code:

1 2 plot(OME1(1:end-1,1), uk) % première colonne plot(OME1(1:end-1,2), uk) % deuxième colonne
?
[EDIT]Plutôt OME que OME1 en abscisses d'après ce que je vois.
30/01/2014, 18h48
doinier

thank jerome

Merciiiiii, Winjerome là je commence à avoir une lieur d'espoir c'était parfait :ave:

1 pièce(s) jointe(s)

choix de courbe

bonsoir !!!
tous d'abord je m'excuse pour le silence. je voudrais ne tracer qu'une seule courbe celle indiquée sur le fichier joint. mais à chaque fois que je

prend n=1,2,3,4,5,6...etc, j'ai comme un mélange de courbe pourtant je n'en veux qu'une seule, celle indiquée sur la courbePièce jointe 155273

voilà mon code.

Code:

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
39
40
41
42
43
44
45
46
47
48
49
C0=540e-12;
 
%*************************************************************
L2=14e-6;
L1=28e-6;
S1=5;
S2=S1*L2/L1;
 
I1=1/(L1*C0);
I2=1/(L2*C0);
 
Ir1=S1/L1;
Ir2=S2/L2;
 
It=I1+I2;
IT=I1+I2;
 
 
j=1;
for k=0:0.001:1.6
%----------------déclaration du polynôme-----------------------------------
Q8=I2.^2.*I1.^2.*(1-cos(2.*k).^2);
Q6=(-2.*Ir1.*Ir2-4.*It).*I2.^2.*I1.^2.*(1-cos(2.*k).^2)-(Ir1+Ir2).^2.*I2.^2.*I1.^2.*cos(2.*k).^2;
Q4=(-2.*I1.^2-2.*I2.^2+2.*It.^2+(-Ir1.*Ir2-2.*It).^2).*I2.^2.*I1.^2.*(1-cos(2.*k).^2)-(2.*(-Ir1.*It-Ir2.*It)).*(Ir1+Ir2).*I2.^2.*I1.^2.*cos(2.*k).^2;
Q2=(2.*(-I1.^2-I2.^2+It.^2)).*(-Ir1.*Ir2-2.*It).*I2.^2.*I1.^2.*(1-cos(2.*k).^2)-(-Ir1.*It-Ir2.*It).^2.*I2.^2.*I1.^2.*cos(2.*k).^2;
Q0=+(-I1.^2-I2.^2+It.^2).^2.*I2.^2.*I1.^2.*(1-(cos(2.*k)).^2);
 
 
pol=[ Q8 0 Q6 0 Q4 0 Q2 0 Q0];
%--------------------------------------------------------------------------
%------------------Résolution du polynôme----------------------------------
eigen=roots(pol);
 
for n=1:numel(eigen)  
 
if  real(eigen(n))>=0 
 
  OMER(j,1)=k ;      % range les éléments dans une matrice 
 
  OME1(j,1)=real(eigen(n)); % range la partie réelle de eigen dans une matrice
 
  j=j+1;
 
end
 
end
end
%---------tracé des courbe----------------
plot(OMER(:,1),OME1(:,1),'.')

et aussi je devrais avoir 8 solutions, mais je n'en n'ai que 6, car pour n=7,8; y a rien. pourtant

Code:

length(eigen)

me donne 8

s'il vous plaît help me.:calim2:

04/09/2014, 11h23
lucaswalker

Bonjour,

Le problème est que pour k = 0 et k = pi/2 la courbe est verticalement asymptotique.

Essayez avec des valeurs de k excluant ces points (for k=0.005:0.001:1.565) par exemple, et on retrouve la jolie courbe. Au passage il y a bien une unique courbe, mais avec les asymptotes tous les points sont concentrés et on comprend plus rien.
05/09/2014, 17h40
doinier

choix de courbe

bonjour lucaswalker!!!

merci pour ton intervention, j’essaie pour voir. mais je me réserve de boucler cette discussion car quelqu'un pourrait entre temps avoir quelque chose à ajouter. merci encore.

choix d'une courbe

bonjour lucaswalker,

j'ai essayé le code suivant pour choisir la courbe sélectionner mais ça n'affiche que le deuxième bloc:
Code:

1 2 3 4 5 6 7 8 9 10 plot(OMER2(:,1),OME1(:,1),'.') % hold on % ylim([0 7130000]) % xlim([0 0.8]) %%%%% deuxième bloc%%%%%% % hold on % plot(OMER2(:,1),OME1(:,1),'.') % ylim([7130000 14110000]) % xlim([0.8 1.57]) % hold off
merci d'y jeter un coup d’œil s'il te plait

le plot ne trace qu'une unique courbe, pour obtenir le plot que tu veux il faut par conséquent sélectionner les bonnes valeurs dans tes vecteurs. Le code suivant peut faire l'affaire:
Code:

1 2 3 4 5 6 %---------tracé des courbe---------------- hold on plot(OMER(2:2:length(OME1(:,1))/2,1),OME1(2:2:length(OME1(:,1))/2,1),'.') plot(OMER(length(OME1(:,1))/2+1:2:length(OME1(:,1)),1),OME1(length(OME1(:,1))/2+1:2:length(OME1(:,1)),1),'.') ylim([0 14000000]) xlim([0 1.6])

Merci lucaswalker, ça fonctionne :D.

mais je voudrais stp que tu m'explique ses lignes . là je pourrais modifier les zone à ma guise c-à-dire
Code:

1 2 plot(OMER2(2:2:length(OME1(:,1))/2,1),OME1(2:2:length(OME1(:,1))/2,1),'.') plot(OMER2(length(OME1(:,1))/2+1:2:length(OME1(:,1)),1),OME1(length(OME1(:,1))/2+1:2:length(OME1(:,1)),1),'.')

18/09/2014, 14h54
lucaswalker
ok,

essaye ces deux lignes:
Code:

1 2 3 plot(OMER2(1:2:length(OME1(:,1)),1),OME1(1:2:length(OME1(:,1)),1),'.') plot(OMER2(2:2:length(OME1(:,1)),1),OME1(2:2:length(OME1(:,1)),1),'.')
l'une plot les éléments dont l'indexe impairs de tes vecteurs, l'autre les éléments pairs.
La courbe que tu désire est la combinaison de la moitié de l'une et de l'autre moitié de l'autre (le code que je t'ai mis précédemment)

Le pourquoi mathématique je ne l'ai pas, a toi de voir.

1:2:length(OME1(:,1)) est un vecteur commençant par 1, allant de 2 en 2, et finissant par le nombre d'élément du vecteur OME1.
ce vecteur est ensuite utilisé comme index dans OMER2 et OME1.
18/09/2014, 17h37
doinier

thanks, lucaswalker,
excellent du moins pour le moment. merci encore