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

这是我写的一段用来通过S参数得到F P E 的matlab代码但是画出来不对 m8V
}E&6  
Rank=4;%滤波器阶数 xE{slDl  
wn=[-2 2];%有限传输零点 yJ="dEn>i"  
RL=20；%回波损耗 I_
AFHrj  
[S11,S21,LB]=nmax(Rank,wn,RL);%返回通带部分的S21 S11数据（nmax是一个多项式综合和求耦合矩阵的子函数） g;Zy3   
w=-1.1:0.11:1.1;%直接使用归一化频率，满足采样点个数大于2N+Nz+2 gT(
8.<h8  
S=1j*w;%根据论文里w变为S 4hxa|f  
S21=diag(S21); "|Kag|(qB  
S11=diag(S11); fRT4>So  
N=length(S); ^#XQ2UN  
Mzero1=zeros(N,Rank+1); s\< @v7A  
Mzero2=zeros(N,length(wn)+1); ]:;dJc'  
[Vn,Vnz]=vd(S,Rank,length(wn));%这是我写的一个求范德蒙行列式的子程序 6a<zZO`Z6+  
X=[Vn Mzero2 -1*S11*Vn; G)q;)n;*=  
   Mzero1 Vnz -1*S21*Vn]; }>EWFE`  
[U,T,V]=svd(X); )}P/xY0  
%取V的最后一列就是 F P E 的多项式系数 E $P?%<o  
V=V(:,13); p1\mjM  
F=V(1:5); FU(2,V
l  
F= fliplr(F); ^wD`sj<Qg  
P=V(6:8); ^X?uAX-RP|  
P=fliplr(P); ;E
,i  
E=V(9:13); I9dX\w}  
E=fliplr(E); 0cDP:EzR;  
F=F/F(1)；P=P/P(1);E=E/E(1);%系数归一化 7?D?s!%\  
w:/3%-  
Fw=poly(roots(F)/1j); cvQAo|  
Pw=poly(roots(P)/1j); tyWDa$u,u  
LB1=polyval(Fw,1); %@:6&  
LB2=polyval(Pw,1); Cg|\UKfy$  
LB=LB1/(LB2*sqrt(10^(RL/10)-1)); rP}[>  
LB=abs(LB); s7sTY   
%WC^aKfY  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

好厉害，支持一下

X;W
0r5T  
F*@2)  
                厉害了

厉害啊，完全看不懂

厉害 厉害