基于MATAB的GPS信号仿真完整源代码123.docVIP

配套毕业设计论文见百度文库 请搜索 《基于MATLAB的GPS信号仿真123》 附录C 仿真程序代码 数据码的产生 function datacode=data(x) y=rand(1,x); for i=1:x if y(i)0.5 datacode(i)=0; else datacode(i)=1; end end y(1)=0; show2(1)=datacode(1); q=2; for i=1:length(datacode) for j=1:100 y(q)=i-1+j*0.01; show2(q)=datacode(i); q=q+1; end end plot(y,show2); axis([0 length(datacode) -0.2 1.2]); C/A码的产生及扩频调制 clc; c=input(请输入数据码的长度：c=); y=rand(1,c); for i=1:c if y(i)0.5 datacode(i)=0; else datacode(i)=1; end end x(1)=0; show(1)=datacode(1); p=2; for i=1:c for j=1:100 x(p)=i-1+j*0.01; show(p)=datacode(i); p=p+1; end end subplot(4,1,1); plot(x,show); title(数据码); axis([0 c -0.2 1.2]); number=input(请输入卫星PRN号码:number=); cacode=CAgenerate(number); temp=cacode(1:100) x(1)=0; show(1)=temp(1); p=2; %下面的循环是为了将结果显示成方波形式 for i=1:length(temp) for j=1:100 x(p)=i-1+j*0.01; show(p)=temp(i); p=p+1; end end %画出仿真结果图 subplot(4,1,2); plot(x,show); title(C/A码); axis([0 100 -0.2 1.2]); %截取CA码的前十个数据进行扩频，每个数据插入5个CA序列 cacode1=cacode(1:10); for i=1:c if datacode(i)==1 datacodek((i-1)*50+1:i*50)=ones(1,50); else datacodek((i-1)*50+1:i*50)=zeros(1,50); end end for i=1:c for j=1:50 addr=rem(((i-1)*50+j),10); if addr==0 addr=10; end kuopindata((i-1)*50+j)=xor(datacodek((i-1)*50+j),cacode1(addr)); end end %下面的循环是为了将结果显示成方波形式 x(1)=0; show(1)=kuopindata(1); p=2; for i=1:length(kuopindata) for j=1:100 x(p)=i-1+j*0.01; show(p)=kuopindata(i); p=p+1; end end subplot(4,1,3); plot(x,show); title(扩频数据); axis([0 length(kuopindata) -0.2 1.2]); %每位数据通过正弦波来调制 Sinwave=sin([0:2*pi/8:2*pi*7/8]); Sinwave=single(Sinwave); GPSsignal=zeros(1,1); Sinwave=[Sinwa