合刚度矩阵,然后用与有限元方法完全同样的方式解一个代数方 程组,得到ψ,其余节点值可以由(1.⑦)确定 如果A是一个可逆矩阵,上述广义特征值问题可以化为普 通特征值问题,以R2左乘(1.12)式得 CAT)IK (AT)-IA Ag Ag 它就是一个普通特征值问题 特征值λ可能是复数,这就需要作复数运算。为避免这 点,我们可以将公式(1.13)胳作变形。特征值与特征向量必定是 成对共轭地出现的,设λ=a±,g=p土i是一对特征值与特 征向量,由(1.11)得 X(t1)=(±iB)(p±iq)。 分离实部与虛帑得 Xp=ap-Bq, Xq=Rp+ag, Y(卩,4)=(,q) 我们替换矩阵T与A中对应的两列,令 不难看出 X=711T 将各对复特征值一一作了以上替换以后,求转移矩阵X的运算就 是实的了。 §2 Fourier方法 上节中給出的方法可适用于任何的相似剖分,其运算量为解