资源简介
gps信号捕获程序,matlab仿真FFT+PMF
代码片段和文件信息
samplesPerCode = round(settings.samplingFreq / ...
(settings.codeFreqBasis / settings.codeLength));
% Create two 1msec vectors of data to correlate with and one with zero DC
signal1 = longSignal(1 : samplesPerCode);
signal2 = longSignal(samplesPerCode+1 : 2*samplesPerCode);
% index=ceil(samplesPerCode*(1 : 40920)/40920);
signal0DC = longSignal - mean(longSignal);
% Find sampling period
ts = 1 / settings.samplingFreq;
% Find phase points of the local carrier wave
phasePoints = (0 : (samplesPerCode-1)) * 2 * pi * ts;
% Number of the frequency bins for the given acquisition band (500Hz steps)
numberOfFrqBins = round(settings.acqSearchBand * 2) + 1;
% Generate all C/A codes and sample them according to the sampling freq.
%--- Initialize arrays to speed up the code -------------------------------
% Search results of all frequency bins and code shifts (for one satellite)
results = zeros(numberOfFrqBins 2046);
% Carrier frequencies of the frequency bins
frqBins = zeros(1 numberOfFrqBins);
%--- Initialize acqResults ------------------------------------------------
% Carrier frequencies of detected signals
acqResults.carrFreq = zeros(1 32);
% C/A code phases of detected signals
acqResults.codePhase = zeros(1 32);
% Correlation peak ratios of the detected signals
acqResults.peakMetric = zeros(1 32);
IQ1_add=zeros(11023*2);
IQ2_add=zeros(11023*2);
fprintf(‘(‘);
% Perform search for all listed PRN numbers ...
for PRN = settings.acqSatelliteList
%% Correlate signals ======================================================
%--- Perform DFT of C/A code ------------------------------------------
caCodes=generateCAcode(PRN);
caCodes1=kron(caCodesones(12));
% caCodes2=[caCodes1 0 0];
caCodeFreqDom = conj(fft(caCodes1));
%--- Make the correlation for whole frequency band (for all freq. bins)
for frqBinIndex = 1:numberOfFrqBins
%--- Generate carrier wave frequency grid (0.5kHz step) -----------
frqBins(frqBinIndex) = settings.IF - ...
(settings.acqSearchBand/2) * 1000 + ...
0.5e3 * (frqBinIndex - 1);
%--- Generate local sine and cosine -------------------------------
sinCarr = sin(frqBins(frqBinIndex) * phasePoints);
cosCarr = cos(frqBins(frqBinIndex) * phasePoints);
%--- “Remove carrier“ from the signal -----------------------------
I1 = sinCarr .* signal1;
Q1 = cosCarr .* signal1;
I2 = sinCarr .* signal2;
Q2 = cosCarr .* signal2;
IQ1=I1 + j*Q1;
IQ2=I2 + j*Q2;
% IQ11=IQ1(:index);
% IQ22=IQ2(:index);
for i=1:1023*2
IQ1_add(i)=sum(IQ1((i-1)*8+1:8*i));
IQ2_add(i)=sum(IQ2((i-1)*8+1:8*i));
end
%--- Convert the baseband signal to frequency domain --------------
- 上一篇:电压型逆变器双闭环仿真
- 下一篇:空气参数计算器
相关资源
- 高灵敏度GPS接收机MATLAB仿真,附捕获
- 基于MATLAB的质点弹道计算与外弹道优
- 阵列天线的matlab仿真
- MATLAB 经典程序源代码大全
- MATLAB小波软阈值去噪代码33473
- 天线阵的波束形成在MATLAB仿真程序及
- 非线性SVM算法-matlab实现
- 《MATLAB 智能算法超级学习手册》-程序
- 组合导航matlab程序
- 读取txt文件内容matlab代码实现
- Matlab实现基于相关的模板匹配程序
- matlab优化工具箱讲解
- 基于MATLAB的快速傅里叶变换
- 光纤传输中的分布傅立叶算法matlab实
- 基于matlab的图像处理源程序
- matlab 椭圆拟合程序
- 算术编码解码matlab源代码
- optical_flow 光流法 matlab 实现程序
- 引导图像滤波器 Matlab实现
- 分形几何中一些经典图形的Matlab画法
- OFDM系统MATLAB仿真代码
- SVM工具箱(matlab中运行)
- 图像小波变换MatLab源代码
- LU分解的MATLAB实现
- 冈萨雷斯数字图像处理matlab版(第三
- 替代数据法的matlab程序
- 用matlab实现的多站定位系统性能仿真
- 通过不同方法进行粗糙集属性约简m
- k近邻算法matlab实现
- matlab识别系统
评论
共有 条评论