请教W-W算法！

hongdan

最近要计算一超越特征值问题，有文章中说可以用W-W算法实现，W-W算法是用来计算频率的个数的吧，怎么用来计算固有频率呀？这个算法最主要的特点是什么？恳请赐教！

风花雪月

袁驷编写的《程序结构力学》一书中有详细的介绍

风花雪月

或者参考第十四届全国结构工程学术会议论文集 2005年中

《杆系结构自由振动精确求解的理论和算法 》一文，也是由袁驷写的，以下摘自该文

WW算法并不直接计算频率，本质上它是一种计数方法

基于WW计数法的频率计算方法有以下一些特点:
1. 方法简单,方便;
2. 结果可以达到计算机允许的任意精度;
3. 不丢根,单根,重根均可处理;
4. 可以直接求第阶频率,或第阶至第kknk+阶的频率;
5. 方法为精确方法,一个杆件一个单元即可;
6. 普通的TLDL分解或Gauss消元法即可(不必选主元);
7. 只须作分解或前消去,无须回代过程;
8. 保留K的带状稀疏性;
9. 由于的计算只需要KJ K对角线上负号元素的个数,而不要数值,因此对数值误差不敏感;万一K奇异,只需稍调整一下频率值即可.
但是,二分法最大的不足是收敛速度慢,只有线性的收敛速率.

hongdan

我看了程序结构力学这本书，看得不是很懂！它的初值是怎么确定的呀？W-W是计数方法，在计算固端频率事要分为轴向和弯曲振动来算，如果不细分这些振动，应该怎么求解呢？