在Simulink - MATLAB和Simulink - M金宝appathWorks中估计功率谱 - 金宝app,下载188bet金宝搏,金宝搏官方网站

估算中功率谱金宝app

时域信号的功率谱(PS)是基于有限数据集的信号中包含的功率在频率上的分布。信号的频域表示通常比时域表示更容易分析。许多信号处理应用，如噪声消除和系统识别，都是基于信号的特定频率修改。功率谱估计的目的是从一个序列的时间样本中估计信号的功率谱。根据对信号的了解，估计技术可以包括参数或非参数方法，并可以基于时域或频域分析。例如，一种常见的参数技术包括将观测值拟合到自回归模型中。一种常见的非参数技术是周期图。采用傅立叶变换方法进行功率谱估计，如韦尔奇法和滤波器组法。对于长度相对较小的信号，滤波器组方法产生的频谱估计具有更高的分辨率，更准确的噪声底噪声，峰值比Welch方法更精确，且频谱泄漏低或没有。这些优势是以增加的计算量和较慢的跟踪速度为代价的。 For more details on these methods, see光谱分析．您还可以使用其他技术，如最大熵方法。

在Si金宝appmulink.^®，您可以使用该动态信号进行实时频谱分析频谱分析仪块。您可以在Spectrum Analyzer中查看光谱数据。要获取进一步处理的最后一个光谱数据，请创建一个SpectrumanalyZerconFiguration对象并运行getspectrumdata.函数。或者，你可以用谱估计块的dspspect3库来计算功率谱，以及数组图块来查看光谱。

用频谱分析仪估算功率谱

您可以使用该信号查看信号的功率谱（PS）频谱分析仪块。PS是实时计算的，并随输入信号的变化和器件特性的变化而变化频谱分析仪块。你可以改变输入信号的动态，并实时观察这些变化对信号频谱的影响。

该模型ex_psd_sa.向该时馈送嘈杂的正弦波信号频谱分析仪块。正弦波信号是两个正弦波的和:一个频率为5000赫兹，另一个频率为10000赫兹。输入处的噪声为高斯分布，均值为零，标准差为0.01。

打开并检查模型

要打开模型，输入ex_psd_sa.在MATLAB^®命令提示符。

下面是模型中块的设置。

块	参数更改	块的目的
正弦波1.	频率到5000年采样时间到1/44100. 样品每帧到1024年	频率为5000hz的正弦信号
正弦波2.	频率到10000. 相抵消(rad)到10 采样时间到1/44100. 样品每帧到1024年	频率为10000hz的正弦信号
随机源	源类型到`高斯` 方差到1E-4 采样时间到1/44100. 样品每帧到1024年	随机源块生成具有通过块对话框指定属性的随机噪声信号
添加	标志名单到`+++`．	添加块给输入信号增加随机噪声
频谱分析仪	单击光谱设置图标．右边出现一个窗格。在里面主要选项窗格中,在类型，选择`权力`．下方法，选择`滤波器组`．在里面跟踪选项窗格中,清除双边频谱复选框。这只显示了光谱的实半部分。如果需要，选择Max-hold跟踪和闵持痕迹复选框。单击配置属性图标和集合y-limits（最低）作为`－100`和Y-limits(最大)作为`40`．	频谱分析仪块显示信号的功率谱密度

播放模式。打开频谱分析仪块来查看正弦波信号的功率谱。有两个频率为5000赫兹和10000赫兹的音调，它们对应于输入端的两个频率。

RBW时，分辨率带宽是频谱分析仪能分辨的最小频率带宽。默认情况下,RBW(赫兹)被设置为汽车．在里面汽车RBW是频率跨度与1024的比值。在双边谱中，这个值是F_年代/ 1024但从单方面的角度来看，事实的确如此（F_年代/ 2）/ 1024．频谱分析仪ex_psd_sa.配置为显示单侧频谱。因此，RBW是(44100/2)/1024或21.53 Hz。

使用该值的RBW，用于计算一个光谱更新的输入样本的数量由N_样品= FS / RBW，即44100/21.53或2048。

在此模式下计算的RBW具有良好的频率分辨率。

