shat.m

本程序的目的是通过变化DCT系数以及寻找可以应用的DCT系数的最小值来进行可扩展的MPEG压缩

function result = shat(originalSI,decodedSI,originalTI,decodedTI)
% m1=rms(abs(5.81*((originalSI-decodedSI)/originalSI)));
% temp=0.108*max([(originalTI-decodedTI) 0]);
% m2=std(conv(temp,[-1,2,-1]));
% m3=max(4.23*log10(decodedTI/originalTI));
% 
% result=4.77-(0.992*m1)-(0.272*m2)-(0.356*m3);


m1=rms(rms(abs(5.81*((originalSI-decodedSI)./originalSI))));
temp=0.108*max([(originalTI-decodedTI) 0]);
m2=std(conv(temp,[-1,2,-1]));
m3=max(4.23*log10(decodedTI./originalTI));

result=4.77-(0.992*m1)-(0.272*m2)-(0.356*m3);