m0_63360141 2021-11-13 19:38 采纳率: 83.3%
浏览 50
已结题

用matlab编写一个方程,输入一个n × n对称正定矩阵A,输出一个n × n矩阵L,计算公式如下

用matlab编写一个方程,L=myfunction(A)实现输入一个n × n对称正定矩阵A,输出一个n × n矩阵L,计算公式如下:

img

  • 写回答

1条回答 默认 最新

  • CSDN专家-Matlab_Fans 2021-11-13 19:56
    关注

    可以直接用循环实现。

    
    function Test()
    
    A = [1 2 3;3 4 5;6 7 8]
    L = myfunction(A)
    
    
    function L = myfunction(A)
    
    n = length(A);
    
    for jj = 1:n
        for kk = 1:n
            if jj < kk
                L(jj,kk) = 0;
            elseif jj == kk
                sumtemp = 0;
                for ii = 1:kk-1
                    sumtemp = sumtemp + L(kk,ii)^2;
                end
                L(jj,kk) = A(kk,kk) - sumtemp;
            else
                sumtemp = 0;
                for ii = 1:kk-1
                    sumtemp = sumtemp + L(jj,ii)*L(kk,ii);
                end
                L(jj,kk) = (A(jj,kk) - sumtemp) / L(kk,kk);
                
            end
        end
    end
    

    结果:

    A =
    
         1     2     3
         3     4     5
         6     7     8
    
    
    L =
    
        1.0000         0         0
        3.0000   -5.0000         0
        6.0000    2.2000  -32.8400
    
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 请采纳用户回复 5月26日
  • 系统已结题 11月21日
  • 已采纳回答 11月13日
  • 创建了问题 11月13日

悬赏问题

  • ¥15 maya的mel里,怎样先选择模型A,然后利用mel脚本自动选择有相同名字的模型B呢。
  • ¥15 Python题,根本不会啊
  • ¥15 Python题,不会啦
  • ¥15 LD衰减结果文件dist的单位?
  • ¥15 Python题,回答一下下啦
  • ¥15 利用Java语言结合数据库实现-个信息管理系统(如学生管理系统、图书管理系统(可以加上借书还书等功能等)或者小游戏
  • ¥15 会会信号与系统和python的来
  • ¥15 关于#python#的问题
  • ¥20 oracle RAC 怎么配置啊,配置
  • ¥15 excel 日常使用中出现问题