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