求瞬时频率的差分算法

‘壹’ 什么是差分算法

在数值计算中，常用差分近似微分。
例如：
向前差分：f'(n)=f(n+1)-f(n)
向后差分：f'(n)=f(n)-f(n-1)

‘贰’ 什么是有限差分算法

有限差分法(FDM)的起源,讨论其在静电场求解中的应用.以铝电解槽物理模型为例,采用FDM对其场域进行离散,使用MATLAB和C求解了各节点的电位.由此,绘制了整个场域的等位线和电场强度矢量分布.同时,讨论了加速收敛因子对超松弛迭代算法迭代速度的影响,以及具有正弦边界条件下的电场分布.
有限差分法
有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。
该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。
分类
对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
构造差分的方法
构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式
时域有限差分法在GIS局部放电检测中的应用
1 前言
GIS由于其占地面积小以及高度的可靠性被广泛应用，但也有因为固定微粒、自由微粒以及绝缘子内部缺陷而发生的绝缘故障。一般发生绝缘故障都伴随有局部放电发生，因而局部放电检测是诊断电力设备绝缘状况的有效方法之一。超高频局部放电检测方法因为具有强的抗干扰能力和故障点定位能力而受到制造厂家和研究部门的普遍关注，并且已有部分产品应用于现场。超高频局部放电检测方法一般直接检测出局部放电脉冲的时域信号或者频谱信号，因为不同的研究者所研制的检测用传感器的带宽和检测系统(内部传感器法和外部传感器法)不同，以及传感器和局部放电源的相对位置对检测结果的影响，检测所得结果存在较大差异，缺乏可比性，因此有必要对局部放电信号的传播规律进行研究。
时域有限差分(Finite-Difference Time-Domain)法最早是由KaneS.Yee在1966年提出的，是一种很有效的电磁场的数值计算方法，不需要用到位函数，是一种在时间域中求解的数值计算方法。这种方法被应用于天线技术、微波器件、RCS计算等方面。
本文借助时域有限差分法对252KV GIS内部局部放电所激发的电磁波传播进行仿真，并用外部传感器超高频局部放电检测方法在实验室对252kV GIS固定高压导体上的固定微粒局部放电信号进行实测，仿真结果和实验结果基本一致，为超高频局部放电检测结果提供了有效的理论依据。
2 时域有限差分法
时域有限差分法是一种在时域中求解的数值计算方法，求解电磁场问题的FDTD方法是基于在时间和空间域中对Maxwell旋度方程的有限差分离散化一以具有两阶精度的中心有限差分格式来近似地代替原来微分形式的方程。FDTD方法模拟空间电磁性质的参数是按空间网格给出的，只需给定相应空间点的媒质参数，就可模拟复杂的电磁结构。时域有限差分法是在适当的边界和初始条件下解有限差分方程，使电磁波的时域特性直接反映出来，直接给出非常丰富的电磁场问题的时域信息，用清晰的图像描述复杂的物理过程。网格剖分是FDTD方法的关键问题，Yee提出采用在空间和时间都差半个步长的网格结构，通过类似蛙步跳跃式的步骤用前一时刻的磁、电场值得到当前时刻的电、磁场值，并在每一时刻上将此过程算遍整个空间，于是可得到整个空间域中随时间变化的电、磁场值的解。这些随时间变化的电、磁场值是再用Fourier变换后变到相应频域中的解。
在各向同性媒质中，Maxwell方程中的两个旋度方程具有以下形式(式(1)~(2))。

式中，ε为媒质的介电常数；μ为媒质的磁导率；σ为媒质的电导率；σ*为媒质的等效磁阻率，它们都是空间和时间变量的函数。
在直角坐标系中，矢量式(1)～(2)可以展开成以下六个标量式。

为了用差分离散的代数式恰当地描述电磁场在空间的传播特性，Yee提出了Yee Cell结构，在这种结构中，每一磁场分量总有四个电场分量环绕，同样每一电场分量总有四个磁场分量环绕，Yee对和分量在网格单位上的分布情况如图1所示。为达到精度，Yee计算和时在时间上错开半个步长，用中心差商展开偏微分方程组，得到x轴方向电场和磁场FDTD迭代公式(式(9)~(10))，Y轴和z轴迭代公式与x轴迭代公式成对称形式(略)。

FDTD方法是Maxwell方程的一种近似求解方法，为了保证计算结果的可靠性，必须考虑差分离散所引起的算法稳定性和数值色散问题，时间步长和空间步长应满足(11)~(12)条件。

