应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » matlab 求解线性方程组Ax=b
111/1返回列表
查看: 3431  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

mathstudy

金虫 (正式写手)

[交流] matlab 求解线性方程组Ax=b

x=inv(A)*b;和x=A\b;和迭代法  有什么区别？ 前面两个求 A逆 的是基于什么？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2014-02-07 13:07:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

tiancanrong

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
如果你要知道inv(A)*b与Gauss 消元法的区别，就要知道计算机运算和你自己用笔算是有区别的。你用笔算可以用分数，理论上是没有误差的；然而计算机运算时受限于数的精度，理论上是有误差的。所以同样的题目，计算机用不同的算法就是误差或效率会有差别。另外说一句，inv（A）也需要很复杂的计算的，Gauss消去法也需要条件的。当然了，由于现在电脑的软硬件都很先进了，一般的方程随便什么方法解都没区别了。
赞一下(1人) 回复此楼
5楼2014-02-08 22:00:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiancanrong

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我大致看了一下matlab help，A\b 根据A的情况用到几种不同的算法，如果A是一个方阵，b是一个列向量，此时Ax=b就是你说的线性方程组，\用的算法就是Gauss消去法；如果A不是方阵，那么Ax=b就是一个矩阵方程，其实这时候精确解不一定存在，matlab中的\就用了最小二乘法的思想求了一个近似的x（也可称为广义解），算法是QR分解。

可以这个简单理解，inv(A)*b是A\b的一种特殊情况（线性方程组时），但是A\b更强大，可以求出广义解。
在计算线性方程组Ax=b时，inv(A)*b和A\b用的数值算法都是Gauss消去法。至于你说的效率，A\b由于考虑的情况比较庞大，调用的时候稍微会慢一点，我不认为有质的区别。
赞一下(1人) 回复此楼
7楼2014-02-09 12:15:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

feixiaolin

荣誉版主 (文坛精英)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
A=[……],  b=[……],  s=rref([A,b])为佳，如果非满秩，也可知道。
赞一下 回复此楼
2楼2014-02-07 13:28:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴钢盔的猪头

木虫 (知名作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
如果A不是方阵，它的逆是不存在的。这时inv(A)无法用，但可以用广义逆pinv(A)。而A\b何时都可用。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
3楼2014-02-07 17:09:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)
A是可逆方阵， 我想请教的的是  A\b 和 inv(A)*b 中 他们求A逆的 原理 有什么区别， 为什么  A\b比 inv（A）*b 的效率高？ 另外  如果是 Gauss 消元法呢？
赞一下 回复此楼
4楼2014-02-07 19:11:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)
引用回帖:
5楼: Originally posted by tiancanrong at 2014-02-08 22:00:36
如果你要知道inv(A)*b与Gauss 消元法的区别，就要知道计算机运算和你自己用笔算是有区别的。你用笔算可以用分数，理论上是没有误差的；然而计算机运算时受限于数的精度，理论上是有误差的。所以同样的题目，计算机用 ...

那  inv(A)*b与 Gauss 和A\b 三者  在matlab 上面的效率为什么不一样？  他们在 matlab 上的运行方式有什么区别？ (比如说 A\b 可能是基于LU 分解而求解方程的)
赞一下 回复此楼
6楼2014-02-09 08:47:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)
引用回帖:
7楼: Originally posted by tiancanrong at 2014-02-09 12:15:12
我大致看了一下matlab help，A\b 根据A的情况用到几种不同的算法，如果A是一个方阵，b是一个列向量，此时Ax=b就是你说的线性方程组，\用的算法就是Gauss消去法；如果A不是方阵，那么Ax=b就是一个矩阵方程，其实这时 ...

谢谢啊.....但是我的运行结果显示  A\b比自己的Gauss程序快很多
赞一下 回复此楼
8楼2014-02-09 19:33:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiancanrong

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
8楼: Originally posted by mathstudy at 2014-02-09 19:33:04
谢谢啊.....但是我的运行结果显示  A\b比自己的Gauss程序快很多...

你自己编的Gauss程序，和matlab自带的程序库里的Gauss算法调用的方式应该有区别的，亲儿子当然快一点。另外问一下，你如何测试程序的快慢的？为什么要测试快慢？

[ Last edited by tiancanrong on 2014-2-10 at 12:26 ]
赞一下 回复此楼
9楼2014-02-10 12:14:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)
引用回帖:
9楼: Originally posted by tiancanrong at 2014-02-10 12:14:12
你自己编的Gauss程序，和matlab自带的程序库里的Gauss算法调用的方式应该有区别的，亲儿子当然快一点。另外问一下，你如何测试程序的快慢的？为什么要测试快慢？
...

好吧  谢过了    测试时间的话 可以用cpu时间啊.  有的时候还是注重算法的快慢的...
赞一下 回复此楼
10楼2014-02-12 13:43:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cced_1998

木虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
It takes almost two and one half times as long to compute the solution with y = inv(A)*b as with z = A\b. Both produce computed solutions with about the same error, 1.e-6, reflecting the condition number of the matrix. But the size of the residuals, obtained by plugging the computed solution back into the original equations, differs by several orders of magnitude. The direct solution produces residuals on the order of the machine accuracy, even though the system is badly conditioned.

The behavior of this example is typical. Using A\b instead of inv(A)*b is two to three times as fast and produces residuals on the order of machine accuracy, relative to the magnitude of the data.
赞一下 回复此楼
11楼2016-03-31 21:37:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mathstudy 的主题更新
111/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 324分 +8 llllkkkhh 2026-03-18 8/400 2026-03-18 21:01 by Catalysis25
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 10/500 2026-03-18 20:59 by 脱颖而出
[考研] 0703化学调剂，求各位老师收留 +9 秋有木北 2026-03-14 9/450 2026-03-18 19:49 by macy2011
[考研] 一志愿吉林大学材料学硕321求调剂 +3 Ymlll 2026-03-18 5/250 2026-03-18 19:32 by Ymlll
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 085601材料工程专硕求调剂 +6 慕寒mio 2026-03-16 6/300 2026-03-18 14:26 by 007_lilei
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
关闭
请填处理意见
关闭