Sau đây là đoạn chương trình của bài toán nghịch:

% ********************************************************************************% Truong Dai hoc Khoa hoc Tu nhien % Truong Dai hoc Khoa hoc Tu nhien

% Bo mon Vat ly Ung dung

% Giang vien huong dan: TS. Le Vu Tuan Hung % Ho ten hoc vien: Phan Trung Vinh

% ********************************************************************************% %

% BAI LAP TRINH 2 - BAI TAP 2 - PROBLEM 11 - BAI TOAN NGHICH - TRUONG HOP 1clc clc

clear all

syms L r1 r2 % Khai bao 3 bien can tim la L, r1, r2 %

% ********************************************************************************% BUOC 1: NHAP VAO CAC SO LIEU DA BIET % BUOC 1: NHAP VAO CAC SO LIEU DA BIET

F1=input('Nhap vao tieu diem thu nhat cua thau kinh cach RP1 (cm): '); F2=input('Nhap vao tieu diem thu hai cua thau kinh cach RP2 (cm): '); H1=input('Nhap vao diem chinh thu nhat cua thau kinh cach RP1 (cm): '); n=input('Nhap vao chiet suat cua tam thuy tinh: ');

n1=input('Nhap vao chiet suat moi truong ben trai cua tam thuy tinh (khong khi): '); n2=input('Nhap vao chiet suat moi truong ben phai cua tam thuy tinh: ');

%

% ********************************************************************************% BUOC 2: VIET BIEU THUC CAC MA TRAN TRUYEN QUA VA MA TRAN KHUC XA % BUOC 2: VIET BIEU THUC CAC MA TRAN TRUYEN QUA VA MA TRAN KHUC XA M1=[1 0;-(n-n1)/r1 1]; % Ma tran khuc xa cua mat cau loi ban kinh r1

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

M3=[1 0;-(n2-n)/r2 1]; % Ma tran khuc xa cua mat cau lom ban kinh r2 M=M3*M2*M1; % Ma tran truyen tia cua ca he quang hoc

%

% ********************************************************************************

% BUOC 3: TIM CAC TIEU CU, CAC TIEU DIEM, CAC DIEM CHINH, CAC DIEM NUTA=M(1,1); % He so A la phan tu dong 1 cot 1 cua ma tran M 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

pt1=n1*D/C-F1; pt2=-n2*A/C-F2; pt3=n1*(D-1)/C-H1;

[L,r1,r2]=solve(pt1,pt2,pt3); disp('****************')

disp('Ban kinh cua mat cong thu nhat (cm):') r1=double(r1)

disp('****************')

disp('Ban kinh cua mat cong thu hai (cm):') r2=double(r2)

disp('****************')

disp('Be day cua thau kinh (cm):'); L=double(L)

Kết quả thu được từ đoạn chương trình trên:

Trường hợp 2: Cho biết vị trí các tiêu điểm F1, F2 (F1 = -18.3333, F2 = 18.6667), vị trí điểm nút thứ nhất L1 (L1 = 5), chiết suất môi trường không khí bên trái thấu kính (n1 = 1), điểm nút thứ nhất L1 (L1 = 5), chiết suất môi trường không khí bên trái thấu kính (n1 = 1), chiết suất thủy tinh (n = 1.5) và chiết suất chất lỏng bên phải thấu kính (n2 = 1.4). Yêu cầu xác định bán kính r1, r2 của từng mặt cầu và bề dày thấu kính.

Bài giải

Sau đây là đoạn chương trình của bài toán nghịch:

% ********************************************************************************% Truong Dai hoc Khoa hoc Tu nhien % Truong Dai hoc Khoa hoc Tu nhien

% Bo mon Vat ly Ung dung

% Giang vien huong dan: TS. Le Vu Tuan Hung % Ho ten hoc vien: Phan Trung Vinh

% ********************************************************************************% %

% BAI LAP TRINH 2 - BAI TAP 2 - PROBLEM 11 - BAI TOAN NGHICH - TRUONG HOP 2clc clc

clear all

syms L r1 r2 % Khai bao 3 bien can tim la L, r1, r2 %

% ********************************************************************************% BUOC 1: NHAP VAO CAC SO LIEU DA BIET % BUOC 1: NHAP VAO CAC SO LIEU DA BIET

F1=input('Nhap vao tieu diem thu nhat cua thau kinh cach RP1 (cm): '); F2=input('Nhap vao tieu diem thu hai cua thau kinh cach RP2 (cm): '); L1=input('Nhap vao diem nut thu nhat cua thau kinh cach RP1 (cm): '); n=input('Nhap vao chiet suat cua tam thuy tinh: ');

n1=input('Nhap vao chiet suat moi truong ben trai cua tam thuy tinh (khong khi): '); n2=input('Nhap vao chiet suat moi truong ben phai cua tam thuy tinh: ');

%

% ********************************************************************************% BUOC 2: VIET BIEU THUC CAC MA TRAN TRUYEN QUA VA MA TRAN KHUC XA % BUOC 2: VIET BIEU THUC CAC MA TRAN TRUYEN QUA VA MA TRAN KHUC XA M1=[1 0;-(n-n1)/r1 1]; % Ma tran khuc xa cua mat cau loi ban kinh r1

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

M3=[1 0;-(n2-n)/r2 1]; % Ma tran khuc xa cua mat cau lom ban kinh r2 M=M3*M2*M1; % Ma tran truyen tia cua ca he quang hoc

%

% ********************************************************************************

% BUOC 3: TIM CAC TIEU CU, CAC TIEU DIEM, CAC DIEM CHINH, CAC DIEM NUTA=M(1,1); % He so A la phan tu dong 1 cot 1 cua ma tran M 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

pt1=n1*D/C-F1; pt2=-n2*A/C-F2; pt3=(D*n1-n2)/C-L1;

[L,r1,r2]=solve(pt1,pt2,pt3); disp('****************')

disp('Ban kinh cua mat cong thu nhat (cm):') r1=double(r1)

disp('****************')

disp('Ban kinh cua mat cong thu hai (cm):') r2=double(r2)

disp('****************')

disp('Be day cua thau kinh (cm):'); L=double(L)