其中，δ=min(△x，△y，△z)；υmax为电磁波在媒质中传播的最大相速；λmin为电磁波在媒质中的最小波长值。
式中△x，△y和△z分别是在x，y和z坐标方向的空间步长，△t是时间步长，ij和k和n是整数。
3 GIS局部放电电磁仿真和超高频检测
SF6气体绝缘的GIS中局部放电的脉冲持续时间极短，其波头时间仅几个ns。为了简化分析，将局部放电电流看成对称脉冲，一般用如下的Gaussian形状的脉冲模型来表示，根据式13和文献6本文仿真用局部放电源高斯脉冲的峰值电流取30mA，脉冲宽度取5ns，波形如图2所示。

GIS局部放电信号频带较宽，用于接收信号的传感器(天线)应该满足检测要求，本文采用超宽带(300MHz~3000MHz)自补结构的双臂平面等角螺旋天线，天线结构如图3所示。

该天线在一定频率范围内可以近似认为具有非频变天线的特性，因为GIS局放信号的频率是在一个范围内变化，对于不同频率的GIS局放信号，该天线的阻抗不随频率变化，可方便实现天线和传输线的阻抗匹配，避免波形畸变。用HP8753D网络分析仪对天线的驻波比进行测试，结果在300MHz~3000MHz的频率范围内驻波比小于2.0，根据电磁理论当驻波比小于2.0时可以不考虑驻波的影响，表明该平面等角螺旋天线在设计频率具有良好的频响特性，所测结果可靠。
超高频法把GIS看作同轴波导(如图4所示)，局部放电产生的短脉冲沿轴向传播，传感器作为接收天线，接收局部放电所激发的电磁波。

本文针对252KV GIS内高压导体上φ0.05×lcm固定突起发生局部放电进行模拟，GIS内部高压导体外直径为10.2cm，外壳内直径为29.4cm，长度为4米。采用1×l×lcm网格进行剖分，边界用完全匹配层(PML)材料吸收边界，其中绝缘子相对介电常数取3.9。采用IMST Empire电磁仿真软件分别对图4的GIS发生局部放电时内部点1和外部点2处的信号进行仿真，仿真结果如图5所示。
图5(a)和(b)的仿真结果表明在GIS内部发生局部放电时，局部放电脉冲可以激发上升沿很陡的信号，由于其内部为不连续波导结构，电磁波在其内部将引起反射和复杂谐振，频率成分可高达GHz。另外，比较内部点1和外部点2处的仿真结果，内部点1处的信号幅值是外部点2处的两倍，表明信号可以从绝缘缝隙泄漏，但由于绝缘子和缝隙的影响幅值将明显发生衰减，并且信号在绝缘缝隙处发生的折射和散射，外部信号比内部信号复杂。图5(c)表明局部放电频带比较宽，可高达GHz，信号成分较为丰富。

采用外部传感器超高频局部放电检测系统对252KV GIS内高压导体φ0.05×1cm固定突起局部放电进行实测。由于局部放电信号比较微弱，加之高频信号传播过程中衰减较大，在测试系统中采用增益不低于20dB的宽带放大器。在实验过程中对空气中的局部放电高频信号进行衰减特性研究发现该检测系统有效检测范围为17米。在外部点2处(距离GIS外壳绝缘缝隙10cm)的检测结果如图6所示。比较图5(b)和图6表明，仿真结果和实测结果基本一致，这个结论为超高频局部放电检测结果提供了理论支持。

超高频局部放电检测方法已经表明是非常有效的局部放电检测方法，本文借用时域有限差分法从信号的时域特征出发来验证局部放电检测结果，但由于不同电压等级的GIS结构存在差异，以及故障微粒的状态不同，对检测结果都有影响，并且目前还没有找出超高频方法和传统检测方法之间的内在关系，有待进一步深入研究。
4 结论
时域有限差分法对GIS局部放电脉冲所激发的电磁波仿真结果表明，局部放电信号上升沿较陡，频率可达GHz；由于绝缘子以及绝缘缝隙的影响，使得同轴波导结构不连续，将产生很复杂的电磁波。
a.由于绝缘子以及绝缘缝隙的影响，使信号幅值发生明显衰减，外部信号的幅值是内部信号幅值的一半。
b.实验结果和仿真结果基本一致，进一步从理论上论证了超高频局部放电检测方法的有效性。

‘叁’ 差分算法是什么

