![]() |
||||||
|
| Script S6_2_12B.m | |||
|
|||
|
%==============================================
%doppler 2B %============================================== % % %light speed in m/sec c=3.0e+08; %negative relative velocities: the two inertial system %are approaching vr=-(0.1:0.1:0.9) %corresponding positive relative velocities v=vr*c %frequency nu of the beam in the xyz system in Hz nu=1.0e+014 %angle of reflection in the xyz system teta_g=30 teta=(pi/180)*teta_g; %gamma arg=1-vr.^2; gamma=1./sqrt(arg) %frequencies nup in the x'y'z' costeta=cos(teta); fatt1=vr*costeta; fatt2=gamma*nu; fatt3=1-fatt1; nup=fatt2.*fatt3 %corresponding angles in degrees fatt4=costeta-vr; tetap=(acos(fatt4./fatt3))*(180/pi) %their plots plot(vr,nup,'r-o'),grid on title('nup varying vr with the assigned nu') figure plot(vr,tetap,'bs-'),grid on title('tetap varying vr with the assigned tetap') %============================================== % |
|||
| Top | |||