%《MATLAB及其在理工課程中的應(yīng)用指南》(十一五規(guī)劃版)第六章例6－2－4程序a
% 碰撞的擴(kuò)展（非正撞情況）
% 西安電子科技大學(xué)出版社出版  陳懷琛編著 2007年7月
% 本程序不夠完善，僅供參考
%
clear
K=100,logspace(-1,1,11);		% 把K設(shè)成自變量數(shù)組,
ur=zeros(9,2);vr=zeros(9,2);phi=zeros(9,2);
for i=1:9			% 對各個K循環(huán)計算
   theta=(i-1)*pi/8;
   ur1=roots([(1+1/K),-2/K*cos(theta),(1/K-1)]);	% 二次方程有兩個解
   ur(i,:)=ur1';%(abs(ur1-1)>0.001);		% 去掉在1鄰近的庸解
   vr(i,:)=(1-ur(i,:)*cos(theta))./K;
   phi(i,:)=asin(ur(i,:)*sin(theta)/K./vr(i,:)); 
end							% 用（5）式求vr,用元素群運算
[[0:pi/8:pi]',ur,vr]									% 顯示輸出數(shù)據(jù)
e=1-ur.^2;
[[0:pi/8:pi]',phi,e]									% 顯示輸出數(shù)據(jù)
plot([0:pi/8:pi]',[ur,vr,e]),grid			%	繪圖
set(gcf,'color''w')


% h=6.626e-34;		普朗克常數(shù)[焦.秒]
% m0=9.11e-31;		電子質(zhì)量[千克]
% c=3e8;				光速[米/秒]
% lx=1e-10;			x射線波長[米]
%?nux=c/lx=3.0000e+018				x射線頻率[1/秒]
%?e=h*nux=  1.9878e-015			x射線能量[焦]
%?mx=e/c/c=  2.2087e-032			x光子質(zhì)量
%?px=mx*c=e/c						x光子動量
%?m0/mx=   41.2466				電子與光子質(zhì)量之比