%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第六章例6-6-2程序
% 多普勒效應的演示
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
x0=500;v=60;y0=30;		% 設定聲源運動參數
c=330;w=1000;				% 音速和頻率
t=0:0.001:30;				%	設定時間數組
r=sqrt((x0-v*t).^2+y0.^2);	% 計算聲源與聽者距離
t1=t-r/c;						% 經距離遲延后聽者的等效時間
u=sin(w*t)+sin(1.1*w*t);	% 聲源發出的信號
u1=sin(w*t1)+sin(1.1*w*t1); 	% 聽者接受到的信號
disp('靜止聲源的聲音')	    % 先后將原信號和接受到的信號恢復為聲音
pause,sound(u);pause(8);		
disp('運動聲源的聲音')
sound(u1);