matlab中怎么追赶法求解矩阵
发布网友
发布时间:2022-05-05 23:28
我来回答
共1个回答
热心网友
时间:2022-06-28 07:31
matlab中怎么追赶法求解矩阵?追赶法是求解三角矩阵方程的一种方法,其特点是计算量少,方法简单,算法比较稳定。
追赶法求解的思路是:
1、预处理,生成方程组的系数u(i)及其除数d(i)
2、追,从y(1)y(2)y(3)。。。y(n)一直往前追
3、赶,从x(n)x(n-1)。。。x(1)一直往后赶
按这思路,就可编写chase()函数。其格式
x=chase(a,b,c,f)
求解例题:
┌ ┐┌ ┐ ┌ ┐
│2 -1 0 0││x1│ │6│
│-1 3 -2 0││x2│=│1│
│0 -1 -2 -1││x3│ │0│
│0 0 -3 5││x4│ │1│
└ ┘└ ┘ └ ┘
a=[0,-1,-1,-3];
b=[2,3,2,5];
c=[-1,-2,-1,0];
f=[6,1,0,1];
x=chase(a,b,c,f)
运行结果为 x=(5,4,23,2)^T