楚列斯基分解法解方程组matlab程序

简介

  • 01

    楚列斯基凤姐发只适用于系数矩阵A是对称正定的情况。

操作方法

  • 01

    先将矩阵A进行楚列斯基分解:A=R'R;

  • 02

    然后接R'y=b;

  • 03

    最后解Rx=y,得到原方程组的接。下面给出利用楚列斯基分解法求解线性方程组Ax=b的函数。

  • 04

    function x=solvebyCHOL(A,b)

  • 05

    % 该函数利用Cholesky分解法求线性方程组Ax=b的解 lambda=eig(A); if lambda>eps&isequal(A,A') [n,n]=size(A); R=chol(A); %解R'y=b y(1)=b(1)/R(1,1); if n>1 for i=2:n y(i)=(b(i)-R(1:i-1,i)'*y(1:i-1)')/R(i,i); end end %解Rx=y x(n)=y(n)/R(n,n); if n>1 for i=n-1:-1:1 x(i)=(y(i)-R(i,i+1:n)*x(i+1:n)')/R(i,i); end end x=x'; else x=[]; disp('该方法只适用于对称正定的系数矩阵!'); end

  • 06

    举例应用

  • 07

    >> A=[3 3 -3;3 5 -2; -3 -2 5]; >> b=[1 2 3]'; >> x=solvebyCHOL(A,b) x = 10/3 -2/3 7/3 >> A*x          %验证解的正确性 ans = 1 2 3

(0)

相关推荐

  • 侠盗猎车手圣安地列斯按键操作详解

    侠盗猎车手圣安地列斯中文版里会有大量的交通工具(从自行车一直到飞机,这已经是GTA的招牌啦)供玩家选择,不过在能使用某种交通工具前CJ必须考取相应的执照才行.游戏中的战斗要比以前更加复杂,近身战斗占到 ...

  • 用Matlab对矩阵进行LU分解法

    简介 01 利用矩阵分解来求先行方程组,可以节省内存,节省计算时间,因此在工程计算中最常用的技术.其中LU分解法是最基本也是最常用的方法. 操作方法 01 现将系数矩阵A进行LU分解,得到LU=PA: ...

  • 如何用matlab求解线性方程组的符号解?用matlab解符号方程组方法介绍

    如何用matlab解符号方程组?在学习线性代数或者其他应用是我们经常遇到要求解线性方程组,这是我们无法跳过的步骤,而matlab给我们提供了一个求解线性方程组的好方法,下面就简单介绍如何用matlab ...

  • 侠盗飞车圣安地列斯主角变成机器人

    侠盗飞车系列游戏因其灵活的操作,高度的自由性,得到了许多玩家的青睐,除此以外,更让人着迷的是它拥有的大量MOD,更是使得游戏的耐玩性得到了大大提高,许多人喜欢将游戏的主人公变成奥特曼,但那只是换了一身 ...

  • 圣安地列斯安卓秘籍/作弊码和安卓CLEO怎么安装

    国外大神 Alexander Blade.开发出了Android平台的CLEO,让手机版圣安地列斯玩家可以体验到PC平台的CLEO.圣安地列斯安卓版是取消了秘籍这一设定的,玩家无法通过键盘软件来激活秘 ...

  • 侠盗飞车圣安地列斯攻略

    <侠盗猎车手:圣安地列斯[3]>将游戏舞台由前作的佛罗里达州移到了加州,并收录了加州两大城市洛杉矶(Los Angeles )(游戏中改称 Los Santos).旧金山(San Fran ...

  • 英雄联盟钻石5补分详解!

    英雄联盟钻石5补分详解!从S5以后只要第一次上了钻石5,不管胜率多少,都会补分,也就是胜一局16分以下,输一局扣24分以上,最多达到赢一局只有10分,输一局扣30分,下面就给你们细说其中的原因. 操作 ...

  • 圣安地列斯任务攻略

    <侠盗猎车手:圣安地列斯[3]>将游戏舞台由前作的佛罗里达州移到了加州,并收录了加州两大城市洛杉矶(Los Angeles )(游戏中改称 Los Santos).旧金山(San Fran ...

  • 圣安地列斯攻略

    这款游戏的全称是<侠盗猎车手:圣安地列斯>,由美国开发商Rockstar North公司开发,2004年上市,是一款较为真实的模拟游戏. 这款游戏有着数以百计的交通工具,高自由度的游戏系统 ...