discrete.c

来自「C语言实现的矩阵的常用计算！包括：矩阵乘法」· C语言 代码 · 共 38 行

void      discrete(
			       Matrix   A,
				   Matrix   B,
			       Matrix   F,
				   Matrix   G,
			       double   T,
			       int      m,
				   int      n)
{
                   Matrix   D,E;
                   int      i,j;

				   for(i=0;i<m;i++)
					   for(j=0;j<m;j++)
						   D[i][j]=A[i][j];
				   for(i=0;i<m;i++)
					   for(j=0;j<n;j++)
						   D[i][j+m]=B[i][j];
                   for(i=m;i<(m+n);i++)
					   for(j=0;j<(m+n);j++)
						   D[i][j]=0;                       //D为(m+n)*(m+n)阶矩阵
                   MNmul(D,T,D,(m+n),(m+n));

				   
                   
				   expm(D,E,(m+n));
				   

                   for(i=0;i<m;i++)
					   for(j=0;j<m;j++)
						   F[i][j]=E[i][j];
				   for(i=0;i<m;i++)
					   for(j=0;j<n;j++)
						   G[i][j]=E[i][j+m];


				   
}