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