function z = prs

% Generate phase reference symbol

%DAB standard, table 39, pp. 148
K = 1536;

% k = k', i, n
t39 = zeros(1536,3);
o = K/2+1;

for k = -768:-737
    t39(k+o:k+o,:,:) = [-768,0,1];
end

for k = -736:-705
    t39(k+o:k+o,:,:) = [-736,1,2];
end

for k = -704:-673
    t39(k+o:k+o,:,:) = [-704,2,0];
end

for k = -672:-641
    t39(k+o:k+o,:,:) = [-672,3,1];
end

for k = -640:-609
    t39(k+o:k+o,:,:) = [-640,0,3];
end

for k = -608:-577
    t39(k+o:k+o,:,:) = [-608,1,2];
end

for k = -576:-545
    t39(k+o:k+o,:,:) = [-576,2,2];
end

for k = -544:-513
    t39(k+o:k+o,:,:) = [-544,3,3];
end

for k = -512:-481
    t39(k+o:k+o,:,:) = [-512,0,2];
end

for k = -480:-449
    t39(k+o:k+o,:,:) = [-480,1,1];
end

for k = -448:-417
    t39(k+o:k+o,:,:) = [-448,2,2];
end

for k = -416:-385
    t39(k+o:k+o,:,:) = [-416,3,3];
end

for k = -384:-353
    t39(k+o:k+o,:,:) = [-384,0,1];
end

for k = -352:-321
    t39(k+o:k+o,:,:) = [-352,1,2];
end

for k = -320:-289
    t39(k+o:k+o,:,:) = [-320,2,3];
end

for k = -288:-257
    t39(k+o:k+o,:,:) = [-288,3,3];
end

for k = -256:-225
    t39(k+o:k+o,:,:) = [-256,0,2];
end

for k = -224:-193
    t39(k+o:k+o,:,:) = [-224,1,2];
end

for k = -192:-161
    t39(k+o:k+o,:,:) = [-192,2,2];
end

for k = -160:-129
    t39(k+o:k+o,:,:) = [-160,3,1];
end

for k = -128:-97
    t39(k+o:k+o,:,:) = [-128,0,1];
end

for k = -96:-65
    t39(k+o:k+o,:,:) = [-96,1,3];
end

for k = -64:-33
    t39(k+o:k+o,:,:) = [-64,2,1];
end

for k = -32:-1
    t39(k+o:k+o,:,:) = [-32,3,2];
end

for k = 0:0
    t39(k+o:k+o,:,:) = [0, 0, 0];
end    

for k = 1:32
    t39(k+o:k+o,:,:) = [1,0,3];
end

for k = 33:64
    t39(k+o:k+o,:,:) = [33,3,1];
end

for k = 65:96
    t39(k+o:k+o,:,:) = [65,2,1];
end

for k = 97:128
    t39(k+o:k+o,:,:) = [97,1,1];
end

for k = 129:160
    t39(k+o:k+o,:,:) = [129,0,2];
end

for k = 161:192
    t39(k+o:k+o,:,:) = [161,3,2];
end

for k = 193:224
    t39(k+o:k+o,:,:) = [193,2,1];
end

for k = 225:256
    t39(k+o:k+o,:,:) = [225,1,0];
end

for k = 257:288
    t39(k+o:k+o,:,:) = [257,0,2];
end

for k = 289:320
    t39(k+o:k+o,:,:) = [289,3,2];
end

for k = 321:352
    t39(k+o:k+o,:,:) = [321,2,3];
end

for k = 353:384
    t39(k+o:k+o,:,:) = [353,1,3];
end

for k = 385:416
    t39(k+o:k+o,:,:) = [385,0,0];
end

for k = 417:448
    t39(k+o:k+o,:,:) = [417,3,2];
end

for k = 449:480
    t39(k+o:k+o,:,:) = [449,2,1];
end

for k = 481:512
    t39(k+o:k+o,:,:) = [481,1,3];
end

for k = 513:544
    t39(k+o:k+o,:,:) = [513,0,3];
end

for k = 545:576
    t39(k+o:k+o,:,:) = [545,3,3];
end

for k = 577:608
    t39(k+o:k+o,:,:) = [577,2,3];
end

for k = 609:640
    t39(k+o:k+o,:,:) = [609,1,0];
end

for k = 641:672
    t39(k+o:k+o,:,:) = [641,0,3];
end

for k = 673:704
    t39(k+o:k+o,:,:) = [673,3,0];
end

for k = 705:736
    t39(k+o:k+o,:,:) = [705,2,1];
end

for k = 737:768
    t39(k+o:k+o,:,:) = [737,1,1];
end

o=1;
% h(i,j)
h(0+o,0+o) = 0;
h(1+o,0+o) = 0;
h(2+o,0+o) = 0;
h(3+o,0+o) = 0;

h(0+o,1+o) = 2;
h(1+o,1+o) = 3;
h(2+o,1+o) = 0;
h(3+o,1+o) = 1;

