Proposed idea to 修文

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