Newton插值法在电力系统采样中应用
牛顿插值法原理及应用
牛顿插值法插值法是利用函数f (x)在某区间中若干点的函数值,作出适当的特定函数,在这些点上取已知值,在区间的其他点上用这特定函数的值作为函数f (x)的近似值。
如果这特定函数是多项式,就称它为插值多项式。
当插值节点增减时全部插值基函数均要随之变化,这在实际计算中很不方便。
为了克服这一缺点,提出了牛顿插值。
牛顿插值通过求各阶差商,递推得到的一个公式:f(x)=f[x0]+f[x0,x1](x-x0)+f[x0,x1,x2](x-x0)(x-x1)+...f[x0,...xn](x-x0 )...(x-xn-1)+Rn(x)。
插值函数插值函数的概念及相关性质[1]定义:设连续函数y-f(x) 在区间[a,b]上有定义,已知在n+1个互异的点x0,x1,…xn上取值分别为y0,y1,…yn (设a≤ x1≤x2……≤xn≤b)。
若在函数类中存在以简单函数P(x) ,使得P(xi)=yi,则称P(x) 为f(x)的插值函数.称x1,x2,…xn 为插值节点,称[a,b]为插值区间。
定理:n次代数插值问题的解存在且唯一。
牛顿插值法C程序程序框图#include<stdio.h>void main(){float x[11],y[11][11],xx,temp,newton;int i,j,n;printf("Newton插值:\n请输入要运算的值:x=");scanf("%f",&xx);printf("请输入插值的次数(n<11):n=");scanf("%d",&n);printf("请输入%d组值:\n",n+1);for(i=0;i<n+1;i++){ printf("x%d=",i);scanf("%f",&x[i]);printf("y%d=",i);scanf("%f",&y[0][i]);}for(i=1;i<n+1;i++)for(j=i;j<n+1;j++){ if(i>1)y[i][j]=(y[i-1][j]-y[i-1][j-1])/(x[j]-x[j-i]);elsey[i][j]=(y[i-1][j]-y[i-1][j-1])/(x[j]-x[j-1]);printf("%f\n",y[i][i]);}temp=1;newton=y[0][0];for(i=1;i<n+1;i++){ temp=temp*(xx-x[i-1]);newton=newton+y[i][i]*temp;}printf("求得的结果为:N(%.4f)=%9f\n",xx,newton);牛顿插值法Matlab程序function f = Newton(x,y,x0)syms t;if(length(x) == length(y))n = length(x);c(1:n) = 0.0;elsedisp('x和y的维数不相等!');return;endf = y(1);y1 = 0;l = 1;for(i=1:n-1)for(j=i+1:n)y1(j) = (y(j)-y(i))/(x(j)-x(i));endc(i) = y1(i+1);l = l*(t-x(i));f = f + c(i)*l;simplify(f);y = y1;if(i==n-1)if(nargin == 3)f = subs(f,'t',x0);elsef = collect(f); %将插值多项式展开f = vpa(f, 6);endend牛顿插值法摘要:值法利用函数f (x)在某区间中若干点的函数值,作出适当的特定函数,在这些点上取已知值,在区间的其他点上用这特定函数的值作为函数f (x)的近似值。
牛顿(newton)插值法
牛顿(newton)插值法牛顿插值法是一种数值分析中的插值方法,它用于找到一个多项式函数,该函数会经过给定的一系列数据点。
该方法最初由英国数学家艾萨克·牛顿(Isaac Newton)发明并称为插值多项式,它也被称作差分插值法。
插值是数学和工程学中的一项重要任务,它是用于在给定数据点之间构建连续函数的一种数值方法。
插值方法通常涉及过渡从观察结果派生出抽象结果的过程,从而使得预测可能的结果取得更加准确。
下面介绍牛顿插值法的基本原理。
插值基础插值基础是插值方法中的一个重要概念。
在这里,我们将对牛顿插值法中用到的插值基础进行简要介绍。
一个插值基础是指一个已知数据点的集合,通常是一个 x 坐标和对应的 y 坐标。
每个插值基础一般定义为一个数据点的函数,该函数包含了给定点的所有信息并将这些信息用于构建连续函数。
在牛顿插值法中,我们使用差分来定义插值基础。
差分是指两个相邻数据点之间 y 坐标的差值。
具体来说,若给定以下节点:x0, y0x1, y1x2, y2...xn, yn我们则通过以下的 "+" 符号所示的不断进行差分的方式来构建一个插值基础:y0y1-y0…yn-yn-1 yn-yn-1 yn-yn-2 ... yn-y0上述图表所展示的差分的值即为定义插值基础的差商(divided difference)。
牛顿插值公式基于上述插值基础和差商,我们现在可以使用牛顿插值公式来实现插值。
具体来说,牛顿插值公式可以表示为:f(x) = y0 + d1*f[x0,x1] + d2*f[x0,x1,x2] + ... + dn*f[x0,x1,...,xn]其中 f(x) 是插值函数,x0, x1, ..., xn 是给定的节点,y0, y1, ..., yn 是对应的 y 值,f[x0,x1] 是差商 f(x0,...,x1) 的值,d1, d2, ..., dn 也是差商。
请注意,插值函数的次数最高为 n - 1,这意味着插值函数与插值基础的次数相同。
newton插值法实验报告
竭诚为您提供优质文档/双击可除newton插值法实验报告篇一:牛顿插值法matlab程序《计算方法》数值实验报告篇二:数值分析实验报告(插值法)武汉理工大学学生实验报告书实验课程名称数值分析开课学院计算机科学与技术学院指导老师姓名学生姓名学生专业班级20XX—20XX学年第一学期实验课程名称:数值分析篇三:c++实现牛顿插值法实验报告数值实验用newton商差公式进行插值姓名:陈辉学号:13349006院系:数据科学与计算机学院专业:计算机科学与技术班级:计科一班日期:20XX-10-11指导老师:纪庆革目录一、二、三、四、五、六、七、八、实验目的................................................. .........................................3实验题目................................................. .........................................3实验原理与基础理论................................................. .......................3实验内容................................................. .........................................6实验结果................................................. .......................................10心得体会................................................. .......................................14参考资料................................................. .......................................14附录(源代码)............................................... (14)一、实验目的编写一个程序,用牛顿差商公式进行插值。
供电系统模拟仿真中的高精度数值计算方法探讨
供电系统模拟仿真中的高精度数值计算方法探讨高精度数值计算是一项在供电系统模拟仿真中至关重要的技术,它能够提供准确和可靠的结果,有效帮助工程师进行系统设计和优化。
本文将探讨在供电系统模拟仿真中的高精度数值计算方法,并分析其在实际应用中的优势和挑战。
供电系统模拟仿真是一个复杂的过程,涉及到电力网络的各个方面,包括发电、输电、配电和用电等环节。
为了准确地模拟系统的运行,需要对电力网络进行建模,并进行数值计算以获取各个参数的精确值。
由于涉及到大量的复杂计算,因此高精度数值计算方法变得至关重要。
在供电系统模拟仿真中,高精度数值计算方法可以提供准确的电流、电压、功率和能量等参数的计算结果。
这些参数对于系统设计和运行至关重要,能够帮助工程师做出合理的决策。
例如,对于输电线路的设计,通过高精度数值计算可以获得准确的负载流量分布,从而确保线路的安全运行。
另外,在电力系统的稳定性分析中,高精度的数值计算可以提供准确的功率平衡和电压稳定性等指标,帮助工程师进行系统优化。
高精度数值计算方法中的一个常见技术是牛顿-拉夫逊方法(Newton-Raphson method)。
该方法能够通过迭代计算的方式,逐步逼近方程的根。
在供电系统模拟仿真中,这一方法可以用于求解电力网络的节点电压和线路电流等参数。
通过不断迭代计算,牛顿-拉夫逊方法可以提供较高的计算精度。
另一个常用的高精度数值计算方法是快速傅里叶变换(Fast Fourier Transform,FFT)。
在供电系统模拟仿真中,FFT可以用于计算频率响应以及特定频段内的信号变换。
通过FFT,工程师可以快速而准确地分析电力系统中的频率成分和谐波等特征。
这对于故障和异常情况的检测以及滤波器的设计都具有重要意义。
除了上述方法之外,高精度数值计算还包括逼近法和插值法等技术。
逼近法可以通过将连续函数逼近为多项式或分段函数来进行数值计算,从而提高计算精度。
插值法则用于根据有限数据点推测未知数据点的数值,从而提供更精确的计算结果。
三种插值方法的应用与比较
样条函数中寻找满足下列插值条件:
c.y(xi)=f(xi() i=1,2,…n); d.一般形如 y"(a)=y"(b)=0 等边界条件;的插值
函数 y(x)的方法.[1,2,3]
特点:三次样条插值函数序列一致收敛于被插
函数,因此可通过增加节点的方法提高插值的精
度.
以函数
f(x)=
1 1+x2
为例,用三种不同的插值方
出版社,2001. 〔2〕陈公宁,沈嘉骥.计算方法[M].北京:高等教育出
版社,2002. 〔3〕 李 有 法 . 数 值 计 算 方 法 [M]. 北 京 : 高 等 教 育 出 版
社,2002. 〔4〕王 沫 然.MATLAB 5.X 与 科 学 计 算[M].北 京 :
清华大学出版社,2000. 〔5〕姜启源,谢金星,叶俊.数学模型(第三版)[M].北
京:高等教育出版社,2005.
-3-
例 题:给出 f(x)=lnx 下面的数值表,用 Larange 插值计算 ln(0.54)的近似值.
x 0.4
0.5
0.6
0.7
0.8
lnx - 0.916291 - 0.693147 - 0.510826 - 0.356675 - 0.223144
在 Matlab 命令窗口中输入: x=[0.4:0.1:0.8]; y= [-0.916291 -0.693147 -0.510826 -0.356675 -0.223144]; lagrange(x,y,0.54)[3] ans=
legend(’拉格朗日插值曲线’,’原曲线’,’分段
插值曲线’)[4]
Runge 现象的产生
作分段线性插值的目的在于克服 Lagrange 插
MPPT算法在光伏发电中的应用
MPPT算法在光伏发电中的应用1Newton插值法1.1最大功率点附近的多项式推导分析光伏电池模型实际是由单个电池串并联组成,其中并联电阻很大,串联电阻很小,两者之比几乎为104倍,因此,在实际应用中忽略了因并联电阻存在的漏电流和因内电阻引起输出电压下降。
那么光伏电池数学模型可以简化为。
光伏发电系统稳定工作点在恒压区,也就是工作点的电压变化很小,因此光伏电池最大功率点附近的工作电压U可近似等于最大功率点处的工作电压Um。
根据简化的模型可以得到在最大功率点附近的输出功率。
1.2Newton插值MPPT算法Newton插值算法利用二次插值的思想,仅通过一步就可以跟踪到最大功率点,在跟踪速度和精准度上有很大优点。
Newton插值MPPT算法原理:采集系统最大功率点附近3个工作点(U0,P(U0)),(U1,P(U1)),(U2,P(U2)),应用Newton插值法构造出光伏电池V-P曲线。
Newton插值法的实现过程如表1所示。
1.3改进型MPPT算法通过对以上两种方法的分析和研究,扰动观察法简单易实现但存在振荡,Newton插值法实现跟踪精度高,所以本文提出了一种将两者相结合的改进型MPPT算法,能够提高系统最大功率跟踪的稳态精度和跟踪的精度。
该算法具体实现的思想是首先用变步长扰动观察法快速跟踪到最大功率点附近。
0时,说明最大功率点在所取的三个点之间,此时用Newton插值法快速拟合出三点之间的曲线,精确的跟踪到最大功率点处。
该算法实现的流程图如图3所示。
情况(1)如图4所示,说明此时在最大功率点的左侧正在上坡,那么需要沿着同方向扰动,扰动量是变步长的,步长扰动量为h=ak1k2,a为扰动速度因子,文中取2,扰动后令电压:uout=u(k-2)+h。
情况(2)如图5所示,说明此时在最大功率点的右侧正在下坡,那么需要向相反方向扰动,扰动后的电压为uout=u(k)-h。
当k1k2<0时,k1>0,k2>0或者k1<0,k2>0,如图6中③所示,说明此时最大功率点在这三点之间,那么用Newton插值法快速跟踪到最大功率点处,用DSP来控制输出电压。
牛顿拉夫逊法计算潮流步骤
牛顿拉夫逊法计算潮流步骤牛顿拉夫逊法(Newton-Raphson method)是一种用于求解非线性方程组的迭代方法,它可以用来计算电力系统潮流的解。
潮流计算是电力系统规划和运行中的重要任务,它的目标是求解电力系统中各节点的电压幅值和相角,以及线路的功率流向等参数,用于分析电力系统的稳定性和安全性,以及进行电力系统规划和调度。
下面是使用牛顿拉夫逊法计算潮流的一般步骤:步骤1:初始化首先,需要对电力系统的各个节点(包括发电机节点和负荷节点)的电压幅值和相角进行初始化,一般可以使用其中一种估计值或者历史数据作为初始值。
步骤2:建立潮流方程根据电力系统的潮流计算模型,可以建立节点电压幅值和相角的平衡方程,一般采用节点注入功率和节点电压的关系来表示。
潮流方程一般是一个非线性方程组,包含了各个节点之间的复杂关系。
步骤3:线性化方程组将潮流方程组进行线性化处理,一般采用泰勒展开的方法,将非线性方程组变为线性方程组。
线性化的过程需要计算雅可比矩阵,即方程组中的系数矩阵。
步骤4:求解线性方程组利用线性方程组的求解方法,比如高斯消元法或LU分解法等,求解线性方程组,得到电压幅值和相角的修正量。
步骤5:更新节点电压根据线性方程组的解,更新各个节点的电压幅值和相角,得到新的节点电压。
步骤6:检查收敛性判断节点电压的修正量是否小于设定的收敛阈值,如果满足收敛条件,则停止迭代;否则,返回步骤3,循环进行线性化方程组和线性方程组的求解。
步骤7:输出结果当潮流计算收敛时,输出最终的节点电压幅值和相角,以及线路的功率流向等参数。
牛顿拉夫逊法是一种高效、快速且收敛性良好的方法,广泛应用于电力系统潮流计算。
在实际应用中,可能会遇到多次迭代或者收敛性不好的情况,此时可以采用退火技术或其他优化算法进行改进。
此外,牛顿拉夫逊法的计算也可以并行化,利用多核处理器或者分布式计算集群来加速计算过程。
总之,牛顿拉夫逊法是一种重要的潮流计算方法,通过迭代计算逼近非线性方程组的解,可以得到电力系统中各节点的电压幅值和相角,用于分析电力系统的稳定性和安全性。
Newton插值法的程序设计与应用
Newton插值法的程序设计与应用Newton插值法是一种常用的数值插值方法,用于通过已知数据点的函数值来估计未知数据点的函数值。
它基于多项式插值的思想,通过构造一个插值多项式来逼近原始函数。
在程序设计中,我们可以使用Newton插值法来实现数据的插值计算。
下面将详细介绍Newton插值法的程序设计与应用。
1. 程序设计思路- 输入:已知数据点的横纵坐标数组,以及待估计的未知数据点的横坐标。
- 输出:通过Newton插值法估计得到的未知数据点的纵坐标。
- 程序设计步骤:1) 根据已知数据点的横纵坐标数组,计算并存储差商表。
2) 根据待估计的未知数据点的横坐标,利用差商表计算插值多项式的系数。
3) 利用插值多项式的系数计算未知数据点的纵坐标。
2. 差商表的计算差商表是Newton插值法的关键部份,用于计算插值多项式的系数。
差商表的计算可以通过递推公式来实现。
- 输入:已知数据点的横纵坐标数组。
- 输出:差商表,存储了各阶差商的值。
- 程序设计步骤:1) 初始化差商表为一个空的二维数组。
2) 将已知数据点的纵坐标作为差商表的第一列。
3) 逐阶计算差商表的其他列,使用递推公式:f[x0, x1, ..., xn] = (f[x1, ..., xn] - f[x0, ..., xn-1]) / (xn - x0)。
4) 返回计算得到的差商表。
3. 插值多项式的系数计算插值多项式的系数可以通过差商表来计算,具体计算方法如下:- 输入:差商表,待估计的未知数据点的横坐标。
- 输出:插值多项式的系数。
- 程序设计步骤:1) 初始化插值多项式的系数为一个空的数组。
2) 逐阶计算插值多项式的系数,使用递推公式:a[i] = f[x0, x1, ..., xi]。
3) 返回计算得到的插值多项式的系数。
4. 未知数据点的纵坐标计算利用插值多项式的系数,可以计算待估计的未知数据点的纵坐标。
- 输入:插值多项式的系数,待估计的未知数据点的横坐标。
几种插值法简介[整理版]
举例来看:可以认为某水文要素T随时间t的变化是连续的,某一个测点的水文要素T可以看作时间的函数T=f(t),这样在实际水文观测中,对测得的(n+1)个有序值进行插值计算来获取任意时间上的要素值。
①平均值法:若求Ti 和Ti+1之间任一点T,则直接取T为Ti和Ti+1的平均值。
插值公式为:T=Ti+Ti+1 2②拉格朗日(Lagrange)插值法:若求Ti 和Ti+1之间任一点T,则可用T i-1、T1、T i+1三个点来求得,也可用T i、T i+1、T i+2这三个点来求得。
前三点内插公式为:T=(t-t i)(t-t i+1)(t i-1-t i)(t i-1-t i+1)T i-1+(t-t i-1)(t-t i+1)(t-t i-1)(t-t i+1)T i+(t-t i)(t-t i-1)(t i+1-t i)(t i+1-t i-1)T i+1后三点内插公式为:T=(t-t i+1)(t-t i+2)(t i-t i+1)(t i-t i+2)T i+(t-t i)(t-t i+2)(ti-t i)(t i-t i+2)T i+1+(t-t i)(t-t i+1)(t i+2-t i)(t i+2-t i+1)T i+2为提高插值结果可靠性,可将前后3点内插值再进一步平均。
③阿基玛(Akima)插值法:对函数T=f(t)的n+1个有序型值中任意两点T i和T i+1满足:f(t i)=T i dfdt|t-ti=k i f’(t i+1)=T’idfdt|t-ti+1=k i+1式中k i,k i+1为曲线f(t)在这两点的斜率,而每点的斜率和周围4个点有关,插值公式为:T=P0+P1(t-t i)+P2(t-t i)2+P3(t-t i)3,来对T i和T i+1之间的一点T进行内差。
④牛顿(Newton)插值法:若求Ti 和Ti+1之间任一点T,插值公式为:T=f(x0)+(x-x0)f(x0,x1)+ (x-x0)(x-x1)f(x0,x1,x2)+…+(x-x0)(x-x1)…(x-x n-2)f(x0,x1,…,x n-1)式中,f(x0,x1),f(x0,x1,x2),…f(x0,x1,…,x n-1)是函数f(x)的1到第n-1阶差商。
插值法及其应用【文献综述】
文献综述信息与计算科学插值法及其应用插值问题是数值计算中基础而又核心的问题. 在许多实际问题及科学研究中, 因素之间往往存在着函数关系, 然而, 这种关系经常很难有明显的解析表达, 通常只是有观察与测试得到一些离散数值.有时即使给出了解析表达式, 却由于表达过于复杂, 不仅使用不便, 而且不易与进行计算与理论分析. 例如在工程实际问题中, 我们也经常会碰到诸如此类的函数计算问题, 被计算的函数有时不容易直接计算, 如表达式过于复杂或者只希望能用一个“简单函数”逼近被计算函数, 然后用该简单函数的函数值近似代替被计算函数的函数值. 这种方法就叫插值逼近或者插值法. 插值法要求给出函数的一个函数表, 然后选定一种简单)(x f 的函数形式, 比如多项式、分段线性函数及三角多项式等, 通过已知的函数表来确定一个简单的函数作为的近似, 概括地说, 就是用简单函数为离散数组建立连续模型. )(x )(x f 插值方法是一类古老的数学方法, 它来自生产实践. 早在数千多年前, 由于经典的牛顿力学尚未诞生, 因而人们无法用解析式描述日月五星的运行规律. 我们的祖先凭借插值方法, 利用对日月五星运行规律的有限个观测值获得了比较完整的日月五星的运行规律. 在一千多年前的隋唐时期, 中华先贤在制定历法的过程中就已经广泛地运用了插值技术. 公元6世纪, 隋朝刘焊已将等距结点的二次插值应用于天文计算. 但插值的基本理论和结果是在微积分产生以后才逐步完善的, 随后其应用也日益增多, 特别是在电子计算机广泛使用以后, 由于航空、造船、精密机械加工等实际问题的需要,使插值法在实践上或理论上显得更为重要, 并得到进一步发展.经典的插值方法以Taylor 插值和Lagrange 插值为代表. Taylor 插值利用函数在定义域内某点处的阶至n 阶导数信息给出复杂函数或未知函数的近似多项式表达式, Lagrange 插0值利用多个离散点的函数信息给出函数的近似多项式的表达式, 进一步根据插值结果对复杂函数或未知函数相关的理论和应用问题做出讨论.因此Taylor 插值和Lagrange 插值有着紧密的联系, Taylor 插值可以看作Lagrange 插值的极限形式;Lagrange 插值则是Taylor 插值的离散化形式.Lagrange 插值的优点是插值多项式特别容易建立, 缺点是增加节点时原有多项式不能利用, 必须重新建立, 即所有基函数都要重新计算, 这就造成计算量的浪费;Newton插值多项式是代数插值的另一种表现形式, 当增加节点时它具有所谓的“承袭性”. .在很多实际应用问题中, 为了保证插值函数能更好的密合原来的函数, 不但要求插值函数“过点”, 即插值函数和被插值函数在节点上具有相同的函数值, 而且要求“相切”, 即两者在节点处还具有相同的导数值, 这类插值称作切触插值, 即Hermite插值.由于Taylor 插值利用的是“一点”的各阶导数信息, Lagrange插值利用的“多点”函数信息, 而Hermite插值即利用函数值信息又利用导数信息, 所以Hermite插值是Taylor插值和Lagrange插值的综合和推广.现在, 插值技术应用越来越广泛了. 当我们尚未认识到某一事物的本质时, 常从其观测点出发, 利用插值技术以加深或拓展对该事物的认识或解决某些特定的问题.密钥共享即是插值法的应用之一. 在现代密码体制中, 数据的加密算法是公开的, 数据的安全性主要取决于对密钥的保护. 现在基于Lagrange插值多项式也研究出了一种密钥共享方法, 解决了密钥保护问题.目前实际中使用的也不仅仅局限与上述的插值方法, 很多都是对经典方法的改进, 例如《空间插值法在地阶梯度场中的分析》一文中介绍了4种空间插值法在房产估值中的应用, 4种空间插值方法都各有优缺点, 作者通过对各种方法研究比较, 最后选择克里金插值法作为住宅用地地价梯度场研究的主要方法.根据其研究成果房地产决策者和规划者可以对城市居住用地的土地利用向最有效使用方向调整, 最大限度的实现土地的最高使用价值.随着计算机的发展以及图像处理的重要性, 插值法在计算机图像处理中也有着重要的作用.图像放大是一种常用的数字图像处理技术, 在航天航天、医学、通讯、多媒体等领域有着广泛的应用, 常用的图像插值算法中, 最临近插值算法的实现最为简单、方便, 但它只是原始像素简单复制到其邻域内, 放大图像会出现明显的方块或锯齿, 即我们平时所说的失真.目前较为好的方法之一即双线性插值算法, 双线性插值算法利用映射点在输入图形的4个邻点的灰度值对映射点进行插值, 即插值点处的数值用待插点最近的4个点的值加权求得. 双线性插值能够较好的消除锯齿, 放大后图像平滑性好.但是其缺点是图像高频信息丢失严重, 即图像细节与轮廓的模糊, 影响了放大图像的清晰度.因此在双线性插值放大技术的基础之上, 加入边缘锐化处理, 增强平滑图像的轮廓, 使放大后的图像有较好的清晰度.插值法的基本理论和结果是在微积分产生以后才逐步完善的, 其应用也日益增多, 特别是在电子计算机广泛使用以后, 由于航空、造船、精密机械加工等实际问题的需要, 使插值法在实践上或理论上显得更为重要, 并得到进一步发展.参考文献[1]李庆扬, 王能超, 易大义.数值分析.第4版[M].北京:清华大学出版社, 2001.[2]黄铎, 陈兰平, 王凤.数值分析[M].北京:科学出版社, 2000.[3]沈燮昌.多项式最佳逼近实现[M].上海:上海科学技术出版社, 1984.[4]Stoer J, Bulirsh R. Introduction to Numerical Analysis[M].New York:Springer-Verlag, 1980 .[5]吴才斌.插值法及其应用[J].湖北大学成人教育学院学报, 1999, 17(5):77-80.[6]杨士俊, 王兴华.Hermite插值多项式的差商表示及其应用[J].高校应用数学学报, 2006,21(1):70-78.[7]姜琴, 周天宏.常见的插值法及其应用[J].郧阳师范高等专科学校学报, 2006, 26(3):6-8.[8]陈文略, 王子羊.三次样条插值在工程拟合中的应用[J].华中师范大学学报, 2004,38(4):418-422.[9]朱春钢.二元线性样条插值[J].应用数学, 2006, 19(3):575-579.[10]李洪杰.关于三次样条插值方法在应用中的一点改进[J].计算机与应用化学, 1991,8(3):187-190.[11]王芳.牛顿插值法在数学中的应用[J].浙江师范大学学报, 1994, 17(4):67-73.[12]张元巨.Hermite插值的一种新形式[J].苏州科技学院学报, 2004, 24(3):27-29.[13]文畅平.埃米尔特插值函数的工程应用[J].人民黄河, 2006, 28(4):69-70.[14]C.R.Selvaraj. Lacunary interpolation by consine polynomials [J].Hungar, 1994, 64(4):361-372.[15]R.D.Riess. Error estimates of hermite interpolation [J].BIT, 1973, 13:338-343.。
牛顿插值法原理及应用
牛顿插值法插值法是利用函数f (x)在某区间中若干点的函数值,作出适当的特定函数,在这些点上取已知值,在区间的其他点上用这特定函数的值作为函数f (x)的近似值。
如果这特定函数是多项式,就称它为插值多项式。
当插值节点增减时全部插值基函数均要随之变化,这在实际计算中很不方便。
为了克服这一缺点,提出了牛顿插值。
牛顿插值通过求各阶差商,递推得到的一个公式:f(x)=f[x0]+f[x0,x1](x-x0)+f[x0,x1,x2](x-x0)(x-x1)+...f[x0,...xn](x-x0 )...(x-xn-1)+Rn(x)。
插值函数插值函数的概念及相关性质[1]定义:设连续函数y-f(x) 在区间[a,b]上有定义,已知在n+1个互异的点x0,x1,…xn上取值分别为y0,y1,…yn (设a≤ x1≤x2……≤xn≤b)。
若在函数类中存在以简单函数P(x) ,使得P(xi)=yi,则称P(x) 为f(x)的插值函数.称x1,x2,…xn 为插值节点,称[a,b]为插值区间。
定理:n次代数插值问题的解存在且唯一。
牛顿插值法C程序程序框图#include<stdio.h>void main(){float x[11],y[11][11],xx,temp,newton;int i,j,n;printf("Newton插值:\n请输入要运算的值:x=");scanf("%f",&xx);printf("请输入插值的次数(n<11):n=");scanf("%d",&n);printf("请输入%d组值:\n",n+1);for(i=0;i<n+1;i++){ printf("x%d=",i);scanf("%f",&x[i]);printf("y%d=",i);scanf("%f",&y[0][i]);}for(i=1;i<n+1;i++)for(j=i;j<n+1;j++){ if(i>1)y[i][j]=(y[i-1][j]-y[i-1][j-1])/(x[j]-x[j-i]);elsey[i][j]=(y[i-1][j]-y[i-1][j-1])/(x[j]-x[j-1]);printf("%f\n",y[i][i]);}temp=1;newton=y[0][0];for(i=1;i<n+1;i++){ temp=temp*(xx-x[i-1]);newton=newton+y[i][i]*temp;}printf("求得的结果为:N(%.4f)=%9f\n",xx,newton);牛顿插值法Matlab程序function f = Newton(x,y,x0)syms t;if(length(x) == length(y))n = length(x);c(1:n) = 0.0;elsedisp('x和y的维数不相等!');return;endf = y(1);y1 = 0;l = 1;for(i=1:n-1)for(j=i+1:n)y1(j) = (y(j)-y(i))/(x(j)-x(i));endc(i) = y1(i+1);l = l*(t-x(i));f = f + c(i)*l;simplify(f);y = y1;if(i==n-1)if(nargin == 3)f = subs(f,'t',x0);elsef = collect(f); %将插值多项式展开f = vpa(f, 6);endend牛顿插值法摘要:值法利用函数f (x)在某区间中若干点的函数值,作出适当的特定函数,在这些点上取已知值,在区间的其他点上用这特定函数的值作为函数f (x)的近似值。
常见的插值法及其应用
见的有分段线性插值和分段三次埃尔米特插值. 但是分段
插值光滑性较差.
2. 6 三次样条插值
工程上常用三次样条插值 ,其基本思想是将插值区间
n 等分 ,在每一个子区间上采用三次 Hermite 插值方法导
出插值函数 S3 ( x) (1) 在每个小区间 [ xi- 1 , xi ] 上是不高于三次的多项
= Ii
该公式我们也称为埃特金 ( Ait ken) 逐次线性插值公式. 这
个算法的优点是适合在计算机上计算 ,且具有自动选节点
并逐步比较精度的特点 ,程序也较简单.
2. 3 Newton 插值法
由表 (1) 构造的 Newton 插值多项式为
N ( x) = f ( x0 ) + f ( x0 , x1 ) ( x - x0 ) + …+ f ( x0 , x1 ,
7 YYSZXB
© 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved.
姜 琴 ,周天宏 : 常见的插值法及其应用
例 2 已知函数 f ( x) ,如表 6 所示 :
2 常见的代数插值公式及其构造
2. 1 Lagrange 插值法
表 (1) 的 n 次 Lagrange 插值多项式 L n ( x) 的数学公
n
∑ 式 :L n ( x) =
f ( xi) li ( x)
i =0
其中 li ( x) ( i = 0 , 1 ,2 , …n) 是插值基函数 ,且 li ( x) =
3 例题
例 1 已知数据 ,如表 3 所示 :
表 3 插值表
牛顿插值法的应用
牛顿插值法在处理磁化曲线和铁损曲线中的应用指导老师:李国霞院系:物理工程学院专业:物理电子学姓名:夏委委学号:201112131526一、牛顿插值法简介在科学研究与其他领域中所遇到的许多实际问题中,经常会出现函数不便于处理或计算的情形。
有时候函数关系没有明显的解析表达式,需要根据实验数据或其他方法来确定与自变量的某些值相对应的函数值;有时候函数虽有明显的解析表达式,但是使用很不方便。
因此,在实际应用中,往往需要对实际使用的函数建立一个简单的便于处理和计算的近似表达式,即用一个简单的函数表达式来近似替代原来复杂的函数。
与用近似数代替准确值一样,这也是计算法中最基本的概念和方法之一。
近似代替又称为逼近。
用多项式逼近列表函数的问题即为多项式插值问题。
根据函数)(x f 已有的数据表格来计算函数)(x f 在一些新的点x 处的函数值,这就是插值法所要解决的问题。
因此,所谓的插值法就是在所给定的函数表格中间在插入一些所需要的新的点上的函数值。
插值法的基本思想:首先设法根据表格中已有的函数值来构造一个简单的函数)(x y 作为)(x f 的近似表达式,然后再用)(x y 来计算新的点上的函数值作为)(x f 的近似值。
通常可以选用多项式函数作为近似函数)(x y ,因为多项式具有各阶的导数,求值比较方便。
用代数多项式作为工具研究插值问题,通常称为代数插值。
代数插值法问题的完整提法如下:设函数)(x f y =在区间[]b a ,上是连续的,且已知)(x f 在区间[]b a ,上1+n 个互异点处的函数值,即n i x f y i i ,......1,0),(== 其中,)(j i x x j i ≠≠。
寻找一个次数不高于n的多项式0111)(a x a x a x a x P n n n n n +++=-- 使满足条件n i x f x P i i n ,,1,0),()( ==称)(x P n 为)(x f 的插值多项式,),,1,0(n i x i =称为插值结点,[]b a ,称为插值区间。
基于Newton插值的光伏发电最大功率点跟踪算法研究
Ab tac sr t I t s ape . e n hi p r h m a h m a i l o o te t te t m del f h ph t ol c ca oov t cel i a alz an d ied。 en h ne ai l s n y ed s d er v t h t e w ma m u xi m p wer o pon ta i m eh b e on it r ckng tod as d Newt n o qu t it r lt i adr i n e poa i ac on s pr os d Th al r h op e e go i m i mor efi en a m o e ut t s e f ci t nd r s i — abe f n e hp l orsigl-c i prce s h L gr nge o s ort an a a qu adrt it polt al r hm. a i ner a i c on go i t Ke wors: y d MPPT,ol o i a alsi, w on p yn m al n y s Ne t qu t it polt adr i ner a i ac on
一
式 中 :h 光 生 电 流 ; - 极 管 反 向 饱 和 电 流 I、 光 伏 电 I一 。 I二 o lU一 池 工 作 电 流 、 作 电 压 ; 一 级 管 特 性 因 子 , 值 在 1 2之 间 ; 工 A 二 其 ~
k B l ma n常 量 ,值 为 13 x 0 / q - ot z n . 1 刊 K; 一电 子 的 电 子 电 荷 8
电力系统负荷预测中的数据采样与处理研究
电力系统负荷预测中的数据采样与处理研究随着电力行业的快速发展和能源需求的不断增长,电力系统负荷预测成为了一个重要的研究方向。
准确预测负荷可以帮助电网管理者合理规划电力资源和制定调度方案,有效降低电力系统的运行成本并确保供电的可靠性。
在该领域的研究中,数据采样与处理是至关重要的环节,决定了负荷预测模型的准确性和可靠性。
数据采样是电力负荷预测的第一步,其目的是收集与负荷变化相关的数据。
一般而言,负荷数据可以通过现场实测或者历史数据获取。
现场实测是指通过电力系统中的安装的传感器设备采集实时数据,如温度、压力、电流、电压等参数。
这种数据采集方式具有直接、实时性强的特点,能够提供高质量的数据样本,但是需要投入较大的人力和物力,并且传感器设备的安装和维护也是一个巨大的工程。
另一种数据采集方式是使用历史数据。
历史数据是过去一段时间内的负荷数据记录,可以通过电力系统的历史档案、仪表记录、数据库等方式获取。
相比于现场实测,历史数据的采集成本更低,但是采样频率和精度较低,可能无法满足高精度负荷预测的要求。
因此,在数据采样时需要权衡各种因素,选择合适的数据采集方式以满足负荷预测模型的需求。
数据处理是负荷预测的关键步骤,其目的是将采集的原始数据转化为可用于建模和分析的形式。
常见的数据处理方法包括数据清洗、数据平滑和特征提取等。
数据清洗是指在数据采集过程中,清除那些不符合要求或者不可靠的数据。
例如,在现场实测过程中,可能会存在传感器故障或者人为操作失误导致的异常数据,需要通过数据清洗方法将其排除。
数据平滑是指对预测数据进行平滑处理,消除其中的噪声和异常点。
常见的数据平滑方法包括移动平均法、指数平滑法和样条插值法等。
这些方法可以有效地消除数据中的尖峰和波动,使得负荷预测更加准确和稳定。
特征提取是指从原始数据中提取出与负荷变化相关的特征参数。
常见的特征参数包括平均值、方差、峰值、谐波等。
特征提取的目的是减少数据维度,提取出能够反映负荷特征的主要参数,并用于负荷预测模型的构建。
matlab编写拉格朗日插值法案例
matlab编写拉格朗日插值法案例拉格朗日插值法是一种常用的插值方法,可以通过已知的离散数据点,构造一个多项式函数来逼近这些数据点。
这种方法的最大优点是能够通过任意多个数据点来构造插值函数,因此适用于各种数据分布情况。
下面将列举10个使用Matlab编写拉格朗日插值法的案例,展示该方法在不同领域的应用。
1. 温度插值:假设我们有几个气象站的测温数据,但是某个地点的测温数据缺失。
可以使用拉格朗日插值法来估计该地点的温度值,从而填补数据缺失的空白。
2. 人口预测:根据已有的人口数据,可以使用拉格朗日插值法来预测未来某个时间点的人口数量。
这对城市规划和资源分配具有重要意义。
3. 股票价格预测:根据过去的股票价格数据,可以使用拉格朗日插值法来预测未来某个时间点的股票价格。
这有助于投资者做出合理的投资决策。
4. 图像处理:在图像处理中,可以使用拉格朗日插值法来放大或缩小图像。
通过对图像中的像素点进行插值,可以得到更高分辨率的图像。
5. 数字信号处理:在数字信号处理中,可以使用拉格朗日插值法来重构缺失的信号样本。
这有助于恢复被噪声污染或丢失的信号。
6. 机器学习:在机器学习中,拉格朗日插值法可以用于处理缺失数据。
通过对已有数据进行插值,可以构造更完整的数据集,提高模型的准确性。
7. 地理信息系统:在地理信息系统中,可以使用拉格朗日插值法来生成等高线图。
通过对不同高度的离散数据点进行插值,可以得到连续的等高线分布。
8. 数值积分:在数值积分中,可以使用拉格朗日插值法来近似计算定积分。
通过将被积函数进行插值,可以将积分转化为对插值多项式的求和。
9. 电力系统分析:在电力系统分析中,可以使用拉格朗日插值法来估计电力负荷曲线。
通过对已知的电力负荷数据进行插值,可以预测未来的负荷变化趋势。
10. 信号滤波:在信号滤波中,可以使用拉格朗日插值法来对信号进行平滑处理。
通过对信号进行插值,可以去除噪声和异常值,得到更干净的信号。
数值计算方法实验报告
数值计算方法实验报告实验目的:本次实验的目的是通过对数值计算方法的实践操作,加深对该方法的理解和掌握。
具体来说,本次实验旨在通过使用 MATLAB 软件对一些常见的数值计算问题进行求解,从而掌握和熟练运用一些数值计算方法,如插值、数值微积分、常微分方程数值解等。
实验过程:1.插值(1) Lagrange 插值法(2) Newton 插值法2.数值微积分(1) 梯形公式(2) Simpson 公式3.常微分方程数值解(1) 古典四步 Runge-Kutta 法(2) 改进四步 Runge-Kutta 法实验结果:本次实验中,我们使用 MATLAB 软件对以上数值计算问题进行了求解,成功得到了相应的数值解,并且通过分析和比较不同的数值计算方法的结果,得出了以下结论:1.在插值问题中,Lagrange 插值法和 Newton 插值法的结果相对较为接近,但是 Newton 插值法的计算速度更快。
2.在数值微积分问题中,梯形公式的结果较为精确,但是 Simpson 公式的精度更高。
3.在常微分方程数值解问题中,古典四步 Runge-Kutta 法和改进四步 Runge-Kutta 法均能得到较为准确的结果,但是改进四步Runge-Kutta 法的精度更高,尤其对于复杂的常微分方程求解有更好的效果。
实验结论:本次实验通过对数值计算方法的实践操作,深入理解了该方法的原理和运用,掌握了一些重要的数值计算方法,如插值、数值微积分、常微分方程数值解等,并且通过实验结果的分析比较,得出了相应的结论。
这些知识和技能对于我们在科研和工程实践中的数值计算问题具有非常重要的意义,具有广泛的应用前景。
解多元高次方程组求解方法与实际应用
解多元高次方程组求解方法与实际应用多元高次方程组是数学中的重要问题之一,它在科学研究和实际应用中具有广泛的应用价值。
本文将介绍一些常见的解多元高次方程组的方法,并探讨其在实际问题中的应用。
一、常见的解多元高次方程组的方法1.1 牛顿法牛顿法是解多元高次方程组的常用方法之一。
它基于泰勒级数展开,通过不断迭代逼近方程组的解。
牛顿法的主要思想是通过线性逼近和迭代求解,直到满足一定的收敛准则。
这种方法适用于近似解和初始解较好的情况下,能够快速求得精确解。
1.2 消元法消元法是解多元高次方程组的另一种常用方法。
它通过逐步消去未知数,将多元方程组化简为一元方程,最终求得解。
消元法主要包括高斯消元法和克拉默法则。
高斯消元法通过矩阵变换和行列式计算,将方程组转化为上三角矩阵,进而求得解。
克拉默法则则是通过行列式的计算,求得每个未知数的值。
消元法适用于方程数较少、未知数较多的情况,但计算过程较为复杂,有时会带来数值误差。
1.3 Lagrange插值法Lagrange插值法是解多元高次方程组的一种常见方法。
它基于拉格朗日插值多项式的思想,通过构造多项式逼近方程组的解。
Lagrange插值法可以精确求解多项式方程组,但计算过程相对复杂,计算量较大,在实际应用中通常用于近似求解。
二、解多元高次方程组的实际应用2.1 机器学习中的参数估计在机器学习中,参数估计是一个重要的问题,常常需要解多元高次方程组来求解最优参数。
例如,在线性回归模型中,需要求解最小二乘法问题,通过解多元高次方程组可以得到最优的回归系数。
解多元高次方程组的方法可以帮助我们快速准确地估计模型的参数,从而提高机器学习模型的性能和效果。
2.2 工程中的优化问题在工程领域中,常常需要优化问题的求解。
解多元高次方程组的方法可以应用于求解最优化问题。
例如,在电力系统中,优化电网输电能力和供电质量是一个重要问题,可以通过建立高次方程组模型,应用解多元高次方程组的方法来求解最优的电网配置和参数。
一种电力大数据的数据清洗方法
一种电力大数据的数据清洗方法电力大数据的数据清洗方法1. 引言电力行业是一个充满数据的行业,各种传感器、监测设备和系统都会产生大量的数据。
然而,这些数据往往存在着噪声、缺失值和异常值等问题,对数据的准确性和可靠性造成了挑战。
因此,开发一种有效的电力大数据清洗方法对于提高数据质量和准确性至关重要。
2. 数据清洗的重要性数据清洗是数据预处理的重要步骤,它可以帮助我们去除数据中的噪声、缺失值和异常值,从而提高数据的质量和可靠性。
清洗后的数据可以为后续的数据分析、建模和决策提供更准确和可靠的基础。
3. 数据清洗方法(1)去除噪声噪声是指数据中的不必要的干扰或错误信息,可能是由于设备故障、传感器误差或数据采集过程中的干扰引起的。
我们可以通过以下步骤去除噪声:- 使用滤波器:应用低通、高通或带通滤波器来去除频域上的噪声。
- 使用平滑算法:例如移动平均法、加权平均法等来平滑数据,减少噪声的影响。
(2)处理缺失值缺失值是指数据中的某些项或属性缺失的情况。
处理缺失值的方法包括:- 删除缺失值:如果缺失值的比例较小,可以直接删除包含缺失值的样本或属性。
- 插值法:使用插值算法(如线性插值、多项式插值、Kriging插值等)来填补缺失值。
- 使用默认值:对于某些属性,可以使用默认值来填补缺失值。
(3)检测和处理异常值异常值是指与其他数据点明显不同的数据点,可能是由于设备故障、测量误差或异常情况引起的。
处理异常值的方法包括:- 基于统计方法的异常值检测:使用统计方法(如均值、标准差、箱线图等)来检测异常值。
- 基于机器学习的异常值检测:使用机器学习算法(如聚类、离群点检测算法等)来检测异常值。
- 替换或删除异常值:可以将异常值替换为均值、中位数或删除包含异常值的样本。
4. 数据清洗流程(1)数据收集:收集电力大数据,包括传感器数据、监测设备数据等。
(2)数据预处理:包括数据清洗、数据集成、数据转换和数据规约等步骤。
(3)数据清洗:按照上述方法对数据进行噪声去除、缺失值处理和异常值检测和处理。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Newton插值法在电力系统采样中的应用
摘要:随着计算机及单片机的普及应用,对实际工程中的变量的测量准确性要求越来越高,由于计算机能够提供高速的计算速度,使一些应用数学中的数值分析法能够应用到实际工作中,本文就是针对newton插值法对电力系统采样中的实际应用进行探讨,介绍了newton插值法的优缺点,讨论newton插值法的应用方法,同时提供了部分计算过程。
关键词:数值分析,newton插值法,电力系统,采样点
1引言
针对实际工作中,对运算过程中的变量精度要求越来越高,如何获得准确的或满足精度要求的变量值越来越重要,随着计算机、单片机及程序化语言的出现,提供了高速计算的硬件环境,推进了应用数学中的数值分析方法的应用速度,通过对部分采样点的分析,能够得出比较准确的变量数值或变量曲线,newton插值法就是其中的一种数值分析方法,本文针对电力系统的采样要求,讨论newton 插值法的应用方法,分析其优缺点,为工程计算提供一种参考。
随着微机保护保护的不断应用,传统的电磁式模拟表记对电力系统中的主要变量电流i、电压u的采集精度已经不能满足精度的要求。
主要因为应用电磁式传统表计采集的电流i、电压u等变量,需要经复杂的转换(如通过采集器处理,再经过规约装换方式将变量转成数字信号)才能输入微机保护,大大降低了精度及可靠性,同时由于微机保护的应用,及计算机处理能力的提高,使用用微机
直接采样已经变为可能。
2newton插值法的应用
2.1 概述
由于电力系统中的主要变量u、i是时间的函数,即随时间变化而变化,想完全采到函数的所有点是不现实的,同时,由于运行监控的要求,需要得到变量u、i的曲线及有效值,这样就需要通过一定的算法,采集一定数量的点,在满足精度要求的前提下,计算出其他个点,绘制出变量u、i的曲线,计算出变量u、i的有效值,由于曲线的绘制,只是将变量的各点连接起来,根据有效值和最大值之间的关系,变量的有效值则是um/ (um是最大值)。
于是问题归为:通过一定的算法,采集一定数量的点,在满足精度要求的前提下,计算出变量的各个点。
2.2 newton插值法的具体应用
本文只讨论通过采集一定数量的点,通过newton插值法,计算出变量u、i的有效值。
要求有效值,必须先求出um,而um则是x= ?%i/2时的函数值,问题转化为:采集一定数量的点,在满足精度要求的前提下,通过newton插值法,计算出变量x= ?%i/2时,u( ?%i/2)的值即可。
只需求出对应参数即可,其对应参数为均差,均差可通过newton 插值表求出。
通过此公式可求变量u在任意点的数值,当然也可通过此公式画出变量u的曲线轨迹。
通过比较变量u的计算数值和实际的采样值,可得出变量的相对误差,也可通过比较变量u的曲线和实际采样中的曲线的重合度观测此公式的准确性,本文列举变量u和实际采样点的相对误差来分析newton插值法的优缺点,同时也为工程应用提供了一种简单实用的误差分析方法,来判别采样的准确性。
3newton插值法的误差分析
一般在实际应用中,变量u用到的都是有效值, 对变量u的有效值的实时监测,可以直接判别处变量u是否在标准范围内波动,从而决定是否对变量u的运行区间进行调整。
鉴于以上原因,本文对变量u的有效值相对误差进行分析。
当然也可对其他任意点的变量u的计算值和采样值进行比较,从而得出相对误差。
由于计算及分析方法相同,本文不再对其他点进行累述。
由于以上说明目的在于介绍newton插值法在采样中应用,所以采样点的个数较少,如果增加采样点的个数,就会将精度提高,此处不再论述。
4newton插值法的探讨
由newton插值法计算公式可以看出,newton插值法优点是计算较简单,尤其是增加节点时,每增加一个采样点,计算只要增加一项,如果采样等距节点(所谓等距节点,是指 x0,x1,...,xn中,相邻两点之间的距离都相等,也就是变量x的变化步长相等)。
这个相等的间距称为步长,记为h。
如果插值节点是等距的,那么整个插值公式将会出现新的规律,
毫无疑问会更加简单,同时公式简单直观,容易编程上机。
由于计算方法相同,此处不对等距newton插值公式进行探讨。
一般在电力系统采用中,一个周期采的点数不少于18个,此时误差将 <<5%,足可以满足精度要求。
5结论
本文对newton插值法在电力系统采样中的应用进行了探讨,对变量u的有效值的相对误差进行了分析,通过以上讨论,可以得出以下结论(1)newton插值法优点是计算较简单,公式直观,容易编程,适合应用于比较简单的采样计算问题,(2)newton插值法精度较高,采样节点越多,精度越高,如果采用等距newton插值法,可以进一步简化计算。
(3)newton插值法提供一种预测变量的手段,通过这种预测,可以判别变量是否在标准范围内运行
参考文献:
[1]张筑生.数学分析新讲[m].北京:北京大学出版社,1990
[2]郑咸义,姚仰新,雷秀仁,陆子强,黄凤辉.应用数值分析[m],广州:华南理工大学出版社,2008
[3]李庆扬,易大义,王能超,现代数值分析[m].北京:高等教育出版社,1995
[4]清华大学应用数学现代应用数学手册编写组.现代应用数学
应用手册(计算方法分册)[m].北京:北京出版社,1990
[5]蔡大用.数值分析[m].北京:清华大学出版社,1987
[6]林成森.数值分析[m].北京:科学出版社,2006
[7]施妙根.顾丽珍.科学和工程计算基础[m].北京:清华大学出版社,1999
[8]陈衍.电力系统稳态分析[m].北京:中国电力出版社,1995
[9]西安交通大学等合编.电力系统计算[m].北京:水利电力出版社,1978
[10]南京工学院主编.电力系统[m].北京:电力工业出版社,1980
[11]core lab, oracle class library for c++
3.00.0.1[db]./ocl/,2005
[12]邱关源.电路[m].北京:高等教育出版社,1997
作者简介:潘爽(1981-),女,吉林省长春市人,本科,助教,研究方向数学与应用数学;伊洋(1981-),男,吉林省长春市人,本科,助工,研究方向电力系统及其自动化、电力系统继电保护。