MATLAB (Analog - Digital)

그냥.

close all;
clear all;

%Llyod_Max quantization를 500번 루프의 반복
r = 500;

%1000개의 랜덤한 변수를 생성
v_m = randn(1,1000);
%생성한 랜덤변수의 sorting
t = sort(v_m);
%1000개의 생성된 랜덤변수에 다시 10level(10개)의 랜덤 변수를 생성
v_s = randn(1,10);
%10개의 랜덤변수를 sorting
v_s_s = sort(v_s);


%Llyod_Max quantization함수의 dody
for l_m = 1:r;
    rst = zeros(1,10);
    buf = zeros(1,9);
    cnt = 0;
    sum = 0;
    l_s = 1;
    l_f = 1;
   
    for i = 1:9;
        buf(i) = (v_s_s(i) + v_s_s(i+1))/2; %Llyod_Max condition 1
    end

    while l_s <1001; %각각의 경계사이에서의 값을 비교해나가면서 대표값을 찾는것
        if l_f <10;
            if t(l_s)<buf(l_f)
                cnt = cnt+1;
                sum = sum + t(l_s);
                l_s=l_s+1;
            else
                if cnt == 0;
                    rst(l_f) = v_s_s(l_f);
                else
                    rst(l_f) = sum/cnt;
                end
                l_f = l_f+1;
                sum = 0;
                cnt = 0;
            end
        else
            sum = sum+t(l_s);
            l_s = l_s+1;
        end
    end

    if cnt == 0;
         rst(10) = v_s_s(10);
    else
         rst(10) = sum/cnt;
    end
   
    for i = 1:10;
        v_s_s(i) = rst(i);
    end
end

%Distortion 함수
ex = sort(v_s)
D=((ex-rst).^2)./10;