解决一个平方线性系统使用bicgstab使用默认设置，然后调整在溶液方法中使用迭代的公差和数量。

创建一个随机稀疏矩阵一种用50％的密度。还创建一个随机向量b对于右手边 $\mathrm{斧头}=\mathit{b}$。

rng默认的A = sprand（400400，0.5）;A = A'* A;B =兰特（400,1）;

解决 $\mathrm{斧头}=\mathit{b}$使用bicgstab。输出显示器包括相对剩余误差的值 $\frac{\Vert \mathit{b}-\mathrm{斧头}\Vert }{\Vert \mathit{b}\Vert }$。

x = bicgstab (A, b);

BICGSTAB停在迭代20而未收敛到所需的公差1E-06，因为达到迭代的最大数量。返回的迭代（数20）具有相对残留0.12。

默认情况下bicgstab使用20次迭代和公差1 e-6，该算法对该矩阵在20次迭代中无法收敛。由于剩余值仍然很大，这是需要更多迭代(或预调节矩阵)的一个好指标。您还可以使用更大的容差来简化算法的收敛。

解决再次用宽容的系统1的军医和100的迭代。

X = BICGSTAB（A，B，1E-4100）;

BICGSTAB停在迭代100而未收敛到所需的公差0.0001，因为达到迭代的最大数量。返回的迭代（号100）具有相对残留0.044。

即使有一个更宽松的容错和更多的迭代，残余误差并没有改善太多。当迭代算法以这种方式停滞时，这很好地表明需要一个预处理矩阵。

计算出的不完全Cholesky分解一种，并使用L '因子如预处理器输入到bicgstab。

L = ichol（A）;X = BICGSTAB（A，B，1E-4100，L'）;

bicgstab在30.5迭代收敛到相对残差为5.3e-05的解。

使用预调节器足够的问题,提高了数值属性bicgstab可以收敛。

检查使用预调剂矩阵的效果bicgstab解一个线性系统。

加载west0479，一个479乘479的实非对称稀疏矩阵。

加载west0479A = west0479;

确定b所以，真正的解决方案 $\mathrm{斧头}=\mathit{b}$是所有1的向量。

b =和(2);

设置宽容和最大迭代次数。

托尔= 1 e-12;麦克斯特= 20;

使用bicgstab发现在所要求的公差和迭代次数的解决方案。指定五个输出返回有关求解过程的信息：

[X，FL0，RR0，IT0，RV0] = BICGSTAB（A，B，TOL，麦克斯特）;FL0

fl0 = 1

rr0

rr0 = 1

IT0

it0 = 0

FL0是1因为bicgstab不收敛到要求的容忍1 e-12在请求的20次迭代中。事实上，的行为bicgstab是如此的可怜以至于最初的猜测x0 = 0(大小(2),1)是最好的解决方案，并返回，由所指示的it0 = 0。

为了帮助缓慢收敛，您可以指定一个预调节矩阵。自一种是不对称，利用ILU生成预条件 $\mathit{中号}=\mathit{大号}\mathit{ü}$。指定落忍忽略非对角项与值小于1 e-6。解决预处理系统 ${\mathit{一种}\mathit{中号}}^{-1}\mathit{中号}\mathit{X}=\mathit{b}$通过指定大号和ü作为输入,bicgstab。

设置=结构(“类型”，'ilutp'，“droptol”1 e-6);[L U] = ilu(一个,设置);(x1, fl1 rr1、it1 rv1] = bicgstab (A, b,托尔,麦克斯特,L, U);fl1

fl1 = 0

rr1

RR1 = 3.8661e-14

IT1

it1 = 3

an的用法ILU预调剂产生的相对残留小于规定的耐受性1 e-12在第三次迭代。输出RV1（1）是规范(b)，以及输出rv1(结束)是规范(b * x1)。

您可以按照进度bicgstab通过在每一次迭代绘制的相对残差。绘制每个解决方案的剩余历史与指定容忍线。

semilogy（0：长度（RV0）-1，RV0 /范数（b）中，'-o')举行在semilogy(0:长度(rv1) 1, rv1 /规范(b),'-o'）YLINE（TOL，'R--'）;传说（“不预条件”，“ILU预调节器”，'公差'，“位置”，“东”)包含(“迭代次数”) ylabel (“相对残存”）

检查供应的效果bicgstab对解决方案有初步的猜测。

创建一个三对角稀疏矩阵。用每一行的和作为的右边的向量 $\mathrm{斧头}=\mathit{b}$所以的期望解 $\mathit{X}$是那些的向量。

n = 900;e =的(n - 1);A = spdiags([e 2*e]，-1:1,n,n);b =和(2);

使用bicgstab解决 $\mathrm{斧头}=\mathit{b}$两次:一次是默认的初始猜测，另一次是对解决方案进行良好的初始猜测。对两种解决方案使用50次迭代和默认容忍度。金宝搏官方网站指定第二个解中的初始猜测为一个所有元素都等于的向量0.99。

麦克斯特= 50;x1 = bicgstab (A, b,[],麦克斯特);

BICGSTAB会聚在迭代20.5与相对残余9.3E-07的溶液。

X0 = 0.99 * E;X2 = BICGSTAB（A，B，[]，麦克斯特，[]，[]，X0）;

BICGSTAB会聚在迭代4与相对残余8.7E-07的溶液。

在这种情况下，提供一个初始猜测将启用bicgstab更快地收敛。

X0 =零（尺寸（A，2），1）;TOL = 1E-8;麦克斯特= 100;对于K = 1：4 [X，旗，relres] = BICGSTAB（A，B，TOL，麦克斯特，[]，[]，X0）;X（：，K）= X;R（k）的= relres;X0 = X;结束

通过提供解一个线性系统bicgstab带有一个可计算的函数句柄* x代替系数矩阵的一种。

其中一个威尔金森测试矩阵生成画廊是一个21通过-21三对角阵。预览矩阵。

A =画廊（“威尔克”, 21)

