求瞬時頻率的差分演算法

『壹』 什麼是差分演算法

在數值計算中，常用差分近似微分。
例如：
向前差分：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格式。

另外，你如果想要解析解的話，估計可能要用特徵線法。或者分離變數法看一下。