exn547c.m

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

%《MATLAB及其在理工课程中的应用指南》(十一五规划版）第五章例5-4-7程序
% 求微分方程的数值解
% 西安电子科技大学出版  陈怀琛著 2007年5月

clear, close all
t0=0;	tf=3*pi; x0=[0; 0];	% 给出初始值
[t,x]=ode23('exn547f',[t0,tf],x0);   % 此处显示结果
y=x(:,1);				% x1为x的第一列
%本题的解析结果为 yt(I)=(1+2/(pi^2))*(1-cos(t(I)))-t(I)^2/(pi^2); 
%在数值积分输出的时间序列点t上计算它的值并画图与数值解作比较，
for I=1:length(t);
  	yt(I)=(1+2/(pi^2))*(1-cos(t(I)))-t(I)^2/(pi^2);	% 解析解计算
end
u=1-(t.^2)/(pi^2);
clf,plot(t,y,'--',t,u,'+',t,yt,'o')
% title('ODE23  ---  Input ++++ Analytical ooo')
legend('数值积分解','输入量','解析解')
xlabel('t'),ylabel('y(t)'),grid