duopule.m

来自「该程序是信号的多普勒现象」· M 代码 · 共 21 行

%多普勒效应
x0=550;
v=150; %车速
y=30;
c=340;%声速
wc=2200;%声的频率
t=0:0.0005:30;
r=sqrt((x0-v*t).^2+y.^2);               %计算声源与听着间的距离
t1=t-r/c;                          %径距离迟延后听者的等效时间
u=1.5*sin(wc*t); %声源发出的信号
u1=1.5*sin(wc*t1);%听者接收到的信号

%sound(u,10000);pause(5),sound(u1,10000);


subplot(211)
plot(u);axis([11500,11700,-2.2,2.2])


subplot(212)
plot(u1,'r');axis([5500,7800,-2.2,2.2])