exn526.m

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

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5-2-6
% 用平行截面分析空间曲面
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月
%
clf,clear
[x,y]=meshgrid(-10:.2:10);		% 确定计算和绘图的定义域网格
z1=(x.^2-2*y.^2)+eps;		% 第一个曲面方程
a=input('a= (-50<a<50) ');
z2=a*ones(size(x));	%第二个曲面方程（平面）
subplot(1,3,1),mesh(x,y,z1);hold on;mesh(x,y,z2);	%分别画出两个曲面
v=[-10 10 -10 10 -100 100];axis(v),grid			%确定第一个分图的坐标系
colormap(gray),hold off,	%取消彩色改为灰度
r0=abs(z1-z2)<=1;			%求两曲面z坐标差小于0.5的网格
zz=r0.*z2;yy=r0.*y;xx=r0.*x;		%求这些网格上的坐标值,即交线坐标值
subplot(1,3,2),plot3(xx,yy,zz,'x');%画出这些点
axis(v),grid						% 使第二个分图取第一个分图的坐标系
pause,subplot(1,3,3),
contour3(x,y,z1,20)   % 用等高线命令将更为方便地求出多根交线
set(gcf,'color','w')