%《MATLAB及其在理工課程中的應用指南》第六章例6－6－2程序
% 多普勒效應的演示
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
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); % 聽者接受到的信號
sound(u);pause(5);			% 先后將原信號和接受到的信號恢復為聲音
sound(u1);