设为首页收藏本站|繁體中文 快速切换版块

 找回密码
 立即加入
搜索
查看: 980|回复: 1

电力系统算法

[复制链接]
  • TA的每日心情
    奋斗
    2018-12-13 09:01
  • 签到天数: 359 天

    连续签到: 2 天

    [LV.8]以坛为家I

    累计签到:359 天
    连续签到:2 天
    发表于 2017-8-7 08:40:42 | 显示全部楼层 |阅读模式
    电子图书
    电子图书名: 列主元消去法
    编者: 电小白
    内容简介: 算法
    所属专业方向: 电力系统
    出版社:
    来源:

    马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!

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

    ×
    function [X]=Gauss(A,B)7 f- v* s  b8 ~9 |7 s- F
    [m,n]=size(A);& O7 F0 Y7 }) r( [# l
    l=length(B);' B/ H6 q) Y$ k- [, A7 V) x
    if n~=m9 M* `/ i+ I) q. Q
        fprintf('wrong');$ _( W4 E$ M1 c& ]
        return;end
    , ^. Y' A3 q  N: Yif m~=l$ I% U1 W( J7 x9 p+ A" R+ a
        fprintf('wrong');+ a! W/ X# h- S* H- N! l
        return;end
    + R6 j& c/ ^& L9 r4 r3 xfprintf('go to calculate');
    - K1 I/ h- e  y0 t" j/ h  K. vfor k=1:n-1* u* F4 F; `# j' W
        max=0;
    , u( I) p% M2 d  `; C2 d* Q    for j=k:n
    & Y: n4 s: Y/ o; X        if abs(A(j,k))>max
    . z  m" D" X/ p4 a! T6 l& v            max=abs(A(j,k));mark=j;9 L. Q/ ^/ g9 t; l% W. C
            end, {2 l( Q: r2 d8 y) _1 N! H0 A
        end
    $ @5 s+ z' c" r- {: `9 m    if abs(max)<1e-5
    0 ^* O+ \4 f8 {& M8 a       fprintf('wrong');9 v! w4 R% V+ ?# s
           return;end+ C. Z, E1 ?+ R6 B* D2 @$ T
           for i=k:n$ }* [- u$ [. J7 Z' }: s9 T
              a=A(k,i);A(k,i)=A(mark,i);A(mark,i)=a;. @6 K* r) x" \4 K- ^
           end' }0 J0 Q+ Q* ?0 |
           b=B(k);B(k)=B(mark);B(mark)=b;7 _+ P, P. x$ z- R, N! w) [1 o
        for j=k+1:n
    0 y9 T5 \, [  F9 `       m=A(j,k)/A(k,k);2 d" ?+ X- e$ D; J( S3 \
           for i=k:n+ ?; B' _- K/ @- x
               A(j,i)=A(j,i)-m*A(k,i);4 ?% w* x: b$ k* e% U
           end! o& k0 s6 f5 M# N2 d1 W+ q+ K3 H% ?
           B(j)=B(j)-m*B(k);+ {% Y, s8 y6 u; ]
       end1 x1 x! [% W9 e9 M' h& |, P: q0 O8 n
    End$ l3 D. S. N: W% C
    for j=n:-1:1) d; e+ C* t& v* v
        for i=j+1:n
    , r4 v  M/ F- w& F1 g    B(j)=B(j)-A(j,i)*X(i);
    / F2 s+ M: H! O  f0 r    end& B) s5 L* d1 o+ f! F% N7 z
        X(j)=B(j)/A(j,j);
    ' `$ f% x! m: b& Eend
    * G4 B( n+ r# i" B! X0 k4 L# b; B* u3 ?) ?3 l" q2 v
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    [发帖际遇]: zibuyu_2017 发帖时在路边捡到威望2 点,偷偷放进了口袋. 幸运榜 / 衰神榜
    楼主热帖
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    奋斗
    2018-12-13 09:01
  • 签到天数: 359 天

    连续签到: 2 天

    [LV.8]以坛为家I

    累计签到:359 天
    连续签到:2 天
     楼主| 发表于 2017-8-7 08:42:58 | 显示全部楼层
    这个帖子不错,大家快来顶起来!
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    [发帖际遇]: zibuyu_2017捐助研学论坛,另外奖励威望1 点. 幸运榜 / 衰神榜
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    回复 推荐 踩下

    使用道具 举报

    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

    小黑屋|手机版|APP下载(beta)|Archiver|电力研学网 ( 赣ICP备12000811号-1|赣公网安备36040302000210号 )|网站地图

    GMT+8, 2025-4-25 02:56

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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