如果已经输入了一个方阵A 的全部数据,那么只须用命令 det(A) 就能求出A的行列式的值。对于一个特殊的n阶行列式的计算必须先分析行列式中数据分布的规律,用算法产生对应的矩阵A,然后用命令 det(A) 计算出具体的行列式的值。 例如对于n阶行列式 分析数据分布的规律可知,对应的矩阵主对角元素是 1、2、3、…、n,而其余的元素均为2。另外,当行列式的阶数不同时数据结果也不一样,所以n扮演参数的角色。请参考下面的程序段计算行列式的值 n=input('input n= ?'); A=2*ones(n,n); for k=1:n A(k,k)=k; end det(A) 这里,第一条语句是要求从键盘上输入一个n的值,第二条语句是创建全部元素为2的n阶方阵,接下来用 for …… end语句将方阵的主对角元修改为1、2、3、…、n,最后用求行列式命令计算并显示出该行列式的值。如果将n的值输入为5,计算结果为 -12;如果将n 的值输入为6,则计算结果为 -48;……,从多次试验的数据结果可以猜测出本题答案为 -2(n-2)!。但是,这样多次使用这段程序计算必须将程序作为一个方件保存在磁盘上,关于MATLAB的文件操作请参考问题 22 。
|