实验概述

假设两个阵元接收到的是单频平面波,时间为TTT,阵元间隔为λ2\frac{\lambda}{2}2λ​,目标与阵列法向方向夹角为θ\thetaθ。利用两个阵元接收到的时域信号分辨平面波的方向。

Method:互谱法

互谱法

  1. 两个阵元接收到的时域信号分别为s1(t)=s(t)s_1(t)=s(t)s1​(t)=s(t)和s2(t)=s(t+τ)s_2(t)=s(t+\tau)s2​(t)=s(t+τ),则两者之间的频谱关系为s1(f)=s(f)s_1(f)=s(f)s1​(f)=s(f)和s2(f)=s(f)ej2πfτs_2(f)=s(f)e^{j2\pi f \tau}s2​(f)=s(f)ej2πfτ。
  2. 在f0f_0f0​处的对两个频谱值共轭相乘shp=s1∗(f0)s2(f0)=∣s(f)∣2ej2πfτshp=s_1^{*}(f_0)s_2(f_0)=|s(f)|^{2}e^{j2\pi f \tau}shp=s1∗​(f0​)s2​(f0​)=∣s(f)∣2ej2πfτ
  3. 求出时间差τ\tauτ τ=12πf0arctan(Im(shp)Re(shp))\tau=\frac{1}{2\pi f_0}arctan(\frac{Im(shp)}{Re(shp)})τ=2πf0​1​arctan(Re(shp)Im(shp)​)
  4. 根据时间差反演出入射角θ\thetaθ θ=arcsin(2cτd)\theta=arcsin(\frac{2c\tau}{d})θ=arcsin(d2cτ​)

代码实现

clc;clear;
theta=pi/6;
fs=1000;
f0=100;
N=1024;%fft点数 考虑了实际工程中 2的幂次方个数的点fft效率高
ts=1;
t=0:1/fs:ts-1/fs;
c0=1500;
lambda=c0/f0;%平面波波长
d=lambda/2;%阵元间距
tao=d*sin(theta)/c0;%阵元接收的信号的时间差
st1=sin(2*pi*f0*t);%阵元1接收到的信号
st2=sin(2*pi*f0*(t+tao));%阵元2接收到的信号 存在时延相差
sf1=fft(st1,N);%补零作fft
sf2=fft(st2,N);
index=fix(f0/(fs/N))+1;%找到f0对应的那根谱线
shp=conj(sf1(index))*sf2(index);%求互谱值
temp=imag(shp)/real(shp);%求相角的正切值
res=atan(temp)/(2*pi*f0);
angel=asin(c0*res/d)/(2*pi)*360;%反演得到平面波的入射方向

