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

 找回密码
 立即加入
搜索
查看: 965|回复: 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)
    ' h" F+ Y2 K; {/ k' B0 U$ ]* a5 F/ h0 J8 J[m,n]=size(A);6 t" G# X) f, U+ I$ j- ^: w
    l=length(B);
    % W1 g8 u4 f9 I  p/ ]if n~=m
    4 c  \/ O% N7 T! u; a% A    fprintf('wrong');0 I* [% V* V$ t8 ~/ l: P
        return;end. s7 A& A5 w% I9 V# z
    if m~=l- d1 r6 g9 d/ k  c* ?( W
        fprintf('wrong');
    ) o, G4 Q% {  s) D    return;end$ H" H7 ?4 n7 S* Z
    fprintf('go to calculate');/ g6 {* c  S8 ~2 O  {& `: i8 f
    for k=1:n-15 M; W1 }/ f: h+ D- U; ~1 n
        max=0;8 c, M# X0 e1 Z0 z' u% U, [
        for j=k:n4 T: b% q) S: E9 d! H- S
            if abs(A(j,k))>max* y+ {0 C& r& M
                max=abs(A(j,k));mark=j;. R# e$ d" D9 |2 S6 l0 h. }4 `
            end
    ! y8 c# {; b0 a, |    end
    8 F8 Z6 l7 l2 z    if abs(max)<1e-5# S3 N/ O4 q# u) \5 u8 d; c
           fprintf('wrong');# C/ i0 J3 g3 N; H) g
           return;end
    . ]7 k4 u4 F- _3 S! K& D. N9 N       for i=k:n$ I: j3 P# f: i& o, Y- I
              a=A(k,i);A(k,i)=A(mark,i);A(mark,i)=a;$ X6 h8 V( M7 y0 h; K9 e( y
           end
    0 i% e2 H# Z# ^3 w' q7 L4 F       b=B(k);B(k)=B(mark);B(mark)=b;
    3 X4 O0 q# }. H3 M* M0 |    for j=k+1:n
    ' A8 ^8 Z( m$ e9 {       m=A(j,k)/A(k,k);
      ?8 [1 \/ m. ^+ Y  M, V7 X       for i=k:n$ q& x2 {* C  e
               A(j,i)=A(j,i)-m*A(k,i);
    ! l5 A/ @7 a3 w1 @; e       end! [0 ^& |5 }- T0 O+ k& U1 ?+ [
           B(j)=B(j)-m*B(k);
    6 g. y0 W$ s" x% A, x, M   end4 A5 b) y; J" o9 J- ^8 X( q
    End. g1 m. `& a& r9 H
    for j=n:-1:18 P. D! T  X! u% a7 X* k2 F. ]' M
        for i=j+1:n  x7 i4 p& A- B3 M0 R( {! Z8 E
        B(j)=B(j)-A(j,i)*X(i);
    7 p& K2 I+ j# [( q( Z    end
    6 j6 E) q3 t; f9 v) O    X(j)=B(j)/A(j,j);$ d0 h0 i: ~! v+ b# O4 _8 ~0 ^" H
    end
    , J! t% Q* T3 X4 Z+ x) j9 `4 Y; z9 v. y; s  |& @# ~( Q
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    [发帖际遇]: 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-4 15:25

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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