js.asv

该代码实现了am

 % d=abs(randn(1,900))              
                           %k=randn(1,1)
k=1
%b=1
 %c=1                      
 d22=abs(randn(1,160));
d1=abs(randn(1,900));
e=abs(randn(1,202));
% preprocess;
 
                            for   j=1:4;                                    %1
                                k1(1,j)=(k+j*40)-41;
                               
                            end  
                            
                    for n=1:39;
                        for L=40:1:120
                            KK=sum(d22(1:n))        
                            b(1,:)=KK+k1(1,1)-L
                        end
                    end
                           
                    for n1=0:39;
                         KK1=sum(d22(1:n1))
                          R0=KK1*b
               end
                  
                  z1=max(R0);                                        %第一帧7比特量化
                   
             
                
             for j=1:4;                             %max  l
                 k1(1,:)=k+j*40;
             end
               for l=40:120;
                 for i=120:1;
                     w(1,:)=k1(1,:)+i;
                     k2(1,:)=(k1(1,:)+i)-l;
                 end
               end
               for i=120:1;
                   for j=1:4;
                 w(1,:)=k2(1,:);
                   
               y(1,:)=max(w(1,:));
                k1(1,:)=k+j*40;
                l=w(1,:)-y(1);
                   end
               end
              
              
              
            for j=1:4;                                                      %第三步      3
                k1(1,:)=k+j*40;
            end
            for i=1:40;
                
       %      e(1,:)=k1(1,:)+i-l
       % end
               
                      f(1,:)=sum(d1(1,k1(1,:)+i-l)^2);
                  end
                      for l=40:120;                                %111111
                        for  j=1:4;
                             for  i=120:1;
                                k1(1,:)=k+j*40;
                                 w(1,:)=k1(1,:)+i;
                       bj(1,:)=w(1,:)/f(1,:);
                   end
               end
           end
                 
                 
             for  j=1:4;
                   for k9=1:40;
                        k1(1,:)=k+j*40;
                        d22(1,k1+k9)=f(1,:)*d22(1,k1(1,:)+k+l);                             %14
                        d1(1,k1+k9)=e(1,k1+k9)+d22(1,k1+k9);;                                 %15
                    end
                end
                

                
                
                  k=1 ;                        %%%%%%%%%% RPE
      %e=rand(1,200)
      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  k=1:40
      for  i=1:11
          q=k+11-i;
      if q>46
       e(1,q)=0;
      elseif q<7
       e(1,q)=0;
      elseif q>=7 & q<=46
      t(1,:)=e(1,q);
      end
     end
      x(1,:)=sum(h(1,i).*t(1,:));
     end
   % end
    for  j=1:4;                            %17
       for i=1:13;
        %for m=1:4;
         k1(1,:)=k+j*40;
                   x10(1,:)=k1(1,:)+3*i  ;                      %  kj+m+3*i
                   x11(1,:)=k1(1,:)+1+3*i;
                   x12(1,:)=k1(1,:)+2+3*i;
                   x13(1,:)=k1(1,:)+3+3*i;
                   
                   x0(1,:)=x(1,:).*x10(1,:);
                   x1(1,:)=x(1,:).*x11(1,:);
                   x2(1,:)=x(1,:).*x12(1,:);
                   x3(1,:)=x(1,:).*x13(1,:);
                   Am1(1,:)=max(x0(1,:),x1(1,:));
                   Am2(1,:)=max(x2(1,:),x3(1,:));
                   Am(1,:)=max(Am1(1,:),Am2(1,:));         
       end
    end
    
    
            for  i=1:13                                                   %  18
                  E10(1,:)=sum(x0(1,:)^2);
                  E11(1,:)=sum(x1(1,:)^2);
                  E12(1,:)=sum(x2(1,:)^2);
                  E13(1,:)=sum(x2(1,:)^2);
                 
                  Em1(1,:)=max(E10(1,:),E11(1,:)) ;
                  Em2(1,:)=max(E12(1,:),E13(1,:));
                   
                  Em(1,:)=max(Em1(1,:),Em2(1,:));
            end
            
                  for  i=1:13;                                        %  19
                     X0(1,:)=x0(1,:)/Am(1,:);
                     X1(1,:)=x1(1,:)/Am(1,:);
                     X2(1,:)=x2(1,:)/Am(1,:);
                     X3(1,:)=x3(1,:)/Am(1,:);
                  end
                  
                       for  i=1:13;                                  %   20
                         Xm0(1,:)=X0(1,:).*Am(1,:);
                         Xm1(1,:)=X1(1,:).*Am(1,:);
                         Xm2(1,:)=X2(1,:).*Am(1,:);
                         Xm3(1,:)=X3(1,:).*Am(1,:);
                         
                       end
                       
             for  j=1:4;
                % for   k9=1:40;
                 for  m=1:4;
                 for  i=1:13;
                     k9=m+3*i;
                     k1(1,:)=k+j*40;
                     if  k9<4
                         e(1,k1+k9)=0;
                     elseif  k9>43
                          e(1,k1+k9)=0;
                     elseif  k9>=4 & k9<=43
                              e(1,k1+k9)=Xm0(1,:);
                              
                     end
                 end
                 end
             end
             
                  for  j=1:4;
                % for   k9=1:40;
                 for  m=1:4;
                 for  i=1:13;
                     k9=m+3*i
                     k1(1,:)=k+j*40;
                     if  k9<4
                         e(1,k1+k9)=0;
                     elseif  k9>43
                          e(1,k1+k9)=0;
                     elseif  k9>=4 & k9<=43
                              e(1,k1+k9)=Xm1(1,:);
                              
                     end
                 end
                 end
             end         
                           
                   for  j=1:4;
                % for   k9=1:40;
                 for  m=1:4;
                 for  i=1:13;
                     k9=m+3*i;
                     k1(1,:)=k+j*40;
                     if  k9<4
                         e(1,k1+k9)=0;
                     elseif  k9>43
                          e(1,k1+k9)=0;
                     elseif  k9>=4 & k9<=43
                              e(1,k1+k9)=Xm2(1,:);
                              
                     end
                 end
                 end
             end                  
                           
                  for  j=1:4;
                % for   k9=1:40;
                 for  m=1:4;
                 for  i=1:13;
                     k9=m+3*i;
                     k1(1,:)=k+j*40;
                     if  k9<4
                         e(1,k1+k9)=0;
                     elseif  k9>43
                          e(1,k1+k9)=0;
                     elseif  k9>=4 & k9<=43
                              e(1,k1+k9)=Xm3(1,:);
                              
                     end
                 end
                 end
             end