蒼翠的路

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 373.8
紅花: 1
帖子: 28
在線: 15小時
蟲號: 3579141
注冊: 2014-12-05
專業(yè): 通信理論與系統(tǒng)

[求助] matlab仿真實現(xiàn)三階鎖相環(huán)路跟蹤高動態(tài)多普勒變化出現(xiàn)問題，急需指點[有源碼] 已有1人參與

RT，最近在忙論文，但仿真出現(xiàn)問題，結(jié)果總是不對。用的模型是JPL的高動態(tài)模型，環(huán)路帶寬根據(jù)跟蹤門限設(shè)為90Hz，積分時間為1ms，但是好像就是一直沒法鎖相，鑒相誤差一直較大，有沒做過這方面的同學或前輩啊，急需指點和討論！！
clc;
clear all;

Chip_rate = 6.3*10^6;       % 擴頻碼速率(chips/s)
Spreadcode_len = 63;          % 擴頻碼長度(chips)
Data_len = 240;             % 數(shù)據(jù)長度(symbols)
Freq_carrier = 4*10^10;       % 射頻載波頻率(Hz)
Freq_sample = 6.3*10^6;       % 接收機采樣頻率(Hz)
Symbol_rate = 1*10^5;       % 符號速率(sym/s)
SNR = -5;                   % 信噪比(dB)
load GoldSequence.mat;
PNcode = 2*GoldSequence(1

ata_len,1:63)-1;    % 偽碼長度：數(shù)據(jù)長度*63
n_chip = 63;

g = 9.8;                      % 重力加速度(m/s^2)
c = 3*10^8;                   % 光速(m/s)
ts = 1/1890;                   % JPL的采樣時間(s),所設(shè)置的值與擴碼序列的長度有關(guān)
V0 = -500;                   % 初始速度(m/s)
Theta0 = 2*pi*randn(1);       % 初始相位(rad)

Time = [3 0.5 2 0.5 2];
Total_length = round(Time/ts);
% Jerk = [zeros(1,3/ts) 100*g*ones(1,0.5/ts) zeros(1,2/ts) -100*g*ones(1,0.5/ts) zeros(1,2/ts)];
Jerk = [zeros(1,Total_length(1)) 100*g*ones(1,Total_length(2)) zeros(1,Total_length(3)) -100*g*ones(1,Total_length(4)) zeros(1,Total_length(5))];
Acceleration = [-25*g cumsum(Jerk(1:end-1)*ts)-25*g];
Velocity = [V0 cumsum(Acceleration(1:end-1)*ts + Jerk(1:end-1)*ts^2/2)+V0];
Freq_Doppler = Freq_carrier / c * Velocity;
Freq_Doppler_der1 = Freq_carrier / c * Acceleration;
Freq_Doppler_der2 = Freq_carrier / c * Jerk;
Theta = Theta0 + 2*pi*cumsum(Freq_Doppler)*ts + 2*pi*cumsum(Freq_Doppler_der1).*ts^2/2 + 2*pi*cumsum(Freq_Doppler_der2).*ts^3/6;

%================== 數(shù)據(jù)序列 ===================%
Data_sequence = randi(2,1,160)-1;

%================== 信道編碼 ===================%
trellis = poly2trellis([5,4],[23 35 0;0 5 13]);
Data_sequence_conv = convenc(Data_sequence,trellis);

%================== BPSK調(diào)制 ===================%
Data_sequence_mod = BPSKMod(Data_sequence_conv);

