矩阵用matlab求解线性方程组,有限差分计算中任意矩阵线性方程组的Matlab求解

我会用

bsxfun

(在最近的Matlab版本中,这个函数可能被实现到解释器中,我不确定):

p = 4;

Q = 8;

dx = 1;

h = -p:p*dx

Qvector = [Q,1:Q-1]'

Matrix = bsxfun(@(Qvector, h)h.^(Qvector)./factorial(Qvector), Qvector, h)

输出:

h =

-4 -3 -2 -1 0 1 2 3 4

Qvector =

8

1

2

3

4

5

6

7

Matrix =

1.6254 0.1627 0.0063 0.0000 0 0.0000 0.0063 0.1627 1.6254

-4.0000 -3.0000 -2.0000 -1.0000 0 1.0000 2.0000 3.0000 4.0000

8.0000 4.5000 2.0000 0.5000 0 0.5000 2.0000 4.5000 8.0000

-10.6667 -4.5000 -1.3333 -0.1667 0 0.1667 1.3333 4.5000 10.6667

10.6667 3.3750 0.6667 0.0417 0 0.0417 0.6667 3.3750 10.6667

-8.5333 -2.0250 -0.2667 -0.0083 0 0.0083 0.2667 2.0250 8.5333

5.6889 1.0125 0.0889 0.0014 0 0.0014 0.0889 1.0125 5.6889

-3.2508 -0.4339 -0.0254 -0.0002 0 0.0002 0.0254 0.4339 3.2508