为了区分显示器中的两个频率，两个频率之间的距离必须至少是RBW。在该示例中，两个峰之间的距离为5000Hz，其大于RBW。因此，您可以清楚地看到峰值。将第二正弦波从10000 Hz变为5015Hz的频率。两个频率之间的差异小于RBW．

放大后，你可以看到峰值是无法区分的。

要提高频率分辨率，就要降低频率RBW为1hz，并运行模拟。

在放大时，两个相隔15赫兹的峰值现在可以区分了

当您增加频率分辨率时，时间分辨率会减小。为了在频率分辨率和时间分辨率之间保持良好的平衡，更改RBW(赫兹)到汽车．

改变输入信号

当在仿真过程中改变输入信号的动态时，信号的功率谱也会实时发生变化。在模拟运行时，更改频率的正弦波1.块8000并点击应用．光谱分析仪中的第二个音调输出转移到8000 Hz，您可以实时查看变化。

更改Spectrum Analyzer设置

中更改设置时频谱分析仪块，可以实时看到对光谱数据的影响。

模型运行时，在痕迹选项窗格的频谱分析仪块,改变规模到日志．PS现在以对数规模显示。

有关如何频谱分析仪设置会影响功率谱数据，请参阅“算法”部分频谱分析仪块参考页面。

转换单位之间的功率

Spectrum Analyzer提供三个单位以指定功率谱密度：瓦特/赫兹，dBm /赫兹，和dbw / hz.．相应的权力单位是瓦特，dBm，和瓦分贝．对于电气工程应用，您也可以查看信号的RMS在VRMS.或者伏特分贝．缺省情况下，spectrum类型为权力在dBm．

将瓦特的电源转换为dbw和dbm

电力瓦分贝是由:

$P_{d B W} ＝ 10. 日志 10. （ p o w e r 我 n w 一个 t t / 1 w 一个 t t ）$

电力dBm是由:

$P_{d B 米} ＝ 10. 日志 10. （ p o w e r 我 n w 一个 t t / 1 米我 l l 我 w 一个 t t ）$

对于振幅为1v的正弦波信号，其单侧频谱的功率为瓦特是由:

$\begin{array}{l} P_{W 一个 t t 年代} ＝ {一个}^{2} / 2 \\ P_{W 一个 t t 年代} ＝ 1 / 2 \end{array}$

在这个例子中，这个功率等于0.5 W。对应的dBm功率为:

$\begin{array}{l} P_{d B 米} ＝ 10. 日志 10. （ p o w e r 我 n w 一个 t t / 1 米我 l l 我 w 一个 t t ） \\ P_{d B 米} ＝ 10. 日志 10. （ 0.5 / {10.}^{- 3.} ） \end{array}$

在这里，功率等于26.9897 dBm。要使用Peak Finder确认此值，请单击工具>测量>峰探测器．

对于白噪声信号，所有频率的频谱都是平坦的。本例中的频谱分析仪在[0 Fs/2]范围内显示单侧频谱。对于方差为1e-4的白噪声信号，单位带宽功率(P_{UnitaBandWidth.})是1的军医。白噪声的总功率美国瓦茨在整个频率范围内为:

$\begin{array}{l} P_{w h 我 t e n o 我年代 e} ＝ P_{u n 我 t b 一个 n d w 我 d t h} ＊ n u 米 b e r o f f r e 问 u e n c y b 我 n 年代， \\ P_{w h 我 t e n o 我年代 e} ＝（ {10.}^{- 4} ）＊（ \frac{F 年代 / 2}{R B W} ）， \\ P_{w h 我 t e n o 我年代 e} ＝（ {10.}^{- 4} ）＊（ \frac{22050}{21.53} ） \end{array}$

频率箱的数量是总带宽与RBW的比率。对于单面频谱，总带宽是采样率的一半。这个例子中的RBW是21.53 Hz。通过这些值，白噪声的总功率美国瓦茨在DBM中是0.1024 W.可以使用白色噪声的功率计算10 * log10（0.1024 / 10 ^ -3），这等于20.103 dBm。

将瓦特的电源转换为dbfs

如果将光谱单位设置为DBFS.并设置满量程（FullScaleSource）到汽车,权力DBFS.计算为：

