基本原理
在实际工程上,我们一般采用正弦波来发射、接收信号。但是由于我们采集到的是离散的点,所以不能通过找信号的最大值来代替信号的幅度值。因此,如何从我们接收到的正弦波序列来估计它的幅度,成为我们研究的关键。
一般我们有两种方法来估计信号的幅度。
第一种是利用能量公式得到目标信号的能量,进而得到目标的功率,最后得到目标的幅度。
第二种是对正弦波序列做离散傅里叶变换,得到该序列的幅度谱。在幅度谱中找到信号的频率所对应的谱峰。然后把谱峰的值除以序列的长度,就得到序列的幅度。
源代码
% 已知正弦波序列,如何估计他的幅度?
close all;
clear all;
clc;
fc=1e3;
fs=1e5;
t=0:1/fs:1/fc-1/fs;
% y=sin(2*pi*fc*t);
y=exp(1i*2*pi*fc*t);
figure(1)
plot(t,y);
% 法一 利用E=sum(x(n)*x(n))计算序列的能量
E1=y*y'; % 计算序列的能量
P1=E1/length(t);
u1=sqrt(P1); % 得到信号的幅度值
fprintf('时域算出来的幅度值为%g\n',u1);
%法二 利用FFT
Y=fftshift(fft(y));
L=length(y);
P=abs(Y)/L; % 幅度谱
f=(0:L-1)*fs/L-fs/2;
figure(2)
plot(f,P);% 幅度谱
index = find(f==fc);
fprintf('FFT算出来的幅度值为%g\n',P(index));
% text(f(index),P(index),num2str(P(index)));
% [mp,mi] = max(P)