exn546a.m

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

%《MATLAB及其在理工课程中的应用指南》(十一五规划版）第五章例5-4-6程序a
% 三重积分的积分区域绘制
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月
%
%本程序给出由x=1,y=x,z=xy三个曲面以及z=0坐标平面围成的积分区域.
clear, close all
[x,y]=meshgrid(0:.05:1);    % 确定计算和绘图的定义域网格
z1=x.*y.*(y-x<0);           % 求z1=xy并去掉定义域外网格区的z1
mesh(x,y,z1);hold on;       	% 画出积分区顶部
% 画出积分区域的几个侧面
x1=[0:0.02:1];y1=x1;sx1=length(x1);		% y=x平面上的自变量x1,y1
zd=[zeros(1,sx1);x1.*y1];			% y=x平面上与曲面交点的数据zd=x1*y1
% 画出y=x上的平行线族,长度为从0到z=x1.*y1,两端点画’*’号
plot3([x1;x1],[y1;y1],zd,'*')		% 画出此平行线族,端点画*
%	画出x=1平面上沿z的平行线族,长度为从0到z=1*y1,两端点画’o’号
line(ones(2,sx1),[y1;y1],[zeros(1,sx1);y1])    %画出x=1平行线族
plot3(ones(2,sx1),[y1;y1],[zeros(1,sx1);y1],'o') %画出x=1平行线族端点o
xlabel('x'),ylabel('y'),zlabel('z'),hold off
set(gcf,'color','w')