阵列信号处理仿真一——延时求和滤波器

问题导入

假设原点处的阵列传感器接收的信号为平面波，表达式如下
$$f(t)=cos(2\pi f_{0}t+{\theta }_0)$$
其中 ${\theta }_0=\frac{\pi }{3}$ ,$f_0=300Hz$
平面波入射的角度$\theta =30^{\circ }$，$\varphi =60^{\circ }$阵列的采样频率$f_s=20000Hz$

假设在Z轴上存在5个阵列关于原点对称，间隔0.1m,试画出$f_1,f_2,f_3,f_4,f_5$接受到的平面波

解

$$f(t,p)=\left[\begin{array}{c}f(t-{\tau }_1)\\ f(t-{\tau }_2)\\ f(t-{\tau }_3)\\ f(t-{\tau }_4)\\ f(t-{\tau }_5)\end{array}\right]$$
不同位置的阵列接受到的波形为原点处的波形的延时
其中$${\tau }_n=\frac{a^T{p}_n}{c}$$
$$a=\left[\begin{array}{c}-sin(\theta )cos(\varphi )\\ -sin(\theta )sin(\varphi )\\ -cos(\theta )\end{array}\right]$$

    t = 0:0.00005:0.007;%采样频率为20000hz,绘制两个周期的长度
    a = [-sin(pi/6)*cos(pi/3);
         -sin(pi/6)*sin(pi/3);
         -cos(pi/6)  ];%初始化a

     Pn = [0,0,0.2;
         0,0,0.1;
         0,0,0;
         0,0,-0.1;
         0,0,-0.2]';%设置位置坐标
         tao = a.'* Pn / 340;%计算延时
         

矩阵的赋值需要维度和大小相同，所以事先初始化新的平面波函数f和时间T

 f = [0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007];
 T = [0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007];

绘制波形

    for i = 1:5
        subplot(5,1,i)
        T(i,:) = tao(i):0.00005:0.007+tao(i); 
        f(i,:) = cos(2*pi*300*T(i,:) + pi/3);
        plot(T(i,:),f(i,:));
    end

效果展示

从上图可以看到，第一幅图和第二幅图在0时刻之前就出现了波形，说明时延${\tau }_n<0$,$f_3$的图形为原点处接收到的波形，第四幅图形和第五幅图形延时${\tau }_n>0$所以波形出现时刻在0时刻之后

对这个阵列实现延时——求和波束形成

解

由上图可知，延时求和波束形成就是将各个阵元接收到的信号的时延通过处理将时延补齐然后求和赋值上取$\frac{1}{N}$后输出
matlab中我们只需将f的各行向量的相位对齐求和即可实现

    delays =[];
    for i = 1:5
        delays(i) = tao(5)-tao(i);%计算时延的差值
    end
    lattic = round(delays/0.00005)%时延的差值除以采样间隔的到采样点的差值
    F =  f(1,21:141) + f(2,16:136)+f(3,11:131)+f(4,6:126)+f(5,1:121);%通过从阵元接受信号的不同时刻进行取样以得到相同的相位
    F = 0.2*F;%1/N = 1/5 =0.2
    figure
    plot(T(5,1:121),F)%时域上的相位和f5的相位对齐，所以取T(5)

lattic的结果为{20，15,10,5,0}，所以不同阵元的取样取值选择了{21，16.11，6，1}

  t = 0:0.00005:0.007;
  
  a = [-sin(pi/6)*cos(pi/3);
        -sin(pi/6)*sin(pi/3);
        -cos(pi/6)  ];
        
    Pn = [0,0,0.2;
        0,0,0.1;
        0,0,0;
        0,0,-0.1;
        0,0,-0.2]';
        
   tao = a.'* Pn / 340;

   f = [0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007];
   T = [0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007;0:0.00005:0.007];
   for i = 1:5
       subplot(5,1,i)
       T(i,:) = tao(i):0.00005:0.007+tao(i); 
       f(i,:) = cos(2*pi*300*T(i,:) + pi/3);
       plot(T(i,:),f(i,:));
   end
   delays =[];
   for i = 1:5
       delays(i) = tao(5)-tao(i);
   end
   lattic = round(delays/0.00005)
   F =  f(1,21:141) + f(2,16:136)+f(3,11:131)+f(4,6:126)+f(5,1:121);
   F = 0.2*F;
   figure
   plot(T(5,1:121),F)