范文一:利用Matlab求线性方程组的通解
篡;Ⅵ渊-酝
利用Matlab求线性方程组的通解
简绍勇杜玲陈勇
(新余高等专科学校教学与信息科学系江西新余338031)
[摘要]讨论线性方程组解的个数及求线性方程组的通解问题是线性代数中的常见问题,介绍利用Ilatlab软件解决这两种问题的方法.[关键词]Matlab线性方程组基础解系通解中围分类号:01-0
文献标识码:A
文章编号:1671--7597(2008)12201卯一01
一、引一
线性代数是数学中的一个重要分支.很多理论问题和实际问题都需要借助于线性代数的理论工具来分析解决。学习线性代数有两大难点:一是概念、理论抽象,二是计算量大。由于这两个难点,初学者往往很难掌握好线性代数的知识理论。若能掌握好用于解数学问题的lIatlab软件,则能轻松快捷的解决很多线性代数问题。
(一)线性方程组有关定理
1
1
毛+毛+而+毛+毛=7
例3求线性方程组3而+2毛+气+‘一3而=—2的通解
毛+勰+2‘+6毛=23
5而+4b+3而+3丘一毛=12
解:
1:3
(1)先判断线性方程组解的情况.输入并运行以下代码:a-[1
2
l
l
l-3;0
1226:543
3-1]:b=[7—2
23
12]’;
对线性方程组Ax=6,其中A=(嘞)。。,x=(耳,而。…,矗y,6=(6I。62。…。k)’
定理l:
(1)若,(一)≠r(A16).则线性方程组无解:
(2)若,(一)=r(A16)=一,则线性方程组存在唯一解;(3)若,(一)=r(A16)<疗,则线性方程组有无穷多个解.
定理2:对线性方程组的导出组AX=b,若,(爿)=,<疗.则导出组的
rl=rank(a),r2=rank(【ab1)运行结果为rl=r2=2<5,说明该线性方程组有无穷多个解。(2)输入以下代码求出线性方程组的一个特解x0=a\b运行结果为:x0=(3.1667,0,0,0.3.8333)’(3)输入代码求出导出组“=o的一个基础解系xx=null(B)运行结果是:
Xx=
0.7530
-0.4167_o.3043-t).3043
0.0176
-0.7464
-t}.0000
基础解系含有刀一,个解向量,(1)若舌,参,…,靠,为导出组的基础解系。
-0.0000
-0.70710.707l--0.0000
则毛毒十如最+…+ko就是导出组朋=0的全部解,也称为通解:
0.45330.4533
-0.1778
(2)若磊是从=6的一个特解,则磊+与毒+屯磊+…+屯一,o,就是AX
=b的通解。
二、曩性方程组有关向■的M玳llab解佳(一)有关№tIab命令
讨论线性方程组解的问题的有关Matlab命令见下表
命令
det
0.2723
(4)线性方程组的通解为
3.1667O
X=
OO3.8333
+c1
0.7530-0.4167-03043-O.30430.2723
+c2
O.0176-41.74640.45330.4533-0.1778
O0
+岛-4).7071
0.707l
0
含义
方阵的行列式方阵的逆矩阵的秩
矩辟的行最简化形求基础解系左除
rank
rref
其中cI,c2,c,为任意常数.
nullt
(=)利甩左除命令给出线性方程组的特懈
对线性方程组允r=6.在Matlab软件里常用矩阵的左除命令。A\b”求出该线性方程组的一个解。若^为方阵,则A\b和inv(^)柚基本一致;若^不为方阵。A\b命令{鲫latlab软件自动选择适当的方法来求解.
若,IX=6无解,则命令“A\b”将给出一个最小二乘意义上的近似解,即使AZ—b的长度达到最小。若允r=b有无穷解,则命令。A\b”将给出一个具有最多零元素的特解。若AX=b有唯一解。则命令“^\b”给出这个唯一解。
(三)讨论线性方程组是否有解的MatIab方法
对线性方程组—tⅣ=b,在Matlab软件里先给系数矩阵^和常数列矩阵b赋值,然后利用命令“rank(A),rank(A。b)”求出系数矩阵和增广矩阵的秩.再根据定理(1)所介绍的数学理论进行判断.
(四)求线性方程组的通解的mtab方法
为求线性方程组AX=b的通解,需先判断方程组是否有解。若有解则用。A\b”求出一个特解。再用命令“null(A)。求出导出组AX=0的一个基础解系得出其通解,然后利用定理(2)的结论给出丘r=6的通解.
作者简介:
2006.
参考文献:
[1】胡良剑、孙晓君。(1latlab数学实验,[1I】,北京:高等教育出版社.
[2]赵静、但琦.出版社,2003.
[3]同济大学应用数学系,出版社.2004.
‘教学建模与数学实验’(第2版)[岫。北京:高等教育
‘高等教学'(第5版)[M].北京t高等教育
简绍勇,男,江西新余人.新余高等专科学校教学与信息科学系助教.
固
万方数据
利用Matlab求线性方程组的通解
作者:作者单位:刊名:英文刊名:年,卷(期):被引用次数:
简绍勇, 杜玲, 陈勇
新余高等专科学校数学与信息科学系,江西,新余,338031硅谷
SILICON VALLEY2008(24)1次
参考文献(3条)
1.胡良剑;孙晓君 Matlab教学实验 20062.赵静;但琦 数学建模与数学实验 20033.同济大学应用数学系 高等数学 2004
本文读者也读过(10条)
1. 姬五胜.Ji Wusheng 线性方程组解法及其MATLAB实践[期刊论文]-天水师范学院学报2009,29(2)2. 温旭东 MATLAB与Mathmatica在线性代数运算中的对比分析[期刊论文]-黑龙江科技信息2008(3)3. 晏林.YAN Lin 多项式环上一次不定方程组的矩阵解法与程序设计[期刊论文]-科技通报2005,21(4)4. 吴静杰.Wu Jingjie 引入Matlab辅助线性代数研究[期刊论文]-太原师范学院学报(自然科学版)2009,8(4)5. 任沂军.REN Yijun 矩阵运算的几种算法实现探讨[期刊论文]-测绘信息与工程2006,31(6)
6. 王春利.黄坚.丁少玲.余文质.袁媛.WANG Chun-li.HUANG Jian.DING Shao-ling.YU Wen-zhi.YUAN Yuan 线性代数计算软件辅助教学的探索与研究[期刊论文]-长春师范学院学报(自然科学版)2010,29(5)7. 朱熙湖 矩阵基本运算的数学实验及思考[期刊论文]-中山大学学报论丛2006,26(6)
8. 赵秉新.郑来运 MATLAB在求解线性方程组中的多种应用[期刊论文]-通化师范学院学报2007,28(12)9. 王力梅 浅谈线性代数教学[期刊论文]-甘肃科技2010,26(13)
10. 陈永胜.刘洋萍.Chen Yongshen.Liu Yangping 基于MATLAB求解非齐次线性方程组[期刊论文]-赤峰学院学报(自然科学版)2009,25(10)
引证文献(1条)
1.谭超.张开仕.贺革.任清华 复杂化学反应体系独立反应数的确定[期刊论文]-宜宾学院学报 2010(6)
本文链接:http://d.wanfangdata.com.cn/Periodical_guig200824128.aspx
范文二:利用Matlab求线性方程组的通解
教育 科学
从摆箭
利用
求线 性 方 程 组 的 通 解
简绍 勇 杜
玲 陈 勇
江西
新余 高等专科 学 校 数 学与信息科学 系
新余
〔 摘
要 〕 论 线性 方 程 组 解 的 个数 及 求 线性 方 程 组 的通 解 问题 是线性代数 中 的常 见 问题 讨 线性方 程 组 基础解系
人
,
介绍 利用
软 件解 决这 两 种 问题 的方 法
。
关健 词 」
通解
中图分 类号
一
文 献标识 码
文
编号
一
了
一
一
、
引,
,
三 很 多理论 问题 和 实 际 问题都需 要 借 学 习 线 性 代 数有 两 大难 点
,
。
、
应 用 实侧
不十毛
毛
十毛
线性代数 是 数 学 中的 一 个重要分支
助 于 线 性 代 数 的理 论 工 具 来 分 析 解 决
念
、
凡
凡
一
凡 凡
一是概
例
求 线性方程组
戈
欢 仁
一 的通解
理 论抽 象
,
二 是 计 算 量大
。 。
。
由 于这两 个 难 点
初学者 往 往很难掌握 好 线 软件
,
毛
石
毛 凡
凡
一
性代数 的知 识 理 论
一
若 能 掌握 好 用 于 解 数 学 问题 的
则能轻松快 解
一
凡
,
捷 的解 决很 多线 性代 数 问题
先判断线性方程组解 的情 况
,
输 入 并 运 行 以 下代 码
一
【
线性方 程 组 有关定理 对 线性 方程 组 其中
…
一
二
〕
二
’
解
,
马
,
,
认
,
凡 … 凡
, ,
, ’
, 气
行 结果 为
,
,
运行 结果为
。
二
,
说 明该线 性方程 运
,
气’
定理
组 有 无 穷多个 解 若 若叹 若
笋
,
输 入 以下 代 码求 出 线性 方程组 的 一 个 特 解
,
则 线性 方程 组无解
。
,
,
,
’
输 入 代 码 求 出 导 出组
则线 性方 程 组 存在唯 一 解 则线 性方程 组 有无穷多个解
,
习的一 个基 础解 系
二
,
二
运行 结果是
,
定理
对 线 性 方程 组 的 导 出 组 注万
一
二
若叹
,
,
则 导 出组 的
,
基础 解 系含 有 则 气 点 气 最 若
二
、
个 解 向量
十
,
帐氛
,
若 点橇 … 就 是 导 出组 月尤
, ,
,
的基 弘 为导 出组 也 础解系 的全 部解 称 为通 解
二
一
司
一
一
氛是
。
的一个 特 解
则 磊十 气 十 气 或 么
嵘豪
,
,
就是 盯
的通解
妞性 方租 组 有关间 有 关物
的
州
解法
命令 见 下 表
一
刃 线性 方程组 的通解为
、 、 ‘ 汽 了门
一
命令
一 内 ,
、
、
产 ‘ 、
讨论线性 方程 组 解 的问题 的有 关物
命令 含义
方 阵 的行列 式 方 阵 的逆 逆 矩 阵 的秩 矩 阵 的行最简 化形 求基 础解系 系 左除
刁 刁 刁
,
…
。
其中
几 马 为 任 愈常数
,
二 利用 左 除 命令给出线 性方 程组 的 特解 对 线性 方程 组 月尤 二 在 软件 里 常用 矩 阵的左 除 命令
,
“
、
”
求 出该 线性 方程 组 的一 个解
不 为方 阵 即 使 注万 唯 一解
。
。
若 为方 阵
“
,
则
和
基 本 一致
。
若
,
,
人
命 令使 无解
,
软件 自动 选择 适当的方法 来求解
人
”
。
若 月尤
一
则命 令
。
将给 出一 个最小 二 乘意义上 的近似解 有无 穷解
有 唯 一解
, ,
的长 度 达 到 最 小
若 注万
则命 令
“
“
、
”
”
将给 出
参考文献
胡 良剑 赵静
出 版社 出版 社
,
、 、
一 个 具 有最 多零元 素 的特 解
若 刁尤
则命令
、
给 出这 个
孙 晓君
,
,
《地
数 学 实验 》 〔 〕 北 京
,
高等 教 育 出 版 社
,
,
三 讨 论线性 方程 组 是否 有解 的山 方法 对 线 性 方 程 组 几丫 在 软件 里 先给系数矩 阵 和 常数 列 矩 阵
,
但琦
《数 学 建模 与 数 学 实验 》
,
第 版
〔〕 北 京
高 等 教育 高等 教育
赋值
,
然后 利用 命 令
,
“
,
,
”
求 出 系数矩 阵和 增广矩
阵 的秩
再 根据 定 理
所 介 绍 的数学理论进行判 断
方法
,
同 济 大 学 应用 数 学 系
,
《高等 数 学 》
第 版
〔
,
北京
四 求 线 性方程 组 的 通解 的 为 求 线 性方程 组 注矛 的通 解
需先 判 断方程组 是 否 有解
“ ”
。
若有解则 的一
。
用
“
”
求 出一 个 特 解
,
,
再 用 命令
求 出 导 出组 月尤
书
作者简介 简绍 勇
,
个 基 础 解系得 出其通解
然后利用定理
的结论给 出 刀万
的 通解
男
,
江西 新余人
,
新余高等 专科 学 校 数学与信 息科学系 助 教
。
口 巧。 口
范文三:【doc】利用Matlab求线性方程组的通解
利用Matlab求线性方程组的通解
-??
?Z?,Matlab?????{???1?????r??
??????????????
(?????????5?????n?????]??????????????????338031)
[????]???????{???1???????m?????????{???1?????r???????O???{?????}???????????)?????Z?,Ilatlab???????????^?????????????]
[???_??]Matlab???{???1???????????r??
?}???b?????r01-0???_?(???-?rA?????x???r1671--7597(2008)12201????01 ??c-?2??
???{?????O?????}?????m?????b?O?]???r?b???????????????????????d ???????{???????b?????????b??????cM???????{?????8?^?{?????r???O?? ??c-?b???????)???O???????{cM?'?????^?m?????)?+???l???????????????? ?{???????????b??cM???????????,??????????????lIatlab?????)?k???????? ???????????r???{????????cM
(??)???{???1???8?????b
?????{???1??Ax=6?)???}A=(?o)cMcM?)x=(???)?)cMf??)??y?)6=(6IcM62cM f?cMk)ek
???bl?r(1)???)(??)f?r(A16)?]?k???{???1???????r
(2)???)(??)=r(A16)=???)?k???{???1?????,?O?????<
(3)???)(??)=r(A16)<???)?k???{???1???8?????r?m???]
???b2?r?????{???1??????????AX=b?)???)(??)=?)<???]?k???????? ?????????4?8?????)?m???6???)(1)???)?)???)f??)???)??????????????????cM ?k??????????+f?+ko?'?O???????1=0?????????)???????r???r (2)?????O??=6?????m?????)?k??+?]??+????+f?+?????)o?)?'?OAX =b???r??cM
??c-???{???1???8???6i???M??llab????
(??)?8??nxtIab????
???????{???1?????????????8??Matlab???????1??
?????4??
det?????????=?I
???????b
rank????????
rref?????????????x??
null??????????
t????
(=)?Z?Z???????????????{???1????????
?????{???1????r=6?]?,Matlab?????????,??????????????cMA??be+ ?????????{???1???????m??cM??^???????)?kA??b??inv(^)?????B?????<?? ^?[??????cMA??b????{??latlab?????????8???M???????????????] ??,IX=6?????)?k????eXA??be+?b???????m???I?????I???0?????????) ????AZe9b?????????c???IcM????r=b?8???????)?k????cMA??be+?b???? ???m???8???r?????e??????cM??AX=b?8?O????cM?k????eX^??be+???????m ?O????cM
(?8)???????{???1???O?v?8????MatIab????
?????{???1??e9tn?=b?)?,Matlab??????????????????^???????=???? b?????)???]?Z?,????eXrank(A)?)rank(AcMb)e+???????????????????? ???????]?????????b(1)?????????????b???????^???]
(??)?????{???1?????r????mtab????
???????{?1??AX=b???r???)?????^???1???O?v?8??cM???8???k ?,cMA??be+???????m????cM???,????eXnull(A)cM??????????AX=0???? ?m???????????????r???)???]?Z?,???b(2)???????????yr=6???r???] ??
??3?????{???1??
??+??+?)+??+??=7
3?)+2??+?9+ey??3?)=e92???r??
??+??+2ey+6??=23
5?)+4b+3?)+3?y????=12
???r(1)???^?????{???1???????????]?????????????1???-?ra-[1l
111?r32ll-3?<01226?r5433-1]?rb=[7e922312]ek?< rl=rank(a)?)r2=rank(cAab1)??????????rl=r2=2<5?)???]?????{???1
???8?????r?m??cM(2)???????1???-???????{???1???????m????x0=a??b??
?????????rx0=(3?]1667?)0?)0?)0?]3?]8333)ek(3)???????-??????????
eX=o?????m????????xx=null(B)?????????O?r Xx=
0?]75300?]0176-t}?]0000
-0?]4167-0?]7464-0?]0000
_o?]30430?]4533-0?]7071
-t)?]30430?]45330?]707l
0?]2723-0?]1778--0?]0000
(4)???{???1?????r????
X=
3?]1667
O
O
O
3?]8333
+c1
0?]7530
-0?]4167
-03043
-O?]3043
0?]2723
???}cI?)c2?)c?)?????I?????]
+c2
O?]0176
-41?]7464
0?]4533
0?]4533
-0?]1778
+??
O
0
-4)?]7071
0?]707l
0
???a???_?r
[1c6??????c-?????<cM(1latlab?????????)[1Ic6?)?=???r???????????????] 2006?]
[2]????c-???v?]ey???????n?]????????ek(?B2??)[??cM?=???r???????? ???????)2003?]
[3]?)???{?????,???????)ey????????'(?B5??)[M]?]?=??t???????? ???????]2004?]
???l?????r
???????)???)???????????]?????????5?????n?????]???????????????]
范文四:求线性方程组的解
求线性方程组的解
高斯消去法是求解线性方程组最学用的方法之一,下面用一个数值例题来说明高斯消去法的计算原理(
求解下列线性方程组的解:
运行教材光盘中VB程序得:
即x=1, x=2, x=-1 123
要求:会应用光盘中VB程序求解即可。
参见光盘中VB程序如下:
高斯消元法解线性方程组
'方程组系数输入窗体
'本窗体负责完成方程组AX=b的系数输入任务,计算任务全部由模块mduCalc完成 '文件 :frmInput.frm
'1999/11/05
'Using Microsoft Visual Basic 6.0
1
'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
'控件说明
'grdA:输入系数矩阵,,由用户输入
'grdB:输入,矩阵,由用户输入
'grdX:显示计算结果,来自mduCalc.Calc过程
''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
Private iCol As Integer '当前所选择的列
Private iRow As Integer '当前所选择的行
Private cGrd As String '当前所选择的矩阵(A或b)
Private Sub cmdCalc_Click()
'使用高斯列主元素消去法解方程组
Dim N As Integer '方程数目
N = CInt(txtN.Text)
ReDim a(N, N + 1) As Double '系数矩阵
Dim i, j As Integer
'设置系数矩阵A(N, N+1)的前N列
For i = 1 To grdA.Rows - 1
For j = 1 To grdA.Cols - 1
If grdA.TextMatrix(i, j) = "" Then
a(i, j) = 0#
Else
a(i, j) = CDbl(grdA.TextMatrix(i, j))
End If
Next j
Next i
'设置第N+1列
For i = 1 To grdB.Rows - 1
If (grdB.TextMatrix(i, 0)) = "" Then
a(i, N + 1) = 0#
Else
a(i, N + 1) = CDbl(grdB.TextMatrix(i, 0))
End If
Next i
'用列主元素消去法计算X
ReDim x(N) As Double '保存计算结果向量
mduCalc.Calc a, x, N
2
'显示计算结果
grdX.ColAlignment(0) = flexAlignLeftTop
For i = 1 To grdX.Rows - 1
grdX.TextMatrix(i, 1) = Format(CStr(x(i)), "0.0000E-00")
Next i
End Sub
Private Sub cmdExit_Click()
Unload Me
End Sub
Private Sub Form_Load()
grdA.ColWidth(0) = TextWidth("10") * 2
grdX.ColWidth(0) = TextWidth("X1") * 2
grdX.ColWidth(1) = grdX.Width - grdX.ColWidth(0) - 100
Call txtN_Change
End Sub
Private Sub grdA_EnterCell()
cGrd = "A"
txtValue.Visible = True
iCol = grdA.Col
iRow = grdA.Row
txtValue.Width = grdA.ColWidth(iCol) - 15
txtValue.Height = grdA.RowHeight(iRow) - 15
txtValue.Text = CStr(grdA.TextMatrix(grdA.Row, grdA.Col))
txtValue.Left = grdA.Left + grdA.ColPos(iCol) + 45
txtValue.Top = grdA.Top + grdA.RowPos(iRow) + 45
txtValue.SetFocus
End Sub
Private Sub grdA_GotFocus()
Call grdA_EnterCell
End Sub
Private Sub grdA_LeaveCell()
cGrd = "NULL"
txtValue.Visible = False
End Sub
Private Sub grdB_EnterCell()
cGrd = "B"
3
txtValue.Visible = True
iCol = grdB.Col
iRow = grdB.Row
txtValue.Width = grdB.ColWidth(iCol) - 15
txtValue.Height = grdB.RowHeight(iRow) - 15
txtValue.Text = CStr(grdB.TextMatrix(grdB.Row, grdB.Col))
txtValue.Left = grdB.Left + grdB.ColPos(iCol) + 45
txtValue.Top = grdB.Top + grdB.RowPos(iRow) + 45
txtValue.SetFocus
End Sub
Private Sub grdB_GotFocus()
Call grdB_EnterCell
End Sub
Private Sub grdB_LeaveCell()
cGrd = "NULL"
txtValue.Visible = False
End Sub
Private Sub grdX_EnterCell()
cGrd = "NULL"
txtValue.Visible = True
iCol = grdX.Col
iRow = grdX.Row
txtValue.Width = grdX.ColWidth(iCol) - 15
txtValue.Height = grdX.RowHeight(iRow) - 15
txtValue.Text = CStr(grdX.TextMatrix(grdX.Row, grdX.Col))
txtValue.Left = grdX.Left + grdX.ColPos(iCol) + 45
txtValue.Top = grdX.Top + grdX.RowPos(iRow) + 45
txtValue.SetFocus
End Sub
Private Sub grdX_LeaveCell()
txtValue.Visible = False
End Sub
Private Sub txtN_Change()
'重新设置Grid控件格式
Dim i, j As Integer
'设置系数矩阵尺寸
grdA.Cols = CLng(txtN.Text) + 1
4
grdA.Rows = CLng(txtN.Text) + 1
grdB.Rows = CLng(txtN.Text) + 1
grdX.Rows = CLng(txtN.Text) + 1
'设置标题行,列
For i = 1 To grdA.Cols - 1
grdA.TextMatrix(0, i) = "X" & CStr(i)
grdA.ColWidth(i) = TextWidth("X1") * 2
Next i
For i = 1 To grdA.Rows - 1
grdA.TextMatrix(i, 0) = CStr(i)
Next i
For i = 1 To grdX.Rows - 1
grdX.TextMatrix(i, 0) = CStr(i)
Next i
End Sub
Private Sub txtValue_Change()
Select Case cGrd
Case "A"
grdA.TextMatrix(iRow, iCol) = txtValue.Text
Case "B"
grdB.TextMatrix(iRow, iCol) = txtValue.Text
End Select
End Sub
5
范文五:线性方程组的通解是否唯一吗?
考研屋 www.kaoyanwu.com
提供各大机构考研、公务员、四六级辅导视频课程
我遇到这样的情况不只一次了,这次又看到,我的答案和书上的答案不一样,看了书上的答案,觉得也对,而我的思路也和书上的差不多,我觉得我的也是正确的(是不是线性方程组的通解不唯一啊,
表达式的前面有一个任意常数“k“,当然就不唯一了呀 [s:2]
通解本来就表示一系列解,基础解系的选择不同,所得出的通解表示形式就不同,但不同基础解系之间可以互相转化。
但我看的不是成倍数等关系,能够看出是可以转化的呢([ 本帖最后由 sweetliwei 于2007-6-14 10:10 PM 编辑 ]
帮我看看对不对啊,
回复 #4 sweetliwei 的帖子由于他的秩是2,是4阶的,所以有两个通解,或许你可以得出其他的解,但那些解一定可以由答案中求的那两个通解线性表出
楼上的,你把概念弄错了吧(不是有两个通解,是通解中两个线性无关的向量
回复 #7 sweetliwei 的帖子是的,好象正规说法是什么 两个线性无关的基础解系,好久没看线代了,那些叫法记不清了
楼主的答案是对的,书上的答案错了。这种题如果给了具体的方程,那你可以把解出来的通解代入原方程进行检验。如果没有具体的方程,那就只有对过程进行检验了。检验两组同解是否可以相互表示的办法之一就是检验一下一组通解中的每个线性无关的向量都可以由另一组通解线性表示。这道题的两个答案只有一个线性无关的向量不一样,其他的都一样,很明显是有一组错了。经仔细查验,是书上的答案错了。
怎么检查的啊,
怎么没人回呢各位帮帮忙哈
书上的第二个解向量是这样得来的2(2a2-a3)-(a1+a2)=2(a2-a3)+(a2-a1)=(2,-1,-40,48)T
我上面的说法错了,你可以检验一下Ax=0的三个解,即(1,1,-5,6)\’,(3,0,-45,54)\’和(2,-1,-40,48)\’这三个解应该是线性相关的,也就是说是可以相互表示的。或者说(2,-1,-40,48)\’是可以由另外两个向量线性表示的,或者说(3,0,-45,54)\’是可以由另外两个向量线性表示的。这样就可以说明两种答案实质上是一样的。另外:这道题中Ax=b的三个解是可以算出来的,α1,(1,1,1,-1)\’,α2,(1,0,-9,11)\’,α3,(0,0,6,-7)\’。你也可以通过这个验证一下。再另外:你的计算明显比书上的答案要简单,所以数上的解法不如你的解法。 专业提供提供各大机构考研、公务员、四六级辅导视频课
考研屋 www.kaoyanwu.com
提供各大机构考研、公务员、四六级辅导视频课程
其实如果按照先算出α1,α2,α3,再计算Ax=0的通解可能更快速一些。虽然计算量稍大一些,但是凑的方法不是很容易看出来的。
如果先把a1,a2,a3求出来,那书上的答案就有点问题了,
a1-a2=(1,0,10,-12)\’a2-a3=(1,0,-15,18)\’a1-a3=(1,1,-5,6)\’后两项是我的答案(与之等价的应该是第一项,但书上的答案(,,,,,,,,,,,),与(1,0,10,-12)\’有出入啊(欢迎来到免费考研网www.freekaoyan.com
大概有点明白了(但考试时,答案就不一定了,这怎是好,
2007-6-16 10:21 AM 发表 如果先把a1,a2,a3求出来,那书上的答案 原帖由 sweetliwei 于
就有点问题了,a1-a2=(1,0,10,-12)\’a2-a3=(1,0,-15,18)\’a1-a3=(1,1,-5,6)\’后两项是我的答案(与之等价的应该是第一项,但书上的答案(,,,,,,,,,,,), ...
你前面都写了2(2a2-a3)-(a1+a2)=2(a2-a3)+(a2-a1)=(2,-1,-40,48)T怎么可能等于a1-a2呢,
原帖由 sweetliwei 于 2007-6-16 10:29 AM 发表 大概有点明白了(但考试时,答案就不一定了,这怎是好,
一般来说考试不会出这种题的。
考研屋www.kaoyanwu.com :提供各大机构考研、公务员、四六级辅导视频课程 专业提供提供各大机构考研、公务员、四六级辅导视频课
转载请注明出处范文大全网 » 利用Matlab求线性方程组