Matlab 1

close all;

clear all;

C1=0.984*10^(-12);

L2=6.438*10^(-9);

C3=3.183*10^(-12);

L4=6.438*10^(-9);

C5=0.984*10^(-12);

S21=[];

S11=[];

F=[];

z=50;

for f=1:10^7:7*10^9

    w=2*pi*f;

    mat1=[1 0;j*w*C1 1];

    mat2=[1 j*w*L2;0 1];

    mat3=[1 0;j*w*C3 1];

    mat4=[1 j*w*L4;0 1];

    mat5=[1 0;j*w*C5 1];

    matf=mat1*mat2*mat3*mat4*mat5;

    a=matf(1,1);

    b=matf(1,2);

    c=matf(2,1);

    d=matf(2,2);

    S21=[S21 ; 2/(a+b/z+c*z+d)];

    S11=[S11 ; (a+b/z-c*z-d)/(a+b/z+c*z+d)];

    F=[F ;f];

end

subplot(2,1,1);

plot(F,abs(S21),'r',F,abs(S11),'g');

subplot(2,1,2);

plot(F,phase(S21),'r',F,phase(S11),'g');