为什么HHT时频谱图的高频区域会出现波浪线,低频区域却没有,只是散点? 5C

function HHT
clear;clc;clf;
N=2048;
%fft默认计算的信号是从0开始的
t=linspace(1,2,N);deta=t(2)-t(1);fs=1/deta;
x=5*sin(2*pi*10*t)+5*sin(2*pi*35*t);
z=x;
c=emd(z);

%计算每个IMF分量及最后一个剩余分量residual与原始信号的相关性
[m,n]=size(c);
for i=1:m;
a=corrcoef(c(i,:),z);
xg(i)=a(1,2);
end
xg;

for i=1:m-1
%--------------------------------------------------------------------
%计算各IMF的方差贡献率
%定义:方差为平方的均值减去均值的平方
%均值的平方
%imfp2=mean(c(i,:),2).^2
%平方的均值
%imf2p=mean(c(i,:).^2,2)
%各个IMF的方差
mse(i)=mean(c(i,:).^2,2)-mean(c(i,:),2).^2;
end;
mmse=sum(mse);
for i=1:m-1
mse(i)=mean(c(i,:).^2,2)-mean(c(i,:),2).^2;
%方差百分比,也就是方差贡献率
mseb(i)=mse(i)/mmse*100;
%显示各个IMF的方差和贡献率
end;
%画出每个IMF分量及最后一个剩余分量residual的图形
figure(1)
for i=1:m-1
disp(['imf',int2str(i)]) ;disp([mse(i) mseb(i)]);
end;
subplot(m+1,1,1)
plot(t,z)
set(gca,'fontname','times New Roman')
set(gca,'fontsize',14.0)
ylabel(['signal','Amplitude'])

for i=1:m-1
subplot(m+1,1,i+1);
set(gcf,'color','w')
plot(t,c(i,:),'k')
set(gca,'fontname','times New Roman')
set(gca,'fontsize',14.0)
ylabel(['imf',int2str(i)])
end
subplot(m+1,1,m+1);
set(gcf,'color','w')
plot(t,c(m,:),'k')
set(gca,'fontname','times New Roman')
set(gca,'fontsize',14.0)
ylabel(['r',int2str(m-1)])

%画出每个IMF分量及剩余分量residual的幅频曲线
figure(2)
subplot(m+1,1,1)
set(gcf,'color','w')
[f,z]=fftfenxi(t,z);
plot(f,z,'k')
set(gca,'fontname','times New Roman')
set(gca,'fontsize',14.0)
ylabel(['initial signal',int2str(m-1),'Amplitude'])

for i=1:m-1
subplot(m+1,1,i+1);
set(gcf,'color','w')
[f,z]=fftfenxi(t,c(i,:));
plot(f,z,'k')
set(gca,'fontname','times New Roman')
set(gca,'fontsize',14.0)
ylabel(['imf',int2str(i),'Amplitude'])
end
subplot(m+1,1,m+1);
set(gcf,'color','w')
[f,z]=fftfenxi(t,c(m,:));
plot(f,z,'k')
set(gca,'fontname','times New Roman')
set(gca,'fontsize',14.0)
ylabel(['r',int2str(m-1),'Amplitude'])

hx=hilbert(z);
xr=real(hx);xi=imag(hx);
%计算瞬时振幅
sz=sqrt(xr.^2+xi.^2);
%计算瞬时相位
sx=angle(hx);
%计算瞬时频率
dt=diff(t);
dx=diff(sx);
sp=dx./dt;
figure(6)
plot(t(1:N-1),sp)
title('瞬时频率')

%计算HHT时频谱和边际谱
[A,fa,tt]=hhspectrum(c);
[E,tt1]=toimage(A,fa,tt,length(tt));
figure(3)
disp_hhs(E,tt1) %二维图显示HHT时频谱,E是求得的HHT谱
pause
figure(4)
for i=1:size(c,1)
faa=fa(i,:);
[FA,TT1]=meshgrid(faa,tt1);%三维图显示HHT时频图
surf(FA,TT1,E)
title('HHT时频谱三维显示')
hold on
end
hold off
E=flipud(E);
for k=1:size(E,1)
bjp(k)=sum(E(k,:))*1/fs;
end
f=(1:N-2)/N*(fs/2);
figure(5)
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('信号幅值');
title('信号边际谱')%要求边际谱必须先对信号进行EMD分解

function [A,f,tt] = hhspectrum(x,t,l,aff)

error(nargchk(1,4,nargin));

if nargin < 2

t=1:size(x,2);

end

if nargin < 3

l=1;

end

if nargin < 4

aff = 0;

end

if min(size(x)) == 1
if size(x,2) == 1
x = x';
if nargin < 2
t = 1:size(x,2);
end
end
Nmodes = 1;
else
Nmodes = size(x,1);
end

lt=length(t);

