色散熵的主要代码

function [Out_DispEn npdf]=DispEn（xmncMAtau）

N=length（x）;
sigma_x=std（x）; %求标准偏差
mu_x=mean（x）;  %求均值

%  nc为归一化相关系数   nc=NC_wmk（AB） A是原始信号，B是提取的水印信息必须是double类型

switch MA
    case   ‘LM‘
        y=mapminmax（x01）;%归一化  归一到0--1
        y（y==1）=1-1e-10;
        y（y==0）=1e-10;
        z=round（y*nc+0.5）;%取整函数
        
    case ‘NCDF‘
        y=normcdf（xmu_xsigma_x）;%正态分布概率密度函数
        y=mapminmax（y01）;
        y（y==1）=1-1e-10;
        y（y==0）=1e-10;
        z=round（y*nc+0.5）;
        
    case ‘LOGSIG‘
        y=logsig（（x-mu_x）/sigma_x）;
        y=mapminmax（y01）;
        y（y==1）=1-1e-10;
        y（y==0）=1e-10;
        z=round（y*nc+0.5）;
        
    case ‘TANSIG‘
        y=tansig（（x-mu_x）/sigma_x）+1;
        y=mapminmax（y01）;
        y（y==1）=1-1e-10;
        y（y==0）=1e-10;