微带功分器的仿真和优化

奇模激励
3端口的奇模电压与2端口的奇模电压大小相同，方向相反：V3o=—V2o。2端口的奇模电压为
S32=V3/V2=（V3e+V3o）/（V2e+V2o）=（V2e—V2o）/（V2e+V2o）
编程画出频率响应：
程序如下：
f0=2.225e9;
f=1.0e9:1e7:3.5e9;
z2e=sqrt(2)*(2+sqrt(2)*j*tan(pi/2*(f./f0)))./(sqrt(2)+2*j*tan(pi/2*(f./f0)));
图（5）
在电子圆图中设定好频率及系统特性阻抗后，可得到电容为2.6pF，电感为12.0nH。
在圆图中的匹配路径上，点3的结点品质因子最大，为 ，所以有载品质因子 ，则可估计带宽为： 。
第二章网络理论
1.用matlab编程对上面第二题中的匹配电路实现ABCD参数矩阵连乘，求匹配电路的带宽，并与由节点品质因子得到的结果进行比较。下面提供3种解法，分别是由功率，S参量及A参量的频率响应特性求解。
C=2.6e-12;
w=logspace(8,10,10000);
Z1=((1-w.^2*L*C)*Z+i*w*L)./(i*w*C*Z+1);
r=(Z1-50)./(Z1+50);
figure(1);
plot(w/(2*pi),abs(r));gridon;
xlabel('频率（Hz）');
ylabel('反射系数');
xlabel('频率(GHz)');
ylabel('|S41|');
%-------------------------------------------------
S41的频率特性曲线如图1所示。理想情况下，S11和S41均为零，此处出现非零值是由于 并不严格等于 的平方所致。
图1
图2是用耦合带状线在ads中的仿真结果。与前面计算结果一致。若用微带线仿真，则曲线形状便不同了，因为其偶模及奇模的传输速度不一致（β不同）。
A=[1 j*w*l1;0 1];％串联电感
B=[1 0;j*w*c2 1];％并联电容
C=[1 0;1/zl 1];％并联负载
D=E*A*B*C;
DA=1/D(1,1);
Hh=subs(DA,ff,f);
semilogx(f,-20*log10(abs(Hh)));
hold on;
y=min(-20*log10(abs(Hh)));
S11=(z1-Z0)./(z1+Z0);
plot(f,-20*log10(abs(S11)));
xlabel('Frequency(GHz)');
ylabel('S11(dB)');
对于S41，
由于： ；
；
其中： ；
； ;
所以：
程序如下：
%参数S41
Pe=(Z0-Z0e)/(Z0+Z0e);
Po=(Z0-Z0o)/(Z0+Z0o);
v2e=2*z2e./(1+z2e);
v2o=(1./(1+1./(j.*sqrt(2).*tan(pi*f./(2*f0)))))./(1./(1+1./(j.*sqrt(2).*tan(pi*f./(2*f0))))+1)*2;
t=(2-sqrt(2))/(2+sqrt(2));
s32=(v2e-v2o)./(v2e+v2o);
图（3）
注意：软件TXLINE 2001—Microstrip中的Electrical Length代表的就是 ，即电长度乘以 ，而不是乘以 。
2.用特性阻抗50 的传输线给天线馈电，其驻波比及驻波相位分别是3.0及10cm，求此天线的输入阻抗。设计匹配电路使天线与传输线达到匹配，并估计其带宽。设工作频率为886MHz。
总之要求带宽越宽，匹配器件越多，不分集总参数还是分布参数。
第三章功率分配器
1.重做等功分器示例，频率为2.2＋0.xxx(GHz)。
解：由辅助设计软件可以得到，特性阻抗为1.414 （即70.7 ）的微带线宽度为1.507mm，长度为17.734mm。如图1所示。
图1
特性阻抗为 （即50 ）的微带线宽度为2.8803mm。如图2所示。
解1： ，其中 是端口1的输入电阻。
功分器传输线长： ， ，由输入电阻计算公式可以得到， 。将此式代入，得到 表达式为：
。（1）
图8是用上式绘得的 的频率响应曲线，图9是第1题中仿真得到的结果。可以看到两者非常一致。
图8
图9
为求得 ,等效电路如下所示，左端为端口2，中间为端口1，右边为端口3.
端口2、3间可以看成两个网络的并联，电阻2 的归一化导纳矩阵为：
解：工作频率为886MHz，则波长为： ，则驻波相位10 cm对应的电长度为0.295。如图（4）所示，点1为波节点，沿驻波比为3的圆逆时针走0.295个点长度到点2，点2即为负载（为天线）阻抗点，即天线的输入阻抗： 。
图（4）
为实现天线与传输线的匹配，可采取并容串感的方案（路线为2—3—4），如图(5)所示。
z=y+3;
plot(f ,z);
运行结果如右图，可以看到
匹配网络是个低通滤波器，
其带宽为1.22G，与节点法
的结果基本一致。
思考题：用2节及3节四分之一波长阻抗变换器实现200欧至50欧的阻抗匹配。要求带宽最宽。
解：先看单节匹配情况。在电子圆图上节点1阻抗为200，节点2为50，单节时匹配路径上最大Q值为0.73
S21=2./(a+b+c+d);
plot(f,20*log10(abs(S21)));
xlabel('频率（MHz）');
ylabel('幅度（dB）');
title('S21的频率响应');
grid;
解3利用圆图匹配时可以看到负载点Q值最大，为1.8，这里的Q是节点品质因数，有载品质因数为0.9，所以带宽BW=0.88e9/0.9=0.98G。
V2e=(1+Pe)*Zine./((Z0+Zine).*(exp(j*bl)+Pe*exp(-j*bl)));
V2o=(1+Po)*Zino./((Z0+Zino).*(exp(j*bl)+Po*exp(-j*bl)));
V4=V2e-V2o;
figure(2);
plot(f,abs(V4));grid;
ze=Z0e1*(Z0+j*Z0e1*tan(pi/2*(f./f0)))./(Z0e1+j*Z0*tan(pi/2*(f./f0)));
zo=Z0o1*(Z0+j*Z0o1*tan(pi/2*(f./f0)))./(Z0o1+j*Z0*tan(pi/2*(f./f0)));
z1=Z0+(2*(ze.*zo-Z0^2)./(ze+zo+2*Z0));
对于两节的情况，关键就是调节各节的特性阻抗值使每一节所经过匹配路径的最大Q值相等。当特性阻抗分别为140及70时，满足条件，Q约为0.34。可见带宽是单节时的2倍多。
同样可得到三节的结果。在电子圆图上的情况如下图
特性阻抗分别为158，100及63.5，满足条件，Q约为0.22。可见带宽是单节时的3倍多。
图2
2.设计耦合度为-10dB平行线定向耦合器。特性阻抗50Ω。使用的板材参数为相对介电常数εr =4.25，介质厚度h=1.45mm,敷铜厚度t＝0.035mm。频率为2.386(GHz)。
解：根据耦合度、特性阻抗、板材参数和频率，利用ADS软件可以得到耦合微带线尺寸，如图3所示。
解：用传输线实现ZL至A的移动，应该并联终端开路传输线，并联短路的也可以，但其电长度会大于0.25。用(如图1)所示电子SMITH圆图，可知并联（终端开路）传输线长度为 (这里的λg不是仅指在波导中的波长，是指在所有传输线中传输的波长)，串联传输线的长度为0.152 。
图（1）
上面得到的是传输线的电尺寸，包括特性阻抗（50欧）及电长度。还需要得到物理尺寸，包括宽度长度（注意不同传输线的物理尺寸不同）。在TXLINE软件中选择微带传输线，把板材参数设置好。接下来由电参数特性阻抗及电长度求物理参数宽度及长度。
figure(2);
VSWR=(1+abs(r))./(1-abs(r));
plot(w/(2*pi),VSWR);gridon;
xlabel('频率（Hz）');
ylabel('驻波比');
%-----------------------------------------------------
若以驻波比小于1.5为带宽的定义，得到的带宽为267.4MHz（如图8）。
图（8）
解2：原匹配电路为串感并容型，得到A参量为
,
再有A矩阵转化为S参量，最后由S21画出频率响应如下图：
由上图可以看出电路带宽为1.2GHz，与节点品质因子计算得到的带宽1.291GHz大致相同。其中计算过程的matlab代码如下：
最优化大作业
——微带功分器的仿真和优化
班级：B1103491
学号：1110349019
姓名：胡亮
日期：2011-11-16
前言：本文主要是基于微带结构的功分器用ADS仿真并用ADS软件进行优化与matlab优化进行了比较：
第一章传输线理论
1.求解下图匹配路径的匹配电路参数。若将并容串线方案在介电常数为4.2介质材料厚度为1.586mm，导体带厚度为0.035mm微带基板上实现，求各段微带线的宽度及长度值。