怎么用matlab赋值圆孔,求一个用MATLAB做圆孔或圆屏衍射的实验,要求程序正确,能出图像,非常感谢...
满意答案
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做圆孔或圆屏衍射的实验,要求程序正确,能出图像,非常感谢...相关推荐
- Matlab三元隐函数求极值,matlab用三重循环求一个三元函数的最大值所对应的x1,x2,x3...
用MATLAB实现for循环 t=2;whileS(t)>Pstrong&&t 求一个MATLAB循环语句表示这个矩阵200分 这样的,i和j是内部虚数变量,避免轻易使用.cle ...
- 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^ ...
- 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 ...
- matlab 测量矩阵,急求一个测量矩阵采用分块多项式矩阵时怎样引用的代码!!!
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 各位大佬们,在压缩感知中,当测量矩阵采用分块多项式矩阵时,怎样写代码,下面是采用高斯矩阵时的代码 clc;clear %% 1. 时域测试信号生成 K=7 ...
- 求拉格朗日多项式matlab,拉格朗日插值多项式积分求圆周率近似Matlab实现
Lagrange 插值多项式积分求圆周率近似 摘要: 公式1:y1=4/(1+x^2) 公式2:y2=4*sqrt(1-x^2) 分别对公式1.公式2求其拉格朗日插值多项式,再对其求0-1上的定积分来 ...
- matlab已知ft求f2t,基于Matlab的信号与系统实验指导2
基于Matlab 的信号与系统实验指导 实验一 连续时间信号在Matlab 中的表示 一.实验目的 1.学会运用Matlab 表示常用连续时间信号的方法 2.观察并熟悉这些信号的波形和特性 二.实验原 ...
- 求一个给定的n阶方阵的鞍点
一个矩阵元素的"鞍点"是指该位置上的元素值在该行上最大.在该列上最小.本题要求编写程序,求一个给定的n阶方阵的鞍点. 输入格式: 输入第一行给出一个正整数n(1≤n≤6).随后n行 ...
- 本题要求编写程序,求一个给定的m×n矩阵各行元素之和。
本题要求编写程序,求一个给定的m×n矩阵各行元素之和. 输入格式: 输入第一行给出两个正整数m和n(1≤m,n≤6).随后m行,每行给出n个整数,其间 以空格分隔. 输出格式: 每行输出对应矩阵行元素 ...
- python求交点坐标_Python - 两圆相交求交点坐标
Python - 两圆相交求交点坐标 Max.Bai 2016-05-16 Python - 两圆相交求交点坐标 三轴机械臂求坐标问题,其实转化为平面问题就是两圆相交求交点问题,交点算出来就可以用反三 ...
最新文章
- 孙立岩 python-basic: 用于学习python基础的课件(161718)
- 网络学习笔记—计算机网络基础
- mysql特殊编码_Mysql 字符编码
- 第八节:EF Core连接MySql和Sqlite数据库
- winfrom软件开发汽车测试_ETci — 全自动软件测试调度(持续集成)平台
- 详细说说ActionScript中function的call()方法和apply()方法
- OpenResty缓存
- JavaScript设计模式入坑 1
- 搜狗拼音皮肤 php文件,搜狗输入法皮肤制作教程(组图版)
- 《Python算法教程_中文版》pdf
- 基于web在线餐饮网站的设计与实现——蛋糕甜品店铺(HTML+CSS+JavaScript)
- 串操作指令---movs,stos,rep
- 人工智能——自然演绎推理
- k8s kubeadm 安装
- vue 移动端音乐(3) amp;amp;gt;热门歌单推荐部分(webpack-dev-conf.js做后端接口代理+scroll插件)
- 509实验室打印机双面打印的方法
- 史上最详细springboot vue UEditor整合(包括遇到的各种坑)
- cogs 1487. 麻球繁衍(概率dp)
- HP大中华区总裁孙振耀撰文谈退休并畅谈人生【全详细】
- 贝壳后台开发面经(22 届春招)