main.m

信号的离散累加

%按照改的坐标系定义 注意几何关系
clear all;
close all;
clc;
f=200;               %频率                  改变频率是没有用处的 因为d=bochang/2固定
c=1500;               %声速
bochang=c/f;          %波长
d=bochang/2;          %阵间距
N=5;                 %阵个数                 阵元个数越多 指向性越尖锐 效果越好 
H=60;                 %基阵到海底的深度
tao=0.1;             %重点：脉宽不能太长
%%%%%%%%%%%%%%%%%%%%%%%根据几何关系
sita=179*pi/180:pi/180:100*pi/180;%基阵主波束与垂直方向夹角
phi=0*pi/180:pi/180:360*pi/180;  %垂直布放,主波束与垂直方向夹角
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
I0=0.001;                 %发射声强
Iref=0.67*10.^(-22);  %参考声强
beam=beam_show(f,N);%显示波束指向性
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
v=0;
zhixiang=zeros(1,91);
%all=0;
for sita=179*pi/180:-pi/180:100*pi/180%指向性也要改
    v=v+1;
    deta=pi-sita;%注意角的意义
    r=H./cos(deta);
    m=(179*pi/180-sita)/(pi/180)+1;
    all=0;
    for i=1:1:(360-0+1)
        B_T=beam(round(m),round(i));%计算发射波束加权值
        all=all+B_T.^2;
    end
    zhixiang(v)=all;
    Ss=reverse(sita);%计算散射系数
    Ib(v)=I0./(r*r*r*r)*all*Ss*r*c*tao./(2*cos(sita-pi/2));
%%%%%%%%%%%%%%%%%%%%%%%%%
end
RL=10*log(Ib/Iref); 
sita=179*pi/180:-pi/180:100*pi/180;
r=H.*tan(pi-sita);
figure(2);
plot(r,RL);
figure(3);
plot(zhixiang);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%造信号
%for sita=179*pi/180:-pi/180:100*pi/180
sita=100*pi/180;%设角度     主波束角度
    rr=H./cos(pi-sita);
    delay=2*rr/c;%此角度回波时间
    mm=(179*pi/180-sita)/(pi/180)+1;
    P=sqrt(zhixiang(mm));
    fs=10*f;
    Ts=1/fs;
    %n=1:0.01*fs;%要求和后面指向性一致
    n=1:0.1*fs;
    signal_T=sqrt(I0)*cos(2*pi*f*n*Ts);
    %signal_T=cos(2*pi*f*n*1/Ts);
    signal_R=P*cos(2*pi*f*(n*Ts+delay));
figure(4);
plot(signal_T);
figure(5);
plot(signal_R);