《通信天线原理Matlab实验.docx》由会员分享，可在线阅读，更多相关《通信天线原理Matlab实验.docx（12页珍藏版）》请在课桌文档上搜索。
1、第三题、稀疏直线阵（选做）在N=20元均匀直线阵（边射状态，d=O.5）的基础上，随机去掉10个阵元（即稀疏率50%,首尾两端阵元保留以保持与原均匀直线阵相同的最大阵列长度）（1）给出一种随机稀疏情形下的10元稀疏阵阵因子，并将其与原20元均匀直线阵阵因子相比较，总结主瓣宽度和副瓣电平的变化规律。（2）以峰值副瓣电平最小化为优化目标，自己尝试提出方法来优化10元稀疏阵的阵元位置（首尾固定，实际上只需优化8个阵元位置），画出优化后的10元稀疏阵阵因子（给出其对应的峰值旁瓣电平），将其与（1）中的10元随机稀疏阵阵因子进行比较。阵元数N=20时阵因子直角坐标函数图Figl2.阵因子归一化极坐标图F
2、ig13.下阵因子归一化直角坐标图Figl4.阵因子归一化极坐标图Fig15.下阵因子归一化直角坐标图Fig16.下阵因子归一化宜角坐标图O5101520253035404550迭代次数Fig17.下阵因子归一化直角坐标图Fig18.下阵因子归一化直角坐标图cpe然喉 跑102030405060708090100迭代次数Fig19.下阵因子归一化宜角坐标图indivlength=30;popsize=30;globalatennum;atennum=20;eranum=50;pc=0.8;pm=0.05;%3fi00,6D;%06fiy%iy6uuIfiy%2ae,AAfiindivlengt
3、h=30;popsize=50;globalatennum;atennum=20;eranum=100;pc=0.8;pm=0.05;%jtaz30%3fi00z6Dj%06a1fiy%iz6iiuty%2aejAAfi%aijAAficlearall;clc;N=input(,CefiafieOoOay,);lamuda=input(,Qefiafie23,);dis=0.2*lamuda;k=2*piIamuda;alfa=O;theta=-pi:O.OOl:pi;w=alfa+k*dis.*cos(theta);Fl=sin(N*w./2);F2=sin(w./2);AF=abs(F1.
4、F2)/N;%AF=abs(sin(N*(w./2)./sin(w./2)/N;%Uiedeta=3*piN;detal=acos(deta(k*dis)%EaAaaOoeciIingdian=2*pi/N;%Eajoefuban=3*piN;SLL=10*logl0(abs(1N)*(sin(N*fuban2)/sin(fuban/2)2)%0e%0e1!Afiz0ecifigure(1);polar(theta,AF);title(O60afiyN=num2str(N)fiOo60e-fiyt,);figure(2);plot(theta*360/(2*pi),N.*AF);title(0
5、o0afiyN=,num2str(N)OoO6C0e-fiyt,);%oe2oe1Afii116%synsx;%bbbb=solve(sin(N*x2)/sin(x2)=sqrt(2)*N2);%banbo=90-acos(real(bbbb(2,1)abs(bbbb(2,1)*360clearall;clc;N=20;%N=input(efiaEeOoOaRy,);Iamuda=Input(,feEaEe23Q);dis=O.5*Iamuda;k=2*pilamuda;alfa=0;a=2:19;K=randperm(1ength(a);count=8;b=a(K(l:count);NSUM
6、=1b20;theta=-pi:0.OOl:piw=alfa+k*dis.*cos(theta);ff=l;fortt=2:10ff=ff+exp(i*NSUM(l,tt).*w);end%AF=abs(sin(N*(w./2).sin(w./2)N;subplot(1,2,1)polar(theta,ff);title(l06afiyN=num2str(N),+O06array(i,3)&m=0array(i,2)=array(i,3);m=l;endarray(i,3)=a;%ifaarray(i,2)array(i,2)=a;endendend%disp(array);fori=l:po
7、psizefitvalue(i,1)=array(i,1)/array(i,2)；%适应度值等于主瓣最大值除以副瓣最大值；endbest,m%保留每代的最优适应值和最优染色体functionbestindivdual,bestfit=best(pop,fitvalue)Maxvalue,indcx=max(fitvalue)将当前种群中最优解保存在MaxValuebestfit=Maxvalue;bestindivdual=pop(index,:);select,m选择%选择更制functionselection=select(fitvalue,pop)m,n=size(pop);Select
8、ion=Zeros(in-1,n);Minvalue,indexl=min(fitvalue);%将当前种群中最次解保存在MinValueMaxvalue,indcx2=max(fitvalue);%将当前种群中最优解保存在MaxValuebestchrom=pop(index2,:)；%最佳染色体保存pop(indexl,:)=bestchrom;index=1:m;index(indexl)=O;index=nonzcros(index);%使最优染色体坐在的位置为零，便于清空，使最优的染色体的原来位置不参与本次的选择、交叉、变异ncwpop=pop(indcx,参与选择交叉，变异的种群f
9、itvalue=fitvalue(indcx,参与选择、交叉、变异的适应度函数值%evo_pOPSiZe=mT;%参与选择、交叉、变异的种群大小fitvalue=cumsum(fitvalue)/sum(fitvalue);rndpos=sort(rand(m-l,1);%产生随即概率fori=l:m-lwhilerndpos(i)fitvalue(j)j=jl;endselection(i,:)=pop(j,:);endcross,m交叉%交叉functionnewpop=cross(selection,pc)m,n=size(selection);newpop=ones(m,n);s=ra
10、nd(l,m);%随机产生与交叉概率相比较的概率yl=find(s=pc);%不参与交叉位lenl=length(yl);iflenl2&mod(lenl,2)=1%如果用来进行交叉的染色体的条数为奇数，将其调整为偶数y2(length(y2)+l)=yl(Ienl);Ienl=Ienl-I;end%lenliflenl=2forp=l:lenl/2crosspos=round(n-2)*rand()+1;a=selection(yl(2*p-l),crosspos:n);b=selection(yl(2*p),crosspos:n);while(length(find(a=l)=length
11、(find(b=l)crosspos=round(n-2)*rand()+1;a=selection(yl(2*p-l),crosspos:n);b=selection(yl(2*p),crosspos:n);end%crossposnewpop(yl(2*pl),:)=selection(yl(2*p-l),Iicrosspos-I),selection(yl(2*p),crosspos:n);newpop(yl(2*p),:)=selection(yl(2*p),Iicrosspos-I),selection(yl(2*p-l),crosspos：n);endendnewpop(y2,:)
12、=Selection(y2,:);mutation,m变异%变异functionnewpop=mutation(newpop,pm)m,n=size(newpop);s=rand(l,m);%随机产生与变异概率相比较的概率fori=l:mifs(i)=pmmutapos=round(n-3)*rand()+2;%随即产生变异位置ifnewpop(i,mutapos)=0a=find(newpop(i,:)=1);a0=round(lcngth(a)-3)*rand()+2;%a(a0)随即找一个为1的位置newpop(i,mutapos)=l;newpop(i,a(a)=0;elseb=fin
13、d(newpop(i,:)=0);b0=round(length(b)-1)*rand()+1;%b(bO)随即找一个为0的位置newpop(i,mutapos)=0;newpop(itb(b)=l;endelsenewpop(i,:)=newpop(i,:);endendmain,m主程序%遗传算法优化稀疏天线阵问题clearallclcformatlonge;indivlength=1000;%个体串长度popsize=50;%初始种群大小globalatennum;atennum=20;%阵元个数eranum=200;%最大迭代次数pc=0.8;%交叉概率pm=0.05;%变异概率%pc
14、min=0.5;%pmmin=0.03;%dc1tapc=(pc-pcmin)*100/eranum;%deltapm=(pm-pmmin)*l00/eranum;Tl=clock;s=sprintfC程序正在运行中，请稍等)；disp(s);%第一步，实数编码%产生初始种群trace=zeros(eranum,indivlength+1);为最优轨迹pop=initpop(popsize,indivlength,atennum);fori=l:eranumfitvalue=Calfitval(popsize,indivlength,pop);Ebcstindivdual,bestfit=be
15、st(pop,fitvalue);selection=select(fitvalue,pop);newpop=cross(selection,pc);ncwpop=mutation(newpop,pm);pop=newpop;bestindivdua1;trace(i,l)=bestfit;trace(i,2:indivlength+l)=bestindivdual;%ifmod(i,100)=0%pc=pc-deltapc;%pm=pc-deltapm;%endendt=l:eranum;plot(t,-20*logl0(trace(:,1);XlabelC迭代次数)；ylabelC每代最优
16、值dB最;figure;%trace(:,1);bestfit,index=max(trace(:,1);%最优适应度值bestchrom=trace(index,2:indivlength+1);%最优染色体strl=sprintf(,进化到%d代，对应染色体为：%s,得本次求解的最优值%.5fdBn,findex,num2str(bestchrom),-20*logl0(bestfit);disp(strl);str2=find(bestchrom=l);str3=sprintf(阵元放置位置为：%s11,num2str(str2);%显示阵元位置disp(str3);i=l；a=zero
17、s(l,1801);%放置方向性函数值fortheta=-pi/2:pi/1800:pi/2fork=l:indivlengtha(i)=a(i)+bestchrom(k)*exp(j*pi*(k-l)*sin(theta);enda(i)=20*1og10(abs(a(i)/atennum);%取dBi=i+l;endtheta=-pi/2:pi/1800:pi/2;plot(theta*180pi,a);XIabe1(角度);ylabel(辐射强度dB);T2=clock;%计算运行时间elapscd_time=T2-Tl;ifelapsedtimc(6)0,elapsedtime(6)=
18、elapsed_time(6)+60;elapsed_time(5)=elapsedtime(5)-1;end%秒ifelapsedtime(5)0,elapsed_time(5)=elapsed_time(5)+60;elapsed_time(4)=elapsedtime(4)-1;end%分ifelapsed_time(4)0,elapsedtime(4)=elapscd_time(4)+24;elapsed_time(3)=elapscd_time(3)-1;end%时str2=sprintf(,程序运行耗时%d小时1分钟%.4f秒,elapsed_time(4),elapsed_time(5),elapsed_timc(6);disp(str2);
