nra.m

来自「国外教材MIMOSignalsAndSystems附带的程序」· M 代码 · 共 52 行

% Program "nra.m"

% Implementation of Node Reduction Algorithm
% Calculates admittance matrix Yr of 
% equivalent two-port network.

% Input Parameters:
% - Number of internal and exterior nodes: P
% - Admittances connecting nodes

P = 1;

while P < 2;
  P = input('Number of Nodes P = ');

   if P < 2
     disp('Specify more nodes!')
   end;
end;

% Input of admittances connecting nodes
disp('Make sure off-diagonal elements have negative sign')
disp(' ')

for k=1:P;
   for n=k:P;

	if k==n;
	 disp('Admittance of'), disp(n), disp('th node');
	 Y(k,n)=input('Y(k, n) = ');
	else	
    disp('Admittances connecting nodes'),disp(k),disp(n);
    Y(k,n)=input('Y(k, n) = ');
    Y(n,k)=Y(k,n);
	end;

    end;
end;
% ===============

Yr=Y;

% Calculation of reduced admittance matrix, Yr
for k=1:P-2;   
   Yr=Yr(1:P-k,1:P-k)-Yr(1:P-k,P+1-k)*Yr(P+1-k,1:P-k)/Yr(P+1-k,P+1-k);
end;
% ===============

disp(' ')
disp('Original(P-by-P) Admittance Matrix:'),Y
disp('Reduced (2-by-2) Admittance Matrix:'),Yr