关于验证Griffiths例3.3中拉普拉斯方程解的边界条件问题
验证Griffiths例3.3中拉普拉斯方程解的边界条件问题
嘿,我来帮你揪出这个问题的根源!
首先,你在写求和式的时候犯了一个小错误:原解里的n是取所有正奇数(1,3,5,...),但你的求和式写成了$\sum_{n=1}^{100} \frac{4 (47)}{(2n+1)\pi} \sin\left(\frac{(2n+1)\pi y}{5}\right)$——这会让你漏掉最关键的第一个奇数项(n=1对应的项)!因为当n从1开始时,2n+1是3,5,7,...,直接跳过了n=1的情况,而这个项是整个求和里贡献最大的正项,漏掉它之后,后面的交替正负项很容易出现负数结果,比如你算y=3时得到的-9.7875就是这么来的。
正确的求和式应该是把奇数项表示为n=2k-1(k从1到100),也就是:
$$V(0,y) = \sum_{k=1}^{100} \frac{4 \times 47}{(2k-1)\pi} \sin\left(\frac{(2k-1)\pi y}{5}\right)$$
这样k=1时对应n=1,k=2对应n=3,以此类推,完整包含了所有奇数项。
其次,咱们得回忆一下傅里叶正弦级数的收敛特性:这个解里的V(0,y)其实是$V_0$在区间$y∈(0,a)$上的奇延拓的傅里叶正弦级数。这意味着:
- 在区间内部($0<y<a$),当求和项数足够多时,级数会收敛到$V_0$(比如你取y=2.5,也就是a=5的中点,用100项求和的话结果会非常接近47)
- 在端点$y=0$和$y=a$处,级数收敛到0,这正好符合题目里给定的边界条件$V(y=0,a)=0$,所以你算这两个点得到0是完全正确的。
给你两个小建议:
- 先修正求和的索引,确保包含n=1这个关键项,再重新计算y=3的结果,应该会趋近于47
- 可以先试前5项、前10项求和,看看结果是不是逐渐向47靠近,直观感受一下傅里叶级数的收敛过程
备注:内容来源于stack exchange,提问作者The Rizzler
相关产品推荐
相关产品推荐

