交叉小波分析(by matlab)

交叉小波分析需先下载matlab小波工具箱，小波分析可实现多种功能，高频的小波谱可研究序列突变或振荡的瞬态，低频的小波谱可研究序列周期性特征，下图为$sinx$（数据1）与$sin(x-\frac{\pi}{4})$（数据2）的交叉小波谱（有数据顺序），交叉小波箭头含义在此不再赘述

clear all;
clc;
%% 
data=readtable('E:\data\站点数据.txt');%各站点溶解氧

time=table2array(data(:,1));%时间
d1=table2array(data(:,2));%数据1
d2=table2array(data(:,3));%数据2
d3=table2array(data(:,5));%数据3
d4=table2array(data(:,6));%数据4
target=table2array(data(:,4));%目标数据

%%
n = length(time);
figure;
%
%set(gcf,'unit','normalized','position',[0.1,0.25,0.8,0.5]); %采用相对值设置，相对屏幕
%set (gca,'position',[0.1,0.1,0.8,0.8] );
%subplot(1,2,1);
xwt(d1,target);%%%调用xwt那个function
%set(gca,'FontName','songti');
set(gca,"xtick",[1:180:n]);%%% 步长
set(gca,'xticklabel',{'2020-01','2020-07','2021-01','2021-07','2022-01','2022-07','2023-01','2023-07'},'fontsize',10.5, 'Fontweight','bold','FontName','Times New Roman');%%% 标题名称+字体格式
title('\fontname{宋体}水温与溶解氧交叉小波功率谱','FontWeight','bold','FontSize',12,'FontName','Times New Roman');%%% 标题名称+字体格式
%xlabel('Date','FontWeight','bold','fontsize',10.5,'FontName','Times New Roman');%%% x标签名称+字体格式
ylabel('周期/日','FontWeight','bold','FontSize',10.5,'FontName','songti');%y轴标题名称+格式
set(gca, 'YDir', 'normal');#reverse：y轴倒置