三点法隐写代码

function hide（orifilemsg1outfilelposrinterab）
%参数说明：
% orifile：原始载体音频信号
% msg1：待隐藏的信息
% outfile:隐藏后的音频
% l：音频信号做DCT变换后的分段长度
% pos：在每段中嵌入的起始位置
% r: 每段选出的系数
% inter：每段中，相邻两个嵌入点的间隔
% a：量化步长
% b：最邻近的量化区间的距离
% msg1：待隐藏的信息
% outfile:隐藏后的音频
fi = fopen（msg1‘r‘）;
    [msgN] = fread（fi‘ubit1‘）;
    fclose（fi）;  
    [sampledataFSnbits]=wavread（orifile）;
    sdmax = max（max（max（sampledata））max（abs（min（sampledata））））;
    x = sampledata/sdmax;
    p=x;
    x1 = x（:1）;
    [Ln] = size（x1）;
    bn = floor（L/l）;
    if N>bn
        error（‘Message is too long‘）;
    end
      audct = zeros（L1）;
    for i=1:bn
        audct（（i-1）*l+1:i*l） = dct（x1（（i-1）*l+1:i*l））;
    end
       audctmax = max（max（max（audct））max（abs（min（audct））））;
    audct = audct/audctmax;
    for i=1:N        
        tmp = audct（（i-1）*l+1:i*l）;
        for j=1:r            
            t = tmp（pos+（j-1）*inter）;
            si = floor（t/a）;
            mid = （