电力系统潮流计算及矩阵分解算法实现.pdfVIP

作业3

问题1

编写一个具有两个输入参数(A,ag)的子程序，该程序能为任何非奇异矩阵A生成输

出(Q,P)，使得若

⚫

flag=0,

一个

=LU,

P

=I

⚫

标志=1,

权限

=LU

其中L是下三角矩阵，U是对角线元素为1的上三角矩阵，且Q

=L+U-I解：

算法如下：

根据上述算法编写

LU.m

文件，具体代码见附录。

问题2

利用问题1的子程序，为以下矩阵A生成Q和P输出：

6−224

A[12−8410]

3−1333

−642−18

解：

在主程序中输入矩阵

A

后，运行[Q,P]LU(A,0)时，结果如下：

3

Problem1

Writeasubroutinewithtwoinputs(A,ag)thatwillgenerateforanynon-singular

matrixA,theoutputs(Q