| 6 | 1/1 | 返回列表 |
| 查看: 1187 | 回復(fù): 5 | ||
[求助]
matlab仿真實(shí)現(xiàn)三階鎖相環(huán)路跟蹤高動態(tài)多普勒變化出現(xiàn)問題,急需指點(diǎn)[有源碼] 已有1人參與
|
|
RT,最近在忙論文,但仿真出現(xiàn)問題,結(jié)果總是不對。用的模型是JPL的高動態(tài)模型,環(huán)路帶寬根據(jù)跟蹤門限設(shè)為90Hz,積分時(shí)間為1ms,但是好像就是一直沒法鎖相,鑒相誤差一直較大,有沒做過這方面的同學(xué)或前輩啊,急需指點(diǎn)和討論。 clc; clear all; Chip_rate = 6.3*10^6; % 擴(kuò)頻碼速率(chips/s) Spreadcode_len = 63; % 擴(kuò)頻碼長度(chips) Data_len = 240; % 數(shù)據(jù)長度(symbols) Freq_carrier = 4*10^10; % 射頻載波頻率(Hz) Freq_sample = 6.3*10^6; % 接收機(jī)采樣頻率(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的采樣時(shí)間(s),所設(shè)置的值與擴(kuò)碼序列的長度有關(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 |
木蟲 (著名寫手)

木蟲 (著名寫手)
|
哪敢說指點(diǎn),我現(xiàn)在被逼著解決一個(gè)satcom下高動態(tài)捕獲跟蹤問題,根本沒有什么資料……嘗試這學(xué)習(xí)下高階鎖相環(huán)方法…… 發(fā)自小木蟲IOS客戶端 |

木蟲 (小有名氣)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +6 | 皮皮皮123456 2026-03-22 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +5 | 上岸6666@ 2026-03-20 | 5/250 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 085600材料與化工 +8 | 安全上岸! 2026-03-16 | 8/400 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|