Cauchy行列式的一种计算方法

称行列式$$\det A=\det\left(\frac{1}{a_{i}+b_{j}}\right)_{n\times n}$$为Cauchy行列式,我们来计算他:

由于$$\det\left(\frac{1}{a_{i}+b_{j}}\right)_{n\times n}=\frac{1}{\prod\limits_{1\leq i,j\leq n}(a_{i}+b_{j})}\det(c_{ij})_{n\times n}$$

如果将$a_{1},\cdots,a_{n},b_{1},\cdots,b_{n}$视作变量,那么每个$c_{ij}$都是一个$n-1$次的多元多项式,因此$\det(c_{ij})_{n\times n}$的结果必然是一个$n(n-1)$次的多元多项式$f(a_{1},\cdots,a_{n},b_{1},\cdots,b_{n})$(简记作$f$).

如果我们将$a_{1}$视作变量,而其余视作常数,那么显然如果存在$a_{k}=a_{1}$,必然有$$f=\det(c_{ij})_{n\times n}=0$$

从而$f$有因子$(a_{1}-a_{k}),k=2,3,\cdots,n$;以此类推可知$f$含有因子$$g=\prod_{1\leq i<j\leq n}(a_{i}-a_{j})(b_{i}-b_{j})$$

注意到$${\rm deg }g=2\binom{2}{n}={\rm deg}f$$

因此$f$和$g$仅相差一个非零常系数.即$$\det A=\frac{C\prod\limits_{1\leq i<j\leq n}(a_{i}-a_{j})(b_{i}-b_{j})}{\prod\limits_{1\leq i,j\leq n}(a_{i}+b_{j})}$$我们来求这个系数$C$,现令$$a_{i}=\frac{1}{2}+ix,b_{j}=\frac{1}{2}-jx$$,那么$$\det A=\det\left(\frac{1}{1+(i-j)x}\right)_{n\times n}$$

注意到上式的结果对于充分大的$x$是连续的,从而我们令$x\to\infty$可知$$\det A=1$$

而此时\begin{align*}\lim_{x\to\infty}\frac{\prod\limits_{1\leq i<j\leq n}(a_{i}-a_{j})(b_{i}-b_{j})}{\prod\limits_{1\leq i,j\leq n}(a_{i}+b_{j})}&=\lim_{x\to\infty}\frac{\prod\limits_{1\leq i<j\leq n}(i-j)x(j-i)x}{\prod\limits_{1\leq i,j\leq n}(1+(i-j)x)}\\&=\frac{-\prod\limits_{1\leq i<j\leq n}(i-j)^2}{\prod\limits_{1\leq i,j\leq n,i\neq j}(i-j)}=1\\\Rightarrow C&=1\end{align*}

综上便知$$\det\left(\frac{1}{a_{i}+b_{j}}\right)_{n\times n}=\frac{\prod\limits_{1\leq i<j\leq n}(a_{i}-a_{j})(b_{i}-b_{j})}{\prod\limits_{1\leq i,j\leq n}(a_{i}+b_{j})}$$

时间: 2024-10-29 13:38:20

Cauchy行列式的一种计算方法的相关文章

求一个数组当中最大(最小)值的两种计算方法

求一个数组当中最大(最小)值的两种计算方法 1.常用方法(以求最大值为例) 1 #include "stdafx.h" 2 #include <string> 3 using namespace std; 4 int _tmain(int argc, _TCHAR* argv[]) 5 { 6 int a[5] = { 1, 23, 2, 6, 7 }; 7 int array_length = sizeof(a) / sizeof(a[0]);//数组大小 8 int m

HR人力资源软件考勤的两种计算方法

说起考勤,最明显的特点莫过于他灵活的规则. 每家企业基本上都各不相同,所以一直以来考勤的计算在劳动密集型企业都是件繁杂的差事.人力部门到月底往往要面对成千上万的考勤明细记录.因此现在企业一般都是采用软件计算的方法.而软件计算的方法也很多,下面通过案例来介绍两种使用比较普遍的考勤计算的方法. 案例: 某月有31天,A员工共出勤18天,公休8天,事假2天,病假1天,法定假2天,平时加班30小时,迟到30分钟.   计算方法一处理步骤: 手工安排员工班次 自动采集员工打卡记录 手工登记事假,病假申请单

