﻿功能:将对称矩阵采用豪斯赫尔德方法化为对称三对角阵

格式:stm(a)

说明:a为对称矩阵

例子:
a =
 [ 5.08534607714291  3.32897863040165  6.17341125205784
   3.32897863040165  4.58814488471865  1.41467333837164
   6.17341125205784  1.41467333837164  6.49441459052936 ]
stm(a)//回车后得到如下结果
ans =
 [ 5.08534607714291        -7.0137796663928       8.88178419700125E-16
  -7.0137796663928          7.24697677231266     -0.0190873774207021
   8.88178419700125E-16    -0.0190873774207025    3.83558270293535     ]
//注意E-16明显那个参数为0,其实上面的结果我们使用命令{zeros<矩阵运算\zeros>}(ans,14)来获取更直接的结果
ans =
 [ 5.08534607714291   -7.0137796663928       0
  -7.0137796663928     7.24697677231266     -0.0190873774207021
   0                  -0.0190873774207025    3.83558270293535  ]