M=[0.0586,-0.0147,-0.2374,-0.0578,0.4314,-0.4385,0;
-0.0147,-0.0810,0.4825,0.3890,0.6585,0.0952,-1.3957;
-0.2374,0.4825,0.2431,-0.0022,0.3243,-0.2075,0.1484;
-0.0578,0.3890,-0.0022,-0.0584,-0.3047,0.4034,-0.0953;
0.4314,0.6585,0.3243,-0.3047,0.0053,-0.5498,-0.1628;
-0.4385,0.0952,-0.2075,0.4034,-0.5498,-0.5848,-0.1813;
0,-1.3957,0.1484,-0.0953,-0.1628,-0.1813,0.0211];
I=eye(7)
SS21=[];
SS11=[];
a=0.7220;
b=2.2354
SS21=[];
SS11=[];
a=0.7220;
b=2.2354;
T=[a,0,0,0,0,0,b];
R=diag(T)
for w=-2:0.001:4;
cx=M+w*I-j*R;
cx2=inv(cx);
ss21=(-2*j*sqrt(a*b)*cx2(7,1));
ss11=1+2*j*sqrt(a*b)*cx2(1,1);
s21=-20*log10(ss21);
s11=-20*log10(ss11);
SS21=[SS21,s21];
SS11=[SS11,s11];
end
w=-2:0.001:4;
plot(w,SS11,'-',w,SS21,'-');
grid on;
axis([-inf,inf,-inf,40])
呼唤英雄啊!! |