K = max(f_saved); q = zeros(1,106); q_index = 1; dq = zeros(1,256); disp(strcat('f(1) (original) =',num2str(f(1)))); %kvantificer den foerste koefficient q(1) = sign(f(1)) * round((abs(f(1))/K)*2^M); disp(strcat('f(1) (kvantificeret)=',num2str(q(1)))); %dekvantificer den foerste koefficient dq(1) = sign(q(1)) * (abs(q(1))/K) * 2^M; disp(strcat('f(1) (dekvantificeret)=',num2str(dq(1)))); % Udregn forskellene mellem hver koeficcient r = zeros(1,106); r(1) = f(1) - sign(q(1))*(abs(q(1))/2^M)*K; disp(strcat('r=',num2str(r(1)))); ------- Giver følgende output på sinus-signalet ------- f(1) (original) =1.3009 f(1) (kvantificeret)=17 f(1) (dekvantificeret)=219.9295 r=-0.013164 f(1) (original) =-8.9344 f(1) (kvantificeret)=-751 f(1) (dekvantificeret)=-63105.3573 r=0.0030196 f(1) (original) =2.5002 f(1) (kvantificeret)=65 f(1) (dekvantificeret)=1678.7769 r=-0.01647 f(1) (original) =9.7312 f(1) (kvantificeret)=73 f(1) (dekvantificeret)=550.9809 r=0.059369 f(1) (original) =-4.5677 f(1) (kvantificeret)=-121 f(1) (dekvantificeret)=-3213.2329 r=-0.011206 f(1) (original) =-8.3084 f(1) (kvantificeret)=-429 f(1) (dekvantificeret)=-22168.6286 r=-0.0065243 f(1) (original) =6.0301 f(1) (kvantificeret)=82 f(1) (dekvantificeret)=1110.8556 r=-0.022914