光学4F系统仿真代码

%这个程序是模拟4f系统滤波处理，在滤波面加菲涅尔滤波片
clc;
clear all;

%模拟4f系统对图像进行滤波处理
V=imread（‘D:\3.bmp‘）;     %读入字母V
V=rgb2gray（V）;
V=im2double（V）;
subplot（121）;imshow（V）;title（‘字母V‘）;
V1=fftshift（fft2（V））;     % 将图像频谱低频部分移到中央，高频部分移到四周
subplot（122）;imshow（abs（V1）*16/max（max（abs（V1））））;title（‘字母V经物面到达频谱面后的频谱图‘）; %字母V经物面到达频谱面后的频谱图


r=16;
a=4*sqrt（r）;
b=0;
%生成菲涅尔滤波片
F=ones（256256）;
for k=1:30
    p=6*sqrt（k*r）;