gps頻域捕獲實(shí)驗(yàn)報(bào)告(編輯修改稿)
2024-09-10 21:16 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 1)*(fc+dopplor)/fs)產(chǎn)生的，其中Len為gpsdata的長(zhǎng)度，fc為信號(hào)中心頻率，dopplor為多普勒頻移，剛開(kāi)始設(shè)置為4000，增量為200。，fs為信號(hào)采樣頻率。然后在加一部信號(hào)的強(qiáng)化： for k=1:floor(length(gpsdata_ddc)/Q) gpsdata_ddc_CIC(k) = sum(gpsdata_ddc((k1)*Q+1:k*Q))。 end接著利用如下程序分出正弦和余弦波： sine=real(gpsdata_ddc_CIC)。 %虛部 cosine=imag(gpsdata_ddc_CIC)。 %實(shí)部 積分和平方串行搜索捕獲方法的最后一部分是與正弦和余弦信號(hào)分別相乘后進(jìn)行積分和平方。平方主要是為了獲得信號(hào)功率。積分則是將對(duì)應(yīng)于處理數(shù)據(jù)長(zhǎng)度的所有38192個(gè)點(diǎn)的值相加。平方在積分過(guò)程之后進(jìn)行。最后一步是將I支路和Q支路的兩個(gè)平方值相加，然后在開(kāi)方。如果本地偽碼與接受信號(hào)偽碼完全對(duì)齊，同時(shí)本地產(chǎn)生載波頻率與接收信號(hào)頻率相匹配，則得到的輸出值比不滿足這些條件的結(jié)果要大得多。 for q=1:1023 for k=1:1023 %積分 Comp1(q)=Comp1(q)+LoCAcode(q1+k)*sine(k)。 Comp2(q)=Comp2(q)+LoCAcode(q1+k)*cosine(k)。 end Comp(q)=sqrt(Comp1(q)*Comp1(q)+Comp2(q)* Comp2(q))。 %平方和，開(kāi)方 end首先程序?qū)⒈镜貍未a序列與載波信號(hào)逐位進(jìn)行相乘并相加，即卷積。然后本地的偽碼序列從第2位開(kāi)始與載波的第一位相乘，以此類推，這就相當(dāng)于本地偽碼序列循環(huán)左移了一位。循環(huán)卷積1023次后就可算出1023個(gè)相關(guān)值。然后令多普勒頻移加200，再進(jìn)行循環(huán)卷積，直到多普勒頻移從4000到3800為止。這樣就可以算出40*1023的相關(guān)性矩陣。利用mesh函數(shù)畫(huà)出圖像即可獲得圖像，見(jiàn)圖3_3和圖3_4。 圖3_3 1號(hào)衛(wèi)星的圖像 圖3_4 9號(hào)衛(wèi)星的圖像從上面兩張圖可以看出，有些衛(wèi)星可以捕獲的信號(hào)，有些衛(wèi)星不能捕獲到信號(hào)。因此，我們還要設(shè)定一些條件來(lái)判斷是否捕獲到了衛(wèi)星信號(hào)。 圖像的刪選由已給程序知，當(dāng)峰值為平均值的4倍的時(shí)候，我們認(rèn)為GPS捕捉到了衛(wèi)星，那么由圖像知，并不是所有衛(wèi)星都被捕捉到了，那么首先我們要算出最大值和平均值。我們用了如下語(yǔ)句來(lái)計(jì)算最大值和平均值：avg= sum(yy(:))/numel(yy)。 [a,l]=max(yy(:))。[
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號(hào)-1