h(0+o,2+o) = 0;
h(1+o,2+o) = 2;
h(2+o,2+o) = 0;
h(3+o,2+o) = 2;

h(0+o,3+o) = 0;
h(1+o,3+o) = 3;
h(2+o,3+o) = 2;
h(3+o,3+o) = 1;

h(0+o,4+o) = 0;
h(1+o,4+o) = 0;
h(2+o,4+o) = 0;
h(3+o,4+o) = 0;

h(0+o,5+o) = 0;
h(1+o,5+o) = 1;
h(2+o,5+o) = 2;
h(3+o,5+o) = 3;

h(0+o,6+o) = 1;
h(1+o,6+o) = 3;
h(2+o,6+o) = 1;
h(3+o,6+o) = 3;

h(0+o,7+o) = 1;
h(1+o,7+o) = 0;
h(2+o,7+o) = 3;
h(3+o,7+o) = 2;

h(0+o,8+o) = 2;
h(1+o,8+o) = 2;
h(2+o,8+o) = 2;
h(3+o,8+o) = 2;

h(0+o,9+o) = 0;
h(1+o,9+o) = 1;
h(2+o,9+o) = 2;
h(3+o,9+o) = 3;

h(0+o,10+o) = 0;
h(1+o,10+o) = 2;
h(2+o,10+o) = 0;
h(3+o,10+o) = 2;

h(0+o,11+o) = 0;
h(1+o,11+o) = 3;
h(2+o,11+o) = 2;
h(3+o,11+o) = 1;

h(0+o,12+o) = 2;
h(1+o,12+o) = 2;
h(2+o,12+o) = 2;
h(3+o,12+o) = 2;

h(0+o,13+o) = 2;
h(1+o,13+o) = 3;
h(2+o,13+o) = 0;
h(3+o,13+o) = 1;

h(0+o,14+o) = 1;
h(1+o,14+o) = 3;
h(2+o,14+o) = 1;
h(3+o,14+o) = 3;

h(0+o,15+o) = 1;
h(1+o,15+o) = 0;
h(2+o,15+o) = 3;
h(3+o,15+o) = 2;

h(0+o,16+o) = 0;
h(1+o,16+o) = 0;
h(2+o,16+o) = 0;
h(3+o,16+o) = 0;

h(0+o,17+o) = 2;
h(1+o,17+o) = 3;
h(2+o,17+o) = 0;
h(3+o,17+o) = 1;

h(0+o,18+o) = 0;
h(1+o,18+o) = 2;
h(2+o,18+o) = 0;
h(3+o,18+o) = 2;

h(0+o,19+o) = 0;
h(1+o,19+o) = 3;
h(2+o,19+o) = 2;
h(3+o,19+o) = 1;

h(0+o,20+o) = 0;
h(1+o,20+o) = 0;
h(2+o,20+o) = 0;
h(3+o,20+o) = 0;

h(0+o,21+o) = 0;
h(1+o,21+o) = 1;
h(2+o,21+o) = 2;
h(3+o,21+o) = 3;

h(0+o,22+o) = 1;
h(1+o,22+o) = 3;
h(2+o,22+o) = 1;
h(3+o,22+o) = 3;

h(0+o,23+o) = 1;
h(1+o,23+o) = 0;
h(2+o,23+o) = 3;
h(3+o,23+o) = 2;

h(0+o,24+o) = 2;
h(1+o,24+o) = 2;
h(2+o,24+o) = 2;
h(3+o,24+o) = 2;

h(0+o,25+o) = 0;
h(1+o,25+o) = 1;
h(2+o,25+o) = 2;
h(3+o,25+o) = 3;

h(0+o,26+o) = 0;
h(1+o,26+o) = 2;
h(2+o,26+o) = 0;
h(3+o,26+o) = 2;

h(0+o,27+o) = 0;
h(1+o,27+o) = 3;
h(2+o,27+o) = 2;
h(3+o,27+o) = 1;

h(0+o,28+o) = 2;
h(1+o,28+o) = 2;
h(2+o,28+o) = 2;
h(3+o,28+o) = 2;

h(0+o,29+o) = 2;
h(1+o,29+o) = 3;
h(2+o,29+o) = 0;
h(3+o,29+o) = 1;

h(0+o,30+o) = 1;
h(1+o,30+o) = 3;
h(2+o,30+o) = 1;
h(3+o,30+o) = 3;

h(0+o,31+o) = 1;
h(1+o,31+o) = 0;
h(2+o,31+o) = 3;
h(3+o,31+o) = 2;

o = K/2+1;
for k = -K/2:1:K/2
    if k~=0
        ks = t39(k+o,1,1);
        i = t39(k+o,2,1);
        n = t39(k+o,3,1);
        phi = pi/2 * ((h(i+1,k-ks+1)+n));
        z(k+o) = exp(j * phi);
    end
    if k==0
        z(k+o) = 0;
    end 
end

end