主要内容

发现使用持久性光谱干扰

想象一个窄带干扰信号嵌入在一个宽带信号。

生成一个啁啾取样1 kHz为500秒。短促声波的频率增加从180赫兹到220赫兹在测量。

fs = 1000;t = (0:1 / fs: 500)”;x =唧唧喳喳(220 t, 180 t(结束),)+ 0.15 * randn(大小(t));

信号还包含一个210 Hz正弦信号。正弦信号的振幅为0.05,目前只对总信号持续时间的1/6。

idx =地板(长度(x) / 6);(1:idx) = x (1: idx) + 0.05 * cos(2 *π* t (1: idx) * 210);

保存信号作为一个MATLAB®时间表。

S =时间表(x秒(t));

开放信号分析仪并拖动的时间表工作空间浏览器显示。单击时频按钮来添加一个谱图视图。在光谱图选项卡,在时间分辨率中,选择指定并输入一个1秒的时间分辨率。设置频率的限制100 Hz, 290 Hz。这两个信号组件是可见的。

回到显示选项卡。单击时间并单击按钮来删除时间视图光谱按钮来添加一个功率谱的观点。的频率范围仍然是100赫兹到290赫兹。疲软的正弦信号被唧唧声。

单击光谱▼按钮来改变光谱为了一个持久性频谱视图。在持久性频谱选项卡,在时间分辨率中,选择指定并输入一个1秒的时间分辨率。指定零重叠的部分。设置权限-50分贝和0分贝密度的限制0.1和4。现在信号组件都清晰可见。

显示选项卡,在分享,点击生成脚本▼并选择持久性频谱脚本。MATLAB脚本出现在编辑器。

%计算持久性频谱%由MATLAB 9.7 (R)和信号处理工具箱8.2。%上生成:26 - 12月- 2018 16:07:45%的参数期限=秒(500 [0]);%秒frequencyLimits = (100 - 290);%赫兹timeResolution = 1;%秒overlapPercent = 0;% %%索引信号时间感兴趣的地区S_x_ROI = S (:,“x”);S_x_ROI = S_x_ROI (timerange(时限(1),期限(2),“关闭”),1);%计算谱估计%运行下面的函数调用不带输出参数绘制结果[P F,压水式反应堆]= pspectrum (S_x_ROI,“坚持不懈”,“FrequencyLimits”frequencyLimits,“TimeResolution”timeResolution,“OverlapPercent”,overlapPercent);

另请参阅

应用程序

功能

相关的例子

更多关于