二次型化标准形的三种方法

二次型化标准形的三种方法

将二次型化为标准形有利于我们了解二次型的简单形式、二次型的各种参数如正负惯性指数、得到二次型的规范形、对称矩阵合同的简单形等等。另外,化标准形也是解析几何化简二次曲线和二次曲面的需要。

下面,我们以两道题目为例说明计算二次型的标准形的三种方法:

  • 配方法
  • 合同变换法
  • 特征值法

第一题

例1 计算二次型XTAX=x12+5x22+5x32+2x1x2−4x1x3\displaystyle X^TAX=x_1^2+5x_2^2+5x_3^2+2x_1x_2-4x_1x_3XTAX=x12+5x22+5x32+2x1x24x1x3的标准形.

方法一:配方法

配方法的要领是:第一次将所有含有x1x_1x1的项集中到一起,进行配方,从而消掉含有x1x_1x1的交叉项,第二次将含有x2x_2x2的项集中到一起进行配方……直到去掉所有的交叉项.

解:将含有x1\displaystyle x_1x1的项集中起来进行配方:

XTAX=x12+2x1(x2−2x3)+5x22+5x32X^TAX=x_1^2+2x_1(x_2-2x_3)+5x_2^2+5x_3^2XTAX=x12+2x1(x22x3)+5x22+5x32

=(x1+x2−2x3)2−(x2−2x3)2+5x22+5x32=(x_1+x_2-2x_3)^2-(x_2-2x_3)^2+5x_2^2+5x_3^2=(x1+x22x3)2(x22x3)2+5x22+5x32

=(x1+x2−2x3)2+4x22+x32+4x2x3=(x_1+x_2-2x_3)^2+4x_2^2+x_3^2+4x_2x_3=(x1+x22x3)2+4x22+x32+4x2x3

=(x1+x2−2x3)2+(2x2+x3)2=(x_1+x_2-2x_3)^2+(2x_2+x_3)^2=(x1+x22x3)2+(2x2+x3)2

=y12+y22,=y_1^2+y_2^2,=y12+y22,

其中,非退化线性替换为:

