Bonjour
J'ai un programme simple en maple 5 qui me bloque.
le programme cosiste a dessiner la reponse en transmission d'un systeme multicouches periodique donc il faut definir les éléments de la matrice elementaire m11, m12, m21 et m22 ( qui sont en fonction de la frequence f) puis calculers les elements de la matrice produit A11, A12, A21 et A22 puis donner l'_expression de la transmission T= abs(2/( A11+ A12+ A21 + A22)) et en fin fait le plot de T en fonction de la frequence f qui donne un graphe qui resemble à un parabole.
J'ai fait l'execusion ligne par ligne mais je ne trouve pas la courbe.
Le programme n'est pas long et il est le suivant:
>
> MLS :
> restart;
########### initialisation
> N:=16;
> fc:=2.3*10^9;
> c:=3*10^8;
> lam0:=c/fc;
> d1:=6*10^(-3);
> d2:=12*10^(-3);
########### definition des indices n1 et n2
> eps:=1+(25/(0.9^2-f^2))+(100/(11.5^2-f^2));
> mu:=1+(9/(0.902^2-f^2));
> n1:=simplify((-sqrt(eps*mu)), sqrt);
> n2:=1;
########### les elements de la matrice elementaire
> m11:=cos(2*Pi*n1*d1/lam0)*cos(2*Pi*n2*d2/lam0)-(n2/n1)*sin(2*Pi*n1*d1/lam0)*sin(2*Pi*n2*d2/lam0);
> m12:=(-I/n2)*sin(2*Pi*n2*d2/lam0)*cos(2*Pi*n1*d1/lam0)+(-I/n1)*sin(2*Pi*n1*d1/lam0)*cos(2*Pi*n2*d2/lam0);
>
> m21:=(-I*n1)*sin(2*Pi*n1*d1/lam0)*cos(2*Pi*n2*d2/lam0)+(-I*n2)*sin(2*Pi*n2*d2/lam0)*cos(2*Pi*n1*d1/lam0);
> m22:=cos(2*Pi*n1*d1/lam0)*cos(2*Pi*n2*d2/lam0)-(n1/n2)*sin(2*Pi*n1*d1/lam0)*sin(2*Pi*n2*d2/lam0);
###########
> X1:=arccos((m11+m22)/2);
########### les elements de la matrice produit
> A11:=(m11*sin(N*X1))/sin(X1)-sin((N-1)*X1)/sin(X1);
> A12:=m12*sin(N*X1)/sin(X1);
> A21:=m21*sin(N*X1)/sin(X1);
> A22:=(m22*sin(N*X1))/sin(X1)-sin((N-1)*X1)/sin(X1);
########### l’expression de la transmission T
> T:=(abs(2/(A11+A12+A21+A22)))^2;
>
########### plot de T =f(f) la frequence f
plot(T,1*10^9..10^10,thickness =2,axes=BOXED,labels=[f,T],labelfont=[TIMES,ITALIC,16],color=black);
>
-----