A =21×2110 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 8 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 7 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 6 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0⋮

威尔金森矩阵具有特殊的结构，这样你就可以显示操作* x带有函数句柄。什么时候一种乘以一个向量，大多数在所得矢量的元素是零。在结果所对应的非零元素具有非零三对角元素一种。而且，只有主对角线上有不等于1的非零。

表达式 $\mathrm{斧头}$变为：

$\mathrm{斧头}=\left[\begin{array}{cccccccccc}10& 1& 0& \cdots & & \cdots & & \cdots & 0& 0\\ 1& 9& 1& 0& & & & & & 0\\ 0& 1& 8& 1& 0& & & & & ⋮\\ ⋮& 0& 1& 7& 1& 0& & & & \\ & & 0& 1& 6& 1& 0& & & ⋮\\ ⋮& & & 0& 1& \mathrm{五}& 1& 0& & \\ & & & & 0& 1& 4& 1& 0& ⋮\\ ⋮& & & & & 0& 1& 3& \ddots & 0\\ 0& & & & & & 0& \ddots & \ddots & 1\\ 0& 0& \cdots & & \cdots & & \cdots & 0& 1& 10\end{array}\right]\left[\begin{array}{c}{\mathit{X}}_1\\ {\mathit{X}}_2\\ {\mathit{X}}_3\\ {\mathit{X}}_4\\ {\mathit{X}}_{\mathrm{五}}\\ ⋮\\ \\ ⋮\\ \\ {\mathit{X}}_{21}\end{array}\right]=\left[\begin{array}{c}10{\mathit{X}}_1+{\mathit{X}}_2\\ {\mathit{X}}_1+9{\mathit{X}}_2+{\mathit{X}}_3\\ {\mathit{X}}_2+8{\mathit{X}}_3+{\mathit{X}}_4\\ ⋮\\ {\mathit{X}}_{19}+9{\mathit{X}}_{20.}+{\mathit{X}}_{21}\\ {\mathit{X}}_{20.}+10{\mathit{X}}_{21}\end{array}\right]$。

将得到的载体可被写为三个矢量的总和：

$\mathrm{斧头}=\left[\begin{array}{c}0+10{\mathit{X}}_1+{\mathit{X}}_2\\ {\mathit{X}}_1+9{\mathit{X}}_2+{\mathit{X}}_3\\ {\mathit{X}}_2+8{\mathit{X}}_3+{\mathit{X}}_4\\ ⋮\\ {\mathit{X}}_{19}+9{\mathit{X}}_{20.}+{\mathit{X}}_{21}\\ {\mathit{X}}_{20.}+10{\mathit{X}}_{21}+0\end{array}\right]$= $\left[\begin{array}{c}0\\ {\mathit{X}}_1\\ ⋮\\ {\mathit{X}}_{20.}\end{array}\right]+\left[\begin{array}{c}10{\mathit{X}}_1\\ 9{\mathit{X}}_2\\ ⋮\\ 10{\mathit{X}}_{21}\end{array}\right]+\left[\begin{array}{c}{\mathit{X}}_2\\ ⋮\\ {\mathit{X}}_{21}\\ 0\end{array}\right]$。

在MATLAB®，编写创建这些矢量并将它们相加，由此得到的值的函数* x：

函数Y = afun（X）Y = [0;X（1:20）] +…[（10：-1：0）';（1:10）']。* X +…[X（2:21）;0];结束

(在示例的最后，该函数被保存为一个本地函数。)

现在，解这个线性方程组 $\mathrm{斧头}=\mathit{b}$通过提供bicgstab用函数句柄计算* x。使用的公差1 e-12和50次迭代。

B =酮（21.1）;托尔= 1 e-12;麦克斯特= 50;X1 = BICGSTAB（@ afun，B，TOL，麦克斯特）

BICGSTAB会聚在迭代11.5与相对残余5.2E-13的溶液。

X1 =21×10.0910 0.0899 0.0999 0.1109 0.1241 0.1443 0.1544 0.2383 0.1309 0.5000⋮

检查afun（X1）生成一个1的向量。

afun（X1）

ANS =21×11.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000⋮

函数Y = afun（X）Y = [0;X（1:20）] +…[（10：-1：0）';（1:10）']。* X +…[X（2:21）;0];结束