满意答案
用C语言设计克莱姆法则的计算程序
【摘要】 本文以线性代数中的克莱姆法则为依据,讨论其计算程序的设计方法,并给出在微机上运行的模拟人工解题的计算程序和应用示例。
关键词 克莱姆法则 程序
一、引 言
我们知道,在代数课的教学中无论是工科院校的线代数还是数学专业的高等代数都离不开克莱姆法则。由于该法则涉及到行列式的计算,在教学中只是挑选一些简单的例子进行教学,作业也比较少,致使学生练习不够充分。现在计算机已普及,但线性代数的教学并没有与算法语言的教学联系起来,而现有文献所研究的都是一些数值解法,从来不涉及象克莱姆法则这样的概念。现在已有按照行列式的定义计算行列式的计算程序[1,2],本文将在该程序的基础上进行扩展,编写出与克莱姆法则相对应的程序。
二、克莱姆法则简述
设A=(aij)是一个n×n矩阵,设D=det A≠0那么方程组
有唯一解
这里Dj是用方程组右边的常数项代替矩阵A中第j列所构成的矩阵的行列式。
三、程序设计步
从克莱姆法则知,首先要把方程组的系数aij及bi( i,j=1,…,n)输入计算机,计算出Dj(j=1,2,…,n)最后算出xj(j=1,2,…,n),具体方法述于下。
1. 方程的个数N及方程组的系数与数组间的关系简述
由于克莱姆法则的计算工作量是很大的,方程组的个数n也不可能取的很大,又由于C编译程序中二维数组的行是从0行开始,每行元素编号也从0开始,这样用二维数组a来代表数学中的(aij),其对应关系为a[i-1][j-1]对应aij (i,j=1,…,n)。
2. 计算D和Dj (j=1,…,n).