該樓層疑似違規已被系統折疊?隱藏此樓查看此樓

附上MATLAB內置的hht繪圖函數hhtPlot.m的代碼：

classdef hhtPlot

properties

hfig

end

methods

function [this] = hhtPlot(insFreq, insEnergy, tvec, FRange, MinThres, FreqLoc, isNormFreq)

isYAxis = strcmpi(FreqLoc,'yaxis');

% patch currently do not support datetime/duration, convert to

% double

if(isduration(tvec) || isdatetime(tvec))

tvec = seconds(tvec);

end

frequnitstrs = getfrequnitstrs;

if isNormFreq

tvec = 2*pi*tvec; % Convert time axis to samples

freqlbl = frequnitstrs{1};

timelbl = getString(message('shared_signalwavelet:hht:general:TimeUnitSamples'));

else

freqlbl = frequnitstrs{2};

timelbl = getString(message('shared_signalwavelet:hht:general:TimeUnitSeconds'));

end

if isYAxis

xlbl = timelbl;

ylbl = freqlbl;

xyrange = [0,tvec(end),FRange(1),FRange(2)];

else

xlbl = freqlbl;

ylbl = timelbl;

xyrange = [FRange(1),FRange(2),0,tvec(end)];

end

this.hfig = newplot;

numIMF = size(insFreq,2);

for i = 1:numIMF

% Plot each IMF

insfi = insFreq(:,i);

insei = insEnergy(:,i);

insei((10*log10(insei))

if(isYAxis)

patch([tvec(1);tvec;tvec(end)], [0;insfi;0], [nan;insei;nan], ...

'EdgeColor','interp','EdgeAlpha','interp',...

'FaceColor', 'none', 'FaceVertexAlphaData',[nan;insei;nan],...

'LineWidth', 2, 'FaceAlpha', 'interp');

else

patch([0;insfi;0], [tvec(1);tvec;tvec(end)], [nan;insei;nan], ...

'EdgeColor','interp','EdgeAlpha','interp',...

'FaceColor', 'none', 'FaceVertexAlphaData',[nan;insei;nan],...

'LineWidth', 2, 'FaceAlpha', 'interp');

end

end

axis(xyrange);

xlabel(xlbl);

ylabel(ylbl);

title(getString(message('shared_signalwavelet:hht:plot:hhtPlotTitle')));

colormap parula

colorbar

end

end

end

總結

以上是生活随笔為你收集整理的matlab算hht,关于MATLAB中HHT谱图的问题【百度知道悬赏100积分】的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。