clear;
clc;

sinfai=zeros(4,8,32);                 %生成4*8*32三維數組
cosfai=zeros(4,8,32);
for d=1:32                            %d為隨機數據bn
    for L=0:3                         %L為象限
        for t=0:7                     %t為8比特抽樣
            m=change(d-1);            %生成bn
            ph=phase(t/8,m)+pi*L/2;   %計算相位
            cosfai(L+1,t+1,d)=cos(ph);
            sinfai(L+1,t+1,d)=sin(ph);%存入數組
            
        end
    end
end

%…………………………
s=1;
for n=1:1:32
    for L=1:1:4
        for j=1:1:8
Cos_table(s,:)=change2(cosfai(L,j,n));
Sin_table(s,:)=change2(sinfai(L,j,n));
s=s+1;
        end
    end
end

coslow(:,1:8)=Cos_table(:,3:10);                   %余弦表低8位
sinlow(:,1:8)=Sin_table(:,3:10);                   %正弦表低8位
csmiddle(:,1:2)=Sin_table(:,1:2);
csmiddle(:,3:4)=0;
csmiddle(:,5:6)=0;
csmiddle(:,7:8)=Cos_table(:,1:2);
coslow2=coslow';                                   %二進制余弦表和正弦表的生成
sinlow2=sinlow'; 
csmiddle2=csmiddle';
Transbincos=flipud(coslow2);
Transbinsin=flipud(sinlow2);
Transbincs=flipud(csmiddle2);
for i=1:1:1024
    TTransbincos(:,8*(i-1)+1)=Transbincos(:,i);
    TTransbinsin(:,8*(i-1)+1)=Transbinsin(:,i);
    TTransbincs(:,8*(i-1)+1)=Transbincs(:,i);
end
fid=fopen('sin.bin','wb');
fwrite(fid,TTransbinsin,'ubit1');                  %生成二進制正弦表文件（低8位）
fid=fopen('cos.bin','wb');
fwrite(fid,TTransbincos,'ubit1');                  %生成二進制余弦表文件（低8位）
fid=fopen('cs.bin','wb');
fwrite(fid,TTransbincs,'ubit1');                   %生成二進制正弦和余弦表文件（高2位）
fclose('all');