眩しさだけは、忘れなかった。|

合同变换法

合同变换法

已知二次型 f=xTAx,求变换 x=Py,使得二次型化为标准型 f=yTΛy,且 PTAP=Λ。该过程的实质是一次合同变换,即

[A,E]A,EA[Λ,PT]

具体的操作看下面几个例子。

一、实对称矩阵 A 对角元素均不为零

【例 1】将二次型 f(x1,x2,x3)=x12+5x22+5x32+2x1x24x1x3 化为标准型。

【解】由合同变换得

[112100150010205001]r2r1[112100042110205001]c2c1[102100042110225001]r3+2r1[102100042110021201]c3+2c1[100100042110021201]()r312r2[10010004211000052121]c312c2[10010004011000052121]

所以标准型为 y12+4y22,所作变换矩阵为 P=[11520112001],使 x=Py

若要求规范型,需对上述继续作合同变换,将 Λ 上的对角元素 a 化为 110,为此需作一次初等倍乘行变换(rn/a),再对应作一次初等倍乘列变换(cn/a)。

r2/4[1001000201212000052121]c2/4[1001000101212000052121]

所以规范型为 z12+z22,所作变换矩阵为 Q=[1125201212001],使 x=Qz

需要注意的是,合同变换的实质仍是配方,但配方法只是用了某种坐标变换,得到标准型的系数,不一定是特征值(不要以为使用该方法得到的 Λ 就是特征值!Λ 只能指示正、负和零特征值的个数,即正、负惯性指数一定是唯一的)。只有进行正交变换得到的系数才是特征值。

由此可知,二次型的标准型并不唯一,但是规范型唯一!如对(*)处还可作如下合同变换

()r3r2[100100021201042110]c3c2[100100012201024110]r32r2[100100012201000512]c32c2[100100010201000512]

所以标准型为 y12+y22,所作变换矩阵为 P=[125001012],使 x=Py

【例 2】(2014 年数二第 14 题)设二次型 f(x1,x2,x3)=x12x22+2ax1x2+4x2x3 的负惯性指数为 1,求 a 的取值范围。

【解】本题可使用配方法,但对于填空题来说比较麻烦。不妨采用合同变换法迅速解决本题。

[10a100012010a20001]r3ar1[10a10001201002a21a01]c3ac1[10010001201002a21a01]r3+2r2[100100012010004a21a21]c3+2c2[100100010010004a21a21]

因为负惯性指数为 1,所以 4a20,解得 2a2

二、实对称矩阵 A 对角元素有零

当发现二次型所对应的实对称矩阵 A 上的对角元素为 0 时,需要先想办法将对角线上的元素变成不为 0 的数,具体看下例。

【例 3】将二次型 f(x1,x2,x3)=2x1x2+2x1x3+2x2x3 化为标准型。

【解】发现对角线上第一个元素为 0,为使其不为 0,可将第二行加到第一行上,相应的就要作一次列变换,将第二列加到第一列上。对角线上其他位置为 0 的元素也是类似的处理方法。

[011100101010110001]r1+r2[112110101010110001]c1+c2[212110101010210001]r212r1[212110012012120210001]c212c1[202110012012120200001]r3r1[202110012012120002111]c3c1[200110012012120002111]

所以标准型为 2y1212y222y32,所作变换矩阵为 P=[11211121001],使 x=Py

若要求规范型,则继续进行变换

r1/2[20012120012012120002111]c1/2[10012120012012120002111]r2/12[10012120012012120002111]c2/12[1001212001012120002111]r3/2[1001212001012120002121212]c3/2[1001212001012120001121212]

所以规范型为 z12z22z32,所作变换矩阵为 Q=[1212121212120012],使 x=Qz

三、实战一道题

【例 4】(2021年数一张宇八套卷卷一第21题)已知实对称矩阵 A=[222a]B=[4331],其中 a 为正整数,求可逆矩阵 C,使得 CTAC=B

【解】由于 A 有未知参数,先对 B 进行合同变换。

[43103101]r234r1[4310054341]c234c1[4010054341]r1/2[20120054341]c1/2[10120054341]r2/52[1012005232525]c2/52[101200132525]

由此可知 B 的正、负惯性指数均为 1,变换矩阵为 C2=[12325025],使得 C2TBC2=[1001]。下面来对 A 进行合同变换。

[22102a01]r2r1[22100a211]c2c1[20100a211]

由于 CTAC=B,即 AB 合同,所以两者的正、负惯性指数相等,于是有 a2<0,又因为 a 为正整数,所以 a=1。继续对 A 进行合同变换得

[20100111]r1/2[201200111]c1/2[101200111]

由此可知变换矩阵为 C1=[12101],使得 C1TAC1=[1001]

因此有 C1TAC1=C2TBC2,即 (C1C21)TA(C1C21)=B,因此所求矩阵为

C=C1C21=[12101][232052]=[232254052]

本文作者:漫舞八月(Mount256)

本文链接:https://www.cnblogs.com/Mount256/p/17808710.html

版权声明:本作品采用CC 4.0 BY-SA许可协议进行许可。

posted @   漫舞八月(Mount256)  阅读(378)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 无需6万激活码!GitHub神秘组织3小时极速复刻Manus,手把手教你使用OpenManus搭建本
· Manus爆火,是硬核还是营销?
· 终于写完轮子一部分:tcp代理 了,记录一下
· 别再用vector<bool>了!Google高级工程师:这可能是STL最大的设计失误
· 单元测试从入门到精通
点击右上角即可分享
微信分享提示
评论
收藏
关注
推荐
深色
回顶
展开
  1. 1 Main Menu Theme Syd Matters
  2. 2 Luminous Memory (Acyanxi Remix) Acyanxi
  3. 3 夏影 麻枝准
  4. 4 潮騒の香り 水月陵
  5. 5 stand still 井口裕香 (いぐち ゆか)
  6. 6 流星雨 麻枝准
  7. 7 Summer Fantasy 傅许
  8. 8 失う 米白
  9. 9 epilogue 霜月はるか
  10. 10 夏に君を待ちながら 小原好美
  11. 11 桜のような恋でした 鹿乃 (かの)
  12. 12 風は微かに、熱を残し… 水月陵
  13. 13 夏凪ぎ 麻枝准/やなぎなぎ
  14. 14 空に光る 戸越まごめ
  15. 15 木漏れ日 riya
  16. 16 Songbirds Homecomings (ホームカミングス)
  17. 17 宝物になった日 麻枝准/やなぎなぎ
  18. 18 夏影~あの飛行機雲を超えた、その先へ~ 雪桜草 (雪樱草)
  19. 19 快晴 Orangestar (蜜柑星P),初音未来 (初音ミク)
  20. 20 永遠 霜月はるか
  21. 21 Sion 天門
  22. 22 遙かな年月-piano- 麻枝准
  23. 23 夏恋慕 kobasolo/春茶
  24. 24 夏凪ぎ-piano ver.- MANYO/麻枝准
  25. 25 Goodbye Seven Seas -piano ver.- MANYO/麻枝准
  26. 26 Light Years 麻枝准/やなぎなぎ
  27. 27 優しさの記憶 鹿乃 (かの)
Luminous Memory (Acyanxi Remix) - Acyanxi
00:00 / 00:00
An audio error has occurred, player will skip forward in 2 seconds.