EWT分解

%% Copyright by Lwcah(公众号：Lwcah) %% 2022-11-19 %% 清空 clc;close all;clear; %% 加载数据 load sineE%横坐标数据 load sortedY%纵坐标数据 x=sortedY;%待分解的数据 %% ewt信号分解（不能控制分解层数） % fs = 250; [mra,cfs] = ewt(x); mra=mra'; cfs=cfs'; %信号重构 d1=mra(1,:); d2=mra(2,:); d3=cfs(1,:); d4=cfs(2,:); % d5=imf(5,:); % d6=imf(6,:); sortedY=d4;%d3+d4 sortedY=sortedY';%要转置哦 %% 画图 figure plot(sineE,x);hold on plot(sineE,sortedY);hold on legend('dSNR','d4');%绘制图例 % 背景颜色 set(gcf,'Color',[1 1 1]) saveas(gcf,'plot1.jpg');%保存当前窗口的图像 %% 画图 imf=[d1;d2]; t=1:length(x); [a b]=size(imf); figure; subplot(a+1,1,1); plot(sineE,x);% the ECG signal is in the first row of the subplot ylabel('dSNR') set(gca,'xtick',[]) % axis tight; for i=2:a subplot(a+1,1,i); plot(sineE,imf(i-1,:)); ylabel (['MRA ' num2str(i-1)]); set(gca,'xtick',[]) % xlim([1 length(x)]) end subplot(a+1,1,a+1) plot(sineE,imf(a,:)) ylabel(['MRA ' num2str(a)]) % 背景颜色 set(gcf,'Color',[1 1 1]) saveas(gcf,'plot1.jpg');%保存当前窗口的图像 %% 箱型图 % figure; % boxplot(its); % % 背景颜色 % set(gcf,'Color',[1 1 1]) % saveas(gcf,'plot3.jpg');%保存当前窗口的图像