首页 > 其他分享 >5.1.2勒让德多项式

5.1.2勒让德多项式

时间:2024-11-14 21:18:53浏览次数:1  
标签:5.1 cos frac 多项式 leq quad theta

勒让德方程

\[\begin{cases} (1-x^2)\frac{d^2y(x)}{dx^2}-2x\frac{dy(x)}{dx}+l(l+1)y(x)=0, \quad -1 \leq x \leq 1 \\ |y(x)| < \infty, \quad -1 \leq x \leq 1 \end{cases}\]

\[y(x) = \sum_{n=0}^{\infty} a_nx^n \qquad a_{k+2} = \frac{(k-l)(k+l+1)}{(k+1)(k+2)}a_k \]

\[y(x) = a_0y_0(x) + a_1y_1(x) \qquad \text{收敛半径} \quad R = 1 \]

如何使 \(y(x)\) 在边界处也满足收敛性条件?\(\longrightarrow\) 将其自然地截断为多项式!

对 \(l\) 的取值作出限制并适当选取待定系数:

\(l\) 为偶数 \(\quad a_{l+2} = 0 \quad y_0(x)\) 自然地截断为多项式 \(\quad y(x)\) 截断为最高幂次为 \(l\)
\(\quad a_1 = 0 \quad\) 去掉依然发散的 \(y_1(x)\) \(\quad\) 次的多项式

\(l\) 为奇数 \(\quad a_{l+2} = 0 \quad y_1(x)\) 自然地截断为多项式 \(\quad y(x)\) 截断为最高幂次为 \(l\)
\(\quad a_0 = 0 \quad\) 去掉依然发散的 \(y_0(x)\) \(\quad\) 次的多项式

令 \(a_l = \frac{(2l)!}{2^l(l!)^2}\) 则 \(\quad y(x) \longrightarrow P_l(x) = \sum_{k=0}^N \frac{(-1)^k(2l-2k)!}{2^lk!(l-k)!(l-2k)!}x^{l-2k},\)

\(l\)-阶勒让德多项式 \(\qquad N = \begin{cases} \frac{l}{2}, & l \text{ is even} \\ \frac{l-1}{2}, & l \text{ is odd} \end{cases}\)

刘维尔本征值问题

\[\begin{cases} (1-x^2)\frac{d^2y(x)}{dx^2}-2x\frac{dy(x)}{dx}+l(l+1)y(x)=0, \quad -1 \leq x \leq 1 \\ |y(x)| < \infty, \quad -1 \leq x \leq 1 \end{cases}\]

本征值 \(\quad l = 0, 1, 2, \cdots, \quad\) 本征函数 \(\quad y_l(x) = P_l(x)\)

常用低阶勒让德多项式

\(P_0(x) = 1, \quad P_1(x) = x,\)

\(P_2(x) = \frac{1}{2}(3x^2 - 1), \quad P_3(x) = \frac{1}{2}(5x^3 - 3x),\)

\(P_4(x) = \frac{1}{8}(35x^4 - 30x^2 + 3),\)

\(x \equiv \cos \theta\)

\(P_0(\cos \theta) = 1, \quad P_1(\cos \theta) = \cos \theta,\)

\(P_2(\cos \theta) = \frac{1}{2}(3\cos^2 \theta - 1) = \frac{1}{2}(2 - 3\sin^2 \theta) = \frac{1}{4}(3 \cos 2\theta + 1),\)

\(P_3(\cos \theta) = \frac{1}{2}(5\cos^3 \theta - 3\cos \theta) = \frac{1}{8}(5\cos 3\theta + 3\cos \theta), \quad \cdots\)

[图表显示了 \(P_0\) 到 \(P_4\) 的函数图像]

标签:5.1,cos,frac,多项式,leq,quad,theta
From: https://www.cnblogs.com/RES-HON/p/18546849

