【单音频率测量】通过DFT、MLE和加权平均进行单音频率测量（估计）（Matlab代码实现）

 👨‍🎓个人主页：研学社的博客    

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

文献来源：

 原文摘要：

Abstract: - The present paper addresses the estimation of the frequency of a single-tone sine wave signal corrupted by additive white Gaussian noise. Brief background information is given about the signal model and the Maximum Likelihood Estimation (MLE) method. Further the newly proposed method is described, involving two-step coarse and fine search of the frequency index that maximized the signal modified periodogram, based on MLE and Weighted Averaging Estimation, respectively. Computer simulations are performed to compare the mean-square-error of the proposed estimator with that of some other frequency estimators along with Cramér-Rao lower bound. Keywords: - single-tone, frequency estimation, method.

本文用于测量单音信号频率。该测量基于信号的离散傅里叶变换（DFT）和两步估计程序，涉及经典的最大似然（ML）粗略估计和作者对频率指数的加权平均（WA）更精细的估计，以最大化修改后的信号周期图。该方法可用于任何类型的稳态周期信号（例如，正弦波、矩形、三角形等）的频率估计。为了阐明函数的用法，给出了一个示例。为方便起见，输入和输出参数在函数的开头给出。

📚2 运行结果

主函数：

clear, clc, close all

% signal parameters
fs = 512;
f0 = 137.365;
T = 0.5;
N = round(T*fs);
SNR = 60;

% signal generation
t = (0:N-1)/fs;
s = sin(2*pi*f0*t);
n = std(s)/(10^(SNR/20))*randn(1, N);
x = s + n;

% signal frequncy estimation
f0tild = frequency(x, fs);

% error of the estimation
mse = (f0tild - f0)/f0*100;

% display the results
f0hatstr = num2str(f0tild);
msestr = num2str(mse);
disp(['Estimated frequency = ' f0hatstr])
disp(['Error of the estimation = ' msestr, ' %'])

commandwindow

运行结果：

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1] H. Zhivomirov, I. Nedelchev, R. Vasilev. A Method for Single-Tone Frequency Estimation. Romanian Journal of Acoustics and Vibration, ISSN: 1584-7284, Vol. XIII, No. 1, pp. 20-24, 2016. 

🌈4 Matlab代码实现