dpcm_function.m

利用MATLAB实现正交滤波器组的设计。

%这个是DPCM编码函数

function [e,r] = DPCM(X)
[m,n]=size(X);

r1=0;r2=0;r3=0;                    %初始化三个初值
for i=1:m
    r1=r1+X(i,1)^2+X(i,2)^2;          %做累加过程
    r2=r2+X(i,1)*X(i,2);              %做累加过程
    for j=3:n
        r1=r1+X(i,j)^2;                %做累加过程
        r2=r2+X(i,j)*X(i,j-1);         %做累加过程
        r3=r3+X(i,j)*X(i,j-2);         %做累加过程
    end
end
r1=r1/(m*n);r2=r2/(m*n);r3=r3/(m*n);   %线性差分编码 运算
a1=r2*(r1-r3)/(r1^2-r2^2);             %线性差分编码 运算
a2=(r3*r1-r2^2)/(r1^2-r2^2);           %线性差分编码 运算
Y=zeros(m,n);
c=zeros(m,n);

for i=1:m
    for j=3:n
        Y(i,j)=a1*X(i,j-1)+a2*X(i,j-2); %数值权值相乘
        c(i,j)=X(i,j)-Y(i,j);           %做差
    end
end
for i=2:m
    Y(i,1)=a1*X(i-1,1)+a2*X(i-1,2);     %数值权值相乘
    c(i,1)=X(i,1)-Y(i,1);               %做差
    Y(i,2)=a1*X(i-1,2)+a2*X(i-1,1);     %数值权值相乘
    c(i,2)=X(i,2)-Y(i,2);               %做差
end
c(1,1)=0;c(1,2)=0;
rc=abs(c);
rc=reshape(rc,m*n,1);
rc=sort(rc);
bit=8;
for i=1:(2^(bit-1))
    q=round(m*n*i/(2^(bit-1)));
    d(i)=rc(q);
end
e=zeros(m,n);
for i=2:(2^(bit-1))
    [k1,k2]=find(abs(c)<=d(i)&abs(c)>d(i-1));                 %对运算后的数值进行筛选
    q1=length(k1);
    for j=1:q1
        if(c(k1(j),k2(j))>=0) e(k1(j),k2(j))=(d(i)+d(i-1))/2; %取均值
        else e(k1(j),k2(j))=-(d(i)+d(i-1))/2;                 %取均值
        end
    end
end
[k1,k2]=find(abs(c)<=d(1));      
q=length(k1);
for j=1:q
    if(c(k1(j),k2(j))>=0) e(k1(j),k2(j))=d(1)/2;
    else e(k1(j),k2(j))=-d(1)/2;
    end
end
e(1,1)=X(1,1);
e(1,2)=X(1,2);
e=round(e*(2^(bit-1)))/(2^(bit-1));


r=1-length(k1)/length(X);