主要内容

sfdr公司

无杂波动态范围

说明

例子

R=sfdr()返回无杂散动态范围(SFDR),R,单位为实际正弦信号的dB,.sfdr公司使用修正的周期图和Kaiser窗口计算功率谱β=38.平均值从在计算功率谱之前。离散傅里叶变换(DFT)计算中使用的点数与信号长度相同,.

R=sfdr(,fs公司)返回时域输入信号的SFDR,,当采样率,fs公司,已指定。的默认值fs公司是1赫兹。

例子

R=sfdr(,fs公司,材料安全数据表)返回SFDR,仅考虑与基本(载波)频率相隔最小杂散距离的杂散,材料安全数据表,以周期/单位时间指定。采样率为fs公司. 如果载波频率是Fc公司,然后中场休息时所有的马刺(Fc物料安全数据表,Fc+msd公司)被忽略。

例子

R=sfdr(sxx公司,F(“电源”)返回实值信号单边功率谱的SFDR,sxx公司.F是与中的功率估计相对应的频率向量sxx公司. 第一个元素F必须等于0。该算法去除了所有从直流仓单调递减的功率。

R=sfdr(sxx公司,F,材料安全数据表(“电源”)返回SFDR,仅考虑与基本(载波)频率相隔最小杂散距离的杂散,材料安全数据表. 如果载波频率是Fc公司,然后中场休息时所有的马刺(Fc物料安全数据表,Fc+msd公司)被忽略。当输入到sfdr公司是功率谱,指定材料安全数据表可以防止高旁瓣电平被识别为杂散。

例子

[R,马刺,杂散频率]=sfdr(___)返回最大杂散的功率和频率。

例子

sfdr公司(___)在没有输出参数的情况下,在当前图形窗口中绘制信号的频谱。它使用不同的颜色来绘制基本成分、DC值和光谱的其余部分。它对SFDR进行着色,并在绘图上方显示其值。它还标明了最基本的和最大的丁坝。

示例

全部折叠

获取10 MHz音调的SFDR,振幅1在100 MHz下采样。在第一次谐波(20 MHz)处有一个振幅为的杂散 3. . 1. 6. × 1. 0 - 4. .

德尔塔=1e-8;fs=1/德尔塔;t=0:deltat:1e-5-deltat;x=cos(2*pi*10e6*t)+3.16e-4*cos(2*pi*20e6*t);r=sfdr(x,fs)
r=70.0063

显示信号的频谱。注释基波、直流值、杂散和SFDR。

sfdr(x,fs);

图中包含一个轴。标题为SFDR:70.01 dB的轴包含9个类型为patch、line和text的对象。这些对象表示SFDR、基本、杂散、DC(不包括)。

获取10 MHz音调的SFDR,振幅1在100 MHz下采样。在第一次谐波(20 MHz)处有一个振幅为的杂散 3. . 1. 6. × 1. 0 - 4. 另一个频率为25兆赫,振幅为 1. 0 - 5. . 通过使用11兆赫的最小杂散距离跳过一次谐波。

德尔塔=1e-8;fs=1/德尔塔;t=0:deltat:1e-5-deltat;x=cos(2*pi*10e6*t)+3.16e-4*cos(2*pi*20e6*t)+...0.1e-5*cos(2*pi*25e6*t);r=sfdr(x,fs,11e6)
r=120.0000

显示信号的频谱。注释基波、直流值、杂散和SFDR。

sfdr(x,fs,11e6);

图中包含一个轴。标题为SFDR:120.00 dB的轴包含9个类型为patch、line和text的对象。这些对象表示SFDR、基本、杂散、DC(不包括)。

获得在100 MHz下采样的振幅为1的10 MHz音调的功率谱。在第一次谐波(20 MHz)处有一个振幅为的杂散 3. . 1. 6. × 1. 0 - 4. . 使用单侧功率谱和相应频率的矢量(单位:Hz)计算SFDR。

德尔塔=1e-8;fs=1/德尔塔;t=0:deltat:1e-6-deltat;x=cos(2*pi*10e6*t)+3.16e-4*cos(2*pi*20e6*t)[sxx,f]=周期图(x,rectwin(长度(x)),长度(x),fs,'电源'); r=sfdr(sxx,f,'电源');

显示信号的频谱。注释基波、DC值、第一个支路和SFDR。

sfdr(sxx,f,'电源');

图中包含一个轴。标题为SFDR:70.01 dB的轴包含9个类型为patch、line和text的对象。这些对象表示SFDR、基本、杂散、DC(不包括)。

确定最大杂散的频率(以MHz为单位)。输入信号是一个10兆赫音调,振幅1在100兆赫采样。在一次谐波(20 MHz)处有一个振幅为的杂散 3. . 1. 6. × 1. 0 - 4. .

德尔塔=1e-8;t=0:deltat:1e-6-deltat;x=cos(2*pi*10e6*t)+3.16e-4*cos(2*pi*20e6*t)[r、 spurpow,spurfreq]=sfdr(x,1/deltat);杂散频率=杂散频率/1e6
杂散频率=20

在白高斯加性噪声中创建频率为9.8、14.7和19.6 kHz的三个正弦叠加。信号采样频率为44.1 kHz。9.8千赫正弦波的振幅为1伏,14.7千赫正弦波的振幅为100微伏,19.6千赫信号的振幅为30微伏。噪声均值为0,方差为0.01微伏。此外,该信号的直流偏移为0.1伏。

rng公司违约Fs=44.1e3;f1=9.8e3;f2=14.7e3;f3=19.6e3;N=900;nT=(0:N-1)/Fs;x=0.1+sin(2*pi*f1*nT)+100e-6*sin(2*pi*f2*nT)...+30e-6*sin(2*pi*f3*nT)+sqrt(1e-8)*randn(1,N);

绘制信号的频谱和SFDR。显示其基本和最大的刺。直流电平不包括在SFDR计算中。

sfdr(x,Fs);

图中包含一个轴。标题为SFDR:79.52 dB的轴包含9个类型为patch、line和text的对象。这些对象表示SFDR、基本、杂散、DC(不包括)。

输入参数

全部折叠

实值正弦信号,指定为行或列向量。从中减去平均数在获得SFDR计算的功率谱之前。

例子:x=cos(π/4*(0:79))+1e-4*cos(π/2*(0:79));

数据类型:双重的

以周期/单位时间为单位的信号采样率,指定为正标量。当时间单位是秒时,fs公司单位为赫兹。

数据类型:双重的

SFDR计算中要忽略的最小离散傅里叶变换(DFT)箱数,指定为正标量。您可以使用此参数来忽略发生在基频附近的杂散或旁瓣。例如,如果载波频率Fc公司,然后所有马刺都在射程内(Fc msd,Fc+msd)被忽略。

数据类型:双重的

用于SFDR计算的单侧功率谱,指定为行或列向量。

功率谱必须用线性单位表示,而不是分贝。使用db2pow公司将分贝值转换为幂次值。

例子:[sxx,f]=周期图(cos(pi./[4;2] *(0:159))'+随机数(160,2),'功率')指定以2πHz采样的噪声双通道正弦波的周期图功率谱估计及其计算频率。

数据类型:双重的

与中的功率估计相对应的频率向量sxx公司,指定为行或列向量。

输出参数

全部折叠

无杂散动态范围(dB),指定为实值标量。无杂散动态范围是峰值频率处的功率与下一最大频率处的功率(杂散)之间的差值,单位为dB。如果输入是时间序列数据,则使用汉明窗从修正的周期图获得功率估计。周期图中使用的DFT长度等于输入信号的长度,. 如果要使用不同的功率谱作为SFDR测量的基础,可以使用'电源'旗帜。

数据类型:双重的

以dB为单位的最大杂散的功率,指定为实值标量。

数据类型:双重的

最大杂散的频率,单位为赫兹,用实值标量表示。如果不提供采样率作为输入参数,sfdr公司假设采样率为1 Hz。

数据类型:双重的

更多关于

全部折叠

失真测量功能

功能总谐波失真,sfdr公司,锡那德,和信噪比测量受正弦激励的弱非线性系统的响应。

当给定时域输入时,sfdr公司使用具有大旁瓣衰减的Kaiser窗口执行周期图。为了找到基频,该算法在周期图中搜索最大的非零谱分量。然后,它计算所有相邻箱子的中心力矩,这些箱子从最大值开始单调减小。为了能够被检测到,基波至少应该在第二个频率槽中。如果一个谐波位于另一个谐波附近的单调递减区域内,则认为其功率属于较大的谐波。这个较大的谐波可能是基波,也可能不是。该算法忽略了所有的功率单调减少远离直流箱。

sfdr公司如果基波不是信号中的最高频谱分量,则会失败。

确保频率分量相距足够远,以适应Kaiser窗口的旁瓣宽度。如果不可行,可以使用'电源'标记并计算具有不同窗口的周期图。

扩展功能

R2013a引入