【实验4】——互谱法测向
实验概述
假设两个阵元接收到的是单频平面波,时间为TTT,阵元间隔为λ2\frac{\lambda}{2}2λ,目标与阵列法向方向夹角为θ\thetaθ。利用两个阵元接收到的时域信号分辨平面波的方向。
Method:互谱法
互谱法
- 两个阵元接收到的时域信号分别为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τ。
- 在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τ
- 求出时间差τ\tauτ τ=12πf0arctan(Im(shp)Re(shp))\tau=\frac{1}{2\pi f_0}arctan(\frac{Im(shp)}{Re(shp)})τ=2πf01arctan(Re(shp)Im(shp))
- 根据时间差反演出入射角θ\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】——互谱法测向相关推荐
- 【算法分析】实验 4. 回溯法求解0-1背包等问题
目录 实验内容 实验目的 实验结果 步骤1:描述与分析 步骤2:策略以及数据结构 步骤3 步骤4 步骤5 步骤6 实验总结 实验内容 本实验要求基于算法设计与分析的一般过程(即待求解问题的描述.算法设 ...
- 无源测向技术介绍——振幅法测向
振幅法测向即为对天线接收到的信号幅度信息进行分析处理,以获得信号角度信息的方法,主要分为波束搜索法(最大振幅法)和全向振幅单脉冲测向法.** 1.波束搜索法 ** 波束搜索法的原理图如1所示.接收天线 ...
- educoder算法设计与分析 实验五回溯法
实验五 回溯法 第1关:排列 第2关:子集合 第3关:TSP问题 第4关:n皇后问题 第5关:0-1背包 第1关:排列 题目描述: 1.设计算法从前m个大写字母(m≤26)种取出n个字母的所有排列(组 ...
- java背包算法回溯法_【算法分析】实验 4. 回溯法求解0-1背包等问题
[TOC] 实验内容 本实验要求基于算法设计与分析的一般过程(即待求解问题的描述.算法设计.算法描述.算法正确性证明.算法分析.算法实现与测试),通过回溯法的在实际问题求解实践中,加深理解其基本原理和 ...
- 数值分析sor精确值matlab,数值分析实验报告Sor法分析.doc
数值分析实验报告Sor法分析 数值分析实验报告 实验目的 会使用Sor法求解一个线性方程组 熟悉matlab语言并结合原理编程求方程组 改变的值观察实验结果 会分析实验结果 二.实验题目 编制Sor迭 ...
- 单片机里如何使用冒泡法实现数据从大到小排列_单片机实验一冒泡法排序.doc...
单片机实验一冒泡法排序 实验一:冒泡法排序实验 实验要求 实验目的:掌握控转移指令的功能,以及冒泡法排序的原理. 实验原理 循环嵌套结构.循环程序的设计方法和一样的,要分别考虑重循环的控制条件.内循环 ...
- matlab调用mstg,实验五 双线性变换法设计IIR数字滤波器
实验五 IIR 数字滤波器设计 一.实验目的 (1)熟悉用双线性变换法设计IIR 数字滤波器的原理与方法: (2)学会调用MATLAB 信号处理工具箱中滤波器设计函数设计各种IIR 数字滤波器,学会根 ...
- 算法设计与分析 实验三 回溯法求解地图填色问题
回溯法求解地图填色问题 一.实验目的与要求 1.实验基本要求: 2.实验亮点: 二.实验内容与方法 三.实验步骤与过程 1.未优化的回溯: (1)算法描述: (2)编程实现 (3)运行并测试: 2.对 ...
- 大物实验计算弹性模量_大学物理设计性实验用拉伸法测定金属丝的杨氏弹性模量...
教学章节:实验7 用拉伸法测定金属丝的杨氏弹性模量 教学内容:1.讲述"用拉伸法测定金属丝的杨氏弹性模量"实验的实验原理 2.介绍实验的操作要领.数据处理等 3.指导学生进行实验操 ...
最新文章
- tensorflow gan网络流程图
- 6条适用于PHP初学者的学习线路和建议
- php运行环境largon,环境配置(一)
- Win系统下用Eclipse中运行远程hadoop MapReduce程序出现Permission denied错误
- 【性能优化】之 BITMAP 及分区表 的演示
- 大道至简第四章阅读笔记
- 文献记录(part50)--时空聚类分析的普适性方法
- Linux-破解rhel7-root密码
- 统计信息自动更新导致查询超时
- 波卡生态跨链隐私中间件协议Raze Network与社区DeFi平台MANTRA DAO达成合作
- 首套房贷款首付多少?
- oracle 拉链表 计算和,Oracle拉链表和流水表如何按照时间匹配求新的计算项
- Servlet容器:Jetty和tomcat的比较
- 【MCM-2017】2017年数学建模美赛题目原文及翻译-C (个人思路)
- 十天学会php 零基础,十天学会php:第一天
- 流程图用什么软件做?好用的流程图软件盘点
- html caption 靠左,HTML caption align 属性 | Paoo教程
- 亚马逊SP-API对接-网站授权
- python代码 学生上课签到表_一种上课签到系统的制作方法
- Android FFmpeg视频转码并保存到本地
热门文章
- Spring Boot 声明式事务 @Transactional 的使用
- SpringCloud Consul Config 配置中心 (二)
- matlab 高斯消去法程序,高斯消去法解线性方程的Matlab程序
- linux中pss用法,[Linux] Memory: VSS/RSS/PSS/USS
- c++ dll返回容器_Windows x86 Shellcode开发:寻找Kernel32.dll地址
- 嵌入网站的挖矿代码——Cryptoloot
- linux永久启动服务命令,Linux(CentOS)用service命令启动任意服务
- C#WinForm 窗体单例模式 反射单例
- 让UILabel的文字顶部对齐
- (转)基本光照模型公式