相关文章

  • 5.1球坐标系下的拉普拉斯方程分离变数法求解
    拉普拉斯方程的球坐标系解法\[\frac{1}{r^2}\frac{\partial}{\partialr}\left(r^2\frac{\partialu}{\partialr}\right)+\frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partialu}{\partial\theta}\right)+\frac{1}{r^2\si......
  • 用两行命令快速搭建深度学习环境(Docker/torch2.5.1+cu118/命令行美化+插件),包含完整的
    深度学习环境的配置过于繁琐,所以我制作了两个基础的镜像,希望可以帮助大家节省时间,你可以选择其中一种进行安装,版本说明:base版本基于pytorch/pytorch:2.5.1-cuda11.8-cudnn9-devel,默认python版本为3.11.10,可以通过condainstallpython==版本号直接修改版本。dl版......
  • 多项式乘幂函数之和 2
    H4.2.1.8.多项式乘幂函数之和2\(n,k\)都是给定数,没什么区别记\(S_k=\sum_{i=1}^ni^kp^i\)\(p=1\)时\(S_k\in\Theta(n^{k+1})\)\(p<1\)时\[\begin{aligned}(1-p)S_k&=\sum_{i=1}^n\left(i^k-(i-1)^k\right)p^i-n^kp^{n+1}\\&=\sum_{i=1}^n\left(......
  • springboot2+mybatis+shardingsphere-5.5.1
    注意:1.druid不能boot-starter方式引入2.snakeyaml需要1.33('voidorg.yaml.snakeyaml.LoaderOptions.setCodePointLimit(int)') #303183.spring.datasource.driverClassName:org.apache.shardingsphere.driver.ShardingSphereDriver4.如果使用了quartz,需要指定独立数据源(Tabl......
  • 多项式板子
    一、数组版本数组版本和poly版本都只涵盖目录中第\(4\sim10\)部分。namespacePoly{intp[maxn],q[maxn],r[maxn],w[maxn];intinum[maxn];intqpow(inta,intk){intres=1;for(;k;a=1ll*a*a%mod,k>>=1)if(k&1)res=1ll*res*a%......
  • Python之函数5.1
    函数:一、什么是函数定义:函数是组织好,可重复使用,用来实现单一,或关联功能的代码段二、pycharm中的表结构项目,包(init)或目录,py文件,py文件包含多个函数或类等三、函数的有哪些优点?1、降低代码冗余2、增加代码的复用性,提高开发效率3、提高程序的拓展性4、封装:就是把代码片段......
  • 生产环境中添加多项式特征实现:将逻辑回归应用于非线性关系
            要将逻辑回归应用于非线性关系,并实现到生产环境中,我们可以通过以下步骤来完成。这里主要使用Python和Scikit-Learn库,因为它们为机器学习任务提供了强大的工具和易于使用的接口。我们将通过添加多项式特征来扩展逻辑回归模型,使其能够处理非线性关系。步骤......
  • macOS15.1及以上系统bug:开发者证书无法打开,钥匙串访问无法打开一直出现图标后立马闪退
    团队紧跟苹果最新系统发现bug:今日设备信息如下,希望能带给遇到这个问题的开发者一点帮助。错误图如下:点击证书文件后,先出现钥匙串访问图标,后立马闪退消失中间试过很多方法,都是一样的表现,最后好在解决了,看网上也没有相关的帖子,这里直接写解决办法和导致原因。&......
  • 最新植物大战僵尸杂交版最新版本2.5.1版,内置触屏+加速+全屏,附PC+安卓+iOS最全安装教程
    植物大战僵尸杂交版链接:https://pan.quark.cn/s/9f12779d1ebd在数字世界的幻境中,总有一些不朽的传奇,它们穿越时光的洪流,持续捕获玩家的眼球与热情。《植物大战僵尸》系列无疑是这股潮流中的佼佼者,以其独创的游戏机制、多彩的角色设计以及幽默诙谐的画风,赢得了无数玩家的倾心......
  • 7.10 已知一组观测数据,如表中7.17.excel(表中第一列为x的值,第二列为y的值)。试用插值方
    importnumpyasnpimportpandasaspdimportmatplotlib.pyplotaspltfromscipy.interpolateimportinterp1d,PchipInterpolator,CubicSplinefromscipy.optimizeimportcurve_fitfromscipy.statsimportnormfile_path='7.17.xlsx'data=pd.rea......