Discussion :

[Benmu]

Salut à tous et à toutes.

Je suis en train de travailler (sur matlab 20008a) avec une fonction définie par morceaux et j'aimerai interpoler cette fonction sur tout l'intervalle d'étude. Cette fonctions est à variable complexe et j'aimerai trouver une fonction définissant la phase et une définissant l'amplitude. Je cherche plutôt une solution graphique mais si je peux en tirer une solution analytique cela me va très bien aussi!!!

Voici mon code:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
Io=1;
D=0.15*10^(-6);
K=0.0002;
l=50*10^(-6);
f1=0.001:0.001:13;
z=0;
f2=13:1*10^5;
Tacthin=(Io.*D./(4.*pi.*K.*l.*f1)).*exp(-i*pi/2);
Tacthick=(Io.*sqrt(D)./(2.*K.*sqrt(2.*pi.*f2))).*exp(sqrt(pi.*f2./D).*z).*exp(i.*(sqrt(pi.*f2./D).*z-pi/4));
f=0.001:10^5;
figure(1)
semilogx(f1,log(abs(Tacthin)),'-r',f2,log(abs(Tacthick)),'-r',f2,unwrap(angle(Tacthick)),'-b',f1,unwrap(angle(Tacthin)),'-b');

Je chercherai donc une fonction qui sur l'intervalle f1 soit le plus proche possible de abs(Tacthin) et sur f2 de abs(Tacthick) (et de même avec unwrap(angle(Tacthin)) sur f1 et unwrap(angle(Tacthick)) sur f2).
EDIT: J'ai regardé les fonctions spline ou convhull mais ce n'est pas vraiment ce dont je recherche...
Avez vous une idée?
Merci d'avance!

[Benmu]

Dans un second temps j'ai essayé de définir cette fonction avec une boucle mais ca ne donne pas vraiment le résultat escompté..

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
Io=1;
D=0.15*10^(-6);
K=0.0002;
l=50*10^(-6);
z=0;
T=0;
hold on
for f=0.01:10^3
    if f<=13
        Tac=(Io.*D./(4.*pi.*K.*l.*f)).*exp(-i*pi/2);
    else
        Tac=(Io.*sqrt(D)./(2.*K.*sqrt(2.*pi.*f))).*exp(sqrt(pi.*f./D).*z).*exp(i.*(sqrt(pi.*f./D).*z-pi/4));
    end
figure(1)
semilogx(f,log(abs(Tac)),'-r',f,unwrap(angle(Tac)),'-b')
end
hold off

[shaiHulud]

Bonjour,

interp2() permet d'interpoler une fonction à valeur dans le plan. Tu peux ensuite passer en coord. polaire pour avoir module et phase.