±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 2633  |  »Ø¸´: 0

shepherd2014

½ð³æ (СÓÐÃûÆø)

[ÇóÖú] ¸µÀïÒ¶¹âѧMATLAB·ÂÕæµÄÎÊÌâ

À´¸ö´óÉñ°É£¡ÎÒ¸ãÁ˺ܶàÌìʵÔڸ㲻³ö²ÅÉÏÕâÀ´Îʵģ¡£¡£¡

ÒªÔÚMATLABÖÐÄ£Äâ ¹âµÄÑÜÉ䣬ÓÃµÄ DavidGeorgeVoelz_Computational_Fourier_Optics_a_MATLAB_tutorial Õâ±¾ÊéµÄ´úÂëºÍ²åͼ
¸µÀïÒ¶¹âѧMATLAB·ÂÕæµÄÎÊÌâ
¸µÀïÒ¶¹âѧMATLAB·ÂÕæµÄÎÊÌâ-1

ʵ¼ÊʵÑéÖУ¬observation plane´¦·ÅÖõÄÊÇÒ»¸öCCDÏà»ú
Ïà»úµÄÏñËØ³ß´çΪdx2 = 29.6 um
ͼƬµÄ´óСΪ       512 x 512£¬     M = 512
L2 = M*dx2    (1.52 cm)
dx2£¬M  £¬ L2 ÕâЩ²ÎÊý¶¼Êǹ̶¨µÄ£¬²»Äܸıä

´«²¥¾àÀë z = 5 cm£¬¹âµÄ²¨³¤ lambda = 675 nm
ÔËÓ÷òÀź̷ÑÑÜÉ乫ʽ
¸µÀïÒ¶¹âѧMATLAB·ÂÕæµÄÎÊÌâ-2
¸µÀïÒ¶¹âѧMATLAB·ÂÕæµÄÎÊÌâ-3

µÃµ½dx1 = lambda*z/L2  (dx1 Ϊsource plane ´¦µÄÏñËØ³ß´ç£¬Îª2.227um)
L1 = M*dx1      (L1Ϊsource planeµÄ±ß³¤ 1.14 mm£¬source plane´¦Ò²ÊÇÒ»¸ö M x M ´óСµÄͼƬ)

ÎҵķÂÕæÊÇÔÚsource plane ´¦·ÅÖÃÒ»¸ö±ß³¤ D = 50 um µÄͼƬ£¨´Ë´¦Óþؿ״úÌæ£©£¬È»ºóµÃµ½ÑÜÉäͼÑù£¬Õâ¸ö D Ò²Êǹ̶¨²»ÄܱäµÄ
µ«ÊÇÕâÑùÒ»À´£¬D/dx1 = 22.4521 ¸öÏñËØ£¬Í¼Æ¬»ù±¾ÉÏ¿´²»³öϸ½Ú
w = D/2
·ÂÕæ½á¹ûÈçÏÂ
¸µÀïÒ¶¹âѧMATLAB·ÂÕæµÄÎÊÌâ-4
¸µÀïÒ¶¹âѧMATLAB·ÂÕæµÄÎÊÌâ-5

ÎÒµÄÎÊÌâÊÇ dx2 £¬M £¬lambda £¬z £¬D ÕâЩ²ÎÊý¶¨ËÀÁË£¬ÄÇôdx1 £¬L1  ÕâЩ²ÎÊý¸ù¾Ý¹«Ê½Ò²¶¨ËÀÁË£¬µ¼ÖÂÎÒ²åÈëµÄͼƬ¿´²»³öϸ½Ú£¨D/dx1 = 22.4521 ¸öÏñËØ£©
ÎÒÏëµÄÊÇÔÚsource plane ´¦µÄ L1 ÄÜ·ñËõС£¬È»ºóM»¹ÊDz»±äµÄ£¬ÕâÑùdx1¾Í±äСÁË£¬D/dx1¾ÍÄܸü´óЩ£¬²åÈëµÄͼƬҲÄÜ¿´µ½¸ü¶àµÄϸ½Ú£¿
Ôڱ߳¤ÎªDµÄ¾ØÐο×Íâ¶¼ÊÇ0£¬ËùÒÔL1µÄËõСÔÚÕæÊµµÄʵÑéÖÐÎ޹ؽôÒª
¾ÍÏñÏÂÃæÕâÑù
¸µÀïÒ¶¹âѧMATLAB·ÂÕæµÄÎÊÌâ-6


%========================================
clear all;clc;

cm = 1e-2;
mm = 1e-3;
um = 1e-6;
nm = 1e-9;
%%
lambda = 675*nm;
k = 2*pi/lambda;

dx2 = 29.6*um;
M = 512;
z = 5*cm;
%%
L2 = M*dx2;
L1 = lambda*z/dx2;
dx1 = L1/M;

w = 25*um;    % Õâ¸ö w = D/2
%%
x1 = -L1/2:dx1:L1/2-dx1;
y1 = x1;            
[X1,Y1] = meshgrid(x1,y1);            
u1 = rect(X1/(2*w)).*rect(Y1/(2*w));

