Discussion :

[doinier]

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

[Invité]

Bonsoir,

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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.

[doinier]

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

[doinier]

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 courbe

voilà mon code.

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

Sélectionner tout - Visualiser dans une fenêtre à part

length(eigen)

me donne 8

s'il vous plaît help me.

[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.

[doinier]

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.

[doinier]

bonjour lucaswalker,

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[lucaswalker]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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])

[doinier]

Merci lucaswalker, ça fonctionne .

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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),'.')

[lucaswalker]

ok,

essaye ces deux lignes:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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.

[doinier]

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