24小时热门版块排行榜    

查看: 3683  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wangzizii

银虫 (初入文坛)

[求助] 请教:关于matlab的菲涅尔衍射FFT计算仿真

菲涅尔衍射可看成 输入函数与相位因子乘积 的傅里叶变换,应用FFT方法可以快速计算输入函数的菲涅尔衍射。
但是在编程过程中遇到采样频率和衍射场的空间频率不知道怎么确定的问题,例如,在以下代码中,保持衍射狭缝宽度不变,改变采样数或者改变观察屏的范围都会导致衍射图样的变化,这是不正确的。
在这里想请教各位大牛,应该如何完善下面的菲涅尔衍射计算的程序才能得到正确的计算结果?
非常感谢

% 本程序用于计算宽度为0.2mm的狭缝的菲涅尔衍射
clear all;close all;
lamda = 0.635e-3;   % 波长
k = 2*pi/lamda;     % 波数
d = 200;            % 传播距离
N = 128;           % 采样数
x = linspace(-1,1,N);   % 衍射屏和观察屏
slit = zeros(1,N);      % 狭缝
slit(abs(x)<=0.1) = 1;  % 宽度为0.2mm的狭缝
coe = exp(1i*k/2/d*x.^2);   % 菲涅尔积分项
Efresnel = fftshift(fft(slit.*coe))*(x(2)-x(1));   % 菲涅尔衍射的傅里叶变换
Ifresnel = abs(Efresnel).^2;
plot(x,Ifresnel);  % 画图
grid on
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangzizii

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by ygjianqiang at 2012-05-30 10:18:02
请问你采用的是菲涅尔衍射的哪个公式?卷积法吗?
没太理解你的意思,如果是想采用卷积法做菲涅尔衍射的话,贴个程序供你参考吧:
clc
clear all
close all
r=512,c=r; a=zeros(r,c);
a(257-1:257+1,257-128: ...

学习了,谢谢。
3楼2012-06-08 22:43:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangzizii 的主题更新
信息提示
请填处理意见