exn537.m

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

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5-3-7程序
% 方波表为多阶正弦波之和
% 西安电子科技大学出版社出版  陈怀琛编著 2007年7月

t = linspace(0,pi,201);						% 设定一个时间数组，有101个点
	y = sin(t);	plot(t,y),figure(gcf),pause     % 频率为w=1(f=1/2π)的正弦基波
	y = sin(t) + sin(3*t)/3; plot(t,y), pause 	% 叠加三次谐波
	%   用1，3，5，7，9次谐波叠加
   y = sin(t) + sin(3*t)/3 + sin(5*t)/5 + sin(7*t)/7 + sin(9*t)/9;plot(t,y)
   %为了绘制三维曲面，要把各次波形数据存为一个三维数组，因此必须
   %重新定义y，重编程.由于拟求至19次谐波，把点取密一些。
y = zeros(10,max(size(t))); x = zeros(size(t));
for k=1:2:19
	x = x + sin(k*t)/k; y((k+1)/2,: ) = x;
end                                     
%      将各波形迭合绘出
pause, plot(t,y(1:9,: ))
%      将各波形绘成三维网格图，看出增加谐波阶次对方波逼近程度的影响
pause, mesh(t,[1:10],y), pause
clc