主要内容

使用持久性频谱查找干扰

可视化嵌入在宽带信号中的干扰窄带信号。

以1 kHz的形式产生一个chirp,持续500秒。在测量过程中,CHIRP的频率从180 Hz增加到220 Hz。

FS = 1000;t =(0:1/fs:500)';x = chirp(t,180,t(end),220) + 0.15*randn(size(t));

该信号还包含210 Hz正弦曲线。正弦曲线的振幅为0.05,仅在总信号持续时间的1/6中存在。

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

将信号保存为MATLAB®时间表。

s =时间表(秒(t),x);

打开信号分析仪并将时间表从工作区浏览器显示。点击时频按钮添加频谱图。在频谱图标签,下时间分辨率, 选择指定并输入1秒钟的时间分辨率。设置频率限制至100 Hz和290 Hz。两个信号组件都是可见的。

回到展示标签。点击时间按钮删除时间视图,然后单击光谱按钮以添加电源谱视图。频率范围继续从100 Hz到290 Hz。弱的正弦曲线被鸣叫掩盖。

单击箭头光谱按钮更改光谱查看持续频谱看法。在持续频谱标签,下时间分辨率, 选择指定并输入1秒钟的时间分辨率。指定相邻段之间的零重叠。设置功率限制到–50 dB和0 dB,密度极限到0.1和4.现在两个信号组件都清晰可见。

展示标签,下分享, 点击生成脚本并选择持久性频谱脚本。该脚本出现在MATLAB编辑器中。

%计算持续频谱由MATLAB(R)9.7和信号处理工具箱8.2生成的%。%生成:26-DEC-2018 16:07:45% 参数timelimits =秒([0 500]);%秒频率limits = [100 290];%Hz计时度= 1;%秒重叠= 0;%%%索引进入信号时间的关注时间区域s_x_roi = s(:,,,'X');s_x_roi = s_x_roi(timerange(timelimits(1),timelimits(2),'关闭'),1);%计算频谱估计%在下面运行函数调用,而无需输出参数以绘制结果[p,f,pwr] = pspectrum(s_x_roi,...“持久性”,,,,...“频率限制”,频率限制,...“计时解决”,计时解决,...“超贴”,夸张);

也可以看看

应用

功能

相关示例

更多关于