%对混响的理解有帮助
%无多普勒频移的海底混响单元散射模型卷积法
clc;clear all;close all;
%参数设置============================================================
fs=200000; %采样频率
f0=30000; %中心频率
%k=1250000; %k =B/t =5kHz/4ms
B=5000;
c=1500; %声速
h=50; %距海底距离
u=-27;
rou=1; %密度
v=20;
azm=pi/6; %声纳水平方位角
startt=0.08; %混响起始时间
endt=0.3; %混响终止时间
T=0.005;%脉冲宽度
ts=0:1/fs:T-1/fs;%脉冲时间序列
s1=exp(j*2.*pi.*f0.*ts); %CW信号
k=B/T; %1000000
s2=exp(j*2.*pi.*(f0-B/2).*ts+j*pi.*k.*ts.*ts); %LFM信号
Ns=length(ts);%Ns 1000个点
k1=2*pi*f0/c;%波数 125.6637
k2=2*pi*(f0-B/2+k*ts)/c; %1*1000
% 信号时域=======================================================================
figure(1);
subplot(2,1,1);
plot(ts,real(s1));%取复信号的实部进行画图,如不取实部,系统自动取实部
axis([0 0.005 -1 1]); %axis([xmin xmax ymin ymax]),单轴xlim([0 5])
xlabel('时间/s');ylabel('幅度');title('CW信号');
subplot(2,1,2);
plot(ts,real(s2));
axis([0 0.005 -1 1]);xlabel('时间/s');ylabel('幅