怎么用matlab作出两个系统的相位差,如何求两个信号的相位差

在信号处理里经常会碰到，用于相位差的非参数估计，我接触到的典型的用处就是在ISAR成像中用于平动相位估计。

假设两个频率相同具有相位差的信号

其时间互相关函数为

相位差

% test programme for phase error estimation

clc; clear;

%% Signal

f0 = 200; % signal frequency (Hz)

fs = 1000; % sampling frequency (Hz)

T = 5 / f0; % sampling time (s)

N = round(T * fs);

N = N + mod(N, 2);

t = [0: N - 1] / fs;

phase_1 = pi * rand(1, 1);

phase_2 = pi * rand(1, 1);

phase_error = phase_1 - phase_2

s1 = exp(1i * 2 * pi * f0 * t + 1i * phase_1);

s2 = exp(1i * 2 * pi * f0 * t + 1i * phase_2);

%% Phase error estimation

corr = xcorr(s1, s2); % cross-correlation function

corr_abs = abs(corr);

phase_error_estimated = angle(corr(corr_abs == max(corr_abs)) / abs(corr(corr_abs == max(corr_abs))))

2.一个简单的办法是构造一个没有lag的，频率相同的标准信号。然后做两者的cross-correlation, 然后找出最大的Lag, 看它对应的时间点是什么。

举例来说

比如你有信号s2, 延时是0.35s(相位差为2pi.35 那么你可以构造一个0相位的s1

t = [0:127]0.02;

f = 1.0;

s1 = sin(2pift);

s2 = sin(2pif*(t-0.35));

做cross-correlation

x = xcorr(s1,s2,‘coeff’);

tx = [-127:127]*0.02;

x的最大值那点所对应的横坐标(tx)就是延时。

两个信号分别与cos(2πft)和sin(2πft)做内积，结果记为X1,Y1；X2,Y2

求出real和imag 使用atan2(real，imag)求得相位

atan2(Y1,X1) 和 atan2(Y2,X2)分别是两个信号的绝对相位，取值范围是 (-π,π]

把绝对相位做差得到相对相位

% H=hilbert(x);%Hilbert变换

% ang=angle(H) %求折叠相位，这个就是瞬时相位

**5.**另有：Z为一个复数时，R=abs(Z)，theta=angle(Z),之后利用Z=R.exp(itheta)可以将Z复原.

怎么用matlab作出两个系统的相位差,如何求两个信号的相位差相关推荐

分别已知两直线上的两点，求两直线交点
分别已知两直线上的两点,求两直线交点求两直线的交点是初中数学的简单问题了,在直角坐标系中直线有很多种表示方式.同时我们知道两点确定一条直线,已知两点坐标自然能求出直线坐标,已知两直线坐标自然能求出两 ...
已知两圆圆心坐标及半径求两圆交点 (C语言|参数方程求解)
已知两圆圆心坐标及半径求两圆交点 (C语言|参数方程求解) 在一个二维平面上给定两个圆的圆心横纵坐标.半径共6个参数, 求交点. 这个问题无非是解二元二次方程组.普通二元二次方程联立消元求解的困难在于 ...
python两个列表的差集_Python求两个list的差集、交集与并集的方法
本文实例讲述了Python求两个list的差集.交集与并集的方法.分享给大家供大家参考.具体如下: list就是指两个数组之间的差集,交集,并集了,这个小学数学时就学过的东西,下面就以实例形式对此加以 ...
python求两数最大公因数_使用辗转相除法求两个数的最大公因数（python实现）
数学背景: 整除的定义: 任给两个整数a,b,其中b≠0,如果存在一个整数q使得等式 a = bq 成立,我们就说是b整除a,记做b|a. 性质1:如果c|a,c|b,且对于任意的整数m,n,则有c| ...
两个圆公切线求法_求两圆的公切线计算
关于两圆的公切线计算我们首先需要搞清楚的是两个圆之间有哪些位置关系,两个圆的位置包括内含.内切.外切.外离几种方式,可以看出位置关系可以根据"两圆心之间的距离"与"两圆半 ...
js寻找两个数组的差集_js求两个数组的交集|并集|差集|去重
let a = [1,2,3], b= [2, 4, 5]; 1.差集 (a-b 差集:属于a但不属于b的集合) a-b = [1,3] (b-a 差集:属于b但不属于a的集合) b-a = [4 ...
c语言有参有类最小公倍数,【C语言】写一个函数，并调用该函数求两个整数的最大公约数和最小公倍数...
程序分析: 在数学中,两个数的最小公倍数=两个数的乘积/两数的最大公约数. 求两个数的最大公约数,运用辗转相除法:已知两个整数M和N,假定M>N,则求M%N. 如果余数为0,则N即为所求:如果余 ...
嵌入式 ARM 汇编编程例题（二维数组按规律求和，求两数 gcd / lcm，求数组 min / max，字符串复制，排序）
文章目录 0x00. 整数加减乘除 0x01. 一维数组按某种规律求和 0x01.1 求 1~100 之和 0x01.2. 求一维数组的和 0x01.2. 求一维数组的所有奇数的和 0x02. 二维数 ...
python编写函数、计算三个数的最大公约数_python 函数求两个数的最大公约数和最小公倍数...
1. 求最小公倍数的算法: 最小公倍数 = 两个整数的乘积 / 最大公约数所以我们首先要求出两个整数的最大公约数, 求两个数的最大公约数思路如下: 2. 求最大公约数算法: 1. 整数A对整数B进行 ...

怎么用matlab作出两个系统的相位差,如何求两个信号的相位差

怎么用matlab作出两个系统的相位差,如何求两个信号的相位差相关推荐

最新文章

热门文章