m0_71277349 2022-05-24 20:21 采纳率: 100%
浏览 61
已结题

有偿解答,MATLAB的

img


利用MATLAB解题,希望带解释。我理解不了。希望有人可以帮助我理解,附带MATLAB的运行图和过程。

  • 写回答

2条回答 默认 最新

  • WaitIKnowYou 2022-05-24 20:30
    关注

    举个栗子 高斯消元法解方程
    列主元消去法

    img

    解方程组

    img

    
    ②LU分解法解方程组
    
    function agui_lu(a,b)
    %求可逆矩阵a的Doolittle分解,1返回单位下三角矩阵,u返回上三角矩阵n=length(a);
    n=length(a);
    u=zeros(n,n);
    l=eye(n,n);
    u(1,:)=a(1,:);
    l(2:n,1)=a(2:n,1)/u(1,1);
        for k=2:n
            u(k,k:n)=a(k,k:n)-l(k,1:k-1)*u(1:k-1,k:n);
            l(k+1:n,k)=(a(k+1:n,k)-l(k+1:n,1:k-1)*u(1:k-1,k))/u(k,k);
        end
    l
    u
    %解Ly=b
    y=zeros(n,1);
    y(1)=b(1);
        for k=2:n
            y(k)=b(k)-l(k,1:k-1)*y(1:k-1);
        end
    %解Ux=y
    x=zeros(n,1);
    x(n)=y(n)/u(n,n);
        for k=(n-1):-1:1
            x(k)=(y(k) -u(k,(k+1):n)*x((k+1):n))/u(k,k);
    end
    x
    
    
    

    运行截图

    img

    方程你自己列嘛

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论 编辑记录
查看更多回答(1条)

报告相同问题?

问题事件

  • 系统已结题 6月5日
  • 已采纳回答 5月28日
  • 创建了问题 5月24日

悬赏问题

  • ¥15 Power query添加列问题
  • ¥50 Kubernetes&Fission&Eleasticsearch
  • ¥15 有没有帮写代码做实验仿真的
  • ¥15 報錯:Person is not mapped,如何解決?
  • ¥30 vmware exsi重置后登不上
  • ¥15 c++头文件不能识别CDialog
  • ¥15 Excel发现不可读取的内容
  • ¥15 关于#stm32#的问题:CANOpen的PDO同步传输问题
  • ¥20 yolov5自定义Prune报错,如何解决?
  • ¥15 电磁场的matlab仿真