23 lines
941 B
Mathematica
23 lines
941 B
Mathematica
|
function [new_G,new_B,GB,Y,r,c,Angle] = modifyadmmatrix(ii,jj,G,B)
|
|||
|
|
%**************************************************************************
|
|||
|
|
% <EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD> : <EFBFBD>Ӻ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>γɽڵ㵼<EFBFBD>ɾ<EFBFBD><EFBFBD><EFBFBD>Y
|
|||
|
|
% <EFBFBD><EFBFBD> <EFBFBD>ߣ<EFBFBD>
|
|||
|
|
% <EFBFBD><EFBFBD><EFBFBD><EFBFBD>ʱ<EFBFBD>䣺2010.12
|
|||
|
|
%**************************************************************************
|
|||
|
|
%% ֧·<EFBFBD><EFBFBD><EFBFBD>ɼ<EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
new_G=G;
|
|||
|
|
new_G(ii,jj)=new_G(ii,jj)-G(ii,jj);
|
|||
|
|
new_G(jj,ii)=new_G(jj,ii)-G(jj,ii);
|
|||
|
|
new_G(ii,ii)=new_G(ii,ii)+G(ii,jj);
|
|||
|
|
new_G(jj,jj)=new_G(jj,jj)+G(ii,jj);
|
|||
|
|
new_B=B;
|
|||
|
|
new_B(ii,jj)=new_B(ii,jj)-B(ii,jj);
|
|||
|
|
new_B(jj,ii)=new_B(jj,ii)-B(jj,ii);
|
|||
|
|
new_B(ii,ii)=new_B(ii,ii)+B(ii,jj);
|
|||
|
|
new_B(jj,jj)=new_B(jj,jj)+B(ii,jj);
|
|||
|
|
|
|||
|
|
%% <EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ʽ
|
|||
|
|
GB = new_G+new_B.*1i; %<EFBFBD><EFBFBD><EFBFBD>絼<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ɺϲ<EFBFBD><EFBFBD><EFBFBD>д<EFBFBD>ɸ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ʽ
|
|||
|
|
Y = abs(GB); %<EFBFBD><EFBFBD><EFBFBD>ڵ㵼<EFBFBD>ɷ<EFBFBD>ֵ
|
|||
|
|
[r,c] = find(Y);
|
|||
|
|
Angle = angle(GB(GB~=0)); %<EFBFBD><EFBFBD><EFBFBD>ڵ㵼<EFBFBD>ɽǶ<EFBFBD>
|