信号的功率是它的时域样本的绝对平方和除以信号长度,或者,等价地,它的均方根水平的平方。这个函数bandpower
让您估计信号功率在一个步骤。
考虑一个单位啁啾嵌入在高斯白噪声中,以1khz采样1.2秒。啁啾的频率在一秒钟内从初始值100 Hz增加到300 Hz。噪声有方差 .重置随机数生成器以获得可重复的结果。
N = 1200;Fs = 1000;t = (0: n - 1) / Fs;σ= 0.01;rng (“默认”) s = chirp(t,100, 1300)+sigma*randn(size(t));
验证由给出的功率估计bandpower
等价于定义。
人口、难民和移民事务局= rms (s) ^ 2
人口、难民和移民事务局= 0.5003
pbp =带功率(s, f,[0 f /2])
powbp = 0.5005
使用obw
函数估计包含99%信号功率的频带宽度,频带的下界和上界,以及频带中的功率。该函数还绘制频谱估计图并标注所占用的带宽。
obw(年代,Fs);
[wd,嗨,权力]= obw(年代,Fs);powtot / 0.99 =权力
powtot = 0.5003
非线性功率放大器的输入为60hz正弦波,输出为三阶失真的噪声信号。信号以3.6 kHz采样2秒。
Fs = 3600;t = 0:1 / Fs: 2 - 1 / f;x =罪(2 *π* 60 * t);Y = polyval(ones(1,4),x) + randn(size(x));
由于放大器引入了三阶失真,输出信号预期有:
一个基本输入频率为60hz。
两个谐波—输入频率的2倍和3倍(120hz和180hz)的频率分量。
使用野
函数可视化信号的频谱。注释基本音和谐波。
(thd (y, Fs);
使用bandpower
确定存储在基波和谐波中的功率。以总功率的百分比和分贝表示每个值。将值显示为表。
pwrTot = bandpower(y, f,[0 f /2]); / /输出功率谐波= {“基本”;“第一”;“第二”};频率= [60 120 180]';Power = 0 ([3 1]);为k = 1:3权力(k) = bandpower (y, Fs,频率(k) + -10 [10]);结束% = / 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