【Java】斐波那契数列(Fibonacci Sequence、兔子数列)的3种计算方法(递归实现、递归值缓存实现、循环实现)

斐波那契数列:0.1.1.2.3.5.8.13………… 他的规律是,第一项是0,第二项是1,第三项开始(含第三项)等于前两项之和. > 递归实现 看到这个规则,第一个想起当然是递归算法去实现了,于是写了以下一段: public class RecursionForFibonacciSequence { public static void main(String[] args) { System.out.println(recursion(10)); } public static double

罗固P5824十二种计算方法(小球计算)

priproblem ++ i + 1n + 1有序集合链接: 如果有:(n \\\\)个球,请将它们放在\\\\(m \\\\)框中,然后找到解决方案的数量. \\\\(\\\\文字{I ..).球是不同的,盒子是不同的. 显然,答案是\\\\(m ^ n :). \\\\(\\\\文字{II ..).球彼此不同,盒子彼此不同,每个盒子最多容纳一个球. 我们将球一个接一个地放在框中,\\\\(i-i + 1 \\\\)在放置球((i \\\\))时有多种选择,因此答案\\\\(m ^:\\\\

n阶行列式

1.首先,在学习n阶行列式定义前,我们先介绍一下排列.逆序数及对换的概念: (1)排列定义:由 1,2,-,n 组成的一个有序数组称为一个 n 级排列(也叫做这 n 个元素的一个全          排列). 例子:所有的 3 级排列: 123     132     213     231     312     321 排列的分类:分为奇排列和偶排列:逆序数是奇数的排列称为奇排列:逆序数是偶数或 0 的排列            称为偶排列 (2)逆序数的定义:在一个排列中,如果一对数的前后

几种压缩算法实现原理详解

gzip .zlib以及图形格式png,使用的压缩算法都是deflate算法.从gzip的源码中,我们了解到了defalte算法的原理和实现.我阅读的gzip版本为 gzip-1.2.4.下面我们将要对deflate算法做一个分析和说明.首先简单介绍一下基本原理,然后详细的介绍实现. 1 gzip 所使用压缩算法的基本原理 gzip 对于要压缩的文件,首先使用LZ77算法的一个变种进行压缩,对得到的结果再使用Huffman编码的方法(实际上gzip根据情况,选择使用静态Huffman编码或者动态

求逆元的四种算法(拓欧费马小线性推欧拉)

求逆元的四种算法 拓展欧几里得算法求逆元 上一篇博客中已经讲过拓展欧几里得算法,并且讲解了求逆元的原理.这里只列出代码 在要求逆元的数与p互质时使用 代码 //扩展欧几里得定理 int ex_gcd(int a,int b,int& x,int& y) { if(b==0) { x=1; y=0; return a; } int ans = ex_gcd(b,a%b,x,y); int tmp = x; x = y; y = tmp-a/b*y; return ans; } int cal

一根电线承载功率准确的计算方法!4平方的铜芯电线能带的功率是5500W,2.5平方能带多少呢?

原文: https://baijiahao.baidu.com/s?id=1601671032391176338&wfr=spider&for=pc 家居杂坛依然首先给出答案:4平方的铜芯电线能带的功率是5500W. 关于这个问题,家居杂坛以前跟大家分析过,在这里再跟大家好好说说一根电线能带多少千瓦,怎么计算?大家可以看看,以后装修会用的上,基础公式是P=UI,P要求的功率,单位是W:U电压,单位是V;I计算电流,单位是A. 以4平方电线为例来进行计算,首先也是需要设定电线的使用条件:家装

高等代数葵花宝典—白皮书

我第一次接触白皮书是和高中同学钟梓源(复旦数学学院 16 级)的交流当中发现的,记得是上半期开学之后,钟梓源给我发了几张他们高等代数"练习册"的照片,还记得是矩阵的 Kronecker 积和摄动法之类的,当时大为惊讶,我虽然对这些东西有所耳闻,但均来源于杂乱无章的各种资料.照片上的几道题均有十足的难度,我当下决定买一本,便开始阅读了起来.现在已经过去大半年,高等代数已经学完了最后一章的内容,白皮书也陪伴我度过了一个半学期.现在,我大概把白皮书上的题目全部做了一遍,不得不惊叹于白皮书乃一