关于微波滤波器柯西法诊断

这是我写的一段用来通过S参数得到F P E 的matlab代码但是画出来不对
Rank=4;%滤波器阶数
wn=[-2 2];%有限传输零点
RL=20；%回波损耗
[S11,S21,LB]=nmax(Rank,wn,RL);%返回通带部分的S21 S11数据（nmax是一个多项式综合和求耦合矩阵的子函数）
w=-1.1:0.11:1.1;%直接使用归一化频率，满足采样点个数大于2N+Nz+2
S=1j*w;%根据论文里w变为S
S21=diag(S21);
S11=diag(S11);
N=length(S);
Mzero1=zeros(N,Rank+1);
Mzero2=zeros(N,length(wn)+1);
[Vn,Vnz]=vd(S,Rank,length(wn));%这是我写的一个求范德蒙行列式的子程序
X=[Vn Mzero2 -1*S11*Vn;
   Mzero1 Vnz -1*S21*Vn];
[U,T,V]=svd(X);
%取V的最后一列就是 F P E 的多项式系数
V=V(:,13);
F=V(1:5);
F= fliplr(F);
P=V(6:8);
P=fliplr(P);
E=V(9:13);
E=fliplr(E);
F=F/F(1)；P=P/P(1);E=E/E(1);%系数归一化
Fw=poly(roots(F)/1j);
Pw=poly(roots(P)/1j);
LB1=polyval(Fw,1);
LB2=polyval(Pw,1);
LB=LB1/(LB2*sqrt(10^(RL/10)-1));
LB=abs(LB);
w=-4:0.001:4;
h=freqs(F,E,w);
g=freqs(P,E,w);
S21=20*log10(abs(g)/LB);
S11=20*log10(abs(h));

好厉害，支持一下

                厉害了

申明：网友回复良莠不齐，仅供参考。如需专业解答，请学习本站推出的微波射频专业培训课程。