www.gusucode.com > matlab编程遗传算法计算匹配电路源码程序 > code1/code/MATLAB源代码/matlab_match/match2.m
function [ObjV,S21_passband_max,root_mean_square,VSWR1_max,VSWR2_max] = match2(Chrom,num_integer,num_parameter) %UNTITLED3 此处显示有关此函数的摘要 % 此处显示详细说明 %ChromCopy=Chrom; column_of_compoment=num_integer+1:num_integer+(num_integer-1)*num_parameter; %代表元件值的基因所在列的列标号 num_structure=1:num_integer; %单元结构的标号 [num_of_individual,num_of_variable]=size(Chrom); %提取出种群的个体数目和变量数 %frequency_1GHz=1e9.*[6 6.1 6.2 6.3 6.4 6.5 6.6 6.7 6.8 6.9 ... % 7.0 7.1 7.2 7.3 7.4 7.5 7.6 7.7 7.8 7.9 ... % 8.0 8.1 8.2 8.3 8.4 8.5 8.6 8.7 8.8 8.9 ... % 9.0 9.1 9.2 9.3 9.4 9.5 9.6 9.7 9.8 9.9 10]; frequency_1GHz=1e9.*[1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 ... 2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 ... 3.0 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 ... 4.0 4.1 4.2 4.3 4.4 4.5 4.6 4.7 4.8 4.9 5.0]; num_of_point=length(frequency_1GHz); %用于计算的频率点的数目 ABCD=zeros(2,2,num_of_individual,num_of_point); %预先给ABCD矩阵分配内存空间 for f=1:num_of_point %频率点循环 w=2*pi*frequency_1GHz(f); for k=1:num_of_individual %个体循环 % ABCD(:,:,k,f)=[1 0;1/(1i*w*L0)+1i*w*C0 1]; ABCD(:,:,k,f)=[1 0;0 1]; for m=2:Chrom(k,1)+1 increment=(m-2)*num_parameter; electrical_length=10/360*2*pi*frequency_1GHz(f)/3e9; Z01=Chrom(k,column_of_compoment(1+increment)); electrical_length1=Chrom(k,column_of_compoment(2+increment))/360*2*pi*frequency_1GHz(f)/3e9; Z02=Chrom(k,column_of_compoment(3+increment)); electrical_length2=Chrom(k,column_of_compoment(4+increment))/360*2*pi*frequency_1GHz(f)/3e9; R=Chrom(k,column_of_compoment(5+increment)); switch Chrom(k,num_structure(m)) %判断单元电路的结构 case 1 %1代表传输线 A=cos(electrical_length1); B=1i*Z01*sin(electrical_length1); C=1i*sin(electrical_length1)/Z01; D=cos(electrical_length1); case 2 %2代表开路并联传输线 A=1; B=0; C=1i*tan(electrical_length1)/Z01; D=1; % case 3 %3代表短路并联传输线 % A=1; % B=0; % C=1/(1i*Z01*tan(electrical_length1)); % D=1; case 3 %4代表两段的开路并联传输线 A=1; B=0; C=1i*(Z02*tan(electrical_length1)+Z01*tan(electrical_length2))/(Z02*(Z01-Z02*tan(electrical_length1)*tan(electrical_length2))); D=1; % case 5 %5代表两段的短路并联传输线 % A=1; % B=0; % C=(Z02-Z01*tan(electrical_length1)*tan(electrical_length2))/(1i*Z02*(Z01*tan(electrical_length1)+Z02*tan(electrical_length2))); % D=1; % case 4 %6代表串联电阻 % A=1; % B=R; % C=0; % D=1; % case 7 %7代表短路并联传输线+电阻 % A=1; % B=0; % C=1/(R+1i*Z01*tan(electrical_length1)); % D=1; case 4 %8代表开路并联传输线+电阻 Yin1=1i*tan(electrical_length1)/(Z01+1i*tan(electrical_length1)*R); Zin1=1/Yin1; Zin=Z01*(Zin1+1i*Z01*tan(electrical_length))/(Z01+1i*Zin1*tan(electrical_length)); A=1; B=0; % C=1i*tan(electrical_length1)/(Z01+1i*tan(electrical_length1)*R); C=1/Zin; D=1; case 5 %9代表两段的开路并联传输线+电阻 Yin1=1/((Z02*(Z01-Z02*tan(electrical_length1)*tan(electrical_length2)))/(1i*(Z02*tan(electrical_length1)+Z01*tan(electrical_length2)))+R); Zin1=1/Yin1; Zin=Z02*(Zin1+1i*Z02*tan(electrical_length))/(Z02+1i*Zin1*tan(electrical_length)); A=1; B=0; % C=1/((Z02*(Z01-Z02*tan(electrical_length1)*tan(electrical_length2)))/(1i*(Z02*tan(electrical_length1)+Z01*tan(electrical_length2)))+R); C=1/Zin; D=1; % case 10 %10代表两段的短路并联传输线+电阻 % A=1; % B=0; % C=1/((1i*Z02*(Z01*tan(electrical_length1)+Z02*tan(electrical_length2)))/(Z02-Z01*tan(electrical_length1)*tan(electrical_length2))+R); % D=1; end ABCD(:,:,k,f)=ABCD(:,:,k,f)*[A B;C D]; end end end S21_all=zeros(num_of_point,num_of_individual); S11_all=zeros(num_of_point,num_of_individual); S22_all=zeros(num_of_point,num_of_individual); VSWR1=zeros(num_of_point,num_of_individual); VSWR2=zeros(num_of_point,num_of_individual); %polar_point=zeros(num_of_point,num_of_individual); for f=1:num_of_point RS=50; XS=0; RL=50; XL=0; ZS=RS+1i*XS; ZL=RL+1i*XL; AA=squeeze(ABCD(1,1,:,f)); BB=squeeze(ABCD(1,2,:,f)); CC=squeeze(ABCD(2,1,:,f)); DD=squeeze(ABCD(2,2,:,f)); Denominator=(AA+BB/50+CC*50+DD); S21_all(f,:)=-20*log10(abs(2./Denominator)); %由ABCD矩阵计算S21 S11_all(f,:)=abs((AA+BB/50-CC*50-DD)./Denominator); VSWR1(f,:)=(1+S11_all(f,:))./(1-S11_all(f,:)); S22_all(f,:)=abs((-AA+BB/50-CC*50+DD)./Denominator); VSWR2(f,:)=(1+S22_all(f,:))./(1-S22_all(f,:)); % S11_all(f,:)=-20*log10(abs((squeeze(ABCD(1,1,:,f))+squeeze(ABCD(1,2,:,f))/50-squeeze(ABCD(2,1,:,f))*50-squeeze(ABCD(2,2,:,f)))./(squeeze(ABCD(1,1,:,f))+squeeze(ABCD(1,2,:,f))/50+squeeze(ABCD(2,1,:,f))*50+squeeze(ABCD(2,2,:,f))))); % polar_point(f,:)=20*log10(abs((squeeze(ABCD(1,1,:,f))+squeeze(ABCD(1,2,:,f))/50+squeeze(ABCD(2,1,:,f))*50+squeeze(ABCD(2,2,:,f))))-2); end ObjV=zeros(num_of_individual,1); S21_all_cor=zeros(num_of_point,num_of_individual); S21_passband_max=zeros(num_of_individual,1); %S21_passband_min=zeros(num_of_individual,1); VSWR1_max=zeros(num_of_individual,1); VSWR2_max=zeros(num_of_individual,1); %average_S21=zeros(num_of_individual,1); root_mean_square=zeros(num_of_individual,1); for k=1:num_of_individual % x=30:88; % y=-0.0129.*x.^2+1.5187.*x-33; f1=6:0.1:10; y=abg(f1); % y=5:-0.1:1; % x=6:0.1:10; % y=-2.712.*x.^2+43.34.*x-161.5; S21_all_cor(:,k)=S21_all(:,k)-y'; S21_passband_max(k)=max(S21_all(:,k)); %通带的最大插损 % S21_passband_min(k)=min(S21_all(:,k)); VSWR1_max(k)=max(VSWR1(:,k)); VSWR2_max(k)=max(VSWR2(:,k)); % average_S21(k)=sum(S21_all_cor(:,k))/length(S21_all_cor(:,k)); % for f=1:num_of_point % root_mean_square(k,1)=root_mean_square(k,1)+(S21_all_cor(f,k)-average_S21(k)).^2; % end % root_mean_square(k,1)=sqrt(root_mean_square(k,1)/length(S21_all_cor(:,k))); % f=1:num_of_point; % x=1-(f-1)./(num_of_point-1); % max_difference_S21=max(abs(S21_all_cor(:,k))); % max_difference_VSWR=VSWR1_max(k)-1+VSWR2_max(k)-1; % ObjV(k,1)=max_difference_S21+max_difference_VSWR; ObjV(k,1)=ObjV(k,1)+sum((S21_all(:,k)-y').^2)+sum((VSWR1(:,k)-1).^2)+sum((VSWR2(:,k)-1).^2); % if S21_passband_max(k)>1.3 % ObjV(k,1)=S21_passband_max(k); % else % ObjV(k,1)=S21_passband_max(k)+log10(root_mean_square(k,1)); % end % ObjV(k,1)=S21_passband_max(k); end end