MATLAB方程求解:1.线性方程组

线性方程组

一.直接解法

例如:

方程组:

x+2y+3z=5

x+4y+9z=-2

x+8y+27z=6

可知:

A=[1 2 3;1 4 9;1 8 27]

b=[5;-2;6]

x=A\b

代码:

matlab 复制代码
A=[1 2 3;1 4 9;1 8 27];
b=[5;-2;6];
x=inv(A)*b;

线性方程组是线性代数研究的主要对象之一。

求解线性方程组的问题,不但在自然科学和工程技术中有所涉及,而且在数值计算方法的其它分支研究中,比如样条插值、最佳平方逼近、微分方程数值解,也往往需要求解线性方程组;

MATLAB软件最早也是起源于对线性方程组的研究;

直接解法

Ax=b

x=A^-1*b

MATLAB矩阵求逆函数:inv(A)

方程的解法为:x=inv(A)*b

从执行时间和数值准确性方面而言,一种更好的方法是使用矩阵反斜杠运算符,即:

x=A\b

这会使用高斯消去法来求解线性方程组,不必求逆矩阵。

二.矩阵分解法

矩阵分解法是指将线性方程组Ax=b分解为矩阵A的一些基本子式,然后利用这些子式来求解线性方程组。

LU分解法

  1. L,U\] = lu(A) LU分解法是将矩阵A分解为**下三角矩阵L和上三角矩阵U**,并使得L的对角线元素为1。注意:这里的矩阵A必须是方阵。

    产生一个上三角矩阵U和一个下三角矩阵L,并使得L的对角线元素为1。P是一个置换矩阵,使得PA=LU。这里的矩阵A必须是方阵。

LU分解:

实现LU分解后,线性方程组Ax=b的解为

x=U(L\b)

x=U(L\Pb)

这样可以大大提高运算速度。

QR分解法

QR分解:对矩阵X进行QR分解,就是把X分解为一个正交矩阵Q和一个上三角阵R的乘积形式;QR分解只能对方阵进行。

  1. MATLAB的函数qr可用于对矩阵进行QR分解,其调用格式为:

    Q,R\]=qr(X)

    Q,R,E\]=qr(X)

    x=R(Q\b)或x=E(R(Q\b))

三.迭代解法

  • 线性方程组的直接解法:精确解;对大型线性方程组的计算效率高:
  • 迭代法是一种不断用变量的旧值递推新值的过程,是用计算机解决问题的一种基本方法。它利用计算机运算速度快、适合做重复性操作的特点,让一组指令重复执行,在每次执行这组指令时,都从变量的原值推出它的新值。
  • 迭代解法非常适合求解大型稀疏矩阵的方程组。在数值分析中,迭代解法主要包括:
    • Jacobi迭代法;
    • Gauss--Serdel迭代法;
    • 超松弛迭代法;

线性方程组的迭代解法一Jacobii迭代法

于是Ax=b转化为:

x=D^-1(L+U)x +D^-1b

与之对应的迭代公式为:

x(k+1)=D-1(L+U)x(k)+D-1b

这就是Jacobi迭代公式。

如果序列{xk+1}收敛于x,则x必是方程Ax=b的解;

matlab 复制代码
function [y,n]=jacobi(A,b,x0,ep)
if nargin==3
    ep=l.0e-6;
elseif nargin<3
    error
    return
end

D=diag(diag(A)):%求A的对角矩阵
L=-tril(A,-1);%求A的下三角阵
U=-triu(A,1);%求A的上三角阵

B=D\(L+U);
f=D\b;
y=B*x0+f:
n=l;%迭代次数
while norm(y-x0)>=ep
    x0=y;
    y=B*x0+f;
end
相关推荐
Json____11 分钟前
使用python的 FastApi框架开发图书管理系统-前后端分离项目分享
开发语言·python·fastapi·图书管理系统·图书·项目练习
人生在勤,不索何获-白大侠32 分钟前
day16——Java集合进阶(Collection、List、Set)
java·开发语言
LIN-JUN-WEI1 小时前
[ESP32]VSCODE+ESP-IDF环境搭建及blink例程尝试(win10 win11均配置成功)
c语言·开发语言·ide·vscode·单片机·学习·编辑器
望获linux1 小时前
【Linux基础知识系列】第四十三篇 - 基础正则表达式与 grep/sed
linux·运维·服务器·开发语言·前端·操作系统·嵌入式软件
不吃肉的羊2 小时前
PHP设置文件上传最大值
后端·php
留不住丨晚霞2 小时前
说说SpringBoot常用的注解?
java·开发语言
hardStudy_h2 小时前
C++——内联函数与Lambda表达式
开发语言·jvm·c++
艾莉丝努力练剑3 小时前
【C语言】学习过程教训与经验杂谈:思想准备、知识回顾(三)
c语言·开发语言·数据结构·学习·算法
witton4 小时前
Go语言网络游戏服务器模块化编程
服务器·开发语言·游戏·golang·origin·模块化·耦合
枯萎穿心攻击5 小时前
ECS由浅入深第三节:进阶?System 的行为与复杂交互模式
开发语言·unity·c#·游戏引擎