满意答案

usrdu

2013.11.22

采纳率:50%    等级:13

已帮助:9126人

%%傅里叶变换法

clear all

z1=20;

r=4;

N=512;

d=180;

lambda=0.0039;

k=2*pi/lambda;

f1=16*pi/lambda;

R=z1+f1^2/z1;

k0=-N/2:N/2-1;

a=sqrt(lambda*d/N);

[x0,y0]=meshgrid(k0*a);

w=4*sqrt(1+(z1/f1)^2);

Ein=100/w*exp(-(x0.^2+y0.^2)/w^2).*exp(-i*(k*(z1+(x0.^2+y0.^2)/2/R)-atan(z1/f1)));

for m=1:N %%%%%%%%%%%%%%%%%%%圆孔

for n=1:N

if (x0(m,n)^2+y0(m,n)^2>(r)^2)

Ein(m,n) = 0;

end

end

end

E=fft2(Ein.*exp(i*k/2/d*(x0.^2+y0.^2))); %%%%傅里叶变换

Eout=fftshift(abs(E));

x=x0;

y=y0;

U=exp(i*k*d)/i/lambda/d.*exp(i*k/2/d*(x.^2+y.^2)).*Eout;

I=U.*conj(U);

figure,mesh(x,y,I)

figure,plot(x,I)

00分享举报

怎么用matlab赋值圆孔,求一个用MATLAB做圆孔或圆屏衍射的实验,要求程序正确,能出图像,非常感谢...相关推荐

  1. Matlab三元隐函数求极值,matlab用三重循环求一个三元函数的最大值所对应的x1,x2,x3...

    用MATLAB实现for循环 t=2;whileS(t)>Pstrong&&t 求一个MATLAB循环语句表示这个矩阵200分 这样的,i和j是内部虚数变量,避免轻易使用.cle ...

  2. matlab怎么对sinx求导,用matlab程序求y=ln(sinx 1)的导数

    哪位高手会画函数f(x,y)=(1+e^y)*cosx-y*e^y的图像啊?求真相,最好是用Matlab并附上程序! clearall;clc;%画函数f(x,y)=(1+e^y)*cosx-y*e^ ...

  3. ceq c matlab,求一个函数MATLAB的编程,求大神指教。

    function f=myfun(x)f=x(1)+x(2)+81;end %%函数,其中x=x(1)y=x(2)function [c,ceq]=mycon(x)c(1)=100000-x(1)*x ...

  4. matlab 测量矩阵,急求一个测量矩阵采用分块多项式矩阵时怎样引用的代码!!!

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 各位大佬们,在压缩感知中,当测量矩阵采用分块多项式矩阵时,怎样写代码,下面是采用高斯矩阵时的代码 clc;clear %% 1. 时域测试信号生成 K=7 ...

  5. 求拉格朗日多项式matlab,拉格朗日插值多项式积分求圆周率近似Matlab实现

    Lagrange 插值多项式积分求圆周率近似 摘要: 公式1:y1=4/(1+x^2) 公式2:y2=4*sqrt(1-x^2) 分别对公式1.公式2求其拉格朗日插值多项式,再对其求0-1上的定积分来 ...

  6. matlab已知ft求f2t,基于Matlab的信号与系统实验指导2

    基于Matlab 的信号与系统实验指导 实验一 连续时间信号在Matlab 中的表示 一.实验目的 1.学会运用Matlab 表示常用连续时间信号的方法 2.观察并熟悉这些信号的波形和特性 二.实验原 ...

  7. 求一个给定的n阶方阵的鞍点

    一个矩阵元素的"鞍点"是指该位置上的元素值在该行上最大.在该列上最小.本题要求编写程序,求一个给定的n阶方阵的鞍点. 输入格式: 输入第一行给出一个正整数n(1≤n≤6).随后n行 ...

  8. 本题要求编写程序,求一个给定的m×n矩阵各行元素之和。

    本题要求编写程序,求一个给定的m×n矩阵各行元素之和. 输入格式: 输入第一行给出两个正整数m和n(1≤m,n≤6).随后m行,每行给出n个整数,其间 以空格分隔. 输出格式: 每行输出对应矩阵行元素 ...

  9. python求交点坐标_Python - 两圆相交求交点坐标

    Python - 两圆相交求交点坐标 Max.Bai 2016-05-16 Python - 两圆相交求交点坐标 三轴机械臂求坐标问题,其实转化为平面问题就是两圆相交求交点问题,交点算出来就可以用反三 ...

最新文章

  1. 孙立岩 python-basic: 用于学习python基础的课件(161718)
  2. 网络学习笔记—计算机网络基础
  3. mysql特殊编码_Mysql 字符编码
  4. 第八节:EF Core连接MySql和Sqlite数据库
  5. winfrom软件开发汽车测试_ETci — 全自动软件测试调度(持续集成)平台
  6. 详细说说ActionScript中function的call()方法和apply()方法
  7. OpenResty缓存
  8. JavaScript设计模式入坑 1
  9. 搜狗拼音皮肤 php文件,搜狗输入法皮肤制作教程(组图版)
  10. 《Python算法教程_中文版》pdf
  11. 基于web在线餐饮网站的设计与实现——蛋糕甜品店铺(HTML+CSS+JavaScript)
  12. 串操作指令---movs,stos,rep
  13. 人工智能——自然演绎推理
  14. k8s kubeadm 安装
  15. vue 移动端音乐(3) amp;amp;gt;热门歌单推荐部分(webpack-dev-conf.js做后端接口代理+scroll插件)
  16. 509实验室打印机双面打印的方法
  17. 史上最详细springboot vue UEditor整合(包括遇到的各种坑)
  18. cogs 1487. 麻球繁衍(概率dp)
  19. HP大中华区总裁孙振耀撰文谈退休并畅谈人生【全详细】
  20. 贝壳后台开发面经(22 届春招)

热门文章

  1. JavaScript中for循环嵌套编程小案例
  2. LeetCode-704
  3. 关于房屋租赁管理系统设计
  4. 比尔盖茨给青年人的十一点忠告
  5. 黑马旅游网项目详细思路和完整代码整理 -附源码
  6. 手机相机变成QQ摄像头 先锋P80W
  7. xss漏洞扫描器开发随想
  8. linux 命令行 time,Linux 命令行日常交互篇
  9. 带圈数字和常用数学符号怎么写(上标、下标)
  10. ROS与V-rep联合仿真实验