主要内容

sinad

信号噪声和失真率

描述

例子

r= sinad (x)返回信号的噪声和失真率(SINAD) dBc实值的正弦信号x。确定了SINAD使用修改后的长度相同的周期图作为输入信号。修改后的周期图使用Kaiser窗β= 38

例子

r= sinad (x,fs)指定采样率fs的输入信号x。如果你不指定fs,然后采样率的缺省值为1。

例子

r= sinad (pxx,f,' psd ')指定的输入pxx作为一个单边功率谱密度(PSD)估计。f是一个向量的频率对应的PSD估计pxx

r= sinad (sxx,f,rbw“权力”)指定输入单边功率谱。rbw的分辨率带宽每个功率估计是集成。

(r,totdistpow)= sinad (___)返回总噪声和谐波失真功率(dB)的信号。

例子

sinad (___)没有输出参数块信号的频谱图窗口,在当前标签的基本组件。它使用不同的颜色来画的基本组件,直流值和噪声。上面的SINAD出现情节。

例子

全部折叠

创建两个信号。两个信号的基频 π / 4 rad /样本与振幅1和第一次谐波频率 π / 2 rad 0.025 /样本与振幅。另外的一个信号加性高斯白噪声的方差 0 0 5 2

创建两个信号。将随机数字生成器设置为默认设置为可再生的效果。确定信号的SINAD没有添加剂噪声理论SINAD和比较结果。

n = 0:159;x = cos(π/ 4 * n) + 0.025 * sin(π/ 2 * n);rng默认的y = cos(π/ 4 * n) + 0.025 *罪(π/ 2 * n) + 0.05 * randn(大小(n));r = sinad (x)
r = 32.0412
powfund = 1;powharm = 0.025 ^ 2;thSINAD = 10 * log10 (powfund / powharm)
thSINAD = 32.0412

确定SINAD为加性噪声的正弦信号。展示包括理论添加剂噪声的方差近似于SINAD。

r = sinad (y)
r = 22.8085
varnoise = 0.05 ^ 2;thSINAD = 10 * log10 (powfund / (powharm + varnoise))
thSINAD = 25.0515

创建一个信号的基频1 kHz和单位振幅,在480千赫采样。第一次谐波振幅的信号还包括0.02和加性高斯白噪声的方差 0 0 1 2

确定与理论SINAD SINAD和比较结果。

fs = 48 e4;t = 0:1 / fs: 1 - 1 / f;rng默认的x = cos(2 *π* 1000 * t) + 0.02 *罪(2 *π* 2000 * t) + 0.01 * randn(大小(t));r = sinad (x, fs)
r = 32.2058
powfund = 1;powharm = 0.02 ^ 2;varnoise = 0.01 ^ 2;thSINAD = 10 * log10 (powfund / (powharm + varnoise * (1 / f)))
thSINAD = 33.9794

创建一个信号的基频1 kHz和单位振幅,在480千赫采样。第一次谐波振幅的信号还包括0.02和加性高斯白噪声标准差0.01。将随机数字生成器设置为默认设置为可再生的效果。

获取信号的周期图和使用周期图作为输入sinad

fs = 48 e4;t = 0:1 / fs: 1 - 1 / f;rng默认的x = cos(2 *π* 1000 * t) + 0.02 *罪(2 *π* 2000 * t) + 0.01 * randn(大小(t));[pxx f] =周期图(x, rectwin(长度(x)),长度(x), fs);r = sinad (pxx, f,psd的)
r = 32.2109

生成一个正弦信号频率2.5 kHz的采样在50千赫。添加高斯白噪声信号标准偏差为0.00005。结果通过一个弱非线性放大器。绘制SINAD。

fs = 5 e4;f0 = 2.5 e3;N = 1024;t = (0: n - 1) / fs;ct = cos(2 *π* f0 * t);cd = ct + 0.00005 * randn(大小(ct));amp = [1 e-5 5 e-6 1 e - 3 6 e-5 1 25 e - 3);胡志明市= polyval (amp、cd);sinad(胡志明市,fs);

图包含一个坐标轴对象。与标题SINAD坐标轴对象:72.10 dB包含7线类型的对象,文本。这些对象代表基本,噪音和失真,直流(排除)。

情节展示了频谱用于计算比例和该地区视为噪音。直流噪声水平和基本被排除在计算。基本是标记。

输入参数

全部折叠

实值正弦输入信号,指定行或列向量。

例子:因为(π/ 4 * (0:159))+ cos(π/ 2 * (0:159))

数据类型:|

采样率,指定为一个积极的标量。采样率是单位时间内样品的数量。如果时间的单位是秒,那么采样率的单位是赫兹。

片面的PSD估计,指定为一个实值,非负列向量。

功率谱密度必须被表达在线性单元,不是分贝。使用db2pow分贝值转换为功率值。

例子:[pxx f] =周期图(cos (pi. / (4; 2) * (0:159)) ' + randn (160 2))指定嘈杂的双通道的周期图PSD估计正弦信号采样2π赫兹的频率计算。

数据类型:|

周期性的频率对应于片面的PSD估计,pxx,指定行或列向量。的第一个元素f必须是0。

数据类型:|

功率谱,指定为一个实值负的行或列向量。

功率谱必须被表达在线性单元,不是分贝。使用db2pow分贝值转换为功率值。

例子:[sxx, w] =周期图(cos (pi. / (4; 2) * (0:159))”+ randn(160 2),“权力”)指定一个双通道的周期图功率谱估计正弦信号嵌入在高斯白噪声和归一化频率的计算。

分辨率带宽,指定为一个积极的标量。分辨率带宽的乘积的频率分辨率离散傅里叶变换和窗口的等效噪声带宽。

输出参数

全部折叠

信号噪声和失真比在dBc,作为实值返回标量。

噪音和总谐波失真的信号,作为实值返回标量用分贝表示。

更多关于

全部折叠

变形测量功能

的函数,sfdr,sinad,信噪比衡量一个弱非线性系统的响应由正弦信号刺激。

当给定时域输入,sinad执行周期图使用Kaiser窗旁瓣衰减大。找到的基本频率,算法搜索最大的非零谱分量的周期图。然后计算中央的所有相邻的垃圾箱,减少单调远离最大。可以被探测到,基本应该至少在第二频率。更高的谐波基频的整数倍。如果一个谐波在单调递减附近的另一个地区,它的力量被认为属于较大的谐波。这个大的谐波可能是也可能不是最基本的。

函数估计噪声级使用的中间力量区域只包含噪声和失真。直流分量是被排除在计算。每个点的噪声估计的水平或点的纵坐标,哪个更小。然后减去噪声从信号的值和谐波。

sinad失败如果根本不是最高的光谱分量信号。

确保频率成分相距足够远的旁瓣宽度来适应凯撒窗口。如果这是不可行的,你可以使用“权力”国旗和计算一个周期图用不同的窗口。

扩展功能

版本历史

介绍了R2013b