Bonjour, monsieur, S.V.P est ce que vous pouvez m'aider et m'expliquer comment on peut calculer la dérivée d’un nombre complexe dépend de la *pul *ce coefficient s’appelle : tsp la partie réelle de ce cœfficient est x et la partie imaginaire est y, x et y dépend de la (pul) moi j’ai fait un programme qui calcule x et y point par point et je voudrai ensuite calculer la dérivée de x et de y par rapport à pul ces dérivées sont nommés xn et yn et enfin je voudrai calculer ce deux termes : Mx= x .* xn et My= y./ yn , est ce que vous pouvez m’aider SVP. Dans l'espoir de recevoir votre réponse je vous prie d'agrée monsieur mes sincères remerciements.
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
 
lmin=1.7*10^15;
lmax=6.8*10^15;
c=3*(10^8);
pul=lmin:0.1:lmax;
A=length(pul);
p=100;
lam0=0.5;
N(1)=1;
for jjj=1:A
    lam(jjj)=(2*pi*c*10^6)/pul(jjj);
 
for m=2:p+2  
      theta(m)= asin((N(m-1)/N(m))*sin(theta(m-1)));
    end
for m=2:p+2 
    phi(m)=((pi/2)*lam0*cos(theta(m)))/lam(jjj);
end
for m=2:p+2
    NumRp(m)=(N(m-1)*cos(theta(m))-N(m)*cos(theta(m-1)));
    NumTp(m)=2*N(m-1)*cos(theta(m-1));
    DenomRp(m)=(N(m-1)*cos(theta(m))+N(m)*cos(theta(m-1)));
    DenomTp(m)=(N(m-1)*cos(theta(m))+N(m)*cos(theta(m-1)));
end
for m=2:p+2
    rp(m)=(NumRp(m)/DenomRp(m));
    tp(m)=(NumTp(m)/DenomTp(m));
end
Sp=[1 0;0 1];
for m=2:p+2
CRp=[exp(-i*phi(m))/tp(m) rp(m)*exp(i*phi(m))/tp(m); rp(m)*exp(-i*phi(m))/tp(m) exp(i*phi(m))/tp(m)];
Sp=Sp*CRp;
end
 
tsp(jjj)=n(p+2)/Sp(1)
x(jjj)=real(tsp(jjj))
y(jjj)=imag(tsp(jjj));
xn(jjj)=diff(tsp(jjj),pul(jjj))% ce ça la derivée de x par rapport à la pul
yn(jjj)=diff(y(jjj),pul(jjj));
end