Discussion :

[Renaud33140]

Bonjour,

Je débute sur MATLAB et je n'ai jamais vraiment programmé non plus, d'où mes soucis à trouver seul les solutions.
Je dois modélisé une fonction mais comme un des paramètre est variable, j'ai décider de compiler ces courbes sur un même graphe avec différentes valeurs du paramètre. Mon code est le suivant :

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
clear all
 
% w= logspace(-1,4);
w= 0:1:10^4;
C=100*10^-9;
R9=33000;
R11=3300;
RVmax=50000;
w2=1/(R11*C)
nCourbes=10;
 
 
for i=1:nCourbes
    RV=(50/nCourbes)*i*10^3
    w1=1./((C*R11*(R9+RV))./(R9+R11+RV))
    H0=(R9+R11+RV)./(R9+RV)
    H=H0*(1+1i*w/w1)./(1+1i*w/w2)
 
    figure(1)
    plot(w, 20*log10(abs(H)))
    hold all
 
     figure (2)
    plot(w, atan(w/w1)-atan(w/w2))
    hold all
 
end

J'ai donc mes courbes, différentiables grâce aux couleurs mais je n'ai aucun moyen de les identifier. J'ai essayé de mettre une légende grace à la commande "LEGEND" mais ça ne me permet pas de mettre la valeur de RV (mon paramètre variable") dans la légende de manière automatique.

Merci d'avance !

R.ROEDERER

[Linterne]

Bonjour,

Est-ce qu'une solution comme celle ci est ce que tu recherches ?

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
clear all
h = figure(1);
h2 = figure(2);
 
% w= logspace(-1,4);
w= 0:1:10^4;
C=100*10^-9;
R9=33000;
R11=3300;
RVmax=50000;
w2=1/(R11*C)
nCourbes=10;
 
for i=1:nCourbes
 
if i==1
    figure(1);
    clf;
    figure(2);
    clf;
end
 
RV=(50/nCourbes)*i*10^3;
w1=1./((C*R11*(R9+RV))./(R9+R11+RV));
H0=(R9+R11+RV)./(R9+RV);
H=H0*(1+1i*w/w1)./(1+1i*w/w2);
 
figure(h);
title(['RV= ' num2str(RV)]);
plot(w, 20*log10(abs(H)))
hold all
S(i) = getframe(h);
 
figure(h2);
title(['RV= ' num2str(RV)]);
plot(w, atan(w/w1)-atan(w/w2))
hold all
S2(i) = getframe(h2);
 
end

[Jerome Briot]

Comme ceci ?

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
clear all
 
w =  0:1:10^4;
C = 100*10^-9;
R9 = 33000;
R11 = 3300;
RVmax = 50000;
w2 = 1/(R11*C);
nCourbes = 10;
 
figure(1)
ax(1) = axes;
hold on
figure(2)
ax(2) = axes;
hold on
 
for i = 1:nCourbes
 
    RV(i) = (50/nCourbes)*i*10^3;
    w1 = 1./((C*R11*(R9+RV(i)))./(R9+R11+RV(i)));
    H0 = (R9+R11+RV(i))./(R9+RV(i));
    H = H0*(1+1i*w/w1)./(1+1i*w/w2);
 
    p(i,1) = plot(ax(1), w, 20*log10(abs(H)));
    p(i,2) = plot(ax(2), w, atan(w/w1)-atan(w/w2));    
 
end
 
legend(ax(1), p(:,1), num2str(RV(:)));
legend(ax(2), p(:,2), num2str(RV(:)));

[Renaud33140]

C'est exactement ça ! Merci beaucoup !!

[Renaud33140]

Jérôme, si je remplace dans ton code "plot(ax[...])" par loglog afin d'avoir une échelle log, ça ne fonctionne pas. Peux-tu me conseiller s'il te plait ?

Merci d'avance.

[Jerome Briot]

Comme ceci :

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
clear all
 
w =  0:1:10^4;
C = 100*10^-9;
R9 = 33000;
R11 = 3300;
RVmax = 50000;
w2 = 1/(R11*C);
nCourbes = 10;
 
figure(1)
clf
ax(1) = axes;
 
figure(2)
clf
ax(2) = axes;
 
for i = 1:nCourbes
 
    RV(i) = (50/nCourbes)*i*10^3;
    w1 = 1./((C*R11*(R9+RV(i)))./(R9+R11+RV(i)));
    H0 = (R9+R11+RV(i))./(R9+RV(i));
    H = H0*(1+1i*w/w1)./(1+1i*w/w2);
 
    p(i,1) = loglog(w, 20*log10(abs(H)), 'parent', ax(1));    
    p(i,2) = loglog(w, atan(w/w1)-atan(w/w2), 'parent', ax(2));
 
    if i==1
        hold(ax(1),'on');
        hold(ax(2),'on');
    end
 
end
 
legend(ax(1), p(:,1), num2str(RV(:)));
legend(ax(2), p(:,2), num2str(RV(:)));

[Linterne]

Ou comme ça ?

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
clear all
 
C = 100*10^-9;
R9 = 33000;
R11 = 3300;
RVmax = 50000;
w2 = 1/(R11*C);
nCourbes = 10;
 
Hsave=[];
Atan=[];
w = logspace(-1,4);
 
for i = 1:nCourbes
 
    RV(i) = (50/nCourbes)*i*10^3;
    w1 = 1./((C*R11*(R9+RV(i)))./(R9+R11+RV(i)));
    H0 = (R9+R11+RV(i))./(R9+RV(i));
    H = H0*(1+1i*w/w1)./(1+1i*w/w2);
 
    Hsave=[Hsave;20*log10(abs(H))];
    Atan=[Atan;atan(w/w1)-atan(w/w2)];       
 
end
 
figure(3)
clf;
loglog(w,Hsave)
grid on
legend(num2str(RV(:)));
 
figure(4)
clf;
loglog(w,Atan)
grid on
legend(num2str(RV(:)));