设为首页收藏本站 |繁體中文

电力研学论坛

 找回密码
 立即加入

QQ登录

只需一步,快速开始

搜索
查看: 171|回复: 1
收起左侧

电力系统算法

[复制链接]
最佳答案
0 
  • TA的每日心情
    调皮
    14 小时前
  • 签到天数: 97 天

    连续签到: 2 天

    [LV.6]常住居民II

    发表于 2017-8-7 08:45:02 | 显示全部楼层 |阅读模式

    快速获得阅读权限、下载权限、提升会员组攻略【新手必看】2016-03-03更新

    电子图书
    电子图书名: 三角分解法
    编者: 电小白
    内容简介: 算法
    所属专业方向: 电力系统
    出版社:
    来源:

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

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

    x
    function [L,U,y,x]=LU(A,B)
    n=length(A);
    L=eye(n);
    U=zeros(n,n);
    for k=1:n
      for j=k:n
          c=0;
          for i=1:k-1
              c=c+L(k,i)*U(i,j);
          end
          U(k,j)=A(k,j)-c;
      end
      if U(k,k)==0
          for i=k+1:n
              z=i;
              if U(k,i)~=0
                  break;end
          end
          for i=1:k
              a=U(i,k);U(i,k)=U(i,z);U(i,z)=a;
          end
          for i=1:n
              a=A(i,k);A(i,k)=A(i,z);A(i,z)=a;
          end
      end
      for i=k+1:n
          c=0;
          for q=1:k-1
              c=c+L(i,q)*U(q,k);
          end
          L(i,k)=(A(i,k)-c)/U(k,k);
      end
    end
      for k=1:n
          c=0;
          for i=1:k-1
              c=c+L(k,i)*y(i);
          end
          y(k)=(B(k)-c)/L(k,k);
      end
      for k=n:-1:1
          c=0;
          for i=k+1:n
              c=c+U(k,i)*x(i);
          end
          x(k)=(y(k)-c)/U(k,k);
      end

    楼主热帖
    真诚赞赏,手留余香
    [发帖际遇]: 天天上研学论坛zibuyu_2017研究有了新进展,老板犒劳了 学分1 点. 幸运榜 / 衰神榜
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    回复 呼我

    使用道具 举报

    最佳答案
    0 
  • TA的每日心情
    调皮
    14 小时前
  • 签到天数: 97 天

    连续签到: 2 天

    [LV.6]常住居民II

     楼主| 发表于 2017-8-7 08:47:56 | 显示全部楼层
    看到这帖子真是高兴!
    真诚赞赏,手留余香
    [发帖际遇]: 一个袋子砸在了 zibuyu_2017 头上,zibuyu_2017 赚了 学分3 点. 幸运榜 / 衰神榜
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

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

    GMT+8, 2017-11-20 23:25

    Powered by Discuz! X3.4

    © 2001-2017 Comsenz Inc.

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