 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1189 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
matlab仿真實(shí)現(xiàn)三階鎖相環(huán)路跟蹤高動(dòng)態(tài)多普勒變化出現(xiàn)問(wèn)題，急需指點(diǎn)[有源碼] 已有1人參與
|
|||
|
RT，最近在忙論文，但仿真出現(xiàn)問(wèn)題，結(jié)果總是不對(duì)。用的模型是JPL的高動(dòng)態(tài)模型，環(huán)路帶寬根據(jù)跟蹤門(mén)限設(shè)為90Hz，積分時(shí)間為1ms，但是好像就是一直沒(méi)法鎖相，鑒相誤差一直較大，有沒(méi)做過(guò)這方面的同學(xué)或前輩啊，急需指點(diǎn)和討論��！ clc; clear all; Chip_rate = 6.3*10^6; % 擴(kuò)頻碼速率(chips/s) Spreadcode_len = 63; % 擴(kuò)頻碼長(zhǎng)度(chips) Data_len = 240; % 數(shù)據(jù)長(zhǎng)度(symbols) Freq_carrier = 4*10^10; % 射頻載波頻率(Hz) Freq_sample = 6.3*10^6; % 接收機(jī)采樣頻率(Hz) Symbol_rate = 1*10^5; % 符號(hào)速率(sym/s) SNR = -5; % 信噪比(dB) load GoldSequence.mat; PNcode = 2*GoldSequence(1  ata_len,1:63)-1; % 偽碼長(zhǎng)度：數(shù)據(jù)長(zhǎng)度*63n_chip = 63; g = 9.8; % 重力加速度(m/s^2) c = 3*10^8; % 光速(m/s) ts = 1/1890; % JPL的采樣時(shí)間(s),所設(shè)置的值與擴(kuò)碼序列的長(zhǎng)度有關(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); %================== 擴(kuò)頻調(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ù)此樓» 猜你喜歡
323求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有4人回復(fù)
-
352求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有8人回復(fù)
-
277材料科學(xué)與工程080500求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有18人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
280分求調(diào)劑 一志愿085802
已經(jīng)有7人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有3人回復(fù)
-
263求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 如何使用matlab仿真信道容量
已經(jīng)有6人回復(fù)
- MATLAB 仿真
已經(jīng)有0人回復(fù)
- adams與matlab聯(lián)合仿真
已經(jīng)有1人回復(fù)
- 基于蒙特卡羅仿真用于多種失效分布底事件的故障樹(shù)分析如何實(shí)現(xiàn)Matlab仿真？
已經(jīng)有2人回復(fù)
- 如何根據(jù)窗口估計(jì)法實(shí)現(xiàn)matlab的濾波仿真
已經(jīng)有0人回復(fù)
- MATLAB-SIMULINK通信系統(tǒng)建模與仿真實(shí)例分析-邵玉斌-清華大學(xué)
已經(jīng)有109人回復(fù)
- 線性調(diào)頻連續(xù)波雷達(dá)的matlab仿真實(shí)現(xiàn)
已經(jīng)有2人回復(fù)
- 怎樣用matlab實(shí)現(xiàn)細(xì)菌覓食算法的仿真
已經(jīng)有0人回復(fù)
- 有關(guān)matlab仿真的問(wèn)題
已經(jīng)有6人回復(fù)
- 求基于稀疏表示的圖像分類的Matlab仿真實(shí)驗(yàn)代碼
已經(jīng)有2人回復(fù)
- 怎么應(yīng)用matlab蒙特卡洛仿真得以實(shí)現(xiàn)
已經(jīng)有1人回復(fù)
- 【已完結(jié)】matlab仿真圖形（交流問(wèn)題即有獎(jiǎng)，呵呵）
已經(jīng)有6人回復(fù)
- matlab 電路仿真
已經(jīng)有40人回復(fù)
- matlAB仿真中l(wèi)ms算法應(yīng)用的一些問(wèn)題
已經(jīng)有0人回復(fù)
- 用matlab仿真wsn
已經(jīng)有4人回復(fù)
- 【求助】如何利用MATLAB仿真實(shí)現(xiàn)服從某分布的時(shí)間序列
已經(jīng)有2人回復(fù)
aknuaa
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 3699.6
- 散金: 485
- 紅花: 4
- 帖子: 1138
- 在線: 211.6小時(shí)
- 蟲(chóng)號(hào): 1709843
- 注冊(cè): 2012-03-22
- 性別: GG
- 專業(yè): 通信理論與系統(tǒng)
|
哪敢說(shuō)指點(diǎn)，我現(xiàn)在被逼著解決一個(gè)satcom下高動(dòng)態(tài)捕獲跟蹤問(wèn)題，根本沒(méi)有什么資料……嘗試這學(xué)習(xí)下高階鎖相環(huán)方法…… 發(fā)自小木蟲(chóng)IOS客戶端 |
4樓2015-12-14 06:50:21
aknuaa
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 3699.6
- 散金: 485
- 紅花: 4
- 帖子: 1138
- 在線: 211.6小時(shí)
- 蟲(chóng)號(hào): 1709843
- 注冊(cè): 2012-03-22
- 性別: GG
- 專業(yè): 通信理論與系統(tǒng)
2樓2015-12-13 12:15:48
|
也不算GPS方面的，但很相似，偽碼用的是Gold碼。由于只研究跟蹤方面的，所以假設(shè)經(jīng)過(guò)捕獲后頻偏大概1k左右，JPL的高動(dòng)態(tài)模型是這樣的：速度變化的過(guò)程持續(xù)8s，其中在3s-3.5s以及5.5s-6s這兩個(gè)時(shí)間段內(nèi)產(chǎn)生100g/s的正加加速度和負(fù)加加速度，其余時(shí)間段則保持25g的常加速狀態(tài)。理論上來(lái)講三階PLL應(yīng)該可以跟蹤上的，但現(xiàn)在無(wú)法鎖相，程序看了幾次都覺(jué)得沒(méi)問(wèn)題，哥們你也是研究這方面的嗎，求指點(diǎn)啊。。 |
3樓2015-12-13 21:10:29
5樓2015-12-14 10:01:11
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 招08考數(shù)學(xué) +4 | laoshidan 2026-03-20 | 8/400 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
|
[考研] 求助 +5 | 夢(mèng)里的無(wú)言 2026-03-21 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+8
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|