非平稳信号的频谱分析方法—(短时傅立叶变换)

非平稳信号的频谱分析方法—(短时傅立叶变换)非平稳信号又称时变信号。对这一类信号,其一阶、二阶统计量和功率谱的估计显然不能简单的使用平稳信号的估计方法,必须考虑它们的时变因素。

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

   非平稳信号又称时变信号。对这一类信号,其一阶、二阶统计量和功率谱的估计显然不能简单的使用平稳信号的估计方法,必须考虑它们的时变因素。

基本原理

   对非平稳信号,人们希望能有一种分析方法把时域分析和频域分析结合起来,即找到一个二维函数,它既能反映信号的频率内容,也能反映出该频率内容随时间变化的规律。研究这一问题的信号处理理论称为信号的联合时频分布。其中最重要的是以Cohen类为代表的双线性时频分布,此分布可表示为

非平稳信号的频谱分析方法---(短时傅立叶变换)   

   式中非平稳信号的频谱分析方法---(短时傅立叶变换)是一个二维的窗函数,给定不同的窗函数可以得到不同的时频分布。在上式中x(t)出现了两次,且是相乘的形式,这一特点称为双线性。

非平稳信号的频谱分析方法---(短时傅立叶变换)       

   式中w是一个一维的窗函数,则(1)式可以简化成如下的谱图

非平稳信号的频谱分析方法---(短时傅立叶变换)          

式中

非平稳信号的频谱分析方法---(短时傅立叶变换)               

称为信号x(t)的短时傅里叶变换,它反映了信号的频谱随时间和频率的分布。

Matlab仿真分析

   生成的信号两个LFM信号,一个频率随时间由小变大,另一个由大变小,中心频率都为150KHz。

Matlab代码为:

f0 = 150e3;            %定义LFM信号的中心频率
d_f = 100e3;            %定义LFM信号的频率变化大小
fs = 5*f0;             %定义采样率
points = 2^15;         %采样点数
SNR = 40;              %信噪比
 
nTs = (0:points-1)/fs;
Delta_f = (0:points-1)*fs/points;
%生成LFM信号
signal = sin(2*pi*(f0 + (d_f*fs)/(2*points)*(nTs - points/(2*fs))).*(nTs - points/(2*fs))) + ...  %频率由小变大
         sin(2*pi*(f0 + (-d_f*fs)/(2*points)*(nTs - points/(2*fs))).*(nTs - points/(2*fs)));      %频率由大变小                                                                    
%加入一定功率的白噪声
signal_noise = signal + 10^(-SNR/20)*randn(1,points);

Jetbrains全家桶1年46,售后保障稳定

N = 2^8;                                          %截取信号片段点数
window = hanning(N);                              %窗函数类型
Noverlap = N - N/2;                               %相邻两段之间重叠的点数
[S,F,T] = specgram(signal_noise,N,fs,window,Noverlap);
abs_S = abs(S);
max_abs_S = max(max(abs_S));
abs_S = abs_S/max_abs_S;
log10_abs_S = 20*log10(abs_S);
[tt,ff] = meshgrid(T,F/1000);
figure;
surf(tt,ff,log10_abs_S);
xlabel('时间/s');
ylabel('频率KHz');
zlabel('归一化功率谱P(w,t)/dB');
string = ['短时傅里叶变换,采样点数=',num2str(points),',LFM信号中心频率f0=',num2str(f0/1000),...
          'KHz,信号频率变化量d_f=',num2str(d_f/1000),'KHz,采样率fs=',num2str(fs/1000),'KHz,信噪比SNR=',...
          num2str(SNR),'dB'];
title(string);


运行一下得到如下图:

非平稳信号的频谱分析方法---(短时傅立叶变换)

调整一下得到如下的俯视图

非平稳信号的频谱分析方法---(短时傅立叶变换)

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/234168.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • ubuntu系统下mysql重置密码和修改密码操作

    ubuntu系统下mysql重置密码和修改密码操作一、忘记密码后想重置密码在介绍修改密码之前,先介绍一个文件/etc/mysql/debian.cnf.其主要内容如下图:里面有一个debian-sys-maint用户,这个用户只有Debian或Ubuntu服务器才有,所以如果您的服务器是Debain或Ubuntu,debian-sys-maint是个Mysql安装之后自带的用户,具体作用是重启及运行mysql服务。所以如果忘了root密码,可以通…

  • spring cloud 入门系列四:使用Hystrix 实现断路器进行服务容错保护「建议收藏」

    在微服务中,我们将系统拆分为很多个服务单元,各单元之间通过服务注册和订阅消费的方式进行相互依赖。但是如果有一些服务出现问题了会怎么样?比如说有三个服务(ABC),A调用B,B调用C。由于网络延迟或C

  • 如何打造一款直播App(方法流程)

    如何打造一款直播App(方法流程)概要分享内容:互联网内容载体变迁历程,文字——图片/声音——视频——VR/AR——…….。从直播1.0秀场时代(YY),2.0游戏直播(斗鱼、虎牙、熊猫)到如今全民直播3.0泛生活娱乐时代(映客、花椒),国外直播app(Meerkat、Periscope),随着VA/AR/MR提出的沉浸式视听体验,直播4.0时代很快就能到来。在这个全民娱乐的时代,直播已经火得不要不要的,

  • goland 激活服务器【2021.10最新】[通俗易懂]

    (goland 激活服务器)本文适用于JetBrains家族所有ide,包括IntelliJidea,phpstorm,webstorm,pycharm,datagrip等。https://javaforall.cn/100143.htmlIntelliJ2021最新激活注册码,破解教程可免费永久激活,亲测有效,上面是详细链接哦~F…

  • 下载网络直播,多个ts文件下载与合并

    下载网络直播,多个ts文件下载与合并一、chrome与迅雷批量下载视频进入录播回放的页面,右键审查元素,选择network,可以知晓目前的文件活动情况。首先出来的是第一个ts文件,如图所示:472.ts,这个是第一个文件。把进度条拉动到最后面,会下载最后一个文件,如1850.ts。接下来可以使用迅雷的批量下载功能。如:http://www.baidu.com/abc001.ts文件,把格式改为http://www.b…

  • npm、cnpm安装「建议收藏」

    npm、cnpm安装「建议收藏」01.npm安装1.node官网https://nodejs.org/zh-cn/2.安装教程https://www.cnblogs.com/goldlong/p/8027997.html01:双击安装02:可以使用默认路径,本例子中自行修改为d:\nodejs03:一路点Next04:点Finish完成05:打开CMD,检查是否正常06:再看看另…

    2022年10月15日

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号