clear;
clc;

%讀入連續四幀圖像
im1=double(imread('img1.bmp'));
im2=double(imread('img2.bmp'));
im3=double(imread('img3.bmp'));
im4=double(imread('img4.bmp'));

[ny,nx]=size(im1);

%對相鄰幀進行差分運算
yn=(im2-im1);

%從殘留圖像中選取大小為100*100的像素塊
yn=yn((ny/2):(ny/2+99),(nx/2):(nx/2+99));

%計算Friedman統計量及概率alpha
[p,table,stats]=friedman(yn(1:100,1:100),1,'off');
s(1)=table(2,5);
a(1)=table(2,6);
[p,table,stats]=friedman(yn(1:50,1:100),1,'off');
s(2)=table(2,5);
a(2)=table(2,6);
[p,table,stats]=friedman(yn(51:100,1:100),1,'off');
s(3)=table(2,5);
a(3)=table(2,6);
[p,table,stats]=friedman(yn(26:75,1:100),1,'off');
s(4)=table(2,5);
a(4)=table(2,6);
[p,table,stats]=friedman(yn(1:30,1:100),1,'off');
s(5)=table(2,5);
a(5)=table(2,6);
[p,table,stats]=friedman(yn(31:60,1:100),1,'off');
s(6)=table(2,5);
a(6)=table(2,6);
[p,table,stats]=friedman(yn(61:90,1:100),1,'off');
s(7)=table(2,5);
a(7)=table(2,6);
[p,table,stats]=friedman(yn(1:20,1:100),1,'off');
s(8)=table(2,5);
a(8)=table(2,6);
[p,table,stats]=friedman(yn(21:60,1:100),1,'off');
s(9)=table(2,5);
a(9)=table(2,6);
[p,table,stats]=friedman(yn(71:90,1:100),1,'off');
s(10)=table(2,5);
a(10)=table(2,6);

%對相鄰幀進行差分運算
yn=(im3-im2);

%從殘留圖像中選取大小為100*100的像素塊
yn=yn((ny/2):(ny/2+99),(nx/2):(nx/2+99));

%計算Friedman統計量及概率alpha
[p,table,stats]=friedman(yn(1:100,1:100),1,'off');
s2(1)=table(2,5);
a2(1)=table(2,6);
[p,table,stats]=friedman(yn(1:50,1:100),1,'off');
s2(2)=table(2,5);
a2(2)=table(2,6);
[p,table,stats]=friedman(yn(51:100,1:100),1,'off');
s2(3)=table(2,5);
a2(3)=table(2,6);
[p,table,stats]=friedman(yn(26:75,1:100),1,'off');
s2(4)=table(2,5);
a2(4)=table(2,6);
[p,table,stats]=friedman(yn(1:30,1:100),1,'off');
s2(5)=table(2,5);
a2(5)=table(2,6);
[p,table,stats]=friedman(yn(31:60,1:100),1,'off');
s2(6)=table(2,5);
a2(6)=table(2,6);
[p,table,stats]=friedman(yn(61:90,1:100),1,'off');
s2(7)=table(2,5);
a2(7)=table(2,6);
[p,table,stats]=friedman(yn(1:20,1:100),1,'off');
s2(8)=table(2,5);
a2(8)=table(2,6);
[p,table,stats]=friedman(yn(21:60,1:100),1,'off');
s2(9)=table(2,5);
a2(9)=table(2,6);
[p,table,stats]=friedman(yn(71:90,1:100),1,'off');
s2(10)=table(2,5);
a2(10)=table(2,6);

%對相鄰幀進行差分運算
yn=(im4-im3);
%從殘留圖像中選取大小為100*100的像素塊
yn=yn((ny/2):(ny/2+99),(nx/2):(nx/2+99));

%計算Friedman統計量及概率alpha
[p,table,stats]=friedman(yn(1:100,1:100),1,'off');
s3(1)=table(2,5);
a3(1)=table(2,6);
[p,table,stats]=friedman(yn(1:50,1:100),1,'off');
s3(2)=table(2,5);
a3(2)=table(2,6);
[p,table,stats]=friedman(yn(51:100,1:100),1,'off');
s3(3)=table(2,5);
a3(3)=table(2,6);
[p,table,stats]=friedman(yn(26:75,1:100),1,'off');
s3(4)=table(2,5);
a3(4)=table(2,6);
[p,table,stats]=friedman(yn(1:30,1:100),1,'off');
s3(5)=table(2,5);
a3(5)=table(2,6);
[p,table,stats]=friedman(yn(31:60,1:100),1,'off');
s3(6)=table(2,5);
a3(6)=table(2,6);
[p,table,stats]=friedman(yn(61:90,1:100),1,'off');
s3(7)=table(2,5);
a3(7)=table(2,6);
[p,table,stats]=friedman(yn(1:20,1:100),1,'off');
s3(8)=table(2,5);
a3(8)=table(2,6);
[p,table,stats]=friedman(yn(21:60,1:100),1,'off');
s3(9)=table(2,5);
a3(9)=table(2,6);
[p,table,stats]=friedman(yn(71:90,1:100),1,'off');
s3(10)=table(2,5);
a3(10)=table(2,6);

%對高斯白噪聲計算Friedman統計量及概率alpha,與上述結果進行比較

yn=randn(100,100);
[p,table,stats]=friedman(yn(1:100,1:100),1,'off');
sn(1)=table(2,5);
an(1)=table(2,6);
[p,table,stats]=friedman(yn(1:50,1:100),1,'off');
sn(2)=table(2,5);
an(2)=table(2,6);
[p,table,stats]=friedman(yn(51:100,1:100),1,'off');
sn(3)=table(2,5);
an(3)=table(2,6);
[p,table,stats]=friedman(yn(26:75,1:100),1,'off');
sn(4)=table(2,5);
an(4)=table(2,6);
[p,table,stats]=friedman(yn(1:30,1:100),1,'off');
sn(5)=table(2,5);
an(5)=table(2,6);
[p,table,stats]=friedman(yn(31:60,1:100),1,'off');
sn(6)=table(2,5);
an(6)=table(2,6);
[p,table,stats]=friedman(yn(61:90,1:100),1,'off');
sn(7)=table(2,5);
an(7)=table(2,6);
[p,table,stats]=friedman(yn(1:20,1:100),1,'off');
sn(8)=table(2,5);
an(8)=table(2,6);
[p,table,stats]=friedman(yn(21:60,1:100),1,'off');
sn(9)=table(2,5);
an(9)=table(2,6);
[p,table,stats]=friedman(yn(71:90,1:100),1,'off');
sn(10)=table(2,5);
an(10)=table(2,6);