10、电路综合-基于简化实频的宽带匹配电路设计方法

clear clc % 使用1GHZ的微带线，最高控制到3GHz f=1; fe=3; we=2*pi*fe; tau=pi/2/we; ele_l=360*tau*f; %使用6个级联微带线进行设计 k=6; %在DC处无零点 q=0; %初始化H的系数 h=[1 1 1 1 1 1 1]; disp(['此处使用在',num2str(f/1e9),'GHz下电长度为',num2str(ele_l),'°的微带线进行实现']); %0-3GHz的点数 Nopt=57; fr_opt=linspace(0,3,Nopt+1); % 定义S21参数的模值 mag=[1,0.952631578947368,0.905263157894737,0.857894736842105,0.810526315789474,0.763157894736842,0.715789473684211... ,0.668421052631579,0.621052631578947,0.573684210526316,0.526315789473684,0.478947368421053,0.431578947368421... ,0.384210526315789,0.336842105263158,0.289473684210526,0.242105263157895,0.194736842105263,0.147368421052631... ,0.100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000... ,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000... ,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000... ,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000... ,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000... ,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000,0.0100000000000000... ,0.0100000000000000,0.0100000000000000,0.0100000000000000]; % 定义相位，此处只设计幅值，相位随意了 phase=linspace(0,0,Nopt+1); % 画出目标电路的幅度 figure (1) plot(fr_opt,mag) xlabel('Frequency') ylabel('MAG-S21') title('要设计的微带线的目标S21参数') % 运行优化算法 n1=length(h); n=n1-1; h(n+1)=0; for i=1:n x0(i)=h(i); end % Call optimization with no transformer OPTIONS=optimset('MaxFunEvals',20000,'MaxIter',50000); x=lsqnonlin('objective',x0,[],[],OPTIONS,fe,q,k,fr_opt,mag, phase); h(n+1)=0; for i=1:n h(i)=x(i); end % 基于优化得到的h计算其他参数 [G,H,F,g]=SRFT_htoG(h,q,k); tau=1/4/fe; N=length(mag); j=sqrt(-1); % 将得到的解析形式画图 for i=1:N teta=2*pi*tau*fr_opt(i); omega=tan(teta); lmbda=j*tan(teta); fval=(-1)^q*(lmbda)^q*(1-lmbda^2)^(k/2); freal=real(fval); fimag=imag(fval); gval=polyval(g,lmbda); greal=real(gval); gimag=imag(gval); S21=fval/gval; MS21(i)=abs(S21); phase_f=atan(fimag/freal); phase_g=atan(gimag/greal); phase_S21(i)=phase_f-phase_g; ph(i)=teta; fr(i)=fr_opt(i); end figure(2) plot(fr,MS21, fr_opt,mag) xlabel('Frequency') ylabel('MS21, modula') title('Arbitraray amplitude approximation via SRFT') %------------------------------------ % 综合得到所需的微带电路 [Z_imp]=UE_sentez(h,g)