$P_{d B F 年代} ＝ 20. \cdot {日志}_{10.} （ \sqrt{P_{w 一个 t t 年代}} / F u l l ＿年代 c 一个 l e ）$

地点:

P_美国瓦茨是瓦特的力量
对于双和浮动信号，full_scale.是输入信号的最大值。
对于定点或整数信号，full_scale.是可表示的最大值。

如果指定手动满量程（设置FullScaleSource到财产),在DBFS.是由:

$P_{F 年代} ＝ 20. \cdot {日志}_{10.} （ \sqrt{P_{w 一个 t t 年代}} / F 年代）$

在哪里FS表格中是否规定了全部比例系数全档财产。

对于振幅为1v的正弦波信号，其单侧频谱的功率为瓦特是由:

$\begin{array}{l} P_{W 一个 t t 年代} ＝ {一个}^{2} / 2 \\ P_{W 一个 t t 年代} ＝ 1 / 2 \end{array}$

在这个例子中，这个功率等于0.5 W，正弦波的最大输入信号是1v。dBFS中相应的幂为:

$P_{F 年代} ＝ 20. \cdot {日志}_{10.} （ \sqrt{1 / 2} / 1 ）$

这里，功率等于-3.0103。要在Spectrum Analyzer中确认此值，请运行以下命令：

Fs = 1000;%采样频率sinef = dsp.SineWave('SampleRate'，Fs，'SamplesPerFrame'，100);范围= dsp.SpectrumAnalyzer(“SampleRate”,Fs,…'SpectrumUnits'，'dBFS'，'PlotAsTwoSidedSpectrum'，false) %% for ii = 1:10万xsin = sinef();范围(xsine)结束

然后，点击工具>测量>峰探测器．

将dBm中的Power转换为vrm中的RMS

电力dBm是由:

$P_{d B 米} ＝ 10. 日志 10. （ p o w e r 我 n w 一个 t t / 1 米我 l l 我 w 一个 t t ）$

RMS中的电压由：

$V_{r 米年代} ＝ {10.}^{P_{d B 米} / 20.} \sqrt{{10.}^{- 3.}}$

从前面的例子，P_dBm= 26.9897 dBm。V_rms.计算为

$V_{r 米年代} ＝ {10.}^{26.9897 / 20.} \sqrt{0.001}$

等于0.7071。

确认此值:

改变类型到rms.．
点击打开峰值查找器工具>测量>峰探测器．

估计功率谱使用谱估计块

或者，您可以使用谱估计块的dspspect3图书馆。您可以获得频谱估计器的输出，并存储数据以作进一步处理。

取代频谱分析仪块ex_psd_sa.与之谱估计块后跟一个数组图块。要查看模型，输入ex_psd_estimatorblock.在MATLAB命令提示符中。此外，要在MATLAB中访问谱估计数据，需要连接到工作区(金宝app模型)块的输出谱估计块。以下是对设置设置的更改谱估计块和数组图块。

块参数更改块的目的

块	参数更改	块的目的
谱估计	频率分辨率的方法到`频段数量`．频率范围到`片面`．	用滤波器组方法计算输入信号的功率谱。
数组图	点击视图和选择`风格`．在样式窗口中，选择情节类型作为`楼梯`．选择`配置属性`．在“配置属性”窗口中，在主要的选项卡,设置样本增量作为`44.1 / 1024.`．在展示选项卡中,改变X-label到`频率(赫兹)`，y-label.到`权力（DBM）`．有关详细信息，请参见“转换”一节`x`- 表示频率'。另外，套装y-limits（最低）到`－100`和Y-limits(最大)到`40`．	显示功率谱数据。

谱估计

频率分辨率的方法到频段数量．
频率范围到片面．

用滤波器组方法计算输入信号的功率谱。

数组图

点击视图和

选择风格．在样式窗口中，选择情节类型作为楼梯．
选择配置属性．在“配置属性”窗口中，在主要的选项卡,设置样本增量作为44.1 / 1024.．在展示选项卡中,改变X-label到频率(赫兹)，y-label.到权力（DBM）．有关详细信息，请参见“转换”一节x- 表示频率'。另外，套装y-limits（最低）到－100和Y-limits(最大)到40．