如何使用Golomb样本信号的尺子吗?

6视图(30天)
我想样品Golomb统治者multi-sine信号的长度为151,如何实现这一目标?
下面是正弦波采样的信号图Golomb长度为151的统治者。
振幅谱上面的样品(可见在期望频率达到24赫兹以上“under-sampling”噪音)
清晰的所有
关闭所有
clc
% %
Fs = 1000;%的采样频率。
T = 1 / f;%采样周期
L = 1024;%长度的信号
t = 2 * (0: l - 1) * t;%的时间向量
x = 0;%的信号开始值
y = 0;%的信号开始值
F = 24;%的信号频率。
seq = [0 4 20 30 57 59 62 76 100 111 123 136 144 145 151];% Golomb Seq。? ?
% x = cos(2 *π* F * seq);
i01 = 1:长度(seq)
u = 0: seq (i01)
x = x + cos(2 *π* F * u +π/ 2);
结束
结束
次要情节(2、2、1)
情节(seq, x)
标题(原始信号的)

接受的答案

Bjorn Gustavsson
Bjorn Gustavsson 2022年1月31日
唯一有意义的想法我抽样这样使用Golomb-ruler是你希望得到的估计及信号的功率谱计算的所有组合complex-conjugated产品的样品。从第一个和第二个样本(对应于时间步的0和4),你会得到一个估计的自相关函数,ACF,滞后4,从第一个和第三个样本(correpsponding时间步0 - 20)你会得到一个估计的ACF的落后20日从第二个和第三个样本(correpsponding样品时间步4和20)你会得到一个ACF-estimate 16的滞后等Golomb-rulers有几个相当不错为这种类型的操作特性。这个建议你应该能够得到ACF的滞后除了46 0和151之间。如果你知道你的信号是实值你也知道真正的ACF的一部分甚至和虚部是奇数。可以设置你Fourier-transforming ACF的功率谱
HTH
5个评论
Bjorn Gustavsson
Bjorn Gustavsson 2022年2月1日
首先,你必须解决 seq - variable上面脚本:
seq = [1 4 20 30 57 59 62 76 100 111 123 136 144 149 151];
不是GR,自从GR的定义应该独特的不可重复的元素之间的差异。在你的 seq 我们有一个区别4和1 3单位,和另一个62和59之间的区别。你必须每个元素加1,不仅改变第一0到1。是这样的:
seq = [0 4 20 30 57 59 62 76 100 111 123 136 144 149 151] + 1;
接下来你应该做什么是计算每个独特的产品从GR-samples一双样品。给你估计的ACF的时间间隔对应于样品之间的时差。如果你做一个双环samples-sequence你从每一步,您会在哪里存储这些组织正确吗?

登录置评。

更多的答案(0)

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!