13.3
close all;home;clc;clear all;
f1=6.65e9-20e6;
f2=6.65e9+20e6;
c=3e8;
a=34.85e-3;
k1=2*pi*f1/c;
k2=2*pi*f2/c;
kc=pi/a;
Lg1=2*pi/sqrt((k1^2)-(kc^2));
Lg2=2*pi/sqrt((k2^2)-(kc^2));
arfa=(Lg1+Lg2)/(3*pi*(Lg1-Lg2))
%=======================================================================
MS1=1.0638;
M12=-0.1513;
M23=0.5328;
M34=-0.5576;
M45=0.4282;
M56=-0.4892;
M14=-0.8887;
M36=-0.7572;
M11=0.0399;
M22=-0.7214;
M33=0.3023;
M44=0.0732;
M55=-0.5516;
M66=0.0399;
B11=M11/arfa;
B22=M22/arfa;
B33=M33/arfa;
B44=M44/arfa;
B55=M55/arfa;
B66=M66/arfa;
KS1=MS1/sqrt(arfa);
BS1=KS1-(1/KS1)
B12=M12/arfa
K23=M23/arfa;
B23=K23-(1/K23)
B34=M34/arfa
K45=M45/arfa;
B45=K45-(1/K45)
B56=M56/arfa
K14=M14/arfa;
B14=K14-(1/K14)
K36=M36/arfa;
B36=K36-(1/K36)
theta11=((3*pi/2)+(1/2)*(acot(BS1/2)-asin(B12-B11)))*180/pi
theta21=((3*pi/2)+(1/2)*(acot(B14/2)-asin(B12-B11)))*180/pi
theta12=((3*pi/2)+(1/2)*(-asin(B12-B22)))*180/pi
theta22=((3*pi/2)+(1/2)*(acot(B23/2)-asin(B12-B22)))*180/pi
theta13=((3*pi/2)+(1/2)*(acot(B23/2)-asin(B34-B33)))*180/pi
theta23=((3*pi/2)+(1/2)*(acot(B36/2)-asin(B34-B33)))*180/pi
arfa =
20.4998
BS1 =
-4.0212
B12 =
-0.0074
B23 =
-38.4496
B34 =
-0.0272
B45 =
-47.8535
B56 =
-0.0239
B14 =
23.0238
B36 =
27.0362
theta11 =
257.0451
theta21 =
272.7495
theta12 =
269.2032
theta22 =
267.7144
theta13 =
269.7132
theta23 =
273.3174
Comments