符号速率估计——小波变换法

@[TOC]符号速率估计------小波变换法

一、原理

1.Haar小波变换

小波变换在信号处理领域被成为数学显微镜,不同于傅里叶变换,小波变换可以观测信号随时间变换的频谱特征,因此,常用于时频分析。

当小波变换前后位置处于同一个码元内或相邻两个码元相同时,其幅度是不变的,当相邻的码元不同时,小波变换的幅度值会响应变化的,其幅值常与相邻码元的幅度、相位和频率的变换有关,Haar小波可以识别出码元跳变位置,因此,可以进一步估计出数字信号的码元速率。Haar小波定义如下:
ϕ a , b ( t ) = { − 1 / a , 0 < = t < a / 2 1 / a , − a / 2 < t < 0 \phi_{a,b}(t)=\{_{-1/\sqrt{a},0<=t<a/2}^{1/\sqrt{a},-a/2<t<0} ϕa,b(t)={−1/a ,0<=t<a/21/a ,−a/2<t<0

小波变换公式如下:
W s ( a , b ) = ∫ − ∞ ∞ s ( t ) ϕ ∗ ( t − b a ) d t W_{s}(a,b)=\int_{-\infin}^{\infin}s(t)\phi^{*}(\frac{t-b}{a})dt Ws(a,b)=∫−∞∞s(t)ϕ∗(at−b)dt

其中,是 s ( t ) s(t) s(t)为平方可积信号。

2.码元速率估计过程

1)产生haar小波,小波长度为一个码元长度;

  1. 小波变换,识别码元跳变位置;

3)傅里叶变换,计算码元速率,频谱中的第一个谱峰值位置即为所求的码元速率。

二、Matlab仿真

1.代码

c 复制代码
%------Function:Estimate Symbol Rate--------
%------Remark:Using Wavelet Transform-------
%------Time:2025.04.11----------------------
%------Author:Clemence----------------------
clc;
close all;
clear all;

%------------------1.Paras----------------------------------------------
S_N = 100;      % Symbol Num
nSamp = 400;     % Samples per Symbol
N = S_N*nSamp;  % Sample Num
M = 2;          % Modulate index
f = 30;        % Signal Frequence
fs = 2000;      % Sample Frequence
ts = 1/fs;      % Sample Interval
t = 0:ts:ts*(N-1);  % Sample time

symbol = randi([0 M-1],S_N,1);  % Symbol
symbol = 2*symbol-1;            % 
symbrate = fs/nSamp;            % Symbol Rate

%------------------2.Create 2FSK Signal------------------------------------
for i=1:S_N     
    s((i-1)*nSamp+1:i*nSamp) = cos((2*pi*f+pi*symbol(i))*t((i-1)*nSamp+1:i*nSamp));
end

%-----------------3.Create Wavelet-----------------------------------------
Tb = nSamp*ts;       % Wavelet length
t = -Tb/2:ts:Tb/2;   % Wavelet Samples

phi_t = (t>-Tb/2 & t<0)-(t>0 & t<-Tb/2);


%-----------------4. Wavelet Transform-------------------------------------
for i = 1:length(s)-length(phi_t)
   x = s(i:i+length(phi_t)-1); 
   y = x.*phi_t;
   sum_y(i) = sum(y);
end

%-----------------5. plot -------------------------------------------------
figure;
subplot(2,2,1)
plot(s);
title('Original Signal')

subplot(2,2,2)
plot(t,phi_t);
title('Haar Wavelet')

subplot(2,2,3)
plot(sum_y);
title('Wavelet Transform Result')

n = length(sum_y)
inx = (0:n-1)*fs/n;

subplot(2,2,4)
plot(inx,abs(fft(abs(sum_y))));
title('Wavelet Transform Frequence')

2.仿真结果


符号速率估计结果

三、总结

从仿真结果可以看出,小波变换法估计出的符号速率为5.00126,与信号的实际符号速率5基本接近,同时,符号速率估计的精度受信号的频率以及Haar小波的参数影响。

相关推荐
Yuroo zhou26 分钟前
无人机在复杂气流中,IMU 如何精准捕捉姿态变化以维持稳定?
单片机·嵌入式硬件·算法·机器人·无人机
xueyongfu31 分钟前
PTX指令集基础以及warp级矩阵乘累加指令介绍
人工智能·线性代数·算法·矩阵
流星白龙35 分钟前
【C++算法】89.多源BFS_01 矩阵
c++·算法·宽度优先
爱代码的小黄人3 小时前
利用劳斯判据分析右半平面极点数量的方法研究
算法·机器学习·平面
今天也好累6 小时前
C 语言基础第16天:指针补充
java·c语言·数据结构·笔记·学习·算法
大千AI助手6 小时前
直接偏好优化(DPO):原理、演进与大模型对齐新范式
人工智能·神经网络·算法·机器学习·dpo·大模型对齐·直接偏好优化
徐小夕8 小时前
再也不怕看不懂 GitHub 代码!这款AI开源项目,一键生成交互架构图
前端·算法·github
风靡晚9 小时前
雷达系统设计学习:自制6GHz FMCW Radar
信息与通信·信号处理·雷达·fmcw
SirLancelot19 小时前
数据结构-Set集合(一)Set集合介绍、优缺点
java·开发语言·数据结构·后端·算法·哈希算法·set
YouQian7729 小时前
label 拓扑排序
数据结构·算法