QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 8542|回复: 28
打印 上一主题 下一主题

【原创】LU分解求解线性方程组Matlab代码

[复制链接]
字体大小: 正常 放大

2620

主题

162

听众

1万

积分

升级  0%

  • TA的每日心情
    开心
    2015-3-12 15:35
  • 签到天数: 207 天

    [LV.7]常住居民III

    社区QQ达人 发帖功臣 新人进步奖 优秀斑竹奖 金点子奖 原创写作奖 最具活力勋章 助人为乐奖 风雨历程奖

    群组第六届国赛赛前冲刺培

    群组国赛讨论

    群组2014美赛讨论

    群组2014研究生数学建模竞

    群组数学中国试看培训视频

    跳转到指定楼层
    1#
    发表于 2015-1-9 10:01 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    • function [x,L,U]=lusolve(A,b)
    • % LU分解求解线性方程组Ax=b
    • %
    • % 参数说明
    • % A,b:方阵A满足前n-1阶顺序主子式不为零,b为非零列向量
    • % x:方程组根
    • % L,U:L为单位下三三角阵,U的上三角阵(Doolittle分解),满足A=LU,
    • %
    • % LU分解的说明:
    • % 方阵A满足前n-1阶顺序主子式不为零
      游客,如果您要查看本帖隐藏内容请回复
    • % 则A经过行变换的得到一个上三角阵U,而初等行变换相当于左乘一个初等矩阵L,故A可分解为A=LU,
    • % 其中L或U为下或上三角矩阵,可以证明LU分解不唯一
    • %
    • % 如果L为单位下三角阵,U为上三角阵,那么称为Doolittle分解
    • % 如果L为下三角阵,U为单位上三角阵,那么称为Crout分解
    • %
    • % 可以证明Doolittle和Crout分解满足唯一性
    • %
    • % LU分解程组的思路:
    • % Ax=(LU)x=L(Ux)=Ly=b,其中x,y,b都是列向量
    • % 先求解方程Ly=b,计算出y,然后通过Ux=y求解x
    • %
    • % 对A进行Doolittle分解
    • [L,U]=Doolittle(A);
    • % 调用trimatrix()求解Ly=b
    • y=trimatrix(L,b,'L');
    • % 调用trimatrix求解Ux=y
    • x=trimatrix(U,y,'U');



    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    WSHXPY        

    62

    主题

    17

    听众

    5901

    积分

    升级  18.02%

  • TA的每日心情
    奋斗
    2017-9-16 20:51
  • 签到天数: 164 天

    [LV.7]常住居民III

    超级版主

    自我介绍
    来自哈工大,喜欢数学建模

    社区QQ达人 新人进步奖 发帖功臣

    群组2015美赛优秀论文解析

    群组2014第三期英语写作

    群组数学中国_哈工大群

    群组数学建模

    群组2016研赛备战群组

    回复

    使用道具 举报

    WSHXPY        

    62

    主题

    17

    听众

    5901

    积分

    升级  18.02%

  • TA的每日心情
    奋斗
    2017-9-16 20:51
  • 签到天数: 164 天

    [LV.7]常住居民III

    超级版主

    自我介绍
    来自哈工大,喜欢数学建模

    社区QQ达人 新人进步奖 发帖功臣

    群组2015美赛优秀论文解析

    群组2014第三期英语写作

    群组数学中国_哈工大群

    群组数学建模

    群组2016研赛备战群组

    回复

    使用道具 举报

    WSHXPY        

    62

    主题

    17

    听众

    5901

    积分

    升级  18.02%

  • TA的每日心情
    奋斗
    2017-9-16 20:51
  • 签到天数: 164 天

    [LV.7]常住居民III

    超级版主

    自我介绍
    来自哈工大,喜欢数学建模

    社区QQ达人 新人进步奖 发帖功臣

    群组2015美赛优秀论文解析

    群组2014第三期英语写作

    群组数学中国_哈工大群

    群组数学建模

    群组2016研赛备战群组

    回复

    使用道具 举报

    WSHXPY        

    62

    主题

    17

    听众

    5901

    积分

    升级  18.02%

  • TA的每日心情
    奋斗
    2017-9-16 20:51
  • 签到天数: 164 天

    [LV.7]常住居民III

    超级版主

    自我介绍
    来自哈工大,喜欢数学建模

    社区QQ达人 新人进步奖 发帖功臣

    群组2015美赛优秀论文解析

    群组2014第三期英语写作

    群组数学中国_哈工大群

    群组数学建模

    群组2016研赛备战群组

    回复

    使用道具 举报

    WSHXPY        

    62

    主题

    17

    听众

    5901

    积分

    升级  18.02%

  • TA的每日心情
    奋斗
    2017-9-16 20:51
  • 签到天数: 164 天

    [LV.7]常住居民III

    超级版主

    自我介绍
    来自哈工大,喜欢数学建模

    社区QQ达人 新人进步奖 发帖功臣

    群组2015美赛优秀论文解析

    群组2014第三期英语写作

    群组数学中国_哈工大群

    群组数学建模

    群组2016研赛备战群组

    回复

    使用道具 举报

    WSHXPY        

    62

    主题

    17

    听众

    5901

    积分

    升级  18.02%

  • TA的每日心情
    奋斗
    2017-9-16 20:51
  • 签到天数: 164 天

    [LV.7]常住居民III

    超级版主

    自我介绍
    来自哈工大,喜欢数学建模

    社区QQ达人 新人进步奖 发帖功臣

    群组2015美赛优秀论文解析

    群组2014第三期英语写作

    群组数学中国_哈工大群

    群组数学建模

    群组2016研赛备战群组

    回复

    使用道具 举报

    WSHXPY        

    62

    主题

    17

    听众

    5901

    积分

    升级  18.02%

  • TA的每日心情
    奋斗
    2017-9-16 20:51
  • 签到天数: 164 天

    [LV.7]常住居民III

    超级版主

    自我介绍
    来自哈工大,喜欢数学建模

    社区QQ达人 新人进步奖 发帖功臣

    群组2015美赛优秀论文解析

    群组2014第三期英语写作

    群组数学中国_哈工大群

    群组数学建模

    群组2016研赛备战群组

    回复

    使用道具 举报

    WSHXPY        

    62

    主题

    17

    听众

    5901

    积分

    升级  18.02%

  • TA的每日心情
    奋斗
    2017-9-16 20:51
  • 签到天数: 164 天

    [LV.7]常住居民III

    超级版主

    自我介绍
    来自哈工大,喜欢数学建模

    社区QQ达人 新人进步奖 发帖功臣

    群组2015美赛优秀论文解析

    群组2014第三期英语写作

    群组数学中国_哈工大群

    群组数学建模

    群组2016研赛备战群组

    回复

    使用道具 举报

    WSHXPY        

    62

    主题

    17

    听众

    5901

    积分

    升级  18.02%

  • TA的每日心情
    奋斗
    2017-9-16 20:51
  • 签到天数: 164 天

    [LV.7]常住居民III

    超级版主

    自我介绍
    来自哈工大,喜欢数学建模

    社区QQ达人 新人进步奖 发帖功臣

    群组2015美赛优秀论文解析

    群组2014第三期英语写作

    群组数学中国_哈工大群

    群组数学建模

    群组2016研赛备战群组

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2025-9-10 10:47 , Processed in 1.810317 second(s), 96 queries .

    回顶部