| 查看: 757 | 回复: 1 | |||
[交流]
增量谐波平衡法(IHB法)迭代初值选择问题
|
|
以增量谐波平衡法(IHB法)为代表的迭代型半解析方法，都要面对一个共同的问题，就是如何选择迭代的初值。 初值选的不好，即便程序是正确的，迭代也不会收敛，原因在于所选取的迭代初值不在收敛域内(选取的初值离准确值太远)。 如何选取迭代初值是一个棘手的问题。在这里介绍一种从数值解中拟合迭代初值的方法，方法的详细介绍可以参见论文“High-precision semi-analytical solution for the quasi-periodic nanobeam system based on the weight time-domain minimum residual method”的附录部分。 步骤： 1：首先通过数值方法(如RK法算出)系统的稳态响应(振动问题一般是周期或准周期解)，记为NUM 2：上述数值解可以展开为傅里叶级数NUM=a_0+∑a_i*cos(iwt)+∑b_i*sin(iwt)=A*J_a,其中A为傅里叶级数的系数向量，J_a为三角函数的矩阵 3：则A=J_a\NUM。倘若系统的频率w也是未知的，可以对数值解NUM做FFT变换，从图像中直接观察出W。至此，获得的系数向量A就可以作为IHB法的迭代初值。 需要注意的是，上述过程获得的初值A只是从数值解中拟合出来的谐波系数，并不是半解析解的系数，上述选择初值A作为迭代初值就可以获得IHB法的半解析解。 上述详细推导过程参见论文“High-precision semi-analytical solution for the quasi-periodic nanobeam system based on the weight time-domain minimum residual method”的附录部分，论文第一页下方也有程序的下载链接。 备注：因域名owenyaa.tk未能及时续费，下载链接需要将owenyaa.tk更换为owenyaa.com就可以下载。 |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有21人回复
-
导师想让我从独立一作变成了共一第一
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有9人回复
-
申请2026年博士
已经有6人回复
2楼2023-11-11 14:54:05