关于Matlab趋势计算program.doc文件的实现，爱文共享的信息具有丰富的相关文档，每天都有成千上万的行业名人在该站点上共享最新信息.

程序代码如下: 读取数据clcclearfilename ='txt'a = textread(filename)n = a(，)pinghengjd = a(，)phjddianya = a(，)jingdu = a(，)b = zeros (，)j = m，n = size(a)fori =: mifa(i，)== j = jb(j)= iendendb矩阵分区a = a(b(): b()b(),: n )a = a(b(): b(),: n)a = a(b(): b(),: n)a = a(b(): b(),: n)a = a( b(): b(),: n)设置初始值vcz = dcz = kmax = k =查找节点导纳矩阵a = zeros(，)fori =: forj =: a(i，j)= a(i ，j)a(，j)= a(，j)endalinei = a(: ，)linej = a(: ，)liner = a(: ，)linex = a(: ，)lineb = a(: ，)branchi = branchj = branchb = G = zeros(，)B = zeros(，)fork =: i = linei(k，)j = linej(k，)r = liner(k，)x = linex(k，)b = GIJ = r(r * rx * x)BIJ = x(r * rx * x)ifk> = lineb(k)〜= k = lineb(k)G(i，j)= GIJkG(j，i)= G (i，j)B(i，j)= BIJkB(j，i)= B(i，j)G(i，i)= G(i，i)GIJkkB(i，i)= B(i，i )BIJkkelseG(j，i)= GIJG(i，j)= G(j，i)B(j，i)= BIJB(i，j)= B(j，i)G(i，i)= G( i，i)GIJb = lineb(k)B(i，i)= B(i，i)BIJbendG(j，j)= G(j，j)GIJB(j，j)= B(j，j)BIJbendGBB = B * iYf = GBY = abs(Yf)alf = angle(Yf)分配雅可比矩阵参数P = zeros(n，)Q = zeros(n，)Pd = zeros(，n)Qd = zeros(，n)dP = zeros(，n)dQ = zeros(，n)PG = a(: ，)PD = a(: ，)QG = a(:, QD = a(: ，)i = a(: ，)forj =: 长度(i)P(i(j))= PG(i(j))PD(i(j))Q(i(j) )= QG(i(j))QD(i(j))enddelt =零(n，)V = ones(n，)V()= V()= ddelt = zer os(n，)dV = zeros( n，)A = zeros(* n，* n)B = zeros(* n，)Jacobian = Jaco(V，delt，n，Ymatlab潮流计算程序，j)= Jaco = HNJLJaco = zeros(* n，* n)fori =: nforj =: nJaco(* i ，* j)= H(i，j)Jaco(* i，* j)= N(i，j)Jaco(* i，* j)= J(i，j)Jaco(* i，* j)= L(i，j)endendA = Jaco数据:

本文来自电脑杂谈，转载请注明本文网址：

http://www.pc-fly.com/a/tongxinshuyu/article-262341-1.html