m0_51649630 2021-06-06 21:13 采纳率: 66.7%
浏览 780
已采纳

用matlab编程判断矩阵是否正定

matlab编写一个脚本,要求使用顺序主子式(不用特征值判断)判断一个矩阵是正定还是负定还是不定型。要在今晚十二点前编写出来。

  • 写回答

3条回答 默认 最新

  • CSDN专家-Matlab_Fans 2021-06-06 21:46
    关注

    直接上代码吧

    function Test()
    
    %% 案例1
    A1= eye(4)
    f1 = Fun_isPD(A1);
    switch f1
        case 1
            disp('正定')
        case -1
            disp('负定')
        otherwise
            disp('不定')
    end
    
    %% 案例2
    A2 = -A1
    f2 = Fun_isPD(A2);
    switch f2
        case 1
            disp('正定')
        case -1
            disp('负定')
        otherwise
            disp('不定')
    end
    
    %% 案例3
    A3 = diag([1 1 1 0])
    f3 = Fun_isPD(A3);
    switch f3
        case 1
            disp('正定')
        case -1
            disp('负定')
        otherwise
            disp('不定')
    end
    
    %% 案例4
    A4 = diag([1 -1 1 1])
    f4 = Fun_isPD(A4);
    switch f4
        case 1
            disp('正定')
        case -1
            disp('负定')
        otherwise
            disp('不定')
    end
    
    end
    
    function f = Fun_isPD(A)
    % f = Fun_isPD(A) 使用顺序主子式方法判断矩阵A的正定性
    % f = 1   A为正定时
    % f = 0   A不定时
    % f = -1  A为负定时
    
    n = length(A);   % 阶数
    f0 = sign(A(1));  % 一阶顺序主子式
    for ii = 2:n
        f = sign( det(A(1:ii,1:ii)) );
        if f <= 0
            if f == (-1)^ii
                f = -1;
                f0 = f;
                continue;
            end
            f = 0;
            break;
        else
            if f0 > 0
                f = 1;
            else
                f = -1;
            end
        end
        f0 = f;
    end
    
    end
    
    

    结果:

    A1 =
    
         1     0     0     0
         0     1     0     0
         0     0     1     0
         0     0     0     1
    
    正定
    
    A2 =
    
        -1     0     0     0
         0    -1     0     0
         0     0    -1     0
         0     0     0    -1
    
    负定
    
    A3 =
    
         1     0     0     0
         0     1     0     0
         0     0     1     0
         0     0     0     0
    
    不定
    
    A4 =
    
         1     0     0     0
         0    -1     0     0
         0     0     1     0
         0     0     0     1
    
    不定
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(2条)

报告相同问题?

悬赏问题

  • ¥30 关于移动Web网页使用TinyMCE富文本编辑器上传图片后的光标定位、压缩等几个问题如何解决:
  • ¥25 activeMq在同服务器centos8下消费很慢
  • ¥20 为什么zynq CAN IP 无法进入config配置模式,XCan_SelfTest函数失效?
  • ¥15 Pycharm中程序直接运行可以但进入调试报错
  • ¥15 MATLAB动图问题
  • ¥15 有段代码不知道怎么理解,const isToken = (config.headers || {}).isToken === false
  • ¥15 我的显卡支持CUDA最高版本是12.3,这个版本也支持VS 2022 17.0这种情况下如果我想下载CUDA11.8,需要下载旧版本的VS2022吗
  • ¥15 pointnet2包安装
  • ¥20 射频功率问题,解答者有酬谢!
  • ¥80 构建降雨和积水的预测模型