rgb=imread('C:\Users\Administrator\Desktop\1.jpg');
subplot(2,2,1),imshow(rgb),title('彩色图像')
%%提取H、S、V分量
HSV=RGB2HSV(rgb);%%RGB颜色空间转换到HSV颜色空间
H=HSV(:,:,1);%%色调
S=HSV(:,:,2);%%饱和度
V=HSV(:,:,3);%%亮度
subplot(2,2,2),imshow(H),title('H分量图像')
subplot(2,2,3),imshow(S),title('S分量图像')
subplot(2,2,4),imshow(V),title('V分量图像')
A=sum(H)
B=sum(A)
%%提取颜色特征参数(i,j,I,J,N,M,H,S,V,)
[I1 J1]=size(H);%%定义矩阵行列数目
N1=I1*J1;
M1H=0;M2H=0;M3H=0;%%初值定义
M1H=sum(sum(H)/N1);
M1H=mean2(H);%%平均值
for i=1:I1
for j=1:J1
M2H=M2H+(((H(i,j)-M1H)*(H(i,j)-M1H))/N1);%%标准差公式
M3H=M3H+(((H(i,j)-M1H)*(H(i,j)-M1H)*(H(i,j)-M1H))/N1);%%方差（斜度）公式
end
end
M2H=sqrt(M2H);%%标准差
M3H=M3H.^(1/3);%%方差
[I2 J2]=size(S);%%定义矩阵行列数目
N2=I2*J2;
M1S=0;M2S=0;M3S=0;%%初值定义
M1S=sum(S)/N2;
M1S=mean2(S);%%平均值
for i=1:I2
for j=1:J2
M2S=M2S+(((S(i,j)-M1S)*(S(i,j)-M1S))/N2);%%标准差公式
M3S=M3S+(((S(i,j)-M1S)*(S(i,j)-M1S)*(S(i,j)-M1S))/N2);%%方差（斜度）公式
end
end
M2S=sqrt(M2S);%%标准差
M3S=M3S.^(1/3);%%方差
[I3 J3]=size(V);%%定义矩阵行列数目
N3=I3*J3;
M1V=0;M2V=0;M3V=0;%%初值定义
M1V=sum(sum(V)/N3);
M1V=mean2(V);%%平均值
for i=1:I3
for j=1:J3
M2V=M2V+(((V(i,j)-M1V)*(V(i,j)-M1V))/N3);%%标准差公式
M3V=M3V+(((V(i,j)-M1V)*(V(i,j)-M1V)*(V(i,j)-M1V))/N3);%%方差（斜度）公式
end
end
M2V=sqrt(M2V);%%标准差
M3V=M3V.^(1/3);%%方差
F=[M2H,M3H,M1H,M1V,M2V,M3V,M1S,M2S,M3S]; %%颜色特征参数
