exn621.m

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

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第六章例6－2－1程序
% 真空中的抛射体运动
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月
%
clear;
y0 = 0;  x0 = 0;         % 初始位置
vMag = input('输入初始速度 (m/s):(书上为50) '); 		% 输入初始速度
vDir = input(' 输入初速方向(度):(书上为40或50) ');
yf = input('输入目标高度（米）:(书上为8)  '); 			% 输入目标高度
vx0 = vMag*cos(vDir* (pi/180));    % 计算x,y方向的初始速度
vy0 = vMag*sin(vDir* (pi/180));    % 
wy = -9.81; wx = 0;           % 重力加速度 (m/s^2)
tf=roots([wy/2,vy0,y0-yf]); % 解方程wy*t^2/2+vy0*t+y0=yh，计算落点tf
tf=max(tf); 			% 去除落点时间tf中的庸解
t=[0:0.1:tf,tf];		% 设定时间数组，因tf不大可能被0.1整除，必须加一个tf点
y = y0 + vy0*t + wy*t.^2/2;		 % 计算轨迹
x = x0 + vx0*t + wx*t.^2/2;
xf = max(x),plot(x,y),grid 		 % 计算射程，画出轨迹
set(gcf,'color','w')