clc;
% 运动声源的接收信号模拟
fs = 4000;
T = 1; % 信号时长
ts = 0: 1/fs: T; % 发射时间序列
fc = 100; % 信号中心频率
st = sin(2* pi * fc * ts); % 信号 s(t)
v = 10; % 声源运动速度(m/s)
L = 100; % 声源和接收的初始距离(m)
c = 1500; % 声速(m/s)
Tr = T * (1 + v/c); % 经过了多普勒后的信号时长
tau0 = L / c; % 信号到达的时间
tr = 0: 1/fs: Tr; % 重采样
trNor = tr / max(tr); % 归一化插值方法
tsNor = ts / max(ts);
sr = interp1(tsNor, st, trNor, 'linear', 'extrap'); % 接收信号
figure;
plot(ts, st, 'DisplayName', '发射信号');
hold on;
plot(tr+tau0, sr, 'DisplayName', '接收信号');
legend();
ylim([-2 2]) 
xlabel('时间 (s)')
ylabel('幅度')