在数值计算中,常用差分近似微分.
最简单的差分格式有向前、向后和中心3种.
向前差分：f'(n)=f(n+1)-f(n)
向后差分：f'(n)=f(n)-f(n-1)
中心差分：f'(n)=[f(n+1)-f(n-1)]/2

‘肆’ 怎样求差分

f(n)的差分定义为：f(n)-f(n-1)
-----------------------------
y=f(x) h=0.05 的差分

= y(x)-y(x-0.05)
-----------------------------
差分是对离散函数的一种运算。相当于连续函数的微分。
-------------------------------------------------
(我的Hi有时不工作。)

‘伍’ 瞬时频率的定义

以解析讯号法定义瞬时频率，直观上，瞬时频率为相位的微分。对于自然界中的实数讯号，如何定义讯号的相位。Gabor提出解析讯号法(Analytic Signal Method)，将实数讯号表示为对应的复数讯号，即可定义复数讯号的大小与相位，将实数讯号的瞬时频率求出。 实数讯号的解析讯号(Analytic Signal)定义为 瞬时相位
瞬时频率
其中虚数项为实数讯号的希尔伯特转换(Hilbert Transform)，将它定义为。称作解析函数的理由是，此型式的复数函数满足柯西-里曼(Cauchy-Riemann)的可微分条件，称之为解析函数(Analytic Function)。因此，解析讯号可以表示为
其中
; 如果是没有任何限制条件的时间讯号，计算出来的瞬时频率可能不是正确的结果。对于平均值为零的局部对称讯号而言，前述定义的瞬时频率才具有物理意义。在1998年，黄锷(Norden E. Huang)博士提出一个有效的算法，将讯号先行分解成具有局部对称之分量，以正确地求得资料的瞬时频率。这个方法称为希尔伯特-黄转换(Hilbert Huang Transform, HHT)。 以下简单的例子来说明，对于平均值为零的讯号，此瞬时频率的定义才具有物理意义。对于一个弦波讯号，
考虑三种情况: (1) β = 0 (2) 0 < β < 1 (3) β > 1
(1) β = 0: 当弦波讯号平均值为零时，在复数平面上的描述是以坐标原点为中心的单位圆，它的相位角θ(t)则是以坐标原点为中心，逆时钟方向呈线性递增，其图形为斜率1的直线，而瞬时频率是一个常数值。
(2) 0 < β < 1: 在复数平面上仍然是一个单位圆，但其圆心从原点偏移了β个单位，其相角θ(t)不再呈现线性递增，瞬时频率出现震荡的现象，而不是应有的常数值。
(3) β > 1: 因为β值超过单位圆的半径，因此的圆心在单位圆之外。如此相位角θ(t)在[ − π/2, π/2]震荡，瞬时频率出现负值，与原讯号的特性有极大的差别。

‘陆’ 帧间差分法的算法描述

（l）、对序列图像进行3×3中值滤波预处理，去掉图像随机噪声。减少以后运算的复杂度，克服噪声对图像处理结果的干扰。
（2）、从视频图像序列中选取出背景图像所阢砂，使其只包含固定的背景图像：
（3）、在视频图像序列中选取连续的两帧图像，其中前一帧图像pk-1(x,y)，当前帧图像pk(x,y)；
（4）、计算当前帧与背景帧的差得FD(x,y)，从 图像中提取出完整的目标；
（5）、计掉当前1帧的差得FG(x,y)，得到目标的变化量；
（6）、求帧差FD（x,y）与，FG(x,y)的交集得到运动目标粗糙的运动区域幽像，
（7）、数学形态学运算使得运动区域封毕、连续、完整，并去掉背持中的噪声。
其中：（略）

‘柒’ matlab中diff函数求差分什么意思

这是我自己翻译的help

Y = diff(X) 对数组的第一维来计算相邻 X的差值（要求长度不能为1） :
（1）如果 X 是一个 m长度的向量, 那么Y = diff(X) 返回一个 m-1长度的向量。 Y 的元素是相邻 X的差值。
Y = [X(2)-X(1) X(3)-X(2) ... X(m)-X(m-1)]
（2）如果X是一个非空，非向量的p*m 矩阵，那么Y = diff(X) 返回(p-1)*m的矩阵，矩阵的元素是X每一行元素间的差值。
Y = [X(2,:)-X(1,:); X(3,:)-X(2,:); ... X(p,:)-X(p-1,:)]（X(2,:)-X(1,:)代表第2行减第1行）
（3）如果X 是一个零矩阵， 那么Y = diff(X) 返回零矩阵。

