rs_base.m

音频文件的RS分析

%RS隐写分析方法中的最核心函数
%在此函数中求解出进行M和-M变换的R组合S组的数目参数
function [R,S,U,MR,MS,MU]=RS_base(a,bufsize)
%l=length(a);
%进行M伪装变换 并计算相应的参数
m=[1,0,1,0];
R=0;
S=0;
U=0;
for i=1:bufsize/4
   for m=1:4
		b(m)=a(4*(i-1)+m);
   end
   r=disc(b);
   %计算未进行函数变换前的噪声值
    %进行函数变换
    c(1)=Fadd(b(1));
    c(2)=Fequal(b(2));
    c(3)=Fadd(b(3));
    c(4)=Fequal(b(4));
    c=[c(1),c(2),c(3),c(4)];
    %计算未进行函数变换后的噪声值
    s=disc(c);
    if(s>r)
        R=R+1;   %计算R（M）值
    elseif(s<r)
        S=S+1;   %计算S（M）值
    else
        U=U+1;
    end
%进行-M伪装变换 并计算相应的参数
    mm=[-1,0,-1,0];
    MR=0;
    MS=0;
    MU=0;
    d(1)=Fminus(b(1));
    d(2)=Fequal(b(2));
    d(3)=Fminus(b(3));
    d(4)=Fequal(b(4));
    d=[d(1),d(2),d(3),d(4)];
    %计算未进行函数变换后的噪声值
    v=disc(d);
    if(v>r)
        MR=MR+1;   %计算R（-M）值
    elseif(v<r)
        MS=MS+1;   %计算S（-M）值
    else
        MU=MU+1;
    end
end