exn547.m

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

%《MATLAB及其在理工课程中的应用指南》(十一五规划版）第五章例5-4-7程序
% 求微分方程的数值解
% 西安电子科技大学出版  陈怀琛著 2007年7月
%
clear, close all
t0=0;	tf=3*pi; x0=[1; 0];	% 给出初始值
[t,x]=ode23('exn547bf',[t0,tf],x0);   % 此处显示结果
y=x(:,1);plot(t,y)			% y为x的第一列
%本题的解析结果为 yt(I)=(1+2/(pi^2))*(1-cos(t(I)))-t(I)^2/(pi^2); 
%在数值积分输出的时间序列点t上计算它的值并画图与数值解作比较，
%u=1-(t.^2)/(pi^2);
xlabel('t'),ylabel('y(t)'),grid
set(gcf,'color','w')