BÀI TẬP LẬP TRÌNH - BÀI TOÁN THUẬN TRANG 108

BÀI TOÁN THUẬNtrang 108Một buồng cộng hưởng quang học bao gồm một gương cầu lồi có độ tụ P1 = 0.25 và một gương phẳng đặt cách nhau 1m, hai mặt phản xạ quay vào nhau.. Xác định tính bền

BÀI TOÁN THUẬN(trang 108)

Một buồng cộng hưởng quang học bao gồm một gương cầu lồi có độ tụ P1 = 0.25 và một gương phẳng đặt cách nhau 1m, hai mặt phản xạ quay vào nhau Xác định tính bền hay không bền của

hệ cộng hưởng và các yếu tố của nó

%%%%

% LE PHUC QUY

%BAI TOAN THUAN (trang 108)

% BUOC 1: NHAP VAO CAC SO LIEU DA BIET

L=input('Nhap vao khoang cach giua guong phan xa cau va guong phang (m): ');

P1=input('Nhap vao do tu cua guong phan xa cau loi (Phai la so am): ');

while P1>0

P1=input('Phai la so am Vui long nhap lai: ');

end

P2=0; % Do tu cua guong phang bang 0

n=input('Nhap vao chiet suat moi truong giua hai guong: ');

lamda=input('Nhap vao buoc song anh sang khao sat (nm): ');

lamda=lamda*10^-9; % Doi don vi tu nm sang m

%

% BUOC 2: VIET BIEU THUC CAC MA TRAN TRUYEN QUA VA MA TRAN PHAN XA

M1=[1 0;-P2 1]; % Ma tran phan xa cua guong phang do tu P2 = 0

M2=[1 L/n;0 1]; % Ma tran truyen qua chieu dai L

M3=[1 0;-P1 1]; % Ma tran phan xa cua guong cau do tu P1

M4=[1 L/n;0 1]; % Ma tran truyen qua chieu dai L

M=M4*M3*M2*M1; % Ma tran truyen tia cua he cong huong

%

% BUOC 3: XAC DINH TINH BEN CUA HE CONG HUONG VA CAC YEU TO CUA HE A=M(1,1); % He so A la phan tu dong 1 cot 1 cua ma tran M

B=M(1,2); % He so B la phan tu dong 1 cot 2 cua ma tran M

C=M(2,1); % He so C la phan tu dong 2 cot 1 cua ma tran M

D=M(2,2); % He so D la phan tu dong 2 cot 2 cua ma tran M

m=(A+D)/2;

if (m>-1)&(m<1)

disp('HE CONG HUONG BEN')

j=(-1)^(1/2);

costheta=(A+D)/2;

sintheta=(1-((A+D)/2)^2)^(1/2);

disp('Tri rieng cua he la:')

lamda1=costheta+j*sintheta

disp('Ty so vector rieng:')

tyso1=(lamda1-D)/C

tyso2=(lamda1-A)/B

disp('Ban kinh cong phuc cua he la:')

q=(A-D)/(2*C)+j*sintheta/C

disp('Ban kinh cong cua chum tia la:')

R=2*B/(D-A)

disp('Do phan ky cua mat song:')

R1=(D-A)/(2*B)

disp('Ban kinh vet cua chum tia:')

omega=(lamda*B/(pi*sintheta))^(1/2)

disp('Vi tri co chum (ben trai so voi mat phang quy chieu):')

z=(A-D)/(2*C)

disp('Ban kinh co chum:')

omega0=(-lamda*sintheta/(pi*C))^(1/2)

disp('Thong so chum dong tieu:')

z0=-sintheta/C

disp('Nua goc truong xa:')

farfield=omega0/z0

else

disp('HE CONG HUONG KHONG BEN')

sinht=(((A+D)/2)^2-1)^(1/2);

if m>1

disp('NHANH DUONG')

cosht=(A+D)/2;

disp('Tri rieng cua he la:')

lamda1=cosht+sinht

disp('Ty so vector rieng:')

tyso1=(lamda1-D)/C

tyso2=(lamda1-A)/B

disp('Ban kinh cong cua he:')

R=(A-D)/(2*C)+sinht/C

else

disp('NHANH AM')

cosht=-(A+D)/2;

disp('Tri rieng cua he la:')

lamda1=-(cosht+sinht)

disp('Ty so vector rieng:')

tyso1=(lamda1-D)/C

tyso2=(lamda1-A)/B

disp('Ban kinh cong cua he:')

R=(A-D)/(2*C)-sinht/C

end

end

BÀI TOÁN NGHỊCH

Nhập vào các dữ kiện: Hệ cộng hưởng không bền nhánh dương; Khoảng cách L (L = 1m) giữa hai gương phản xạ; Chiết suất môi trường n (n = 1) giữa hai gương và Bán kính cong R (R = 3m) Xác định bán kính r1 của gương cầu phản xạ lồi

% LE PHUC QUY

%BAI TOAN NGHICH (trang 108)

% BUOC 1: NHAP VAO CAC SO LIEU DA BIET

L=input('Nhap vao khoang cach giua guong phan xa cau va guong phang (m): ');

R=input('Nhap vao ban kinh cong cua he (m): ');

P2=0; % Do tu cua guong phang bang 0

n=input('Nhap vao chiet suat moi truong giua hai guong: ');

%

% BUOC 2: VIET BIEU THUC CAC MA TRAN TRUYEN QUA VA MA TRAN PHAN XA

M1=[1 0;-P2 1]; % Ma tran phan xa cua guong phang do tu P2 = 0 M2=[1 L/n;0 1]; % Ma tran truyen qua chieu dai L

M3=[1 0;-P1 1]; % Ma tran phan xa cua guong cau do tu P1 M4=[1 L/n;0 1]; % Ma tran truyen qua chieu dai L

M=M4*M3*M2*M1; % Ma tran truyen tia cua he cong huong

%

% BUOC 3: XAC DINH BAN KINH CUA GUONG PHAN XA LOI A=M(1,1); % He so A la phan tu dong 1 cot 1 cua ma tran M B=M(1,2); % He so B la phan tu dong 1 cot 2 cua ma tran M C=M(2,1); % He so C la phan tu dong 2 cot 1 cua ma tran M D=M(2,2); % He so D la phan tu dong 2 cot 2 cua ma tran M disp('HE CONG HUONG KHONG BEN NHANH DUONG')

sinht=(((A+D)/2)^2-1)^(1/2);

cosht=(A+D)/2;

lamda1=cosht+sinht;

R0=R;

R=(lamda1-D)/C;

pt=(lamda1-D)/C-R0;

disp('Do tu cua guong phan xa loi la:')

P1=solve(pt);

P1=double(P1) % Chuyen ket qua sang so thap phan

disp('Ban kinh cua guong phan xa loi la (m):')

r1=-2*n/P1