formangle.cpp

电力系统发输电系统基于非序贯的蒙特卡罗抽样的风险评估程序

#include<math.h>
#include<stdlib.h>
#include<stdio.h>

#include"datastruct.h"
//给定相乘矩阵的维数a(m,n)*b(M,N)=c(m,N)
void formAngle(int m,int n,int M,int N,double *a,double *b,double *c,int Busnum,
			   struct Busdata *Busdata,int *balanceno)
{
     int i,j,h;
	 double *d;
	 d=new double[m*N];
     if(n!=M)
	 {
	    printf("error!");
	    exit(-1);
	 }
     //矩阵相乘
	 for(i=0;i<m;i++)
		 for(j=0;j<N;j++)
		 {
		    d[i*N+j]=0;
		    for(h=0;h<n;h++)
				d[i*N+j]+=a[i*n+h]*b[h*N+j];
		 }

      for(i=0;i<Busnum;i++)
	  {
	      if(i<*balanceno)
			  c[i]=d[i];
		  else
		  if(i==*balanceno)
			  c[i]=0.0;
	      else
			  c[i]=d[i-1];
          //将结果赋给节点数据
		  Busdata[i].angle=c[i];
	  }
     delete []d;
	 d=NULL;

}