%================== 擴頻調(diào)制 ===================%
Spread_sequence = Data_sequence_mod(1) * PNcode(1,

;
for i=2:1

ata_len
Spread_sequence = cat(2,Spread_sequence,Data_sequence_mod(i) * PNcode(i,

);
end

%================== 加頻偏和噪聲 ===================%
Spread_sequence_offset = Spread_sequence .* exp(1i*Theta);
Spread_sequence_noise = Noise_Generate(Spread_sequence_offset,SNR);

PNcode1 = reshape(PNcode.',1,[]);
Signal_despread = Spread_sequence_noise.*PNcode1;

w_nF_PLL = 90;
T_IC = 1e-3;
a3 = 1.1;
b3 = 2.4;
Simulation_len = length(Signal_despread);
Signal_PLL = zeros(1,Simulation_len);
NCO_Phase = zeros(1,Simulation_len);
Discriminator_Out = zeros(1,Simulation_len);
phase_third_derivate_temp = 0;
phase_second_derivate1_temp = 0;

for i=2:Simulation_len

Signal_PLL(i) = Signal_despread(i)*exp(-1i*mod(NCO_Phase(i-1),2*pi));
I_PLL(i) = real(Signal_PLL(i));
Q_PLL(i) = imag(Signal_PLL(i));
Discriminator_Out(i)=sign(I_PLL(i))*Q_PLL(i)/abs(Signal_PLL(i));
% Discriminator_Out(i) = atan(Q_PLL(i)/I_PLL(i));

phase_third_derivate = T_IC*w_nF_PLL^3*Discriminator_Out(i);
phase_third_derivate1 = phase_third_derivate + phase_third_derivate_temp;
phase_third_derivate1_temp = phase_third_derivate1 + phase_third_derivate_temp;
phase_second_derivate = a3*w_nF_PLL^2*Discriminator_Out(i) + 1/2*phase_third_derivate1_temp;
phase_second_derivate1 = phase_second_derivate*T_IC + phase_second_derivate1_temp;
phase_second_derivate2 = phase_second_derivate1 + phase_second_derivate1_temp;
phase_derivate = 1/2*phase_second_derivate2 + b3*w_nF_PLL*Discriminator_Out(i);

NCO_Phase(i)=NCO_Phase(i-1) + T_IC*phase_derivate;

phase_third_derivate_temp = phase_third_derivate1;
phase_second_derivate1_temp = phase_second_derivate1;

end

回復(fù)此樓

» 猜你喜歡

332求調(diào)劑已經(jīng)有5人回復(fù)
085701求調(diào)劑初試286分已經(jīng)有3人回復(fù)
材料工程294求調(diào)劑已經(jīng)有4人回復(fù)
本科新能源科學與工程，一志愿華理能動285求調(diào)劑已經(jīng)有10人回復(fù)
070305高分子化學與物理 304分求調(diào)劑已經(jīng)有5人回復(fù)
083000學碩274求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京化工大學材料與化工（085600）296求調(diào)劑已經(jīng)有10人回復(fù)
求調(diào)劑一志愿海大，0703化學學碩304分，有大創(chuàng)項目，四級已過已經(jīng)有11人回復(fù)
071000生物學求調(diào)劑，初試成績343 已經(jīng)有7人回復(fù)
0703化學調(diào)劑，求導(dǎo)師收已經(jīng)有8人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

如何使用matlab仿真信道容量已經(jīng)有6人回復(fù)
MATLAB 仿真已經(jīng)有0人回復(fù)
adams與matlab聯(lián)合仿真已經(jīng)有1人回復(fù)
基于蒙特卡羅仿真用于多種失效分布底事件的故障樹分析如何實現(xiàn)Matlab仿真？已經(jīng)有2人回復(fù)
如何根據(jù)窗口估計法實現(xiàn)matlab的濾波仿真已經(jīng)有0人回復(fù)
MATLAB-SIMULINK通信系統(tǒng)建模與仿真實例分析-邵玉斌-清華大學已經(jīng)有109人回復(fù)
線性調(diào)頻連續(xù)波雷達的matlab仿真實現(xiàn) 已經(jīng)有2人回復(fù)
怎樣用matlab實現(xiàn)細菌覓食算法的仿真已經(jīng)有0人回復(fù)
有關(guān)matlab仿真的問題已經(jīng)有6人回復(fù)
求基于稀疏表示的圖像分類的Matlab仿真實驗代碼已經(jīng)有2人回復(fù)
怎么應(yīng)用matlab蒙特卡洛仿真得以實現(xiàn) 已經(jīng)有1人回復(fù)
【已完結(jié)】matlab仿真圖形（交流問題即有獎，呵呵）已經(jīng)有6人回復(fù)
matlab 電路仿真已經(jīng)有40人回復(fù)
matlAB仿真中l(wèi)ms算法應(yīng)用的一些問題已經(jīng)有0人回復(fù)
用matlab仿真wsn 已經(jīng)有4人回復(fù)
【求助】如何利用MATLAB仿真實現(xiàn)服從某分布的時間序列已經(jīng)有2人回復(fù)

1樓 2015-12-11 10:56:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

蒼翠的路

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 373.8
紅花: 1
帖子: 28
在線: 15小時
蟲號: 3579141
注冊: 2014-12-05
專業(yè): 通信理論與系統(tǒng)

引用回帖:

4樓: Originally posted by aknuaa at 2015-12-14 06:50:21
哪敢說指點，我現(xiàn)在被逼著解決一個satcom下高動態(tài)捕獲跟蹤問題，根本沒有什么資料……嘗試這學習下高階鎖相環(huán)方法……
...

額額，好吧。。。

回復(fù)此樓

5樓2015-12-14 10:01:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個回答

aknuaa

木蟲 (著名寫手)

應(yīng)助: 17 (小學生)
金幣: 3699.6
散金: 485
紅花: 4
帖子: 1138
在線: 211.6小時
蟲號: 1709843
注冊: 2012-03-22
性別: GG
專業(yè): 通信理論與系統(tǒng)

GPS方面的？跟蹤頻偏大嗎？JPL的高動態(tài)模型具體指什么？

贊一下

回復(fù)此樓

wantingdrenched

2樓2015-12-13 12:15:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

蒼翠的路

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 373.8
紅花: 1
帖子: 28
在線: 15小時
蟲號: 3579141
注冊: 2014-12-05
專業(yè): 通信理論與系統(tǒng)

引用回帖:

2樓: Originally posted by aknuaa at 2015-12-13 12:15:48
GPS方面的？跟蹤頻偏大嗎？JPL的高動態(tài)模型具體指什么？

也不算GPS方面的，但很相似，偽碼用的是Gold碼。由于只研究跟蹤方面的，所以假設(shè)經(jīng)過捕獲后頻偏大概1k左右，JPL的高動態(tài)模型是這樣的：速度變化的過程持續(xù)8s，其中在3s-3.5s以及5.5s-6s這兩個時間段內(nèi)產(chǎn)生100g/s的正加加速度和負加加速度，其余時間段則保持25g的常加速狀態(tài)。理論上來講三階PLL應(yīng)該可以跟蹤上的，但現(xiàn)在無法鎖相，程序看了幾次都覺得沒問題，哥們你也是研究這方面的嗎，求指點啊。。

贊一下

回復(fù)此樓

3樓2015-12-13 21:10:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

aknuaa

木蟲 (著名寫手)

應(yīng)助: 17 (小學生)
金幣: 3699.6
散金: 485
紅花: 4
帖子: 1138
在線: 211.6小時
蟲號: 1709843
注冊: 2012-03-22
性別: GG
專業(yè): 通信理論與系統(tǒng)

引用回帖:

3樓: Originally posted by 蒼翠的路 at 2015-12-13 21:10:29
也不算GPS方面的，但很相似，偽碼用的是Gold碼。由于只研究跟蹤方面的，所以假設(shè)經(jīng)過捕獲后頻偏大概1k左右，JPL的高動態(tài)模型是這樣的：速度變化的過程持續(xù)8s，其中在3s-3.5s以及5.5s-6s這兩個時間段內(nèi)產(chǎn)生100g/s的 ...

哪敢說指點，我現(xiàn)在被逼著解決一個satcom下高動態(tài)捕獲跟蹤問題，根本沒有什么資料……嘗試這學習下高階鎖相環(huán)方法……

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

wantingdrenched

4樓2015-12-14 06:50:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 343求調(diào)劑 +4	愛羈絆 2026-03-28	4/200	2026-03-28 17:44 by wxiongid
[考研] 學碩274求調(diào)劑 +8	Li李魚 2026-03-26	8/400	2026-03-28 16:53 by xuxiang
[考研] 346求調(diào)劑一志愿070303有機化學 +3	蘿卜燉青菜 2026-03-28	3/150	2026-03-28 14:11 by 唐沐兒
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[考研] 277跪求調(diào)劑 +5	1915668 2026-03-27	9/450	2026-03-28 09:58 by zhshch
[考研] 291求調(diào)劑 +7	孅華 2026-03-22	7/350	2026-03-28 04:02 by fmesaito
[考研] 265求調(diào)劑11408 +3	劉小鹿lu 2026-03-27	3/150	2026-03-27 20:53 by nihaoar
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚的貓
[論文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陸小果畫大餅
[考研] 286求調(diào)劑 +4	lim0922 2026-03-26	4/200	2026-03-27 10:28 by guoweigw
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[碩博家園] 北京林業(yè)大學碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 294分080500材料科學與工程求調(diào)劑 +4	柳溪邊 2026-03-26	4/200	2026-03-26 21:14 by XPU李慶
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 281求調(diào)劑 +6	Koxui 2026-03-24	7/350	2026-03-26 15:37 by 無際的草原
[考研] 一志愿河工大 081700 276求調(diào)劑 +4	地球繞著太陽轉(zhuǎn) 2026-03-23	4/200	2026-03-26 14:27 by zzll406
[考研] 總分293求調(diào)劑 +6	加一一九 2026-03-25	8/400	2026-03-26 13:30 by yujianx
[考研] 材料與化工304求B區(qū)調(diào)劑 +3	邱gl 2026-03-25	3/150	2026-03-25 19:03 by Ainin_
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch