求教这一段是什么意思，怎么改数据

q591388

x0=150; %m  
v=50; %车速m/s
y0=30; %m  
c=330; %声音在空气中的速度m/s
f0=2000; %声音的频率Hz
fs=8000; %采样频率
t=0:1/fs:6; %规定t的范围和步进
x_t=-x0+v.*t; %火车与观察者的水平距离
y=y0; %火车与观察者的垂直距离
r=sqrt(x_t.^2+y.^2); %火车与观察者间距离
costheta=x_t./r; %声速矢量与火车和听者间连线夹角的余弦
f=f0./(1+v.*costheta./c); %听者接收到的频率
fmin=min(f);
fmax=max(f);  
voltage=(f-fmin)./(fmax-fmin)*2-1; %归一化调频电压在-1 ~ +1之间
signal=0.5*vco(voltage,[fmin fmax],fs); %压控震荡器，由输入电压控制输出信号频率
u=sin(2*pi*f0*t); %声源发出的信号
sound(u,fs); %播放声源发出的信号
pause(5); %暂停
sound(signal,fs); %播放听者接收到的信号
wavwrite(signal,fs,'dopp.wav'); %音频数据存盘
figure(1);
plot(t,f);
xlabel('Time/s');
ylabel('Freq/Hz'); %做出听者接收到的信号的频率变化曲线
figure(2);  
spectrogram(signal,kaiser(256,5),220,512,fs,'yaxis'); %做出时间频率图

TestGuru

多普勒效应？