%光學系統的傳遞函數程序
close?all
clear?all							%清除workspace中的變量
clc
addpath?D:\MATLAB\Myfunction
Na?=?127;
wl?=?5e-7;
for?ii?=?2:5
k?=?2*pi./wl;
zr?=?ii*zernike（3Na）;
pupil=Pupil（Na）;?			
pupil?=?pupil.*exp（1i*2*pi*zr）;	?		%瞳孔函數

pupil_1=fftshift（fft2（pupil））;
%?pupil_2=fftshift（pupil_1）;
psf=abs（pupil_1）;				%點擴散函數psf
mtf_1=fft2（psf）;
mtf_2=fftshift（mtf_1）;
mtf_2=abs（mtf_2）;
mtf_2=mtf_2/max（max（mtf_2））;	%調制傳遞函數mtf
title（‘光學傳遞函數‘）
x=1:1:61;
plot（xmtf_2（301301:361））		%繪圖
hold?on
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????547??2020-11-16?20:23??MTF.m