exn519a.m

基天子matlab的双缝干涉从单缝(线光源)发出的单色光射到与之平行的双缝上

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5-1-9a程序
% 求曲线的极值（直接方法）
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月
%
clear,clf,x=linspace(0,pi,10000);
y=2*sin(2*x).^2+5/2*x.*cos(x/2).^2; 
plot(x,y),grid on,hold on
ymax1=max(y),pause
n01=find(abs(x-1)<2e-4); 	% 找两个边界处的y的下标
n02=find(abs(x-2)<2e-4);	
ymin=min(y([n01:n02]))		% 求x=1~2之间y的最小值
ymax2=max(y([n02:end])) 	% 求x=2~pi之间y的最大值
pause
nm2=find(y([n01:n02])==ymin)+n01-1	% 求y的最小值处的下标
nm3=find(y([n02:end])==ymax2)+n02-1	% 求y的最大值处的下标
xm2=x(nm2),xm3=x(nm3)		% 求y的最小最大值处的x值
pause
[xm1,mym1]=fminbnd('-(2*sin(2*x).^2+5/2*x.*cos(x/2).^2)',0,1);
[xm2,ym2]=fminbnd('2*sin(2*x).^2+5/2*x.*cos(x/2).^2',1,2);
[xm3,mym3]=fminbnd('-(2*sin(2*x).^2+5/2*x.*cos(x/2).^2)',2,3);
xm=[xm1,xm2,xm3],ym=[-mym1,ym2,-mym3]
set(gcf,'color','w')