找回密码
 注册
关于网站域名变更的通知
查看: 463|回复: 1
打印 上一主题 下一主题

方程AX=b的解的讨论(特解、通解、零空间向量等概念)及其MATLAB实现

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-5-28 10:49 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
求矩阵形式线代方程组,讨论AX=b的解是最基本的一项内容。
+ v8 `% z9 ~  Z  S) D
; N9 G& x5 K" _- @  eAX=b的解 = 特解 + 矩阵零空间向量
# h& \$ @& d0 l7 k: M* n
0 N, Q0 R5 P/ L8 Y. k+ r% A特解:AX=b的自由变量都=0时x的解。3 r! N% V; n+ z" z4 g! J) Q% N

4 ~; o9 ^5 F. v  y& M矩阵零空间向量:AX=0时x的解空间。矩阵零空间向量又牵扯到了零空间的概念,就不赘述了。我们可以简单记为:+ m* R8 j# L$ f

' d; P2 i# ]4 n5 ]. M& WX = X* + + H- P: }# U  t! J- Z# E
- c( n6 f: b2 L. u% r

* s' Y+ K7 b; I/ Y
9 u8 n# _' o' q( z: J1 o3 L4 p- o零空间向量:+ X5 X4 m8 Q8 h, ?! T8 P
# S7 M) y4 i, H9 `& ]4 i
/ R9 O  [1 j' u3 w4 v' t# C, p, G8 [2 t% J, O
; f0 }# P# V4 r' {$ m
关于可解性:( |  l. p  v4 ^. v) @

, ^/ F2 y; n4 d. q' u. c: i* e
/ Y! u; ?0 ~, _$ V2 s* T/ d/ O. K% N0 M" m# u# v0 W! h
通解、特解:
5 t+ v/ C/ k0 f! { + F  ~6 j6 ?8 G" D: J" J

1 b- O5 j2 T9 T1 u2 }. Y6 g- X3 z/ Z6 s% @5 l5 r8 S2 G$ D, R
对上述例子,写了个简单的MATLAB程序,用以求AX=b的解。更全面的代码,可以参考文末的参考文献。
# i$ V! s4 E2 O  G3 j- i* L
( f) L5 [9 h( g  J8 M; ?# AA = [ 1 2 2 2;
, l$ J$ {1 e: s9 V    2 4 6 8;( V6 ^+ [; o4 B6 s$ r+ D
    3 6 8 10];' F) f  g& N$ r/ {8 z: ?/ J( m
b = [1;
5 ?& q+ T% m' x! W+ f, {     5;' [% y+ D) x% f5 T# a, O
     6];
& Y2 ^" X2 G" C  O2 z
  }: p2 q6 d# ?. `) F- gformat rat;
6 `  f# K4 T/ D' `syms n1 n2;
# x* c! z4 r' e% aX0 = A\b %零空间向量,即AX=0时X的解
1 n5 R7 _# L7 W" K0 wC = null(A,'r');
3 k" x- D' P7 |+ bX = C(:,1)*n1 + C(:,2)*n2 + X0  %X通解& W" |3 ^; t( D  J" G; E6 Y' A
/ ?6 s% a. N$ m3 Y7 G
  d/ P3 \- C6 f+ t9 o" n; ~
4 _( [( N2 ~2 N0 Q9 }  K) ~, Q
- V/ G2 s$ h, Y9 v" U9 f. n
) s; v9 Z/ V/ V( `( h8 x) i9 Y

该用户从未签到

2#
发表于 2020-5-28 13:16 | 只看该作者
方程AX=b的解的讨论(特解、通解、零空间向量等概念)及其MATLAB实现
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-24 11:40 , Processed in 0.156250 second(s), 26 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

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