- 博客(6)
- 收藏
- 关注
原创 Krylov子空间 FOM GMRES CG算法
文章目录krylov子空间代码算法1:基于Gram-Schmidt正交化的Arnoldi过程算法2:基于MGS的Arnoldi过程算法3:FOM算法4:实用FOM算法5:实用GMRES算法6:共轭梯度法算例结果krylov子空间代码投影方法给定初值 x(0)∈Rnx^{(0)}\in \mathbb{R}^nx(0)∈Rn,采用仿射空间 x(0)+Kx^{(0)}+\mathcal{K}x(0)+K,其中 K\mathcal{K}K 是krylov子空间。找到一个 x~∈x(0)+K\tilde{x
2021-12-11 22:50:11 3676 1
原创 哈密顿算子和拉普拉斯算子 格林公式高斯公式和斯托克斯公式 多重积分的分部积分公式
哈密顿算子(Hamiltionian) ∇\nabla∇标量性质u(x,y)u(x,y)u(x,y)是一个标量函数,∇u=(∂u∂x∂u∂y)\nabla u=\begin{pmatrix}\frac {\partial u}{\partial x} \\\frac {\partial u}{\partial y} \\\end{pmatrix}∇u=(∂x∂u∂y∂u) 表示梯度矢量性质u=(u1u2)u=\begin{pmatrix}u_1 \\u_2\\\end{pmatr
2021-11-22 21:40:30 12524 2
原创 偏微分方程基本概念
一般形式D(x,y,...,u,∂u∂x...∂2u∂x2,∂2u∂x∂y,....)=0D(x,y,...,u, \frac {\partial u}{\partial x}...\frac {\partial ^2u}{\partial x^2},\frac {\partial^2 u}{\partial x \partial y},....)=0D(x,y,...,u,∂x∂u...∂x2∂2u,∂x∂y∂2u,....)=0自变量 x y未知量 U阶数 U的最高次偏导数是几阶导次数
2021-11-22 15:27:43 1190
原创 线性抛物方程的通用有限差分格式及代码(未完)
目录题目通用差分格式代码题目通用差分格式现在是type1 dirichlet边界条件将显示格式和隐式格式写成一个格式,当θ=1\theta=1θ=1是隐式格式,θ=0\theta=0θ=0是显式格式,θ=1/2\theta=1/2θ=1/2是CN格式。离散区域 M=(L−l)/Δx,N=T/Δt,k=1...M+1M=(L-l)/\Delta x,N=T/\Delta t,k=1...M+1M=(L−l)/Δx,N=T/Δt,k=1...M+1,n=1…N+1离散初值 uk1=g(k∗Δx),
2021-11-17 22:51:53 1053
原创 RL分数阶微积分
文章目录将积分和导数统一RL分数阶积分RL分数阶导数先积分后求导先求导后积分(t−α)β(t-\alpha)^\beta(t−α)β的分数阶导数和GL的联系将积分和导数统一重新表示n重积分为下面这个函数,n为整数f(−n)(t)=∫atdτ1∫aτ1dτ2...∫aτn−1f(τ)dτ(1) f^{(-n)}(t)=\int^t_a\mathrm{d}\tau_1\int_a^{\tau_1}\mathrm{d}\tau_2...\int_a^{\tau_{n-1}}f(\tau)\mathrm{d}
2021-11-16 22:08:33 1605
原创 GL分数阶微积分
预备公式1.Γ(z+1)=zΓ(z)\Gamma(z+1)=z\Gamma(z)Γ(z+1)=zΓ(z)2.Γ(z)=∫0∞e−ttz−1dt=limn→∞n!nzz(z+1)...(z+n)=limn→∞nz[pr](p+r)\Gamma(z)=\int_{0}^{\infty} e^{-t}t^{z-1}\mathrm{d}t=\lim_{n\to\infty}\frac{n!n^{z}}{z(z+1)...(z+n)}=\lim_{n\to\infty}\frac{n^{z}}{\begin{
2021-11-14 17:17:13 2128 1
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人