如何将计算出来的矩阵,提取出来.和提取矩阵元素是一样的吗?
我编了一个程序,计算得出L,C矩阵均为13个,我要将L C矩阵逐个提取出来,进行 eig(C*L)运算,该如何提取?谢谢大家!
程序如下:
clear
a=4*pi*1E-7;b=8.85419*1E-12;% 介质常数
h=17;
h1=h;h2=h; %导线距地面高度
r1=0.3215; %导线半径
R2=1.46; %绞线之间的间距
r2=0.73; %绝缘半径
c=0:pi/6:2*pi;
%计算电感矩阵
SIN=sin(c);COS=cos(c);
L11=a/(8*pi)+(a/(2*pi))*log(2*h/r1);
L12=(a/(4*pi))*log(1+4*h*(h-r2*SIN)./((h-r2*COS).^2+(r2*SIN).^2));
L13=a*log(1+4*h.*(h+r2*SIN)./((h+r2*COS).^2+(r2*SIN).^2))./(4*pi);
L23=(a/(4*pi))*log(1+4*(h-r2*SIN).*(h+r2*SIN)/(R2^2));
L22=a/(8*pi)+a.*log(2*(h-r2*SIN)/r1)/(2*pi);
L33=a/(8*pi)+(a/(2*pi))*log(2*(h+r2*SIN)./r1);
L21=L12;L31=L13;L32=L23;
for i=1:13
L=[L11 L12(i) L13(i);
L21(i) L22(i) L23(i);
L31(i) L32(i) L33(i)];
C=a*b*inv(L);
end clear
a=4*pi*1E-7;b=8.85419*1E-12;% 介质常数
h=17;
h1=h;h2=h; %导线距地面高度
r1=0.3215; %导线半径
R2=1.46; %绞线之间的间距
r2=0.73; %绝缘半径
c=0:pi/6:2*pi;
%计算电感矩阵
SIN=sin(c);COS=cos(c);
L11=a/(8*pi)+(a/(2*pi))*log(2*h/r1);
L12=(a/(4*pi))*log(1+4*h*(h-r2*SIN)./((h-r2*COS).^2+(r2*SIN).^2));
L13=a*log(1+4*h.*(h+r2*SIN)./((h+r2*COS).^2+(r2*SIN).^2))./(4*pi);
L23=(a/(4*pi))*log(1+4*(h-r2*SIN).*(h+r2*SIN)/(R2^2));
L22=a/(8*pi)+a.*log(2*(h-r2*SIN)/r1)/(2*pi);
L33=a/(8*pi)+(a/(2*pi))*log(2*(h+r2*SIN)./r1);
L21=L12;L31=L13;L32=L23;
for i=1:13
L=[L11 L12(i) L13(i);
L21(i) L22(i) L23(i);
L31(i) L32(i) L33(i)];
C=a*b*inv(L);
e=eig(C*L);ans(i)=e(1);
end
eig是求特征值吧，好像每个C*L的三个特征值都相同（所以我只取了其中一个），而且好像每一个矩阵答案都一样啊。
[ 本帖最后由 ericlin 于 2006-8-10 15:11 编辑 ]
页:
[1]