tt=t((l+1):(lt-l));

for i=1:Nmodes

an(i,:)=hilbert(x(i,:)')';
f(i,:)=instfreq(an(i,:)',tt,l)';
A=abs(an(:,l+1:end-l));

if aff
disprog(i,Nmodes,max(Nmodes,100))
end

end

function disp_hhs(im,t,inf)

% DISP_HHS(im,t,inf)
% displays in a new figure the spectrum contained in matrix "im"
% (amplitudes in log).
%
% inputs : - im : image matrix (e.g., output of "toimage")
% - t (optional) : time instants (e.g., output of "toimage")
% - inf (optional) : -dynamic range in dB (wrt max)
% default : inf = -20
%
% utilisation : disp_hhs(im) ; disp_hhs(im,t) ; disp_hhs(im,inf)
% disp_hhs(im,t,inf)

figure
colormap(bone)
colormap(1-colormap);

if nargin==1
inf=-20;
t = 1:size(im,2);

end

if nargin == 2
if length(t) == 1
inf = t;
t = 1:size(im,2);
else
inf = -20;
end
end

if inf >= 0
error('inf doit etre < 0')
end

M=max(max(im));

im = log10(im/M+1e-300);

inf=inf/10;

imagesc(t,fliplr((1:size(im,1))/(2*size(im,1))),im,[inf,0]);
set(gca,'YDir','normal')
xlabel(['time'])
ylabel(['normalized frequency'])
title('Hilbert-Huang spectrum')
function [f,z]=fftfenxi(t,y)
L=length(t);N=2^nextpow2(L);
%fft默认计算的信号是从0开始的
t=linspace(t(1),t(L),N);deta=t(2)-t(1);
m=0:N-1;
f=1./(N*deta)*m;
%下面计算的Y就是x(t)的傅里叶变换数值
%Y=exp(i*4*pi*f).*fft(y)%将计算出来的频谱乘以exp(i*4*pi*f)得到频移后[-2,2]之间的频谱值
Y=fft(y);
z=sqrt(Y.*conj(Y));![图片说明](https://img-ask.csdn.net/upload/201804/08/1523194608_477078.jpg)图片说明
为什么同样的程序我画出来的时频图,高频区域会出现波浪线,低频区域却没有,只是散点?求大神帮忙

1个回答

Csdn user default icon
上传中...
上传图片
插入图片
抄袭、复制答案,以达到刷声望分或其他目的的行为,在CSDN问答是严格禁止的,一经发现立刻封号。是时候展现真正的技术了!
其他相关推荐
图像低频高频区域分离
简介   本篇整理记录利用小波分离图像的高频、低频部分信息。 具体实现   实现代码参考资料:小波变换 C++ opencv 实现。 小波变换   小波生成和参考资料中一致。小波变换中,首先抽取出原图像的每一行,进行小波分解,获得水平方向的高、低频信息。 接着抽取出原图像的每一列,进行小波分解,获得垂直方向的高、低频信息。   生成的对应结果如下:
如何找出BMP图像的高频与低频区域
如题
HHT时频谱+HHT边际谱
基于HHT的matlab实现,通过EMD对信号进行分解,利用HHT变化得到时频谱
区域图
本课程课程内容突出实战。
opencv 随机选取区域并保存区域图
#include   #include   #include   #include using namespace cv;   using namespace std;   char filename[100]; char filename1[100]; void main() {     for(int i=1;i { sprintf(filename,"C:/Us
点与封闭区域的关系
看了好长时间的书!<>清华大学 周陪德。。编著!rn1.1rn判断点是否在一个封闭的多边形中!rn看了好长时间,仍不解!思路也没理清!rn故来此请教各位!rn看是否有有人研究过这个问题?rn给点思路,当然更希望各位能给点代码!rn
区域点内个数
区域内点的个数 Time Limit: 1000MS Memory Limit: 65536KB Problem Description X晚上睡不着的时候不喜欢玩手机,也不喜欢打游戏,他喜欢数星星。 Input 多组输入。 每组先输入一个整数N(N <= 10000),接着输入两个点代表矩形的左下点B(x,y)和右上点T(x,y),然后输入N个(X,Y)代表N
济南市区域图
山东省济南市边界面图层,非行政区划图层,可用于裁减其他要素的济南市区域等操作。
VB+EXCESS制作区域图
初学VB代码,VB操作excess,然后制作区域图的一个代码。
中国的城市区域图
shp数据,中国的城市区域图,可以用在项目中。
车间区域图
车间CAD图电子图
广州区域图
可使用的广州shp文件,利用配准后的图像所获得的地图画出来的广州区域边界
JFreeChart创建区域图
前言   只是一个demo,效果:运行程序后保存生成的图表。 过程   1.需要的jar包可以从官网上下载    jcommon-1.0.23.jar    jfreechart-1.0.19-experimental.jar    jfreechart-1.0.19.jar    mysql-connector-java-5.1.34.jar,数据库可以取对应数据库的jar
echart折线区域图
在引入echart区域折线图时,没有出现对应的区域图 当发现引入下面代码到自己的代码中并没有对应的区域图 option = { xAxis: { type: 'category', boundaryGap: false, data: mydata.xname(改成自己的数据) }, yAxis: { ...
Highcharts 曲线区域图
一 代码 &amp;lt;html&amp;gt; &amp;lt;head&amp;gt; &amp;lt;meta charset=&quot;UTF-8&quot; /&amp;gt; &amp;lt;title&amp;gt;Highcharts 曲线区域图&amp;lt;/title&amp;gt; &amp;lt;script src=&quot;http://apps.bdimg.com/libs/jquery/2.1.4/jquery.min.js&quot;&amp;gt;&amp;l
中国自然地理区域图
很好用的地图,分部简单明了,区域准确,值得推荐。
Highcharts 区间区域图
一 代码 &amp;lt;html&amp;gt; &amp;lt;head&amp;gt; &amp;lt;meta charset=&quot;UTF-8&quot; /&amp;gt; &amp;lt;title&amp;gt;Highcharts 区间区域图&amp;lt;/title&amp;gt; &amp;lt;script src=&quot;http://apps.bdimg.com/libs/jquery/2.1.4/jquery.min.js&quot;&amp;gt;&am
Highcharts 堆叠区域图
一 代码 &amp;lt;html&amp;gt; &amp;lt;head&amp;gt; &amp;lt;meta charset=&quot;UTF-8&quot; /&amp;gt; &amp;lt;title&amp;gt;Highcharts 堆叠区域图&amp;lt;/title&amp;gt; &amp;lt;script src=&quot;http://apps.bdimg.com/libs/jquery/2.1.4/jquery.min.js&quot;&amp;gt;&amp;l
低频和高频的概念
在数字图像处理中,通常说得的低频信号,高频信号,究竟是什么意思?rn看到不同的一些说法,那个达人清楚,出来讲讲明白啊!100分!rn
微软bug~为什么会出现点击的区域偏移现象??help~~~~~~~
有时候 。。aspx页面中,会出现控件操作区域偏移,比如一个button ,点在上面点不了,要将鼠标点击在button底部才能点击。。其他控件也有类似问题。。how to solve it ????
低频,中频,高频的划分???
低频,中频,高频好像挺模糊的,大体该如何划分????
STVD出现红色区域
在STVD中,当一行的字符长度大于某个值的时候,一行的后面就会出现红色区域,而且也可能使中文字符乱码。 这本来是STVD的一项功能,用于提示开发者一行的输入不要太长。 可以通过以下方法取消: 菜单“Tools”--“Options” 在“Edit”标签下,把下面的勾去掉就行了。表示不限制一行的长度。
有效区域与无效区域
我们到底该怎么理解这个无效区域与有效区域??
区域
区域(zone):包括开区域(开域)和闭区域(闭域) 注:无特殊说明,区域指开区域 开区域:满足下列条件的平面点集(记为D) (1)D是一个开集 (2)D是连通的,即D中的任意两点都可以用一条折线连接起来,且折线上的点全部在D内 闭区域:开区域及其边界点构成的平面点集 有界区域:如果一个区域D可以被包含在一个以原点为中心的圆里面,那么称区域D有界,即存在正数M,使区域D内的每个点都满足...
谷歌地球为什么出现规则白块区域
在使用谷歌地球的时候出现这种地图,看不清楚啊怎么办,有一个边界明显的区域是为什么啊 怎么解决啊 万能的吧友。。。。。rn[img=https://img-bbs.csdn.net/upload/201312/22/1387722694_589304.png][/img]
无效区域与有效区域
绘图时,是在无效区域上绘图,还是在有效区域上绘图,也就是说是在无效区域上才能绘图,还是在有效区域上才能绘图??
如何判断区域包含点
 如题,对包含弧线等的对边形封闭区域,如何保存下来,并判断某个坐标是否在多边形内部,用什么样的结构能保存这样的一个区域?望高手指点!
高德地图通过点连成区域
最近要用高德地图做一个需求,通过在地图进行标记,然后根据标记圈出一块区域。我的思路是先把每次点击marker的坐标存入一个数组,再用数组里面的坐标生成一个多边形。(每次生成多边形要把上个清除掉,不然多边形会相互覆盖。) 先上效果图: &lt;!doctype html&gt; &lt;html&gt; &lt;head&gt; &lt;meta charset="utf-8"&gt...
flash和flex交互 区域图
flash和flex交互 区域图 说明文档 步骤文档
利用区域图分析报表
区域图用于强调数量随时间而变化的程度,也可用于引起人们对总值趋势的注意。还可以显示部分与整体关系。本实例使用区域图显示了某个杀毒软件的全年客户满意度调查表,实例中生成图表的数据是随机生成的,所以每次刷新页面,这个区域都会改变。rn技术要点:rn本实例使用了新的JFreeChart数据集对象和创建区域图的静态工厂方法。下面分别介绍。rn1.TimeSeries类rnTimeSeries类是描述时间序列的一个数据集对象,该对象的add()方法可以为数据集添加新的数据项。其语法格式如下:rnpublic void add(RegularTimePeriod period double value)rn2.创建JFreeChart对象rn本实例使用ChartFactory工厂类的createXYAreaChart()静态方法创建JFreeChart对象。ChartFactory工厂类提供了很多创建不同图表的静态方法,本实例主要使用的是创建区域图的方法,其语法格式如下:rnpublic static JFreeChart createXYAreaChart(java.lang.String title,String xAxisLabel,String yAxisLabel,XYDataset dataset,PlotOrientation orientation,boolean legend,boolean tooltips,boolean urls) rn文章转自http://www.syym.org
matplot之区域画多图
机器学习中,比较热门的语言Python。其中以Numpy、Pandas、Matplotlib为常用。rn本套课程就给大家讲解下Matplotlib的相关用法.rnAnacoda正确的使用姿势rnNotebook基本使用rnmatplot之基本框绘制rnmatplot之折线图优化rnmatplot之区域画多图rnmatplot之其他操作
Highcharts 丢失值区域图
一 代码 &amp;lt;html&amp;gt; &amp;lt;head&amp;gt; &amp;lt;meta charset=&quot;UTF-8&quot; /&amp;gt; &amp;lt;title&amp;gt;Highcharts 教程 | 菜鸟教程(runoob.com)&amp;lt;/title&amp;gt; &amp;lt;script src=&quot;http://apps.bdimg.com/libs/jquery/2.1.4/jquery.min.js&quo
Echart——饼图和堆叠区域图
主要记录下饼图和堆叠区域图的使用,方便工作中快速生成图表代码<!DOCTYPE html> <html> <head> <script src="https://cdn.bootcss.com/jquery/3.2.1/jquery.min.js"></script> <script src="https://cdn.bootcss.com/echarts/3.8.5/echarts.min
Highcharts 使用负数区域图
一 代码 &amp;lt;html&amp;gt; &amp;lt;head&amp;gt; &amp;lt;meta charset=&quot;UTF-8&quot; /&amp;gt; &amp;lt;title&amp;gt;Highcharts 使用负数区域图&amp;lt;/title&amp;gt; &amp;lt;script src=&quot;http://apps.bdimg.com/libs/jquery/2.1.4/jquery.min.js&quot;&amp;gt;&amp
整张图点击区域问题
我现在用一个div的background做一个图层,这里面有好多个区域需要点击,应该怎么做?rn先说明不能用div去划分,也不能用onXXX事件去做单击事件,只能用做单击rnrn
Highcharts 百分比堆叠区域图
一 代码 &amp;lt;html&amp;gt; &amp;lt;head&amp;gt; &amp;lt;meta charset=&quot;UTF-8&quot; /&amp;gt; &amp;lt;title&amp;gt;Highcharts 百分比堆叠区域图&amp;lt;/title&amp;gt; &amp;lt;script src=&quot;http://apps.bdimg.com/libs/jquery/2.1.4/jquery.min.js&quot;&amp;gt;&am
为什么在打印输出时,剪切区域无效呢?
使用GDI绘制图形时使用了SelectClipRgn,rn后来想打印输出此图形,发现SelectClipRgn无效,请问有没有解决办法。rn//打印模式为:rnPrintDC.SetMapMode(MM_ANISOTROPIC);rn
为什么在打印输出时,剪切区域效呢?
使用GDI绘制图形时使用了SelectClipRgn,rn后来想打印输出此图形,发现SelectClipRgn无效,请问有没有解决办法。rn//打印模式为:rnPrintDC.SetMapMode(MM_ANISOTROPIC);rn
区域销售情况细分图
本课程全面介绍 Excel商业图表的制作,循着由常规操作到高端综合应用的思路来安排案例。适合各类层次的学员学习;
Echarts去掉叠堆折线区域图的区域颜色
一开始是这样的,折线图设置了区域阴影的颜色,现在不想要echarts折线区域图的区域颜色,需要去掉。 只需要修改series里面的一系列相关的参数,属性即可。 series : [ { name:'订单流入总数', type:'line', stack: '总量', areaStyle: { normal: { color: '#8cd5c...
相关热词 c# stream 复制 android c# c#监测窗口句柄 c# md5 引用 c# 判断tabtip 自己写个浏览器程序c# c# 字符串变成整数数组 c#语言编程写出一个方法 c# 转盘抽奖 c#选中treeview