exjobb/scilab/plot s21/test1.sce

clear
close

function A=angle(C)
    A = atan(imag(C),real(C))/%pi*180;
endfunction

M = csvRead('test1.csv', ",", ".", "double", [], [], [], 2);
M = M(2:size(M)(1),1:size(M)(2)-1);

f = M(:,1);
t1 = M(:,2) + M(:,3)*%i;
t2 = M(:,4) + M(:,5)*%i;


t3 = M(:,6) + M(:,7)*%i;
t4 = M(:,8) + M(:,9)*%i;

//plot2d("ll", f, [abs(t1), abs(t2), abs(t3), abs(t4)]);
plot2d("ln", f, [angle(t1), angle(t2));



M = csvRead('test2.csv', ",", ".", "double", [], [], [], 2);
M = M(2:size(M)(1),1:size(M)(2)-1);

f = M(:,1);
t1 = M(:,2) + M(:,3)*%i;
t2 = M(:,4) + M(:,5)*%i;