mk突變檢驗代碼

A. 什麼是mk檢驗

mk檢驗是曼－肯德爾法，又稱Mann—Kenddall 檢驗法，是一種氣候診斷與預測技術，應用Mann-Kendall檢驗法可以判斷氣候序列中是否存在氣候突變，如果存在，可確定出突變發生的時間。Mann-Kendall檢驗法也經常用於氣候變化影響下的降水、乾旱頻次趨勢檢測。

由於最初由曼(H.B.Mann)和肯德爾(M.G.Kendall)提出了原理並發展了這一方法，故稱其為曼—肯德爾 (Man-Kendall)法。

(1)mk突變檢驗代碼擴展閱讀：

計算方法：

對於具有n個樣本量的時間序列X，構造一秩序列：

UFi為標准正態分布，它是按時間序列x順序x1,x2,…,xn計算出的統計量序列，給定顯著性水平α，查正態分布表，若|UFi|>Ua，則表明序列存在明顯的趨勢變化。

按時間序列x逆序xn,xn-1,…,x1，再重復上述過程，同時使UBk=–UFk，k=n,n–1,…,1)，UB1=0。這一方法的優點在於不僅計算簡便，而且可以明確突變開始的時間，並指出突變區域。因此，是一種常用的突變檢測方法。

B. origin可以做mk突變檢驗嗎

不可以。
可以用MATLAB軟體做，輸入代碼就可以實現。
1、Mann-Kendall突變檢測方法的簡要計算步驟：計算順序時間序列的秩序列，按照上述公式計算UFk。
2、計算逆序時間序列的秩序列，按照上述公式計算UBr
3、給定顯著性水平，如a=0.05，對於臨界值為Uo.os=加1.96，將UFr與UBx兩個統計量序列曲線與Uo.os=土1.96兩條直線繪制在-一個平面直角坐標繫上。另外，a=0.10對應U0.10=加1.28，u=0.01對應U0.01=士2.32。分析繪制出的UFr與UBk曲線圖，若UFr或UBk的值大於0，則表明序列呈上升趨勢，小於0則呈下降趨勢。當它們超出臨界直線時，表明上升或下降趨勢顯著。超過臨界線的范圍確定為出現突變的時間區域。若UFx與UB，兩條曲線出現交叉點，且交叉點在臨界線之間，它們交叉點對應的時刻便是突變開始的時間。

C. MK突變檢驗（Matlab）

function [ UF,UB ] = MannKendall( x,y,p )
% x表示時間如1982-2015
% y表示對應時間的結果
% p表示顯著水平
N = length(y);
UF = SMK(y);
yy = reshape(y,1,length(y));
yy = fliplr(yy);
UB = -fliplr(SMK(yy));

function U = SMK( Y )
N = length(Y);
s = zeros(1,N);
U(1) = 0;
for k=2:N
r = 0;
s(k) = 0;
for j=1:k-1
if Y(k)>Y(j)
r = r+1;
end
s(k) = s(k-1)+r;
end
E = k (k-1)/4;
VAR = k (k-1) (2 k+5)/72;
U(k) = (s(k)-E)/sqrt(VAR);
end

D. 關於M-K突變檢驗的一些問題，現在用fortran軟體編程，簡單程序看不懂，求大神幫忙，急急急！！！！

1. 我建議你靜下心來認真學習一下Fortran，一本基礎的教材是你需要的。
   

2. MK突變檢驗，反正我是沒聽說過。專業不同，其他行業的人很難理解你的問題。

3. 我認為讓你理解Fortran，和讓別人理解MK突變檢驗。前者更容易。

4. 你既不貼出代碼，又不解釋MK突變檢驗。直接來求助，我確實不知道你期望得到什麼答案？
   

E. mk突變分析圖怎麼看

mk突變分析圖看：mk檢驗是曼－肯德爾法，又稱Mann—Kenddall 檢驗法，是一種氣候診斷與預測技術，應用Mann-Kendall檢驗法可以判斷氣候序列中是否存在氣候突變，如果存在，可確定出突變發生的時間。

1、看mk的外包裝，開封的時候有mk標志的透明袋，並且用mk自己水印的防潮包裝紙包裹完整，正品的防塵袋還非常的精緻工整，假貨的包裝可能就會很粗糙，包裝紙包裹不完整。

2、正品的五金每個字母都很清晰和整齊，而且字母的邊緣也沒有鋸齒，假貨的五金字體傾斜嚴重，做工也不夠精緻。

通關劇情：

當她在地球上伺候她的主人時，她的SHOKAN種族正在外部世界遭殃。紹康現在更喜歡MOTARO的CENTAURS種族，並且幫助他們打SHOKAN族。知道這件事後，SHEEVA開始反叛紹康。她打敗了MOTARO並且在狂怒之下打死紹康。

在拯救了地球和外部世界後，她返回家鄉開始恢復她的種族的榮耀。Sigh，原來也是一女魔頭，在shokan族裡，也許她算漂亮的吧。

F. 用spss怎麼實現mk突變檢驗

主要是看你分了幾組，每組數據的分布如何，而不在於是什麼突變

G. 怎樣在matlab中輸出mk檢驗結果求具體步驟

UF,UB]=MannKendall(Y);
% % Y即為待檢測的時間序列
figure,hold on;
plot(UF,'-b');
plot(UB,'-r');

function [UF,UB]=MannKendall(timeseries)
N=length(timeseries);
UF=SMK(timeseries,N);
for i=1:N
YY(i)=timeseries(N+1-i);
end
u_res=SMK(YY,N);
for i=1:N
UB(i)=-u_res(N+1-i);
end

function u_res=SMK(Y,N)
m_res=zeros(N,1);md_res=zeros(N,1);u_res=zeros(N,1);
m_res(1)=0;
for i=2:N
m_res(i)=0;
md_res(i)=0;
for j=1:i-1
if Y(i)<Y(j)
m_res(i)=m_res(i)+0;
else
m_res(i)=m_res(i)+1;
end
md_res(i)=md_res(i-1)+m_res(i);
end
end
u_res(1)=0;
for i=2:N
E=i*(i-1)/4;
VAR=i*(i-1)*(2*i+5)/72;
u_res(i)=(md_res(i)-E)/sqrt(VAR);
end

H. 關於m-k突變性檢驗的問題： 下面為網上下載到的程序，應該沒有問題，我又向別人要了一個，和這個是一樣的

test.txt中數據應去掉第一列（年份），年份應該在運行後再輸入，這樣出來的圖就不是鋸齒狀的了

I. mk突變檢驗用什麼軟體

mk是氣象學中常用來進行突變檢驗的一種方法，一般會用dps就可以。