 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1188 | 回復(fù): 5 | ||
[求助]
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ù)長度*63n_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_lenSpread_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ù)此樓» 猜你喜歡
280分求調(diào)劑 一志愿085802
已經(jīng)有6人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有3人回復(fù)
-
263求調(diào)劑
已經(jīng)有4人回復(fù)
-
石河子大學（211、雙一流）碩博研究生長期招生公告
已經(jīng)有3人回復(fù)
-
一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有13人回復(fù)
-
北科281學碩材料求調(diào)劑
已經(jīng)有6人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有10人回復(fù)
-
298求調(diào)劑
已經(jīng)有6人回復(fù)
-
310求調(diào)劑
已經(jīng)有4人回復(fù)
-
311求調(diào)劑
已經(jīng)有6人回復(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ù)
aknuaa
木蟲 (著名寫手)
- 應(yīng)助: 17 (小學生)
- 金幣: 3699.6
- 散金: 485
- 紅花: 4
- 帖子: 1138
- 在線: 211.6小時
- 蟲號: 1709843
- 注冊: 2012-03-22
- 性別: GG
- 專業(yè): 通信理論與系統(tǒng)
2樓2015-12-13 12:15:48
3樓2015-12-13 21:10:29
aknuaa
木蟲 (著名寫手)
- 應(yīng)助: 17 (小學生)
- 金幣: 3699.6
- 散金: 485
- 紅花: 4
- 帖子: 1138
- 在線: 211.6小時
- 蟲號: 1709843
- 注冊: 2012-03-22
- 性別: GG
- 專業(yè): 通信理論與系統(tǒng)
|
哪敢說指點，我現(xiàn)在被逼著解決一個satcom下高動態(tài)捕獲跟蹤問題，根本沒有什么資料……嘗試這學習下高階鎖相環(huán)方法…… 發(fā)自小木蟲IOS客戶端 |
4樓2015-12-14 06:50:21
5樓2015-12-14 10:01:11
hwz1002
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2768.4
- 紅花: 2
- 帖子: 160
- 在線: 40.7小時
- 蟲號: 2274025
- 注冊: 2013-02-04
- 性別: GG
- 專業(yè): 信息處理方法與技術(shù)
6樓2015-12-26 21:09:16
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 蒼翠的路 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 310求調(diào)劑 +4 | baibai1314 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研]
 一志愿東華大學化學070300，求調(diào)劑
+6
|
2117205181 2026-03-21 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | 來好運來來來 2026-03-21 | 4/200 |
|
|
|
[考研] 一志愿華中科技大學071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研] 268求調(diào)劑 +9 | 簡單點0 2026-03-17 | 9/450 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調(diào)劑） +3 | 123阿標 2026-03-17 | 3/150 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|