求助大神帮忙写一个matlab程序，用来仿真单缝夫琅禾费衍射的问题，
用 surf画更好看一点（个人偏好）
%单缝的夫琅禾费衍射频谱图x=-2:.05:0.2;y=x;[X,Y]=meshgrid(x,y);F = abs(fftshift(fft2(rect(x))));F1=repmat(F,length(x),1);%h = mesh(X,Y,F1);h = surf(X,Y,F1)title('单缝夫琅禾费衍射图样三维强度绘图')xlabel('x'), ylabel('y'), zlabel('z');set(h,'EdgeColor','none');shading interp;
a=1e-3;D=1;
ym=3*lam*D/a;
ny=51;
ys=linspace(-ym,ym,ny);
np=51;
yp=linspace(0,a,np);
for i=1:ny
sinphi=ys(i)/D;
alpha=2*pi*yp*sinphi/lam;
sumcos=sum(cos(alpha));
sumsin=sum(sin(alpha));
B(i,:)=(sumcos^2+sumsin^2)/np^2;
end
N=255;
Br=(B/max(B))*N;
subplot(1,2,1)
image(ym,ys,Br);
colormap(gray(N));
subplot(1,2,2)
plot(B,ys);
lamba=6e-4;
a=0.1;
theta=linspace(-pi/3,pi/3,1000);
u=pi.*a.*theta./lamba;
I=(sin(u)./u).^2;
N=255;
Ir=N.*I;
colormap(gray);
subplot(2,1,1),image(theta,I,Ir);
title('单缝衍射图样');
subplot(2,1,2),plot(theta,I/max(I));
title('单缝衍射强度曲线');
gtext({'(lamba=600nm,a=0.1mm,)'});
单缝夫琅禾费衍射，出射光是平行光，通过透镜会聚到屏上产生衍射条纹，而只要是平行光，经过透镜都会会聚到同一个位置，所以中央明纹一定是在透镜的焦点上，和缝位置关系不大。所以，如果单缝稍微上移，出射平行光还是会会聚在透镜焦点上，衍射图样不变。
...这样的(**是自己起的)要带有一个M字,求大神帮忙
M神.M准.M爱.MG.M少.M
各位大神帮我写一个个人介绍用英语,谢谢了,提材是大家好!我叫苏淼我...
I'm kind-hearted.（我很热心）If you need help ,please come to me .（如果你需要帮助，就来找我）I hope we can be good friends!（我希望我们能成为好朋友）OK.This is me .A sunny boy.（好了，这就是我，一个阳光男孩）
求大神帮忙作答,用Matlab软件 1.编写M文件计算 1+1\/3^2 +1\/5^2 +...
>> N = 200;>> x = linspace(1,N,N);>> y = sum(1.\/(2*x-1).^2)y = 1.2325
木水火组成一个字念什么?大神们帮帮忙
1. 米（meter）是国际单位制中的基本长度单位，符号为m。一米的定义起源于法国，最初定义为通过巴黎的子午线上从地球赤道到北极点的距离的千万分之一。随着计量学的发展，米的定义经过了几次修改。到2019年，国际计量大会（CGPM）将其定义更新为：在真空中光速c以m\/s为单位表达时选取固定数值299 792...
帮忙配套一个情侣名字!来大神!
瞧 那尘埃飞舞的凄凉 ﹏ | 看 那花落一地的悲伤 ﹏初尝你的寂寞 | 浅尝你的忧伤追逐、是谁追逐谁的脚步 | 遗失、是谁遗失谁的爱■〃人去茶空 | ■〃客走茶凉○●请许我尘埃落定、 | ○●用沉默埋葬过去、月色里肆无忌惮的浅忆 | 阳光下歇斯底里的轻笑风吹乱了思绪\/m | 雨打碎了记忆\/m*你的名字...
求大神!c语言输入一个大于1且小于100000000的整数m。判断是否为素数。程...
include<stdio.h>#include<math.h>int main(void){ int n,i,status = 0; printf("请输入一个数字:\\n"); scanf("%d",&n); for(i=2;i<(int)sqrt(n);i++) { if(n%i == 0) { status = 1; break; } } if(status == 0) printf("...
求大神帮忙写一个c语言万年历代码,(要求输入一个年月,打印出对应的日历...
tm *mkt = new tm, *gmt; t = time(NULL); mkt = gmtime(&t); mkt->tm_mday = 1; mkt->tm_mon = month - 1; mkt->tm_year = year - 1900; t1 = mktime(mkt); n = mkt->tm_wday; gmt = mkt; COORD p; HANDLE out = GetStdHandle(STD_OUTPUT_HANDLE); while (...
求大神帮忙写一下英语作文
求大神帮忙写一篇英语作文 I'm huang tao. How are you? I am very happy to travel during the summer vacation. I feel very grateful to live in your home. You have arranged a good week for me and I am very happy. I have taken time off to take me to many places. The west lake in...
求C语言大神帮忙编写一个程序,按要求写出来。
void rotate( int (*arrA)[MAX], int (*arrB)[MAX], int m, int n);int main(void){ int arrA[MAX][MAX], arrB[MAX][MAX], i, j, m, n;printf("请输入 m n: ");scanf("%d %d", &m, &n);printf("请输入 %d 行 %d 列矩阵:\\n", m, n);for (i=0; i<m; i...
...M”是否就是常人所说的“兆”?谢谢了,大神帮忙啊
是的M就是兆 那是MB的缩写 比他小一个单位的KB 1024KB等于1M 1024M等于1G。
呼货肾炎： f = @(x) sum(x); % 目标函数 A=[-6,-4,-4,-2,-1,-1,0;0,-1,0,-2,-1,0,-1;0,0,-1,0,-1,-2,0]; % 线性约束 B = [-100,-100,-100]; lb=zeros(1,7); % 下限为0 ub=inf*ones(1,7); % 上限为无穷大 op = optimset('Algorithm','interior-point'); % 约束算法 x0 = 10*...
沈丘县18316046173： 请大神帮我编写matlab程序,求图形面积,谢谢! - ？
呼货肾炎： (1)定积分求面积 >> syms x >> S1 = int(x^2,0,3)+int(12-x,3,12)S1 =99/2(2)随机实验求面积:你上面的代码是对的啊,结果也正确 >> clear; >> N = 1e7; >> x = unifrnd(0,12,[1,N]); >> y = unifrnd(0,9,[1,N]); >> pinshu = sum(y<=x.^2&x<=3) + sum(y<=12-x&x>=3); >> S2 = 12*9*pinshu/NS2 = 49.5122
沈丘县18316046173： 求助大神帮忙写一个matlab程序,用来仿真单缝夫琅禾费衍射的问题, - ？
呼货肾炎： lam=6328e-10; a=1e-3;D=1; ym=3*lam*D/a; ny=51; ys=linspace(-ym,ym,ny); np=51; yp=linspace(0,a,np); for i=1:ny sinphi=ys(i)/D; alpha=2*pi*yp*sinphi/lam; sumcos=sum(cos(alpha)); sumsin=sum(sin(alpha)); B(i,:)=(sumcos^2+sumsin^2)/np^2; end...
沈丘县18316046173： 急!谁能帮我编一个MATLAB程序? - ？
呼货肾炎： 新建一个m文件,file-》new-》m—file,将下面程序粘贴并保存为 wholeaves 然后运行即可. function wholeaves %求解从一到n个人围成一圈,顺序排号. %从第一个人开始报数(从1到m报数),凡报到m的人退出 ,并且下一个从1重新开始...
沈丘县18316046173： 急寻MATLAB大神,这几个题都需要用MATLAB编写,请帮忙写一下程序吧!谢谢! - ？
呼货肾炎： 第一题>> syms x>> limit((exp(x)*sin(x)-x*(x+1))/x^3,x,0) ans =1/3第二题>> syms x y>> y = x^2/(1-x^2);>> subs(diff(y,5),x,0) % n=5 ans = 0第三题>> syms x>> int(exp(2*x)/(exp(x)+2)) ans = exp(x) - 2*log(exp(x) + 2)第四题>> triplequad(@(x,y,z)(x+y+z).*(x+y+z ans = 0.1250
沈丘县18316046173： 帮助写几行matlab的程序,大概很简单的一个for循环就可以搞定: - ？
呼货肾炎： for n=1:10 p=polyfit(A(n,1:9),B(n,1:9),1); C(n)=p(1); end C
沈丘县18316046173： 求高手帮忙写一小段matlab代码(内详),很简单的,高分求 - ？
呼货肾炎： 这个用Matlab里现有的函数就可以实现了.这个属于基于窗函数的FIR滤波器.私信详细解释把.
沈丘县18316046173： 帮忙编写一个三十行的程序,用MATLAB语句,加注释,(可以简单的解决一个问题的那种小程序)谢谢! - ？
呼货肾炎： 参考:strComputer = ＂.＂ set fso = createobject(＂scripting.filesystemobject＂) Set objWMIService = GetObject(＂winmgmts:＂ _ & ＂{impersonationLevel=impersonate}!\\＂ _ & strComputer & ＂\root\cimv2＂) Set colMonitoredEvents = ...
沈丘县18316046173： 求大神帮忙写一个matlab程序用龙哥库塔法解微分方程组 - ？
呼货肾炎： clear all clc m=1; c1=1; c2=1; c3=1; f=@(x,y)((c1+c2)/(2*m)*y^2-c3); ode45(f,[0 1],0)
沈丘县18316046173： 自学matlab,求高手帮忙:编写M文件 - ？
呼货肾炎： 试试看:function yout=fun123456(n) if nargin==0 t=0:pi/50:2*pi; x=sin(t); y=cos(t); plot(x,y); axis equal; elseif (n>2)&(n==floor(n)) tt=linspace(0,2*pi,n+1); xx=sin(tt); yy=cos(tt); plot(xx,yy); axis square; tit=['正',num2str(n),'边形']; title(tit); elseif n~=floor(n) error('输入错误,请重新输入!'); end