{x1+x2−2x3=y12x2+x3=y2x3=y3,\begin{cases} x_1+x_2-2x_3&=&y_1\\ \quad \quad \quad 2x_2+x_3&=&y_2\\\quad \quad \quad \quad \quad \quad x_3&=&y_3\end{cases},x1+x22x32x2+x3x3===y1y2y3,

即,

{x1=y1−12y2+52y3x2=12y2−12y3x3=y3□\begin{cases}x_1&=&y_1-\frac{1}{2}y_2+\frac{5}{2}y_3\\x_2&=& \frac{1}{2}y_2-\frac{1}{2}y_3\\ x_3&=&y_3\end{cases}\quad\quad\squarex1x2x3===y121y2+25y321y221y3y3

方法二:合同变换法

解:
[11−2150−205100010001]→−r1+r2[11−2042−205100010001]\begin{bmatrix}1&1&-2\\1&5&0\\-2&0&5\\1&0&0\\0&1&0\\0&0&1\end{bmatrix}\xrightarrow{-r_1+r_2}\begin{bmatrix}1&1&-2\\0&4&2\\-2&0&5\\1&0&0\\0&1&0\\0&0&1\end{bmatrix}112100150010205001r1+r2102100140010225001

→−c1+c2[10−2042−2251−10010001]→2r1+r3[10−20420211−10010001]\xrightarrow{-c_1+c_2}\begin{bmatrix}1&0&-2\\0&4&2\\-2&2&5\\1&-1&0\\0&1&0\\0&0&1\end{bmatrix}\xrightarrow{2r_1+r_3}\begin{bmatrix}1&0&-2\\0&4&2\\0&2&1\\1&-1&0\\0&1&0\\0&0&1\end{bmatrix}c1+c21021000421102250012r1+r3100100042110221001

→2c1+c3[1000420211−12010001]→−12r2+r3[1000420001−12010001]\xrightarrow{2c_1+c_3}\begin{bmatrix}1&0&0\\0&4&2\\0&2&1\\1&-1&2\\0&1&0\\0&0&1\end{bmatrix}\xrightarrow{-\frac{1}{2}r_2+r_3}\begin{bmatrix}1&0&0\\0&4&2\\0&0&0\\1&-1&2\\0&1&0\\0&0&1\end{bmatrix}2c1+c310010004211002120121r2+r3100100040110020201

→−12c2+c3[1000400001−15201−12001]\xrightarrow{-\frac{1}{2}c_2+c_3}\begin{bmatrix}1&0&0\\0&4&0\\0&0&0\\1&-1&\frac{5}{2}\\0&1&-\frac{1}{2}\\0&0&1\end{bmatrix}21c2+c310010004011000025211

所以,标准形为,

XTAX=y12+4y22.X^TAX=y_1^2+4y_2^2.XTAX=y12+4y22.

非退化线性替换矩阵为,

C=[1−15201−12001].□C=\begin{bmatrix}1&-1&\frac{5}{2}\\0&1&-\frac{1}{2}\\0&0&1\end{bmatrix}.\quad\squareC=10011025211.

方法三:特征值法

解:令

A=[11−2150−205].A=\begin{bmatrix}1&1&-2\\1&5&0\\-2&0&5\end{bmatrix}.A=112150205.

由矩阵的特征多项式

∣λE−A∣=∣λ−1−12−1λ−5020λ−5∣|\lambda E-A|=\begin{vmatrix}\lambda-1&-1&2\\-1&\lambda-5&0\\2&0&\lambda-5\end{vmatrix}λEA=λ1121λ5020λ5

=2r2+r3∣λ−1−12−1λ−5002(λ−5)λ−5∣\overset{2r_2+r_3}{=}\begin{vmatrix}\lambda-1&-1&2\\-1&\lambda-5&0\\0&2(\lambda-5)&\lambda-5\end{vmatrix}=2r2+r3λ1101λ52(λ5)20λ5

=−2c3+c2∣λ−1−52−1λ−5000λ−5∣\overset{-2c_3+c_2}{=}\begin{vmatrix}\lambda-1&-5&2\\-1&\lambda-5&0\\0&0&\lambda-5\end{vmatrix}=2c3+c2λ1105λ5020λ5

按照第三行展开得,

=(λ−5)(λ2−6λ)=(\lambda-5)(\lambda^2-6\lambda)=(λ5)(λ26λ)

得到A\displaystyle AA的特征值为0,5,6.

此时,得到二次型XTAXX^TAXXTAX的标准形为,

XTAX=5y12+6y22.□X^TAX=5y_1^2+6y_2^2. \quad\quad \squareXTAX=5y12+6y22.

注:(1)由于二次型的标准形不唯一,所以三种方法得到的结果都不一样。
(2)由于只要求标准形,所以第三种方法所使用的正交变换我们省去不求了,到时专门就“使用正交变换化二次型为标准形”写一篇文章介绍。
在这里插入图片描述

第二题

例2 计算二次型XTAX=−4x1x2+2x1x3+2x2x3X^TAX=-4x_1x_2+2x_1x_3+2x_2x_3XTAX=4x1x2+2x1x3+2x2x3的标准形.

由于这道题目没有平方项,所以在使用配方法和合同变换法的时候需要先处理一下,以便化为例1的情况进行求解。特征值法与例1相同,所以不再赘述。

简要解答:

方法一:配方法

分析:为了得到一个平方项,可以将x1x_1x1设为两数之和,将x2x_2x2设为这两数之差.

解: 作非退化线性替换,

{x1=y1+y2x2=y1−y2x3=y3,\begin{cases}x_1=y_1+y_2\\x_2=y_1-y_2\\x_3=y_3\end{cases},x1=y1+y2x2=y1y2x3=y3,

XTAX=−4y12+4y22+4y1y3,X^TAX=-4y_1^2+4y_2^2+4y_1y_3,XTAX=4y12+4y22+4y1y3,

下面的步骤与例1同.

方法二:合同变换法

分析:为了将左上角变出一个非零元,可以考虑将第二行(或者第三行)加到第一行,同时将第二列(或者第三列)加到第一列.

解:

[0−21−201110100010001]→r3+r1[1−11−201110100010001]\begin{bmatrix}0&-2&1\\-2&0&1\\1&1&0\\1&0&0\\0&1&0\\0&0&1\end{bmatrix}\xrightarrow{r_3+r_1}\begin{bmatrix}1&-1&1\\-2&0&1\\1&1&0\\1&0&0\\0&1&0\\0&0&1\end{bmatrix}021100201010110001r3+r1121100101010110001

→c3+c1[2−11−101110100010101]\xrightarrow{c_3+c_1}\begin{bmatrix}2&-1&1\\-1&0&1\\1&1&0\\1&0&0\\0&1&0\\1&0&1\end{bmatrix}c3+c1211101101010110001

下面的步骤与例1同.


欢迎扫描下图中的二维码关注微信公众号:大哉数学之为用

欢迎关注微信公众号:大哉数学之为用
如果您觉得本文对您有帮助,欢迎赞赏!您的支持是作者继续下去的动力!在这里插入图片描述

### 将二次型化为标准型的 MATLAB 实现 在 MATLAB 中,可以通过特征值分解的方法二次型化为标准型。以下是具体的过程以及相应的代码示例: #### 1. 构造二次型矩阵 给定一个二次型 \( Q(x) = X^TAX \),其中 \( A \) 是对称矩阵,\( X \) 是变量向量。首先需要构建对应的对称矩阵 \( A \)[^1]。 ```matlab % 定义二次型系数矩阵A A = [a11, a12, ..., a1n; a21, a22, ..., a2n; ... an1, an2, ..., ann]; ``` #### 2. 特征值分解 通过 `eig` 函数计算矩阵 \( A \) 的特征值和特征向量[^2]。这一步骤用于寻找正交变换矩阵 \( P \) 和对角阵 \( \Lambda \),使得 \( A = P\Lambda P^{-1} \) 或者更简单地表示为 \( A = PDP' \)(因为 \( P \) 正交)。 ```matlab [V, D] = eig(A); ``` 这里,`V` 表示特征向量组成的矩阵,而 `D` 则是对角矩阵,其对角元即为 \( A \) 的特征值。 #### 3. 变换至标准型 令新的坐标系下的变量为 \( Y = V'^{-1}X \),则原二次型可以写成如下式的标准型: \[ Q(Y) = Y^TDY = d_1y_1^2 + d_2y_2^2 + ... + d_ny_n^2 \] 其中 \( d_i \) 即为矩阵 \( A \) 的第 \( i \) 个特征值。 ```matlab syms y1 y2 yn real; % 定义新变量 Y = [y1; y2; ... ; yn]; % 新的二次型表达式 Q_standard = Y'*D*Y; simplify(Q_standard); % 简得到最终结果 disp('标准后的二次型:'); disp(Q_standard); ``` 以上就是完整的流程说明及其对应的部分代码片段展示。 #### 注意事项 - 如果输入矩阵并非严格意义上的实对称矩阵,则需先将其转置部分平均处理后再继续后续操作。 - 对于数值稳定性较高的情况建议采用奇异值分解(SVD)替代传统意义上基于特征值的方式完成转换过程。
评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值