codemain3.m

该代码实现了am

% function code_f(r)% preprocess_f(r)

k=1
preprocess;

                           for   j=1:4;                                    %1
                                k1(1,j)=(k+j*40)-41;
                               
                            end  
                            
                    for na=1:40;                                       %第一帧
                        for L=40:1:120
                            KKa=sum(d22(1:na));        
                            ba(1,:)=KKa+k1(1,1)-L;
                        end
                    end
                           
                    for n1=2:40;
                         KK1=sum(d22(1:n1));
                          R1=KK1*ba;
               end
                  
                  zz1=max(R1);                                        %第一帧7比特量化
                  zz10=round(zz1);
                
            bb1=zz1/(KKa+k1(1,1)+zz1/KK1)/(KKa+k1(1,1)+zz1/KK1);      %第一帧2比特量化
                      bb10=round(bb1);
                      
                      
                       for nb=41:80;                                       %第二帧
                        for L=40:1:120
                            KKb=sum(d22(1:nb));        
                            bb(1,:)=KKb+k1(1,2)-L;
                        end
                    end
                           
                    for n2=42:80;
                         KK2=sum(d22(1:n2));
                          R2=KK2*bb;
               end
                  
                  zz2=max(R2);                                        %第二帧7比特量化
                   zz20=round(zz2);
                
             bb2=zz2/(KKb+k1(1,2)+zz2/KK2)/(KKb+k1(1,2)+zz2/KK2);     %第二帧2比特量化
                      bb20=round(bb2);
                      
                       for nc=81:120;                                       %第三帧
                        for L=40:1:120
                            KKc=sum(d22(1:nc));        
                            bc(1,:)=KKc+k1(1,3)-L;
                        end
                    end
                           
                    for n3=82:120;
                         KK3=sum(d22(1:n3));
                          R3=KK3*bc;
               end
                  
                  zz3=max(R3);                                        %第三帧7比特量化
                   zz30=round(zz3);
                
          bb3=zz3/(KKc+k1(1,3)+zz3/KK3)/(KKc+k1(1,3)+zz3/KK3);         %第三帧2比特量化
                      bb30=round(bb3);
                      
                      
                        for nd=121:160;                                       %第四帧
                        for L=40:1:120
                            KKd=sum(d22(1:nd));        
                            bd(1,:)=KKd+k1(1,4)-L;
                        end
                    end
                           
                    for n4=122:160;
                         KK4=sum(d22(1:n4));
                          R4=KK4*bd;
               end
                  
                  zz4=max(R4);                                        %第四帧7比特量化
                   zz40=round(zz4);
                
                bb4=zz4/(KKd+k1(1,4)+zz4/KK4)/(KKd+k1(1,4)+zz4/KK4);  %第四帧2比特量化
                      bb40=round(bb4);
                  
                        
                   dz1(1,:)=bb10*(K(1,1)+KKa-zz1/KK1);
                   
                   dz2(1,:)=bb20*(K(1,2)+KKb-zz2/KK2);
                  
                   dz3(1,:)=bb30*(K(1,3)+KKc-zz3/KK3);
                   
                   dz4(1,:)=bb40*(K(1,4)+KKd-zz4/KK4);
                  
                   ee1(1,:)=KK1-dz1;
                  
                   ee2(1,:)=KK2-dz2;
                   
                   ee3(1,:)=KK3-dz3;
                   
                   ee4(1,:)=KK4-dz4;
                   
                 
                   h=zeros(1,11);
                   h(1,1)=-134/(2^13);
                   h(1,2)=-374/(2^13);
                   h(1,3)=0;
                   h(1,4)=2054/(2^13);
                   h(1,5)=5741/(2^13);
                   h(1,6)=8192/(2^13);
                   h(1,7)=5741/(2^13);
                   h(1,8)=2054/(2^13);
                   h(1,9)=0;
                   h(1,10)=-374/(2^13);
                   h(1,11)=-134/(2^13);
                   
                  
                   for  na=1:40
              for  ia=1:11
                      q=na+5-ia;
                         if q>39
                          q=0;
                         elseif q<0
                             q=0;
                        elseif q>=0 & q<=39
                         end
                       end
                  x(1,:)=sum(h(1,ia).*q);
                   end
                  
               
                   for nn=1:13
                       xm1(1,nn)=k1(1,1)+1+3*nn;
                       XM1(1,nn)=xm1(1,nn)^2;
                  
                       xm2(1,nn)=k1(1,2)+2+3*nn;
                       XM2(1,nn)=xm2(1,nn)^2;
                   
                       xm3(1,nn)=k1(1,3)+3+3*nn;
                       XM3(1,nn)=xm3(1,nn)^2;
                   
                       xm4(1,nn)=k1(1,4)+4+3*nn;
                       XM4(1,nn)=xm4(1,nn)^2;
                   
                   end
                     
                      XM10=max(XM1);                               %  第一帧2 bits量化
                
                      XM20=max(XM2);                               %  第二帧2 bits量化
                  
                      XM30=max(XM3);                               %  第三帧2 bits量化
                  
                      XM40=max(XM4);                               %  第四帧2 bits量化
                  
                  
                      xm10=max(xm1);                               %  第一帧6 bits量化
                      
                      xm20=max(xm2);                               %  第二帧6 bits量化
                      
                      xm30=max(xm3);                               %  第三帧6 bits量化                              
                      
                      xm40=max(xm4);                               %  第四帧6 bits量化
                      
                      for nn=1:13
                          XX1(1,nn)=xm1(1,nn)/xm10*32767;
                          XX2(1,nn)=xm2(1,nn)/xm20*32767;
                          XX3(1,nn)=xm3(1,nn)/xm30*32767;
                          XX4(1,nn)=xm4(1,nn)/xm40*32767;
                          
                          XX10=round(XX1);                         %  第一帧3 bits量化
                          XX20=round(XX2);                         %  第二帧3 bits量化
                          XX30=round(XX3);                         %  第三帧3 bits量化
                          XX40=round(XX4);                         %  第四帧3 bits量化
                      end        
    rptltppcm;                  
    frame1pcm;                  
    frame2pcm;
    frame3pcm;
    frame4pcm;                           
    rerank;
    juanjicode;