|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
根据破凰的边际谱程序,改了一下,做瞬间能量的图。
但是由于本人信号处理不熟悉,无法判断结果正确与否。
先把程序和图贴上来,大家帮忙看看啊。
%Hilbert瞬时能量谱图
function [N,IE]=ie(x,n,dt,fs)
%x 原始信号
%n采样点数
%dt采样周期
%fs=1/dt;%采样频率
imf=emd(x);
%[S,freq]=hhte(imf,2) %Y.Lzhang
[A,fa,tt]=hhspectrum(imf);
[E,tt1]=toimage(A,fa,tt,length(tt));
for k=1:size(E,2)
IE(k)=sum(E(k,:))^2*2*pi*fs; %瞬时能量谱
end
N=(0:n-3)/n
clear all
n=1000;
dt=0.002;
%n=250;
%dt=0.002;
fs=1/dt;%采样频率
for i=1:n
t(i)=i*dt;
x(i)=exp(-0.5*t(i))*(sin(10*pi*t(i))+t(i))+0.2*cos(2*pi*t(i));
%x(i)=cos(2*pi*t(i)+0.5*sin(2*pi*t(i)))*exp(-0.2*t(i))+0.05*sin(30*i*t(i));
%%任春论文
end
[N,IE]=ie(x,1000,0.002,1/0.002);
%plot(N(1:25),IE(1:25),'k');
plot(N(1:50),IE(1:50),'k');
xlabel('时间/ Hz','Fontsize',10);
set(gca,'XTick',0:0.01:0.05)
ylabel('瞬时能量','Fontsize',10);
title('时间','FontSize',10) |
|