mainsimulation.m

标准GPS星座的仿真程序。注：采用24颗卫星

                  
              end
              
           elseif n == 5 
              p2=k;
              [a2 b2]=size(p2);
              
              for i7=1:a2
                  
                  for i8=1:b2
                      w2(i7,i8)=g(p2(i7,i8));
                  end
                  
              end
              
           elseif n == 6
              p3=k;
             [a3 b3]=size(p3);
              
             for i7=1:a3
                  
                 for i8=1:b3
                      w3(i7,i8)=g(p3(i7,i8));
                 end
                  
             end
              
           elseif n == 7 
              p4=k;
              [a4 b4]=size(p4);
             
              for i7=1:a4
                  
                  for i8=1:b4
                      w4(i7,i8)=g(p4(i7,i8));
                  end
                  
              end
              
           elseif n == 8 
              p5=k;
              [a5 b5]=size(p5);
              
              for i7=1:a5
                  
                  for i8=1:b5
                      w5(i7,i8)=g(p5(i7,i8));
                  end
                  
              end
              
            elseif n == 9
                w6=[g(1);g(2);g(3);g(4);g(5);g(6);g(7);g(8);g(9)];
              
           end 
          
         for i5=1:b
               
               for i6=1:a
                   x5=X(1,g(k(i6,i5)))-XJ;
                   y5=X(2,g(k(i6,i5)))-YJ;
                   z5=X(3,g(k(i6,i5)))-ZJ;
                   x6=x5/ S(k(i6,i5));
                   y6=y5/ S(k(i6,i5));
                   z6=z5/ S(k(i6,i5));
                   B(:,i6)=[x6,y6,z6,-1]';
               end
               
               A=B';
               Q=inv(A'*A);
               r=r+1;
               x7=Q(1,1);
               y7=Q(2,2);
               z7=Q(3,3);
               PDOP(r)=sqrt(x7+y7+z7);
               TDOP(r)=sqrt(Q(4,4));
               GDOP(r)=sqrt(PDOP(r)^2+TDOP(r)^2);
           end
          
      end
      
  case 10
       for i4=1:c
          sa(i4)=i4;
       end
       
       for n=4:10
           k=combinesingle(sa,n);
           [a b]=size(k);
           
           if n == 4
              p1=k;
              [a1 b1]=size(p1);
              
              for i7=1:a1
                  
                  for i8=1:b1
                      w1(i7,i8)=g(p1(i7,i8));
                  end
                  
              end
              
           elseif n == 5 
              p2=k;
              [a2 b2]=size(p2);
             
              for i7=1:a2
                 
                  for i8=1:b2
                      w2(i7,i8)=g(p2(i7,i8));
                  end
                  
              end
              
           elseif n == 6
              p3=k;
             [a3 b3]=size(p3);
              
             for i7=1:a3
                  
                 for i8=1:b3
                      w3(i7,i8)=g(p3(i7,i8));
                 end
                  
             end
              
           elseif n == 7 
              p4=k;
              [a4 b4]=size(p4);
              
              for i7=1:a4
                  
                  for i8=1:b4
                      w4(i7,i8)=g(p4(i7,i8));
                  end
                  
              end
              
           elseif n == 8 
              p5=k;
              [a5 b5]=size(p5);
              
              for i7=1:a5
                  
                  for i8=1:b5
                      w5(i7,i8)=g(p5(i7,i8));
                  end
                  
              end
              
           elseif n == 9
              p6=k;
              [a6 b6]=size(p6);
              
              for i7=1:a6
                  
                  for i8=1:b6
                      w6(i7,i8)=g(p6(i7,i8));
                  end
                  
              end
              
            elseif n == 10
                w7=[g(1);g(2);g(3);g(4);g(5);g(6);g(7);g(8);g(9);g(10)];
              
           end 
           
          for i5=1:b
               
               for i6=1:a
                   x5=X(1,g(k(i6,i5)))-XJ;
                   y5=X(2,g(k(i6,i5)))-YJ;
                   z5=X(3,g(k(i6,i5)))-ZJ;
                   x6=x5/ S(k(i6,i5));
                   y6=y5/ S(k(i6,i5));
                   z6=z5/ S(k(i6,i5));
                   B(:,i6)=[x6,y6,z6,-1]';
               end
               
               A=B';
               Q=inv(A'*A);
               r=r+1;
               x7=Q(1,1);
               y7=Q(2,2);
               z7=Q(3,3);
              PDOP(r)=sqrt(x7+y7+z7);
              TDOP(r)=sqrt(Q(4,4));
              GDOP(r)=sqrt(PDOP(r)^2+TDOP(r)^2);
           end
           
       end
     
   case 11
       for i4=1:c
           sa(i4)=i4;
       end
       
       for n=4:11
           k=combinesingle(sa,n);
           [a b]=size(k);
           
           if n == 4
              p1=k;
              [a1 b1]=size(p1);
              
              for i7=1:a1
                  
                  for i8=1:b1
                      w1(i7,i8)=g(p1(i7,i8));
                  end
                  
              end
              
           elseif n == 5 
              p2=k;
              [a2 b2]=size(p2);
              
              for i7=1:a2
                 
                  for i8=1:b2
                      w2(i7,i8)=g(p2(i7,i8));
                  end
                  
              end
              
           elseif n == 6
             p3=k;
             [a3 b3]=size(p3);
              
             for i7=1:a3
                  
                 for i8=1:b3
                      w3(i7,i8)=g(p3(i7,i8));
                 end
                  
             end
              
           elseif n == 7 
              p4=k;
              [a4 b4]=size(p4);
              
              for i7=1:a4
                  
                  for i8=1:b4
                      w4(i7,i8)=g(p4(i7,i8));
                  end
                  
              end
              
           elseif n == 8 
              p5=k;
              [a5 b5]=size(p5);
             
              for i7=1:a5
                  
                  for i8=1:b5
                      w5(i7,i8)=g(p5(i7,i8));
                  end
                  
              end
              
           elseif n == 9
              p6=k;
              [a6 b6]=size(p6);
              
              for i7=1:a6
                  
                  for i8=1:b6
                      w6(i7,i8)=g(p6(i7,i8));
                  end
                  
              end
              
           elseif n == 10
              p7=k;
              [a7 b7]=size(p7);
              
              for i7=1:a7
                 
                  for i8=1:b7
                      w7(i7,i8)=g(p7(i7,i8));
                  end
                  
              end
              
            elseif n ==11
                w8=[g(1);g(2);g(3);g(4);g(5);g(6);g(7);g(8);g(9);g(10);g(11)];
              
           end 
           
           for i5=1:b
               
               for i6=1:a
                   x5=X(1,g(k(i6,i5)))-XJ;
                   y5=X(2,g(k(i6,i5)))-YJ;
                   z5=X(3,g(k(i6,i5)))-ZJ;
                   x6=x5/ S(k(i6,i5));
                   y6=y5/ S(k(i6,i5));
                   z6=z5/ S(k(i6,i5));
                   B(:,i6)=[x6,y6,z6,-1]';
               end
               
               A=B';
               Q=inv(A'*A);
               r=r+1;
               x7=Q(1,1);
               y7=Q(2,2);
               z7=Q(3,3);
               PDOP(r)=sqrt(x7+y7+z7);
               TDOP(r)=sqrt(Q(4,4));
               GDOP(r)=sqrt(PDOP(r)^2+TDOP(r)^2);
           end
           
      end

  case 12
       for i4=1:c
          sa(i4)=i4;
       end
       
       for n=4:12
           k=combinesingle(sa,n);
           [a b]=size(k);
           
           if n == 4
              p1=k;
              [a1 b1]=size(p1);
              
              for i7=1:a1
                 
                  for i8=1:b1
                      w1(i7,i8)=g(p1(i7,i8));
                  end
                  
              end
              
           elseif n == 5 
              p2=k;
              [a2 b2]=size(p2);
              
              for i7=1:a2
                  
                  for i8=1:b2
                      w2(i7,i8)=g(p2(i7,i8));
                  end
                  
              end
              
           elseif n == 6
              p3=k;
             [a3 b3]=size(p3);
              
             for i7=1:a3
                 
                 for i8=1:b3
                      w3(i7,i8)=g(p3(i7,i8));
                 end
                  
             end
              
           elseif n == 7 
              p4=k;
              [a4 b4]=size(p4);
              
              for i7=1:a4
                 
                  for i8=1:b4
                      w4(i7,i8)=g(p4(i7,i8));
                  end
                  
              end
              
           elseif n == 8 
              p5=k;
              [a5 b5]=size(p5);
              
              for i7=1:a5
                  
                  for i8=1:b5
                      w5(i7,i8)=g(p5(i7,i8));
                  end
                  
              end
              
           elseif n == 9
              p6=k;
              [a6 b6]=size(p6);
              
              for i7=1:a6
                 
                  for i8=1:b6
                      w6(i7,i8)=g(p6(i7,i8));
                  end
                  
              end
              
           elseif n == 10
              p7=k;
              [a7 b7]=size(p7);
              
              for i7=1:a7
                 
                  for i8=1:b7
                      w7(i7,i8)=g(p7(i7,i8));
                  end