| 
				Doppler Effects(多普勒效应)     Apparent change in the 
				wave frequency (or in the period of periodic events) due to 
				motion of wave source and/or observer is known as Doppler 
				effect. In the 17th century (well before Doppler effect was 
				predicted), Roemer (Danish astronomer) estimated the speed of 
				light from the change in the revolution period of one of 
				Jupiter's moons. When the earth is approaching Jupiter, the 
				period is shortened, while receding, the period is elongated. 
				The first animation shows wave patterns (surfaces of constant 
				phase) emitted by a moving wave source. Wavelengths in front of 
				the moving source are squeezed (shortened) according to and 
				those behind are elongated,where cw is the wave (sound) velocity and Vs is the velocity of 
				the wave source. Therefore, the frequency of a wave source 
				approaching (receding from) an observer is higher (lower) than 
				the "true" frequency,Similar effect occurs when the observer is 
				moving. For sound waves, the velocities of source and observer 
				are well defined relative to air, the wave medium. For 
				stationary sound source and moving observer, the sound speed 
				appears to change (the wavelength remains the same). Therefore, 
				the Doppler shifted frequency for a stationary sound source 
				heard by a moving observer is given by where Vo is the velocity 
				of the observer relative to air, + sign for the case of 
				approaching and ?sign for the case of receding.
 
 For electromagnetic waves in vacuum, the velocity with 
				respect to the wave medium (vacuum) is ill defined. In Doppler 
				effect of electromagnetic waves, only the relative velocity 
				between wave source and observer matters, where is the relative 
				velocity between the wave source and observer normalized by the 
				velocity of light c. The factor in the numerator, is due to 
				intrinsic Doppler shift caused by time dilation effect.
 
 Some objects (e.g., supersonic planes) can move faster than 
				sound waves. Animation below shows wave fronts emitted by a 
				supersonic object when the Mach number is 2. No waves can exist 
				outside the cone created by the object. When the cone surface 
				hits an observer, sudden arrival of a sharp pulse is detected as 
				shock wave. In vacuum, nothing can travel faster than light. 
				However, in material media (glass, plastic, etc.), 
				electromagnetic waves propagate in general at a speed slower 
				than c. Therefore, if a relativistic charged particle with a 
				velocity close to c enters, say, a plastic block, 
				electromagnetic shock wave is emitted. This phenomenon is called 
				Cerenkov radiation and is in wide use for detecting energetic 
				particles.
 
 Doppler effect. Volocity of wave source = one half of the 
				wave velocity. Wavelengths in front of the source are squeezed 
				while those behind are elongated.
 
				animate({[u*sin(t),u*cos(t),t=-Pi..Pi], 
				[(u-1)*sin(t)*Heaviside(u-1)+.5,(u-1)*cos(t)*Heaviside(u-1),t=-Pi..Pi],
				[(u-2)*sin(t)*Heaviside(u-2)+1.,(u-2)*cos(t)*Heaviside(u-2),t=-Pi..Pi],
 [(u-3)*sin(t)*Heaviside(u-3)+1.5,(u-3)*cos(t)*Heaviside(u-3),t=-Pi..Pi],
 [(u-4)*sin(t)*Heaviside(u-4)+2.,(u-4)*cos(t)*Heaviside(u-4),t=-Pi..Pi],
 [(u-5)*sin(t)*Heaviside(u-5)+2.5,(u-5)*cos(t)*Heaviside(u-5),t=-Pi..Pi],
 [(u-6)*sin(t)*Heaviside(u-6)+3.,(u-6)*cos(t)*Heaviside(u-6),t=-Pi..Pi],
 [(u-7)*sin(t)*Heaviside(u-7)+3.5,(u-7)*cos(t)*Heaviside(u-7),t=-Pi..Pi]},
 u=0..8,view=[-8..8,-8..8],frames=40,color=red);
 
				  Velocity of wave source = 
				twice of the wave velocity. Formation of shock wave. 
				animate({[u*sin(t),u*cos(t),t=-Pi..Pi], 
				[(u-1)*sin(t)*Heaviside(u-1)+2,(u-1)*cos(t)*Heaviside(u-1),t=-Pi..Pi],
				[(u-2)*sin(t)*Heaviside(u-2)+4,(u-2)*cos(t)*Heaviside(u-2),t=-Pi..Pi],
 [(u-3)*sin(t)*Heaviside(u-3)+6,(u-3)*cos(t)*Heaviside(u-3),t=-Pi..Pi],
 [(u-4)*sin(t)*Heaviside(u-4)+8,(u-4)*cos(t)*Heaviside(u-4),t=-Pi..Pi],
 [(u-5)*sin(t)*Heaviside(u-5)+10,(u-5)*cos(t)*Heaviside(u-5),t=-Pi..Pi],
 [(u-6)*sin(t)*Heaviside(u-6)+12,(u-6)*cos(t)*Heaviside(u-6),t=-Pi..Pi],
 [(u-7)*sin(t)*Heaviside(u-7)+14,(u-7)*cos(t)*Heaviside(u-7),t=-Pi..Pi]},
 u=0..8,view=[-8..16,-8..8],frames=40,color=red);
 
				  > 
				interface(plotdevice=inline): with(plots):
 implicitplot({x^2+y^2=4.5^2,(x-2)^2+y^2=3.5^2,(x-4)^2+y^2=2.5^2,
 (x-6)^2+y^2=1.5^2,(x-8)^2+y^2=0.5^2,x+1.732*y=9,x-1.732*y=9},x=-5..9,y=-5..5,numpoints=1000);
 
				 |