声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1323|回复: 1

[编程技巧] 如何将计算出来的矩阵,提取出来.和提取矩阵元素是一样的吗?

[复制链接]
发表于 2006-8-6 16:24 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我编了一个程序,计算得出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
回复
分享到:

使用道具 举报

发表于 2006-8-10 15:06 | 显示全部楼层
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 编辑 ]
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-14 07:58 , Processed in 0.058109 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表