Y = diff(X,n) 通过迭代计算diff(X) n次来计算第n次的差值。事实上，这就意味着diff(X,2) 等价于diff(diff(X))。

Y = diff(X,n,dim) 对 dim所指定的维来计算n次差值。 这个dim参数是一个正整数标量。

‘捌’ 怎么用matlab求差分

调用filter函数解差分方程。

1）yn=filter(B,A,xn)是计算输入向量xn的零状态响应输出信号yn，yn与xn长度相等，
其中B=[b0,b1，……bn], A=[a0,a1,……an]。其中a0=1。

2）yn=filter(B,A,xn,xi)是计算全响应的函数。xi是等效初始条件的输入序列，xi能由初始条件确定。此时需要调用filtic函数。xi=filtic(B,A,ys,xs)。其中ys,xs是初始条件向量。
其中ys=[y(-1),y(-2)……y(-N)]，xs=[x(-1),x(-2),……x(-M)]
另外若xn为因果序列xs=0可缺省。

举个例子若y(n)-0.8y(n-1)=x(n)，初始条件y(-1)=1。
a=0.8,ys=1;
xn=[1,zeros(1,30)];
B=1, A=[1,-a];
xi=filtic(B,A,ys);
yn=filter(B,A,xn,xi);
%以下是解出来的yn图像
n=0:length(yn)-1;
stem(n,yn,'.');

‘玖’ 求 隐式差分 具体公式 急用 最好能对公式做下解释

fenbufa
分步法
method of fractional steps
把复杂的问题的每个时间步分解成若干个中间步,例如把多维问题按坐标分解成几个一维问题,然后用差分法解这些比较简单的各中间步,最后得到原始问题的近似解,这类方法叫作分步法.交替方向隐式法、预测校正法、局部一维方法、时间分裂法等都属此类.
1955年D.W.毕斯曼与H.H.瑞契福特在(,)平面上用交替方向隐式法（简称ADI方法）,解二维热传导问题
[207-8] (1)时,对[207-9]与[207-10]进行不同处理,一个取成显式(显式差分方法),一个取成隐式,并依次交替以保持对称性.取＝＝时,可得出如下格式[207-11]格式(2)用了两步合成一个循环,一般称之为P-R格式由于P-R格式交替地沿各个空间方向作一维隐式计算,也称为交替方向隐式法,(2)的每个方程组都是系数矩阵为三对直线矩阵的线性方程组,从(2)中消去[207-14]经整理可得
[208-1]把方程(1)的光滑解代入上式,其截断误差为(+),这表明P-R格式具有二阶精度.格式(2)的增长因子是
[208-2]式中[208-3](=1,2).由于对任何[208-4]都有││≤1 因此P-R格式(2)是无条件稳定的.P-R格式不宜向三维问题推广,J.道格拉斯和瑞契福特又提出了一种三维问题的交替方向隐式法,也称D-R方法.考虑三维热传导方程
[208-5] (3)取空间步长[208-6]D-R方法就是
[208-7] (4)在(4)中消去[208-8],[208-9],可得等价格式[208-10]这可说明(4)与微分方程(3)相容,(5)的增长因子是
[208-11]式中[208-12] (=1,2,3).对于一切[208-13],││≤1,因此 D-格式(4)是无条件稳定的.交替方向隐式格式除上述两种外,还有其他各种变形格式,ADI方法从计算(要分几步完成,中间要计算[208-14]或[208-15],[208-22]等.
对于热传导方程(3),H.H.亚年科1959年还提出了更简单的格式
[208-21] (6)消去[208-8],[208-9]之后,得等价格式
[208-16]展开成的幂次式,得
[208-17]这说明(6)与微分方程(3)相容,(6)的增长因子是
[208-18]所以对于一切[208-19],它是稳定的.通常称(6)是局部一维方法,它也是一种分步方法.上述方法的另一特点是把差分算子分解成为较简单的差分算子的积,因而又称算子分解法.
具体公式见下面连接

‘拾’ 偏微分方程数值解法

稳定性分析是针对某一特定的差分算法来说的。而并不是对偏微分方程来说的。一般是用Fouier分析的办法来做。
你可以看一下
余德浩,汤华中编的科学出版社出版的“微分方程数值解法”里面216页有一些相关的东西。
比较常用的差分算法有Lax_Wendroff格式以及MacCormack格式。

另外，你如果想要解析解的话，估计可能要用特征线法。或者分离变量法看一下。