首页 > 其他分享 >Poisson 方程有限差分(一维+二维)

Poisson 方程有限差分(一维+二维)

时间:2023-09-10 21:35:19浏览次数:29  
标签:phi right frac Poisson epsilon nabla 差分 一维 Delta

Poisson equation can be writtern as follows:

\[\nabla\cdot[\epsilon(r)\nabla\phi(r)] = -q(p-n+N_D-N_A)\\ \nabla\epsilon(r)\cdot\nabla\phi(r) + \epsilon(r)\nabla^2\phi(r) = -q(p-n+N_D-N_A) \]

Since Poisson equation is based on the finite difference method (FDM), discretization is required for implementation. In 1D example, the discretization with grid size can be written as,

\[\left(\frac{\epsilon_{n+1}-\epsilon_{n}}{\Delta}\right)\left(\frac{\phi_{n+1}-\phi_{n}}{\Delta}\right)+\epsilon_n\frac{\phi_{n+1}-2\phi_{n}+\phi_{n-1}}{\Delta^2}= -q(p-n+N_D-N_A)\\ \frac{\epsilon_{n+1}}{\Delta^2}\phi_{n+1}-\frac{\epsilon_{n+1}+\epsilon_{n}}{\Delta^2}\phi_n +\frac{\epsilon_{n}}{\Delta^2}\phi_{n-1}= -q(p-n+N_D-N_A) \]

For 2D- examples,

\[\left(\frac{\epsilon_{i+1,j}-\epsilon_{i,j}}{\Delta_i}+\frac{\epsilon_{i,j+1}-\epsilon_{i,j}}{\Delta_j}\right)\left(\frac{\phi_{i+1,j}-\phi_{i,j}}{\Delta_i}+\frac{\phi_{i,j+1}-\phi_{i,j}}{\Delta_j}\right)+\epsilon_{i,j}\left(\frac{\phi_{i+1,j}-2\phi_{i,j}+\phi_{i-1,j}}{\Delta_i^2}+\frac{\phi_{i,j+1}-2\phi_{i,j}+\phi_{i,j-1}}{\Delta_j^2}\right)\\= -q(p-n+N_D-N_A)\\ \frac{\epsilon_{n+1}}{\Delta^2}\phi_{n+1}-\frac{\epsilon_{n+1}+\epsilon_{n}}{\Delta^2}\phi_n +\frac{\epsilon_{n}}{\Delta^2}\phi_{n-1}= -q(p-n+N_D-N_A) \]

标签:phi,right,frac,Poisson,epsilon,nabla,差分,一维,Delta
From: https://www.cnblogs.com/ghzhan/p/17691987.html

相关文章

  • 树上差分
    树上差分与线性差分差不多,只不过是在树上进行差分,每次将两个点x和y的标志加1,将lca(x,y)和fa(lca(x,y))的标志减1,最后来一次深搜求和,就可以得到值了下面给出几道例题1.P3128[USACO15DEC]MaxFlowP解析:树上差分板子题,直接套班子,求完值后,求最大值即可代码:#include<bits/std......
  • 通过指针变量存取一维数组元素
    通过指针变量存取一维数组元素下面展示一下。#include<stdio.h>intmain(){ inta[10],*p; for(p=a;p<a+10;p++) { scanf("%d",p); }for(p=a;p<a+10;p++) { printf("%d",*p); } printf("\n"); return0;}测试输入......
  • 在方差分析摘要中,”F“、”P值“、”P值摘要“、 ”除手段非常显著性差异 (P < 0.00)
    在方差分析摘要中,“F”、“P值”、“P值摘要”、“除手段非常显著性差异(P<0.00)吗?”、"R平方"分别代表以下内容:“F”:F值是用来衡量组间差异与组内差异之比的统计量。F值越大,说明组间差异相对于组内差异越大,也就意味着不同组之间的差异更加显著。“P值”:P值是用来衡量观察......
  • 白盒AES和SM4实现的差分故障分析
    DFA攻击背景介绍传统的密码安全性分析环境被称为黑盒攻击环境,攻击者只能访问密码系统的输入与输出,但随着密码系统部署环境的多样化,该分析模型已经不能够反映实际应用中攻击者的能力。2002年,Chow等人[1]提出了白盒攻击环境的概念,该攻击环境中的攻击者对算法运行环境具备完全的控制......
  • 差分
    目录差分例题综合运用差分例题综合运用思维+差分牛客小白77C小Why的商品归位......
  • 前缀和与差分
    前缀和一维前缀和公式:\[s[i]=s[i-1]+a[i]\]模板:constintN=10000+10;intn,m;inta[N],s[N];intmain(){ scanf("%d%d",&n,&m);for(inti=1;i<=n;i++){scanf("%d",&a[i]);s[i]=s[i-1]......
  • 学习笔记413—python实现BP神经网络进行预测和误差分析(附源代码)
    python实现BP神经网络进行预测和误差分析(附源代码)反向传播算法也称为BP神经网络,是一种带有反馈的神经网络反向学习方法,它可以对神经网络的各层上的各个神经元的各个神经元之间的连接权重进行不断迭代修改,使神经网络将输入数据转换成期望的输出数据 BP神经网络的学习过程由正向......
  • c 语言 数组(一维)做函数参数
    @TOC前言函数参数:函数参数是函数内外连接的接口,可以互通数据。一、传递一维数组函数调用时,实参是给形参初始化,所以,实参传递什么类型的数据,形参就以什么类型去接住。比如一维数组,如下:函数fun1传递a,因为数组名就是数组的首地址,所以用int*p形参。函数fun2传递&a,是一维数组......
  • 一维数组java练习
    1、打印下列图形*****************************************图形一:publicclassHomeWork8_24{publicstatic......
  • 10.Acwing基础课第797题-简单-差分
    10.Acwing基础课第797题-简单-差分题目描述输入一个长度为n的整数序列。接下来输入m个操作,每个操作包含三个整数l,r,c,表示将序列中[l,r]之间的每个数加上c。请你输出进行完所有操作后的序列。输入格式第一行包含两个整数n和m。第二行包含n个整数,表示整数序列......