
上面图片中是线性方程组的样子,以及一些条件
,求解上面的线性方程组。A是非方阵,并且里面含有待求的参数a

关注P_n-1,n=1,2,...15,那P_0是啥?
这是用matlab算的,我令n=1的时候,a=P1了,还有为什么b的元素全都是-aT_1呢?P的元素全部相同的时候是这样的,但其他的情况就不一样了吧,比如:
clc,clear,close all;
P = [2,4,2,8,4,2,4,2,2,4,8,4,2,8,2]';
a = 0; % 初始化
X = zeros(1,15)';
for i = 1:15
X(i) = sum(P(1:i));
if i>1
a = (P(i-1)-P(i))/P(i);
end
end
disp('X =')
disp(X)
disp('P =')
disp(P)
A = [-(2+a),1,zeros(1,13);
1,-(2+a),1,zeros(1,12);
0,1,-(2+a),1,zeros(1,11);
zeros(1,2),1,-(2+a),1,zeros(1,10);
zeros(1,3),1,-(2+a),1,zeros(1,9);
zeros(1,4),1,-(2+a),1,zeros(1,8);
zeros(1,5),1,-(2+a),1,zeros(1,7);
zeros(1,6),1,-(2+a),1,zeros(1,6);
zeros(1,7),1,-(2+a),1,zeros(1,5);
zeros(1,8),1,-(2+a),1,zeros(1,4);
zeros(1,9),1,-(2+a),1,zeros(1,3);
zeros(1,10),1,-(2+a),1,zeros(1,2);
zeros(1,11),1,-(2+a),1,zeros(1,1);
zeros(1,12),1,-(2+a),1;
]
b = A*X
disp('a =')
disp(a)
打印结果:
X =
2
6
8
16
20
22
26
28
30
34
42
46
48
56
58
P =
2
4
2
8
4
2
4
2
2
4
8
4
2
8
2
A =
-5 1 0 0 0 0 0 0 0 0 0 0 0 0 0
1 -5 1 0 0 0 0 0 0 0 0 0 0 0 0
0 1 -5 1 0 0 0 0 0 0 0 0 0 0 0
0 0 1 -5 1 0 0 0 0 0 0 0 0 0 0
0 0 0 1 -5 1 0 0 0 0 0 0 0 0 0
0 0 0 0 1 -5 1 0 0 0 0 0 0 0 0
0 0 0 0 0 1 -5 1 0 0 0 0 0 0 0
0 0 0 0 0 0 1 -5 1 0 0 0 0 0 0
0 0 0 0 0 0 0 1 -5 1 0 0 0 0 0
0 0 0 0 0 0 0 0 1 -5 1 0 0 0 0
0 0 0 0 0 0 0 0 0 1 -5 1 0 0 0
0 0 0 0 0 0 0 0 0 0 1 -5 1 0 0
0 0 0 0 0 0 0 0 0 0 0 1 -5 1 0
0 0 0 0 0 0 0 0 0 0 0 0 1 -5 1
b =
-4
-20
-18
-52
-62
-64
-80
-84
-88
-98
-130
-140
-138
-174
a =
3