飛行目標(biāo)的多基雷達探測及攻擊問題數(shù)模論文-資料下載頁
2025-06-27 17:33本頁面
【正文】 ,39。z039。)。 %記錄三點確定的坐標(biāo)值 u(l)=x0(1)。 v(l)=y0(1)。 w(l)=abs(z0(1))。 l=l+1。 end end end end end enddiary(39。39。)。% % 將數(shù)據(jù)寫入excel文件% xlswrite(39。39。,u,39。sheet139。,39。A1:A383639。)。 % A就是待寫的數(shù)據(jù)% xlswrite(39。39。,v,39。sheet139。,39。B1:B383639。)。 % A就是待寫的數(shù)據(jù)% xlswrite(39。39。,w,39。sheet139。,39。C1:C383639。)。 % A就是待寫的數(shù)據(jù)%計算坐標(biāo)的平均值Tips=sprintf(39。目標(biāo)最優(yōu)坐標(biāo)如下::\n39。)Aver_x=mean(u)Aver_y=mean(v)Aver_z=mean(w)%方差計算相關(guān)M=[Aver_x,Aver_y,Aver_z]。M2=[11422,964,14033]。%甲組數(shù)據(jù)殘差最優(yōu)解% M2=[88294,1140,68289]。%乙組數(shù)據(jù)殘差最優(yōu)解% M2=[85190,100860,71060]。%丙組數(shù)據(jù)殘差最優(yōu)解u_s2=zeros(1,counter)。v_s2=zeros(1,counter)。w_s2=zeros(1,counter)。%模型一方差計算for i=1:1:counteru_s2(i)=(u(i)M(1))^2。 v_s2(i)=(v(i)M(2))^2。w_s2(i)=(w(i)M(3))^2。endTips=sprintf(39。模型一的方差如下:\n39。)x_s2=mean(u_s2)y_s2=mean(v_s2)z_s2=mean(w_s2)%模型二方差計算for i=1:1:counteru_s2(i)=(u(i)M2(1))^2。 v_s2(i)=(v(i)M2(2))^2。w_s2(i)=(w(i)M2(3))^2。endTips=sprintf(39。模型二的方差如下:\n39。)x_s2=mean(u_s2)y_s2=mean(v_s2)z_s2=mean(w_s2)diary off。%%繪出給定三個雷達坐標(biāo)的散點圖% i=1。j=7。k=13。% Read_data = xlsread(39。39。, 39。sheet139。)。% x = Read_data(: , 1)。% y = Read_data(: , 2)。% z = Read_data(: , 3)。% r = Read_data(: , 4)。% x_x=[x(i) x(j) x(k)]。% y_y=[y(i) y(j) y(k)]。% plot(x_x,y_y,39。r*39。)附件三：%文件名：%功能: 用于求解各組數(shù)據(jù)的殘差最小值function f=fun_Location(k)% 定義全局變量syms x y z r sum_e = 0。%function f=fun_Coordinate(x)。% sheet1, sheet2, sheet3分別對應(yīng)甲乙丙的三組數(shù)據(jù)%Read_data = xlsread(39。39。,39。sheet139。,39。F3:I3439。)%在同一個sheet中讀取需要的數(shù)據(jù)區(qū)域%飛機甲的數(shù)據(jù):A3:D32。飛機乙的數(shù)據(jù):F3:I34。飛機丙的數(shù)據(jù):K3:N14Read_data = xlsread(39。39。,39。sheet139。)。% 讀取矩陣尺寸[rows, cols] = size(Read_data)。%將讀取到的坐標(biāo)數(shù)據(jù)分別存放到相應(yīng)的數(shù)組中x = Read_data(: , 1)。y = Read_data(: , 2)。z = Read_data(: , 3)。r = Read_data(: , 4)。%計算殘差的平方和for i = 1: 1: rows sum_e = sum_e + r(i)^2+(k(1)x(i))^2+k(2)y(i)^2 + k(3)^2 2*r(i)*sqrt((k(1)x(i))^2+k(2)y(i)^2+k(3)^2)。endf =sum_e。% clear。% clc。% syms x0 y0 z0% Read_data = xlsread(39。39。,39。sheet139。)。% x = Read_data(: , 1)。% y = Read_data(: , 2)。% z = Read_data(: , 3)。% r = Read_data(: , 4)。% % 讀取矩陣尺寸% [rows, cols] = size(Read_data)。% % 計算x0, y0, z0初始值% eq1 = (x0 x(1)) ^ 2 + (y0 y(1)) ^ 2 + z0 ^ 2 r(1) ^ 2。% eq2 = (x0 x(rows/2)) ^ 2 + (y0 y(rows/2)) ^ 2 + z0 ^ 2 r(rows/2) ^ 2。% eq3 = (x0 x(rows)) ^ 2 + (y0 y(rows)) ^ 2 + z0 ^ 2 r(rows) ^ 2。% [x0, y0, z0] = solve(eq1, eq2, eq3)。% k0=[x0(1),y0(1),abs(z0(1))]。% k0=eval(k0)。% options=optimset(39。LargeSCale39。,39。on39。)。% diary(39。39。)。% [x,fval,exitflag,output]=fminunc(39。fun_Location39。, k0,options)% diary off。附件四：%文件名：%功能: 用于求解導(dǎo)彈擊中目標(biāo)的飛行軌跡clear。clc。Xa=11771。%x甲Ya=。%y甲Xb=100000。%x乙Yb=0。%y乙k=。diary(39。39。)。w0=(YaYb)/(XaXb)。c1=(w0+sqrt(1+w0^2))/((XbXa)^(k))c2=Ya((XbXa)^(1+k))/(2*c1*(1+k))+(c1*((XbXa)^(1k)))diary off。x=Xa:1:Xb。y=((Xbx).^(1+k))/(2*c1*(1+k))(c1*((Xbx).^(1k)))/(2*(1k))+c2。plot(x,y,39。r39。)[Y,I]=max(y)。% text(x(I),Y,num2str(Y))hold on。plot(x(I),Y,39。mpentagram39。)text(x(I)+800,Y,39。擊中位置39。)hold on。plot(Xa,Ya,39。 39。)text(Xa,Ya+2000,39。導(dǎo)彈發(fā)射點39。)hold on。a=0:Xb。b=Xb*ones(1,Xb+1)。line(b,a)axis([1,120000 0 65000])。 31
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1