Pastebin

  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