飛行目標(biāo)的多基雷達(dá)探測(cè)及攻擊問(wèn)題數(shù)模論文-wenkub.com
2025-06-24 17:33 本頁(yè)面
【正文】 line(b,a)axis([1,120000 0 65000])。導(dǎo)彈發(fā)射點(diǎn)39。)hold on。plot(x(I),Y,39。plot(x,y,39。w0=(YaYb)/(XaXb)。%y乙k=。Xa=11771。fun_Location39。% diary(39。LargeSCale39。% [x0, y0, z0] = solve(eq1, eq2, eq3)。% % 讀取矩陣尺寸% [rows, cols] = size(Read_data)。% x = Read_data(: , 1)。39。endf =sum_e。y = Read_data(: , 2)。sheet139。飛機(jī)乙的數(shù)據(jù):F3:I34。sheet139。%function f=fun_Coordinate(x)。% y_y=[y(i) y(j) y(k)]。% y = Read_data(: , 2)。, 39。j=7。endTips=sprintf(39。模型一的方差如下:\n39。%模型一方差計(jì)算for i=1:1:counteru_s2(i)=(u(i)M(1))^2。%乙組數(shù)據(jù)殘差最優(yōu)解% M2=[85190,100860,71060]。目標(biāo)最優(yōu)坐標(biāo)如下::\n39。,39。 % A就是待寫(xiě)的數(shù)據(jù)% xlswrite(39。sheet139。)。,u,39。39。 v(l)=y0(1)。,39。 [x0,y0,z0]=solve(eq1, eq2, eq3,39。(y(j)==y(k))))) if(((y(i)y(j))*(x(j)x(k)))~=((y(j)y(k))*(x(i)x(j)))) %求解三點(diǎn)確定的坐標(biāo)值 syms x0 y0 z0 。(x(j)==x(k)))||((y(i)==y(j))amp。if(~flag) for i=1:1:(rows2) for j=i+1:1:(rows1) for k=j+1:1:rows if(~(((x(i)==x(j))amp。 end end end end end counter l=1:1:counter。(y(j)==y(k))))) if(((y(i)y(j))*(x(j)x(k)))~=((y(j)y(k))*(x(i)x(j)))) %打印出滿足判定條件的坐標(biāo)下標(biāo)% sprintf(39。(x(j)==x(k)))||((y(i)==y(j))amp。 [x0, y0, z0] = solve(eq1, eq2, eq3)end%flag為0時(shí)進(jìn)入條件判斷求解飛機(jī)坐標(biāo)的平均值(計(jì)算坐標(biāo)均值)if(~flag) for i=1:1:(rows2) for j=i+1:1:(rows1) for k=j+1:1:rows if(~(((x(i)==x(j))amp。 k=7。r = Read_data(: , 4)。)。%讀取表格中的雷達(dá)坐標(biāo)數(shù)據(jù)Read_data = xlsread(39。附件二：%文件名：%功能: 用于平均值法求解飛機(jī)坐標(biāo)的最優(yōu)值% 根據(jù)求解的結(jié)果計(jì)算方差clear。乙組39。g+39。y = Read_data(: , 2)。, 39。b*39。y = Read_data(: , 2)。, 39。ro39。y = Read_data(: , 2)。, 39。 模型的提高首先，為了使定位更加準(zhǔn)備，在時(shí)間允許的條件下，有必要對(duì)定位精度對(duì)雷達(dá)各坐標(biāo)及距離誤差的變化的靈敏度進(jìn)行分析；另外還可以通過(guò)求模型一、模型二的估計(jì)結(jié)果的均方誤差，來(lái)判斷哪個(gè)模型更優(yōu)。六、模型優(yōu)點(diǎn)與提高 模型的優(yōu)點(diǎn)（1）問(wèn)題1中確定定位飛行物的最少雷達(dá)數(shù)目時(shí)，我們通過(guò)討論的方式論證了1個(gè)，2個(gè)及在同一條直線上的3個(gè)雷達(dá)是不滿足要求的；另外在分析距離誤差和雷達(dá)坐標(biāo)誤差與定位精度的關(guān)系時(shí)，使用了控制變量法，并且使用一階泰勒展開(kāi)式進(jìn)行降次處理。于是，把三維空間的求解問(wèn)題轉(zhuǎn)化到二維平面中。 導(dǎo)彈的追蹤敵機(jī)的軌跡從圖中可以得到敵機(jī)被擊毀的坐標(biāo)為，假設(shè)導(dǎo)彈發(fā)射時(shí)刻時(shí)間，則擊毀的時(shí)刻。 2）若型導(dǎo)彈從任意點(diǎn)出發(fā)，要保證該導(dǎo)彈一定在敵機(jī)到達(dá)安全區(qū)的途中命中敵機(jī)，則導(dǎo)彈的速度必須不小于敵機(jī)速度的倍。（1）當(dāng)以該導(dǎo)彈基地上空（即）為空對(duì)空導(dǎo)彈發(fā)射點(diǎn)且以敵機(jī)的作為坐標(biāo)的原點(diǎn)（即），那么微分方程的定解條件為：，則可解得： ()該方程滿足臨界條件：：解得：，由于，為大于0的數(shù)，因此?。? 由上式可得以下結(jié)論： 1）若型導(dǎo)彈從基地出發(fā)，并在臨界點(diǎn)擊中敵機(jī)，則導(dǎo)彈的速度需要為敵機(jī)速度的倍。因題中給出當(dāng)敵機(jī)向正北方向飛行公里后就到達(dá)安全區(qū)，由于敵方電子干擾的作用,敵機(jī)一旦進(jìn)入安全區(qū)后,導(dǎo)彈將失去追蹤目標(biāo)，無(wú)法將其擊毀。，陣列雷達(dá)利用數(shù)量較多的雷達(dá)來(lái)捕捉信號(hào)較弱的目標(biāo)，達(dá)到及時(shí)預(yù)警的功能。模型一估算結(jié)果的方差小于模型二,這是由于在正態(tài)分布的誤差下，樣本均值本身與CRLB重合[2]，其本身就是最小方差無(wú)偏估計(jì)；然后因?yàn)閷?shí)際情況中，定位誤差不可能?chē)?yán)格服從正態(tài)分布，因此估算出來(lái)的小部分結(jié)果的方差要略大于模型二結(jié)果的方差。開(kāi)始進(jìn)行迭代運(yùn)算，選擇保留最小的殘差平方和，其具體求解過(guò)程參見(jiàn)附錄1?？紤]到第三組所給的數(shù)據(jù)量相對(duì)第一組和第二組而言較少所以我們先計(jì)算第三組數(shù)據(jù)。 模型二: 基于最小殘差平方和的最優(yōu)定位模型（1）模型建立過(guò)程假設(shè)組網(wǎng)的個(gè)雷達(dá)之間的距離不大，可以忽略地球表面是曲面的影響，認(rèn)為它們同在一個(gè)平面上，第部雷達(dá)的坐標(biāo)為，目標(biāo)的坐標(biāo)為，是第部雷達(dá)的測(cè)得的目標(biāo)距離，即有 （）這樣就得到類(lèi)似于上式的個(gè)方程，如果每個(gè)雷達(dá)測(cè)得距離是精確的，那么任意選擇 3 個(gè)方程組成方程組可解出目標(biāo)的精確位置。 三組雷達(dá)坐標(biāo)散點(diǎn)圖由圖可知，定位雷達(dá)呈現(xiàn)陣列分布。那么，任何無(wú)偏估計(jì)量的方差必定滿足 ()并且，對(duì)于某個(gè)函數(shù)g和I，當(dāng)且僅當(dāng) ()時(shí)，對(duì)所有的下限無(wú)偏估計(jì)就可以求得。因而我們一般要求三個(gè)雷達(dá)在其所確定的圓周上分布均勻且這三個(gè)雷達(dá)所確定的圓周半徑不能太小（否則也可能會(huì)產(chǎn)生較大的誤差），這樣飛行物誤差會(huì)比較小。（3）影響結(jié)果比較通過(guò)上面的計(jì)算我們知道不論是距離誤差還是坐標(biāo)誤差，它們對(duì)空中目標(biāo)定位的影響都與系數(shù)矩陣,有密切的聯(lián)系（）。進(jìn)一步推導(dǎo)得 ()注意是飛行物的實(shí)際位置向量,而是前述方法測(cè)量解算獲得的飛行物三維位置向量。（2）雷達(dá)自身坐標(biāo)誤差對(duì)定位精度的影響模型設(shè)地面三個(gè)站點(diǎn)橫坐標(biāo)真實(shí)值為、坐標(biāo)誤差分別為、令，觀測(cè)值為 則有： ()同理可得 () ()令, 分別代表三維坐標(biāo)觀測(cè)值向量矩陣及三維坐標(biāo)值
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1