一个图像融合函数

[zhangyang]

这是一个图像融合程序，我看不太懂，也没有运行出结果，:请各位达人帮我分析一下，谢谢啦！
%IHS变换融合
rgb_image1=imread('D:\My Documents\MATLAB\R2007a\toolbox\images\imdemos\mul.tif');
rgb_image1=im2double(rgb_image1);
R=rgb_image1(:,:,1);
G=rgb_image1(:,:,2);
B=rgb_image1(:,:,3);
I=(R+G+B);
H=(G-B)/(3*(I-B));
S=1-B/I;
rgb_image2=imread('D:\My Documents\MATLAB\R2007a\toolbox\images\imdemos\pan.tif');
rgb_image1=im2double(rgb_image1);
rgb_image2=im2double(rgb_image2);
I=rgb_image2;
R=I*(1+2*S-3*H*S);
G=I*(1-S+3*H*S);
B=I*(1-S);
rgb_image=cat(3,R,G,B); %将R,G,B组成一个三维数组
subplot(1,3,1);image(rgb_image1);
subplot(1,3,2);image(rgb_image2);
subplot(1,3,3);image(rgb_image);