- 积分
- 2670
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-12-21
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2018-12-17 16:50:08
|
显示全部楼层
第二个:
clc
clear
%Pettitt突变分析
x=xlsread('D:\APP-X\SI_tem_mon_eemd.xlsx','01','A2:A35');%时间数据
data=xlsread('D:\APP-X\SI_tem_mon_eemd.xlsx','01','AM2:AM35');%SAL数据,剔除奇异点
n=length(data);
s=0;
Ut=zeros(n-1,1);
for i=1:n-1;
for j=i+1;n;
s=s+sign(data(i)-data(j));
end
Ut(i)=s;
end
Kmax=max(abs(Ut));
chg_point=find(abs(Ut)==Kmax); %找出最大值对应的序列号
coef=2*exp(-6*Kmax^2/((n-1)^3+(n-1)^2));
figure(1);
hold on
plot (x(1:n-1),Ut,'K-','linewidth',1.5); %画出整个曲线的变化趋势
plot([x(chg_point),x(chg_point)],[-Kmax,Kmax],'r--','linewidth',1.5);
Klim=Kmax-Kmax*0.1;
plot([x(1),x(n)],[-Klim,-Klim],'b-');
plot([x(1),x(n)],[Klim,Klim],'b-');
figure(2);
plot(x,data); |
|