【实验4】——互谱法测向相关推荐

  1. 【算法分析】实验 4. 回溯法求解0-1背包等问题

    目录 实验内容 实验目的 实验结果 步骤1:描述与分析 步骤2:策略以及数据结构 步骤3 步骤4 步骤5 步骤6 实验总结 实验内容 本实验要求基于算法设计与分析的一般过程(即待求解问题的描述.算法设 ...

  2. 无源测向技术介绍——振幅法测向

    振幅法测向即为对天线接收到的信号幅度信息进行分析处理,以获得信号角度信息的方法,主要分为波束搜索法(最大振幅法)和全向振幅单脉冲测向法.** 1.波束搜索法 ** 波束搜索法的原理图如1所示.接收天线 ...

  3. educoder算法设计与分析 实验五回溯法

    实验五 回溯法 第1关:排列 第2关:子集合 第3关:TSP问题 第4关:n皇后问题 第5关:0-1背包 第1关:排列 题目描述: 1.设计算法从前m个大写字母(m≤26)种取出n个字母的所有排列(组 ...

  4. java背包算法回溯法_【算法分析】实验 4. 回溯法求解0-1背包等问题

    [TOC] 实验内容 本实验要求基于算法设计与分析的一般过程(即待求解问题的描述.算法设计.算法描述.算法正确性证明.算法分析.算法实现与测试),通过回溯法的在实际问题求解实践中,加深理解其基本原理和 ...

  5. 数值分析sor精确值matlab,数值分析实验报告Sor法分析.doc

    数值分析实验报告Sor法分析 数值分析实验报告 实验目的 会使用Sor法求解一个线性方程组 熟悉matlab语言并结合原理编程求方程组 改变的值观察实验结果 会分析实验结果 二.实验题目 编制Sor迭 ...

  6. 单片机里如何使用冒泡法实现数据从大到小排列_单片机实验一冒泡法排序.doc...

    单片机实验一冒泡法排序 实验一:冒泡法排序实验 实验要求 实验目的:掌握控转移指令的功能,以及冒泡法排序的原理. 实验原理 循环嵌套结构.循环程序的设计方法和一样的,要分别考虑重循环的控制条件.内循环 ...

  7. matlab调用mstg,实验五 双线性变换法设计IIR数字滤波器

    实验五 IIR 数字滤波器设计 一.实验目的 (1)熟悉用双线性变换法设计IIR 数字滤波器的原理与方法: (2)学会调用MATLAB 信号处理工具箱中滤波器设计函数设计各种IIR 数字滤波器,学会根 ...

  8. 算法设计与分析 实验三 回溯法求解地图填色问题

    回溯法求解地图填色问题 一.实验目的与要求 1.实验基本要求: 2.实验亮点: 二.实验内容与方法 三.实验步骤与过程 1.未优化的回溯: (1)算法描述: (2)编程实现 (3)运行并测试: 2.对 ...

  9. 大物实验计算弹性模量_大学物理设计性实验用拉伸法测定金属丝的杨氏弹性模量...

    教学章节:实验7 用拉伸法测定金属丝的杨氏弹性模量 教学内容:1.讲述"用拉伸法测定金属丝的杨氏弹性模量"实验的实验原理 2.介绍实验的操作要领.数据处理等 3.指导学生进行实验操 ...

最新文章

  1. tensorflow gan网络流程图
  2. 6条适用于PHP初学者的学习线路和建议
  3. php运行环境largon,环境配置(一)
  4. Win系统下用Eclipse中运行远程hadoop MapReduce程序出现Permission denied错误
  5. 【性能优化】之 BITMAP 及分区表 的演示
  6. 大道至简第四章阅读笔记
  7. 文献记录(part50)--时空聚类分析的普适性方法
  8. Linux-破解rhel7-root密码
  9. 统计信息自动更新导致查询超时
  10. 波卡生态跨链隐私中间件协议Raze Network与社区DeFi平台MANTRA DAO达成合作
  11. 首套房贷款首付多少?
  12. oracle 拉链表 计算和,Oracle拉链表和流水表如何按照时间匹配求新的计算项
  13. Servlet容器:Jetty和tomcat的比较
  14. 【MCM-2017】2017年数学建模美赛题目原文及翻译-C (个人思路)
  15. 十天学会php 零基础,十天学会php:第一天
  16. 流程图用什么软件做?好用的流程图软件盘点
  17. html caption 靠左,HTML caption align 属性 | Paoo教程
  18. 亚马逊SP-API对接-网站授权
  19. python代码 学生上课签到表_一种上课签到系统的制作方法
  20. Android FFmpeg视频转码并保存到本地

热门文章

  1. Spring Boot 声明式事务 @Transactional 的使用
  2. SpringCloud Consul Config 配置中心 (二)
  3. matlab 高斯消去法程序,高斯消去法解线性方程的Matlab程序
  4. linux中pss用法,[Linux] Memory: VSS/RSS/PSS/USS
  5. c++ dll返回容器_Windows x86 Shellcode开发:寻找Kernel32.dll地址
  6. 嵌入网站的挖矿代码——Cryptoloot
  7. linux永久启动服务命令,Linux(CentOS)用service命令启动任意服务
  8. C#WinForm 窗体单例模式 反射单例
  9. 让UILabel的文字顶部对齐
  10. (转)基本光照模型公式