速率信号法

源自师姐论文《卫星信号频谱监测系统关键技术研究与实现-董雪》

速率信号法的实现

思想

算法思想看论文，这里主要就是注意index那里的细节，记得减1.

没有用czt细化频谱，用了zoomfft来细化，这个方法更加容易工程实现。所有的频谱变换全部只用移动数组就可以，不用去做指数计算。减少计算量。zoomfft原理可以参考了这里

问题：为什么细化后得到的结果不对，和估计的结果值有所偏差，目前还没想明白。可能在计算一些地方的时候的近似运算会导致这个结果。

源代码

% 测试速率信号法
clc
clear all
close allsr = 1000;
fs = 8e3;N = sr;srcData = randsrc(1, N, [0:3]);
modData = pskmod(srcData, 4);
% 脉冲成型滤波器
ipoint = fs / sr;
span = 21;
pointFir = span * ipoint;
upData = upsample(modData, ipoint);
sendPulser = rcosdesign(0.5, span, ipoint, 'sqrt');
sendData = conv(upData, sendPulser);
% 形成发射信号，这里主要减去卷积带来的时延部分，详细请看数字信号处理
sendData = sendData(pointFir/2 + 1 : end - pointFir/2);% 计算速率信号
for i = 1:length(sendData) %求得速率信号if(i<length(sendData))res(i) = (norm(sendData(i+1)-sendData(i)))*fs; %norm(x)求x的2范数  end
end% 寻找小于长度的最大的2的幂次方长度
n=nextpow2(length(res)); %求出 2^n >= length(res) 的n值
n=n-1;
len=2^n;
% 数据长度改变，分辨率发生变化
df = fs / len;res=res(1:len);
freResponse = fft(res);
ampResponse = abs(freResponse);
% 显示频谱， 去除0频后，最大的频率分量就是码元速率的粗估计
plot((0:length(ampResponse)-1)*df, ampResponse);
% 去除0频
a = ampResponse(1);
ampResponse(1) = 0;
[~, index] = max(ampResponse);
ampResponse(1) = a;
% 这里有个细节，注意matlab是从1到index开始取得，这样的话就会错过0频，
% 所以我们的index转频率的时候，要记得-1
coarseSymbolRate = (index-1) * df% 细化频谱
% 设置细化倍数，需要设置为2的n次幂
m = 64;
n0 = len / m;% 移动右侧频谱到0频的右边
n02 = n0 / 2;
right = index + n02;
tmp = right + 1 - index;
assert(tmp == n02+1);
freResponse(1:tmp) = freResponse(index:right);
% 左侧频谱到负频率
left = index - n02 + 1;
tmp = n02 - 1;
assert(tmp == (index - left));
freResponse((length(freResponse) - tmp+1):end) = freResponse(left:(index-1));
freResponse(right+1:end-tmp) = 0;% ifft
timeSignal = ifft(freResponse);
% 下采样
downTimeSignal = downsample(timeSignal, m);
downTimeSignalFreqResponse = fft(downTimeSignal, len);
leftFreq = (left-1) * df;
rightFreq = (right-1) * df;
downTimeSignalFreqResponse = circshift(downTimeSignalFreqResponse, floor(len/2));
f = linspace(leftFreq, rightFreq, len);
ampDownTimeSignalFreqResponse = abs(downTimeSignalFreqResponse);
figure()
plot(f, ampDownTimeSignalFreqResponse);
[~, index] = max(ampDownTimeSignalFreqResponse);
df = (rightFreq - leftFreq) / len;
accurateSr = (index - 1) * df + leftFreq
% 频谱图中最大点就是结果

码元速率估计-速率信号法相关推荐

码元速率 matlab,[转载]基于小波变换的移相键控信号符号速率估计（matlab仿真）...
西安电子科技大学西电大宝在电子技术迅猛发展的当代社会,空间中充满了各种各样不同频率.不同调制类型的通信信号.正常通信条件下,发送方和接受方进行的是合作通信,即接受方预先知道发送信号的频率.调制类 ...
【计算机网络】物理层 : 相关参数 ( 码元 | 码元进制 | 速率 | 码元传输速率 | 波特 | 信息传输速率 | 带宽 | 码元速率计算示例 )★
文章目录一.码元概念二.码元进制三.速率四.码元传输速率与波特五.信息传输速率六.带宽七.码元速率计算一.码元概念 "码元" 概念 : 使用固定时长的信号 ...
wcdma码片速率_WCDMA中的码片速率,符号速率,信息速率(bit rate)之间的关系
很久前我写过码元(码片)速率和信息速率的关系.这里在WCDMA中对应一下它们之间的关系,从网上抄了一段: 比特(Bit),符号(Symbol) ,码片(Chip)矩阵通信技术论坛--3G论坛| NGN ...
国家电网一二次融合配电终端(FTU)发展新方向馈线自动化:具备集中型馈线自动化/就地型馈线自动化包括电压时间型、电压电流型、自适应综合型)及零序电流，零序电压或外施信号法的单相接地故障选线功能
近日国家电网一种新型一二次融合馈线终端DAF-810馈线终端产成功上线挂网运行.DAF-810新型一二次融合终端馈线自动化终端具备集中型馈线自动化.就地型馈线自动化(包括电压时间型.电压电流型.自适应 ...
码元速率的盲估计-包络谱法
码元速率的盲估计实现自论文主要实现的是包络谱的检测方法. 算法的大致原理请看论文,本文这里是帮助实现了论文代码.分享出来. 源代码 clc clear close all % 这是一个对码元速率的 ...
基于代价函数小波脊相位的MFSK信号符号速率估计MATLAB仿真及代码（2020.12.14更新）
算法来源王勇, 王李福, 邹辉,等. 一种小波脊相位提取方法: 中国专利. 仿真结果引言当前,脊点的选取较准确也较经典的方法是Liu等提出的里程碑式的基于代价函数的小波脊相位提取方法(即代价函数 ...
关于视频中的速率问题，海思视频速率（高速信号与高频信号区分与解释）
写在前面,因这篇文档是很早之前整理的,忘记在哪里看到的别人总结好的,我稍作修改整理.后面关于海思视频处理过程及速率是我自己写的. 高速信号与高频信号区分与解释 Candence公司对高速信号的判断: ...
matlab 多速率,多速率信号处理系统设计与实现.PDF
多速率信号处理系统设计与实现.PDF 第26卷第18期电子设计工程 2018年9月 Vol.26 No.18 ElectronicDesignEngineering Sep.2018 多速率信号处 ...
电路原理 | 非线性电阻电路的小信号法，动态电路电感
非线性电路的小信号分析方法解题步骤 (1) 求解非线性电路的静态工作点; (2) 求解非线性电阻元件在静态工作点处的动态电导或动态电阻; (3) 作出给定的非线性电阻在静态工作点处的小信号等效电路; ...

码元速率估计-速率信号法

速率信号法

思想

源代码

码元速率估计-速率信号法相关推荐

最新文章

热门文章