1.問題描述：

陸基制導(dǎo)系統(tǒng)地面站布設(shè)策略

2.部分程序：

?

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

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

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

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

? ? gen ? = 0;

? ? %計算對應(yīng)的目標(biāo)值
? ? 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,:);
? ? ? ? ? ? ? %計算對應(yīng)的目標(biāo)值
? ? ? ? ? ? ? [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;?
? ? ? ? ? %保存參數(shù)收斂過程和誤差收斂過程以及函數(shù)值擬合結(jié)論
? ? ? ? ? Error(gen) = mean(JJ);
? ? end?

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

? ? %根據(jù)優(yōu)化結(jié)果輸出部署結(jié)果
? ? 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.仿真結(jié)論：

A06-63

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。