【正文】 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(1)+rb(8),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==8。 %第八顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(2)+rb(9),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(2)+rb(9),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==9。 %第九顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(3)+rb(10),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(3)+rb(10),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==10。 %第十顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(2)+rb(3),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(2)+rb(3),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==11。 %第十一顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(3)+rb(4),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(3)+rb(4),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==12。 %第十二顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(5)+rb(6),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(5)+rb(6),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==13。 %第十三顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(6)+rb(7),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(6)+rb(7),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。end if a==14。 %第十四顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(7)+rb(8),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(7)+rb(8),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==15。 %第十五顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(8)+rb(9),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(8)+rb(9),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==16。 %第十六顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(9)+rb(10),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(9)+rb(10),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==17。 %第十七顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(1)+rb(4),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(1)+rb(4),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==18。 %第十八顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(2)+rb(5),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(2)+rb(5),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==19。 %第十九顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(3)+rb(6),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(3)+rb(6),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==20。 %第二十顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(4)+rb(7),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(4)+rb(7),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==21。 %第二十一顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(5)+rb(8),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(5)+rb(8),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==22。 %第二十二顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(6)+rb(9),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(6)+rb(9),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==23。 %第二十三顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(1)+rb(3),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。 rc(1)=[mod(rb(1)+rb(3),2)]。endrd=rd39。%ra=ra39。z=rd。plot(z)save ca_code z。endif a==24。 %第二十四顆衛(wèi)星ra=ones(1,10)。rb=ones(1,10)。m=ones(1,3)。rc(1)=[mod(rb(4)+rb(6),2)]。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。