| Date 06/09/2010 |
|
||||||
|
| Script S6_2_13.m | |||
|
|||
|
%==============================================
%doppler 3 %============================================== % %light speed in m/sec c=3.0e+08; %negative relative and absolute velocities % vr=-0.7 % vr=-0.9 vr=-0.80; v=vr*c %frequencies nup observed in the x'y'z' system nup1=(3.0)*10^+014 nup2=(1.67)*10^+014 nup=linspace(nup1,nup2,10) %corresponding angles in degrees tetapg1=0; tetapg2=36.9 tetapg=linspace(tetapg1,tetapg2,10) %and in radians tetap=(pi/180)*tetapg; %gamma arg=1-vr^2; gamma=1/sqrt(arg); % values of nu in xyz system fatt1=vr*cos(tetap); fatt3=1+fatt1; fatt2=gamma*nup; nu=fatt2.*fatt3 %corresponding angles in degrees fatt4=cos(tetap)+vr; arg1=fatt4./fatt3; tetag=(acos(arg1))*(180/pi) %plots plot(nup,nu,'r-o'),grid on title('nu varying nup between the assigned values') figure plot(tetapg,tetag,'bs-'),grid on title('tetag varying tetapg between the assigned values') |
|||
| Top | |||
