基于matlab的gps信号仿真完整源代码123

基于matlab的gps信号仿真完整源代码123

ID:3704428

大小:91.00 KB

页数:18页

时间:2017-11-23

基于matlab的gps信号仿真完整源代码123_第1页
基于matlab的gps信号仿真完整源代码123_第2页
基于matlab的gps信号仿真完整源代码123_第3页
基于matlab的gps信号仿真完整源代码123_第4页
基于matlab的gps信号仿真完整源代码123_第5页
资源描述:

《基于matlab的gps信号仿真完整源代码123》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库

1、配套毕业设计论文见百度文库请搜索《基于MATLAB的GPS信号仿真123》附录C仿真程序代码1、数据码的产生functiondatacode=data(x)y=rand(1,x);fori=1:xify(i)<0.5datacode(i)=0;elsedatacode(i)=1;endendy(1)=0;show2(1)=datacode(1);q=2;fori=1:length(datacode)forj=1:100y(q)=i-1+j*0.01;show2(q)=datacode(i);q=q+1;endendplot(y,show2)

2、;axis([0length(datacode)-0.21.2]);1、C/A码的产生及扩频调制clc;c=input('请输入数据码的长度:c=');y=rand(1,c);fori=1:cify(i)<0.5datacode(i)=0;elsedatacode(i)=1;endendx(1)=0;show(1)=datacode(1);p=2;fori=1:cforj=1:100x(p)=i-1+j*0.01;show(p)=datacode(i);p=p+1;endendsubplot(4,1,1);plot(x,show);titl

3、e('数据码');axis([0c-0.21.2]);number=input('请输入卫星PRN号码:number=');cacode=CAgenerate(number);temp=cacode(1:100)x(1)=0;show(1)=temp(1);p=2;%下面的循环是为了将结果显示成方波形式fori=1:length(temp)forj=1:100x(p)=i-1+j*0.01;show(p)=temp(i);p=p+1;endend%画出仿真结果图subplot(4,1,2);plot(x,show);title('C/A码'

4、);axis([0100-0.21.2]);%截取CA码的前十个数据进行扩频,每个数据插入5个CA序列cacode1=cacode(1:10);fori=1:cifdatacode(i)==1datacodek((i-1)*50+1:i*50)=ones(1,50);elsedatacodek((i-1)*50+1:i*50)=zeros(1,50);endendfori=1:cforj=1:50addr=rem(((i-1)*50+j),10);ifaddr==0addr=10;endkuopindata((i-1)*50+j)=xor(

5、datacodek((i-1)*50+j),cacode1(addr));endend%下面的循环是为了将结果显示成方波形式x(1)=0;show(1)=kuopindata(1);p=2;fori=1:length(kuopindata)forj=1:100x(p)=i-1+j*0.01;show(p)=kuopindata(i);p=p+1;endendsubplot(4,1,3);plot(x,show);title('扩频数据');axis([0length(kuopindata)-0.21.2]);%每位数据通过正弦波来调制Sin

6、wave=sin([0:2*pi/8:2*pi*7/8]);Sinwave=single(Sinwave);GPSsignal=zeros(1,1);Sinwave=[SinwaveSinwaveSinwaveSinwaveSinwave];fori=1:length(kuopindata)GPSsignal=[GPSsignalkuopindata(i)*Sinwave];endGPSsignal=GPSsignal(2:length(GPSsignal));subplot(4,1,4);plot(GPSsignal(1:500));ti

7、tle('调制后数据');C/A码产生的子程序CAgenerate:functioncacode=CAgenerate(number)if(number<1)

8、(number>37)disp('输入参数必须在1~37之间取值');returnendCACode=zeros(1,1023);%生成一个1*1023的零矩阵%设置寄存器初相Reg1=[1,1,1,1,1,1,1,1,1,1];Reg2=[1,1,1,1,1,1,1,1,1,1];%设置反馈点,1表示需要反馈gp1=[0,0,1,0,0,0,0,0,0,1];gp2=[0,1,1,

9、0,0,1,0,1,1,1];%抽头G2Table=[2,3,4,5,1,2,1,2,3,2,3,5,6,7,8,9,1,2,3,4,5,6,1,4,5,6,7,8,1,2,3,

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。