找回密码
 欢迎注册
查看: 94|回复: 10

[提问] 求一个 n 阶行列式的值 2024 - 10 - 13

[复制链接]
发表于 3 天前 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
\(D_n = \begin{vmatrix}  
1+x & x & x & \cdots & x & x \\
x & 2+x & x & \cdots & x & x \\
\vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\
x & x & x & \cdots & x & n+x \end{vmatrix}\)

这个行列式的值为 \(D_n=a x+n!\) ,其中\(a\)  是一个正整数,其值如何求出?
经实际计算知:
\(D_2 = 3x + 2\)
\(D_3 = 11x + 6\)
\(D_4 = 50x + 24\)
\(D_5 = 274x + 120\)
\(D_6 = 1764x + 720\)
\(D_7 = 13068x + 5040\)
\(D_8 = 109584x + 40320\)
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 3 天前 来自手机 | 显示全部楼层
$n!H_n$

点评

nyy
思路是什么?归纳?还是?  发表于 前天 09:00
nyy
Hn是调和级数求和吗?  发表于 前天 08:55
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
回复

使用道具 举报

发表于 前天 08:57 | 显示全部楼层
D4中
a=4!*(1/1+1/2+1/3+1/4)=50,根据mathe的意思,应该是这样
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 前天 08:59 | 显示全部楼层
  1. n = 8; Sum[1/k, {k, 1, n}]*n!
复制代码

得到结果
109584
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 前天 11:48 | 显示全部楼层
  1. Table[{n,Det[Table[x+k*Boole[k==c],{k,n},{c,n}]]},{n,20}]
复制代码


1        1+x
2        2+3 x
3        6+11 x
4        24+50 x
5        120+274 x
6        720+1764 x
7        5040+13068 x
8        40320+109584 x
9        362880+1026576 x
10        3628800+10628640 x
11        39916800+120543840 x
12        479001600+1486442880 x
13        6227020800+19802759040 x
14        87178291200+283465647360 x
15        1307674368000+4339163001600 x
16        20922789888000+70734282393600 x
17        355687428096000+1223405590579200 x
18        6402373705728000+22376988058521600 x
19        121645100408832000+431565146817638400 x
20        2432902008176640000+8752948036761600000 x

点评

nyy
你的代码比我的好  发表于 前天 13:20
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 前天 13:18 | 显示全部楼层
northwolves 发表于 2024-10-14 11:48
1        1+x
2        2+3 x
3        6+11 x

  1. Clear["Global`*"];(*mathematica11.2,win7(64bit)Clear all variables*)
  2. (*求出行列式的值*)
  3. aa=Table[{n,Det[Table[x+k*Boole[k==c],{k,n},{c,n}]]},{n,20}]
  4. (*求出行列式值的x项的系数*)
  5. bb=Coefficient[#[[2]],x]&/@aa
  6. (*求出通项公式函数*)
  7. cc=FindSequenceFunction[bb,n]//FullSimplify
复制代码


借你的锅下面

求出结果
行列式的值是
{{1, 1 + x}, {2, 2 + 3 x}, {3, 6 + 11 x}, {4, 24 + 50 x}, {5,
  120 + 274 x}, {6, 720 + 1764 x}, {7, 5040 + 13068 x}, {8,
  40320 + 109584 x}, {9, 362880 + 1026576 x}, {10,
  3628800 + 10628640 x}, {11, 39916800 + 120543840 x}, {12,
  17280 (27720 + 86021 x)}, {13, 17280 (360360 + 1145993 x)}, {14,
  241920 (360360 + 1171733 x)}, {15,
  3628800 (360360 + 1195757 x)}, {16,
  29030400 (720720 + 2436559 x)}, {17,
  29030400 (12252240 + 42142223 x)}, {18,
  1567641600 (4084080 + 14274301 x)}, {19,
  1567641600 (77597520 + 275295799 x)}, {20,
  156764160000 (15519504 + 55835135 x)}}

x的系数是
{1, 3, 11, 50, 274, 1764, 13068, 109584, 1026576, 10628640, \
120543840, 1486442880, 19802759040, 283465647360, 4339163001600, \
70734282393600, 1223405590579200, 22376988058521600, \
431565146817638400, 8752948036761600000}

通项表达式是n! HarmonicNumber[n]
\[n! H_n\]

点评

nyy
mathematica太智能了,连这通项表达式都能找到,真不知道是如何找到的!  发表于 前天 13:21
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 前天 17:04 | 显示全部楼层
i) 增加一行一列,第一行为1 x x ... x
第一列除了第一行为1,其余都是0,所以行列式值不变。
ii)  每行减去第一行,那么除了第一行,第一列,主对角线后其它元素都是0.
    主对角线全部常数,第一行除了第一格全x,第一列除了第一个全-1
iii) 将第k行 x/(k-1)加到第一行 (k>=2), 矩阵变为对角阵,第一格1+H_n *x

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-10-16 16:19 , Processed in 0.024365 second(s), 18 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表