crosscouping_synthesis 微波滤波器耦合矩阵综合

微波滤波器耦合矩阵综合， 1 根据传输零点位置，回波损耗求出耦合矩阵； 2 从S21 S11多项式出发 综合出耦合矩阵； 3 已知传输零点位置，滤波器阶数，回波损耗画出滤波器的S参数，分别以4阶和八阶滤波器为例。
(

%滤波器（交叉）的耦合矩阵的综合  %%从S21 S11多项式出发 综合出耦合矩阵
%五阶滤波器
%通带频率 3.4—3.8GHz，反射系数 S11=20dB，零点位于 3.3GHz 和 3.9GHz 
clear;
wp=[inf inf inf 3.3e9 3.9e9 ]; %给定零点位置
RL=20;f1=3.4e9;f2=3.8e9;N=length（wp）; %通带内的回波损耗、滤波器通带边缘频率 N=5;
f0=sqrt（f1*f2）;bw=（f2-f1）/f0;dw=f2-f1;%中心频率、带宽、相对带宽
wp=（1/bw）*（wp./f0-f0./wp）;%归一化零点位置
syms w
w_yipi=sqrt（w^2-1）;% 表示的是w’
U=w-1/wp（1）;% 初始U1（w）
V=w_yipi*sqrt（1-1/wp（1）^2）;%初始V1（w）
for n=2:N
    Utemp=U;  Vtemp=V;
c=w-1/wp（n）;d=w_yipi*sqrt（1-1/wp（n）^2）; 
U=Utemp*c+Vtemp*d; V=Vtemp*c+Utemp*d;
end
P=1-w/wp（1）;
for n=2:N
     P=P*（1-w/wp（n））;%多项式P可以根据零点直接写出来
end
F=sym2poly（U）;%sym2poly的作用是提取多项式的系数，并按照降次排列。
P=sym2poly（P）;
F=F/F（1）; P=P/P（1）; %归一化
epsl=abs（polyval（P1）/polyval（F1））/sqrt（10^（RL/10）-1）;
%polyval的作用是计算多项式的函数值。
F2=conv（FF）;%conv的作用是计算多项式的乘积，计算结果是新的

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        3348  2020-10-31 06:49  crosscouping_synthesis\crosscouping_S_parameter.m
     文件        2801  2020-10-31 06:49  crosscouping_synthesis\crossmatrix_synthesis.m
     文件         763  2020-10-31 06:49  crosscouping_synthesis\cross_4_order.m
     文件         826  2020-10-31 06:49  crosscouping_synthesis\cross_8_order.m
     文件         570  2020-10-31 06:49  crosscouping_synthesis\M_rotation.m