您所在位置：主题：[讨论]如何提高精度计算
liutao1214
[专家分：30] 发布于 2011-11-21 21:57:00
现有三个向量，
Y1=(-9.5241308E-02,-3.4517303E-02,0.2090452),
Y2=(-8.0674507E-02,-4.2283848E-02,0.1591480)
Y3=(0.2652879,0.8108509,0.3997021)
如何对于这三个向量进行高精度的Gram-Schmidt正交化？
本人进行尝试过，效果十分差，因为其中涉及到浮点数的计算问题。
更加一般的问题是，有两个向量正交向量X,Y，其中Y的范数（模、长度）很小，比如说次数为-7.理论上计算X*Y=0,实际上计算X*Y的次数为-8次。如果再次对向量X,Y进行规范化(即标准化，使每个向量的模为1)，这时再看X*Y的值可以发现次数很高，比如说-2次。这个结果很不符合正交化的特征，如何解决这一问题。
本人在利用wolf计算Lyapunov指数谱，如果有相同兴趣的朋友可以一起讨论。
Y1=(-9.5241308E-02,-3.4517303E-02,0.2090452),
Y2=(-8.0674507E-02,-4.2283848E-02,0.1591480)
Y3=(0.2652879,0.8108509,0.3997021)
如何对于这三个向量进行高精度的Gram-Schmidt正交化？
本人进行尝试过，效果十分差，因为其中涉及到浮点数的计算问题。
更加一般的问题是，有两个向量正交向量X,Y，其中Y的范数（模、长度）很小，比如说次数为-7.理论上计算X*Y=0,实际上计算X*Y的次数为-8次。如果再次对向量X,Y进行规范化(即标准化，使每个向量的模为1)，这时再看X*Y的值可以发现次数很高，比如说-2次。这个结果很不符合正交化的特征，如何解决这一问题。
本人在利用wolf计算Lyapunov指数谱，如果有相同兴趣的朋友可以一起讨论。