figure(1)  %irradiance image
imagesc(x1./mm,y1./mm,u1);
xlabel(\'x (mm)\'); ylabel(\'y (mm)\'); title(\'source plane z = 0\');
colormap(\'gray\');
axis square;
axis xy;
%%
[u2,L2]=propFF(u1,L1,lambda,z);
dx2 = L2/M;
x2 = -L2/2:dx2:L2/2-dx2; %obs ords
y2 = x2;
I2 = abs(u2.^2);

figure(2)  %irradiance image
imagesc(x2./cm,y2./cm,nthroot(I2,3));
xlabel(\'x (cm)\'); ylabel(\'y (cm)\'); title(strcat(\'observation plane z = \',num2str(z/cm),\' cm\'));
colormap(\'gray\');
axis square;
axis xy;
%=================================================
function [out] = rect(x)

out=abs(x)<=1/2;
end
%=================================================
function[u2,L2]=propFF(u1,L1,lambda,z)
% propagation - Fraunhofer pattern  
% assumes uniform sampling
% u1 - source plane field
% L1 - source plane side length
% lambda - wavelength
% z - propagation distance
% L2 - observation plane side length
% u2 - observation plane field
%
[M,N] = size(u1); %get input field array size
dx1 = L1/M; %source sample interval
k = 2*pi/lambda; %wavenumber
%
L2 = lambda*z/dx1; %obs sidelength
dx2 = lambda*z/L1; %obs sample interval
x2 = -L2/2:dx2:L2/2-dx2; %obs coords
[X2,Y2] = meshgrid(x2,x2);
%
c = 1/(j*lambda*z)*exp(j*k/(2*z)*(X2.^2+Y2.^2));
u2 = c.*ifftshift(fft2(fftshift(u1)))*dx1^2;
end
%====================================================

[ Last edited by shepherd2014 on 2015-3-19 at 10:31 ]
»Ø¸´´ËÂ¥

» ÊÕ¼±¾ÌûµÄÌÔÌûר¼­ÍƼö

source matlab

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÖÇÄÜ»úÆ÷ÈË

Robot (super robot)

ÎÒÃǶ¼°®Ð¡Ä¾³æ

ÕÒµ½Ò»Ð©Ïà¹ØµÄ¾«»ªÌû×Ó£¬Ï£ÍûÓÐÓÃŶ~

¿ÆÑдÓСľ³æ¿ªÊ¼£¬ÈËÈËΪÎÒ£¬ÎÒΪÈËÈË
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ shepherd2014 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +5 Ò»ÑùYWY 2026-04-06 5/250 2026-04-06 13:08 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 338Çóµ÷¼Á +3 ÎÒÏëÉϰ¶ii 2026-04-05 3/150 2026-04-05 19:59 by nepu_uu
[¿¼ÑÐ] 282Çóµ÷¼Á +7 aaa³µÁ¾ 2026-04-02 11/550 2026-04-05 17:24 by yulian1987
[¿¼ÑÐ] 306·Ö²ÄÁÏÓ뻯¹¤Çóµ÷¼Á +7 Àè°ÉÀ²À²ÄãºÜÓÐà 2026-04-03 7/350 2026-04-05 17:18 by Hdyxbekcb
[¿¼ÑÐ] 313Çóµ÷¼Á +5 º£ÈÕº£ÈÕ 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[¿¼ÑÐ] 283·ÖÇóµ÷¼Á +9 ÊÔÊÔ¿´ß 2026-04-04 9/450 2026-04-05 10:27 by ¹û¶³´óÍõ
[¿¼ÑÐ] 083200 333Çóµ÷¼Á +3 Ê®¶þ£¡£¡ 2026-04-04 3/150 2026-04-05 08:28 by barlinike
[¿¼ÑÐ] 283Çóµ÷¼Á +4 mcbbc 2026-04-03 5/250 2026-04-04 20:51 by imissbao
[¿¼ÑÐ] 348·Ö»·¾³¹¤³Ì¡¤µ÷¼Á +10 ÎâÑå׿24k 2026-04-03 11/550 2026-04-04 14:19 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] Çóµ÷¼Á +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[¿¼ÑÐ] 282Çóµ÷¼Á +5 ºôÎü¶¼ÊǼõ·Ê 2026-03-31 5/250 2026-04-03 12:03 by 1753564080
[¿¼ÑÐ] Ò»Ö¾Ô¸Öйú¿ÆÑ§Ôº´óѧ265Çóµ÷¼Á +9 Ìñµ­ye 2026-03-31 10/500 2026-04-03 11:10 by txp1986
[¿¼ÑÐ] ±§Ç¸ +5 ÌïºéÓÐ 2026-03-30 5/250 2026-04-03 10:24 by linyelide
[¿¼ÑÐ] 312 »¯¹¤»òÖÆÒ©µ÷¼Á +8 ССī123 2026-04-02 9/450 2026-04-03 09:12 by zhouxiaoyu
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÄÏʦ·¶´óѧ-22408¼ÆËã»ú-292·Ö-Çó»ªÄÏʦ·¶´óѧµ÷¼Á +4 °®¶ÁÊéµÄСöùÓã 2026-04-02 4/200 2026-04-02 18:35 by Çóµ÷¼Ázz
[¿¼ÑÐ] 354Çóµ÷¼Á +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[¿¼ÑÐ] Ò»Ö¾Ô¸±±½»´ó²ÄÁϹ¤³Ì£¬×Ü·Ö358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by ÉÐË®¸óÖ÷
[¿¼²©] 26ÄêÉ격 +3 staryer 2026-03-30 4/200 2026-04-01 23:21 by ai4pharm
[¿¼ÑÐ] °²È«¹¤³Ì 285 Çóµ÷¼Á +3 Xinyu56 2026-04-01 4/200 2026-04-01 21:50 by ¾²¾²¾²¾²¾²¾²¾²¾
[¿¼ÑÐ] Çó0861½»Í¨ÔËÊäר˶or²ÄÁÏר˶µ÷¼Á +4 ÀÕ²¼ÀÊ@ 2026-03-31 4/200 2026-04-01 09:54 by Ò»Ö»ºÃ¹û×Ó?
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û