- 注册时间
- 2022-8-7
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 105
- 在线时间
- 小时
|
发表于 2024-1-4 23:32:51
|
显示全部楼层
Sowya 提供了解线性方程组的求解函数 solve().
运行 Sowya.exe, 在提示符 >> 后键入矩阵 A 和 b, 然后输入 solve(A*x==b) 回车.
>> A=[1,1,-5,-1,1;8,5,-1,0,4;3,-7,-5,-1,1;2,1,8,3,-1;0,6,2,7,-2]
input> [1,1,-5,-1,1;8,5,-1,0,4;3,-7,-5,-1,1;2,1,8,3,-1;0,6,2,7,-2]
--------------------
1 1 -5 -1 1
8 5 -1 0 4
3 -7 -5 -1 1
2 1 8 3 -1
0 6 2 7 -2
--------------------
>> b=[9;2;4;-8;15]
input> [9;2;4;-8;15]
--------------------
9
2
4
-8
15
--------------------
>> solve(A*x==b)
>> 1 1 -5 -1 1 9
8 5 -1 0 4 2
3 -7 -5 -1 1 4
2 1 8 3 -1 -8
0 6 2 7 -2 15
The solution is:
x= E_0 + C_1*E_1 + C_2*E_2 + C_3*E_3 + C_4*E_4 + C_5*E_5
where E_0 is the special solution and the others form the base of the solution:
---------------
E_0 E_1 E_2 E_3 E_4 E_5
------------------------
2509|1670 0 0 0 0 0
1671|1670 0 0 0 0 0
-3853|1670 0 0 0 0 0
118|167 0 0 0 0 0
-68009|15698 0 0 0 0 0
------------------------
|
|