1.問題描述：

陸基制導系統(tǒng)地面站布設策略

2.部分程序：

?

clc;
clear;
close all;
warning off;
addpath 'func\'
addpath 'GA_toolbox\'
global SCALE
SCALE = 10e3;
global STEPS
STEPS = 1e3;
?
?

for tt=1:1%為了仿真速度快點，這里只顯示某一個時刻的部署結果
? ? %輸出每一時刻的無人機位置 ? ?
? ? if mod(tt,1000) == 1%每1000s更新一次部署，這個1000可以自己修改
? ? ? ?target_pos= [1e3, 10e3, 0]';
? ? end

? ? N ? ? ?= 3;%機動站個數，自己可以調整，調整后重新運行即可。?
? ? MAXGEN = 100;
? ? NIND ? = 50;
? ? Nums ? = 3*N; %
? ? Chrom ?= crtbp(NIND,Nums*10);

? ? %sh
? ? Areas = [];
? ? %設置地圖的最大范圍
? ? for i = 1:2*N
? ? ? ? Areas = [Areas,[-SCALE;SCALE]];%雷達最大部署范圍距離
? ? end
? ? for i = 1:N
? ? ? ? Areas = [Areas,[0;100]];%雷達最大部署范圍距離
? ? end
? ??
? ? FieldD = [rep([10],[1,Nums]);Areas;rep([0;0;0;0],[1,Nums])];

? ? gen ? = 0;

? ? %計算對應的目標值
? ? tmps ? ? ? ? = repmat([400],[1,Nums]);
? ? X ? ? ? ? ? ?= tmps;%初始值
? ? [fobj,Vgdop] = func_obj(X,target_pos);
? ? E ? ? ? ? ? ?= fobj;
? ? Js ? ? ? ? ? = E*ones(NIND,1);
? ? Objv ? ? ? ? = (Js+eps);
? ? gen ? ? ? ? ?= 0;?

? ? while gen < MAXGEN; ??
? ? ? ? ? gen
? ? ? ? ? rng(gen);
? ? ? ? ? Pe0 = 0.995;
? ? ? ? ? pe1 = 0.005;?

? ? ? ? ? FitnV=ranking(Objv); ? ?
? ? ? ? ? Selch=select('sus',Chrom,FitnV); ? ?
? ? ? ? ? Selch=recombin('xovsp', Selch,Pe0); ??
? ? ? ? ? Selch=mut( Selch,pe1); ??
? ? ? ? ? phen1=bs2rv(Selch,FieldD); ??

? ? ? ? ? for a=1:1:NIND ?
? ? ? ? ? ? ? X1 ? ? ? ? ?= phen1(a,:);
? ? ? ? ? ? ? %計算對應的目標值
? ? ? ? ? ? ? [fobj,Vgdop]= func_obj(X1,target_pos);
? ? ? ? ? ? ? JJ(a,1) ? ? = fobj;
? ? ? ? ? ? ? XYp{a} ? ? ?= X1;
? ? ? ? ? ? ? Vp{a} ? ? ? = Vgdop;
? ? ? ? ? end?
? ? ? ? ? Objvsel=(JJ); ? ?
? ? ? ? ? [Chrom,Objv]=reins(Chrom,Selch,1,1,Objv,Objvsel); ??
? ? ? ? ? gen=gen+1;?
? ? ? ? ? %保存參數收斂過程和誤差收斂過程以及函數值擬合結論
? ? ? ? ? Error(gen) = mean(JJ);
? ? end?

? ? figure;
? ? plot(Error,'linewidth',2);
? ? grid on
? ? xlabel('迭代次數');
? ? ylabel('遺傳算法優(yōu)化過程');

? ? %根據優(yōu)化結果輸出部署結果
? ? figure; ? ?
? ? [V,I] = min(JJ);
? ? XY ? ?= XYp{I};
? ? Gdop ?= Vp{a};
? ? X ? ? = XY(1:length(XY)/3);
? ? Y ? ? = XY(1+length(XY)/3:2*length(XY)/3);
? ? Z ? ? = XY(1+2*length(XY)/3:end);
? ? JZ_pos= [X',Y',Z'];
? ? x ? ? =-2*SCALE:STEPS:2*SCALE;
? ? y ? ? =-2*SCALE:STEPS:2*SCALE;
? ? xx ? ?= [0:0.5:20];
? ? [c,handle]=contour(x/1e3,y/1e3,Gdop,xx);
? ? clabel(c,handle);
? ? hold on;
? ? plot(JZ_pos(:,1)/1e3,JZ_pos(:,2)/1e3,'r.');
? ? xlabel('x(km)')
? ? ylabel('y(km)')
? ? title('GDOP圖')

end

3.仿真結論：

A06-63

總結

以上是生活随笔為你收集整理的陆基制导系统地面站布设策略的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。