麦克雷 Mavom.cn

标题: 【求助】解线性方程组 [打印本页]

作者: Bbephenzese    时间: 5 小时前
标题: 【求助】解线性方程组
两个矩阵如下,其中第一个是多层矩阵:

(, 下载次数: 0)

(, 下载次数: 0)
我是想解一个方程组Ax=b,其实x的元素一共有450个,相当于第一个eeta矩阵有450层,但x(1)~x(15),x(88)~x(450)已经事先求出,所以我要求的是x(16)~x(87)。如图所示,第二个number矩阵里的元素其实可以说是x的元素的下标,相当于方程组有16~87共72行,方程组的系数是eeta矩阵中每层第2行,相当于是:
①:eeta(2,1,16)*x(16)+eeta(2,2,16)*x(17)+eeta(2,3,16)*x(18)+eeta(2,4,16)*x(86)+...+eeta(2,21,16)*x(84)=0
②:eeta(2,1,17)*x(17)+eeta(2,2,17)*x(16)+eeta(2,3,17)*x(18)+eeta(2,4,17)*x(19)+...+eeta(2,21,17)*x(124)=0
③:eeta(2,1,18)*x(18)+eeta(2,2,18)*x(19)+eeta(2,3,18)*x(17)+eeta(2,4,18)*x(16)+...+eeta(2,21,18)*x(25)=0
......
由于x(1)~x(15),x(88)~x(450)已经求出来,所以只要把含有x(1)~x(15),x(88)~x(450)的项移项到等号右边,就构成了向量b,也就是说,x=[x(16),x(17),x(18),...x(87)]T。我现在绞尽脑汁也想不出代码如何写,有谁能提供下代码?在下无比感谢。
作者: ertwrgtewrt    时间: 5 小时前
这活还没人接??




欢迎光临 麦克雷 Mavom.cn (https://mavom.cn/) Powered by Discuz! X3.5