主要内容

测量信号的功率

信号的功率是其时域样本的绝对平方和除以信号长度,或等效地,它的均方根水平的平方。这个函数bandpower允许您一步估计信号功率。

考虑一个嵌入在高斯白噪声中的单位啁啾,并以1千赫采样1.2秒。啁啾的频率在一秒钟内从初始值100hz增加到300hz。噪声有方差 0 0 1 2 .重置随机数发生器可重现的结果。

N = 1200;Fs = 1000;t = (0:N-1)/Fs;σ = 0.01;rng (“默认”) s = chirp(t,100, 1300)+sigma*randn(大小(t));

验证所给出的功率估计bandpower等于定义。

pRMS = rms(s)^2
pRMS = 0.5003
powbp =带功率(s,Fs,[0 Fs/2])
Powbp = 0.5005

使用obw函数来估计包含信号99%功率的频带宽度、频带的下界和上界以及频带中的功率。该函数还绘制了频谱估计,并标注了占用的带宽。

obw(年代,Fs);

图中包含一个轴对象。标题为99%已占用带宽:241.242 Hz的axis对象包含4个类型为line, patch的对象。

[wd,lo,hi,power] = obw(s,Fs);功率=功率/0.99
波托= 0.5003

一个非线性功率放大器输入一个60hz的正弦波,输出一个三阶失真的噪声信号。信号在3.6千赫采样2秒。

Fs = 3600;t = 0:1/Fs:2-1/Fs;X = sin(2* *60*t);Y = polyval(ones(1,4),x) + randn(size(x));

由于放大器引入了三阶失真,因此输出信号预期具有:

  • 一个基本与输入频率相同的组件,60hz。

  • 两个谐波-两倍和三倍于输入频率的频率组件,120和180 Hz。

使用函数可视化信号的频谱。标注基波和谐波。

(thd (y, Fs);

图中包含一个轴对象。标题为THD: -9.56 dB的axes对象包含16个类型为line, text的对象。这些对象代表基本,谐波,直流和噪声(不包括)。

使用bandpower来确定存储在基波和谐波中的功率。表示每个值占总功率的百分比,单位为分贝。将值显示为表。

pwrTot = bandpower(y,Fs,[0 Fs/2]);谐波= {“基本”“第一”“第二”};频率= [60 120 180]';功率= 0 ([3 1]);k = 1:3权力(k) = bandpower (y, Fs,频率(k) + -10 [10]);结束百分比= Power/pwrTot*100;inDB = pow2db(功率);T =表(频率,功率,百分比,inDB,“RowNames”谐波)
T =3×4表频率功率百分比inDB _____ _______ _______ _______基本60 1.5777 31.788 1.9804第一120 0.13141 2.6476 -8.8137第二180 0.04672 0.9413 -13.305

另请参阅

|||

相关的话题