正定二次型

定义:对于任意非零向量x=[x1,x2...xn]T,恒有

f(x1,x2...xn)=\sum_{i=1}^{n}\sum_{j=1}^{n}a_{ij}x_{i}x_{j}>0

则称二次型f为正定二次型

f正定

\LeftrightarrowA的正惯性指数p=r=n,r为矩阵的秩,n为未知数个数

\LeftrightarrowA的每一个特征值都大于0

\Leftrightarrow A\simeq E

\Leftrightarrow A=D^{T}D,D为可逆矩阵

\LeftrightarrowA的全部顺序主子式大于0,即

此外,若C=\begin{vmatrix} A & 0\\ 0 & B \end{vmatrix}为正定矩阵,则A^{T}=A,B^{T}=B,AB均为正定矩阵

  • 2
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
要判断一个二次是否为正定二次,可以使用以下步骤: 1. 将二次化为矩阵形式:设二次为 $Q(x_1,x_2,\cdots,x_n)=\sum_{i=1}^n\sum_{j=1}^na_{ij}x_ix_j$,则对应的矩阵为 $A=(a_{ij})$。 2. 判断矩阵 $A$ 是否为对称矩阵,若不是,则该二次不是正定二次。 3. 计算 $A$ 的所有顺序主子式,即 $D_k=\begin{vmatrix} a_{11} & a_{12} & \cdots & a_{1k} \\ a_{21} & a_{22} & \cdots & a_{2k} \\ \vdots & \vdots & \ddots & \vdots \\ a_{k1} & a_{k2} & \cdots & a_{kk}\end{vmatrix}$,其中 $k=1,2,\cdots,n$。 4. 如果所有顺序主子式 $D_k$ 的符号均为正,则该二次正定二次;如果 $D_k$ 的符号交替出现,则该二次为不定二次;如果存在某个 $D_k$ 的符号为负,则该二次为负定二次。 在MATLAB中,可以使用以下代码来实现上述步骤: ```matlab % 定义二次的系数矩阵 A = [1 2 3; 2 5 6; 3 6 9]; % 判断矩阵是否为对称矩阵 if ~isequal(A, A') disp('该二次不是正定二次'); return; end % 计算顺序主子式 n = size(A, 1); D = zeros(1, n); for k = 1:n D(k) = det(A(1:k, 1:k)); end % 判断二次 if all(D > 0) disp('该二次正定二次'); elseif all(mod(1:n, 2) == 1 & D < 0) || all(mod(1:n, 2) == 0 & D > 0) disp('该二次为不定二次'); elseif any(D < 0) disp('该二次为负定二次'); end ``` 其中,`isequal(A, A')`用于判断矩阵 $A$ 是否为对称矩阵,`det(A(1:k, 1:k))`用于计算顺序主子式 $D_k$,`all(D > 0)`表示所有顺序主子式的符号均为正,`all(mod(1:n, 2) == 1 & D < 0) || all(mod(1:n, 2) == 0 & D > 0)`表示顺序主子式的符号交替出现,`any(D < 0)`表示存在某个顺序主子式的符号为负。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值