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

@[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)={\int }_{-\infty }^{\infty }s(t){\varphi }^{\ast }(\frac{t-b}{a})dt$$
  其中，是$s(t)$为平方可积信号。

2.码元速率估计过程

1）产生haar小波，小波长度为一个码元长度;
2) 小波变换，识别码元跳变位置；
3）傅里叶变换，计算码元速率，频谱中的第一个谱峰值位置即为所求的码元速率。

二、Matlab仿真

1.代码

%------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 timesymbol = 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 Samplesphi_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小波的参数影响。