matlab高斯消元法列主元,matlab用列主元高斯消去法解方程组,
一、
实验名称:
线性方程组高斯消去法.
二、
实验目的:
进一步熟悉理解
Guass
消元法解法思路,提高
matlab
编程能力.
三、
实验要求:
已知线性方程矩阵,
利用软件求解线性方程组的解.
四、
实验原理
:
消元过程:
设
0
)
0
(
11
a
,令乘数
)
0
(
11
)
0
(
1
1
/
a
a
m
i
i
,做(消去第
i
个方程组的
i
x
)
操作
1
i
m
×第
1
个方程
+
第
i
个方程(
i=2
,
3
,
.n
)
则第
i
个方程变为
1
)
1
(
2
)
1
(
2
...
i
n
in
i
b
x
a
x
a
这样消去第
2
,
3
,
.
.
.
,
n
个方程的变元
i
x
后.原线性方程组变
为:
)
1
(
)
1
(
2
)
1
(
2
)
1
(
2
)
1
(
2
2
)
1
(
22
)
0
(
1
)
0
(
1
1
)
0
(
11
...
.
.
...
...
n
n
nn
n
n
n
n
n
b
x
a
x
a
b
x
a
x
a
b
x
a
x
a
这样就完成了第
1
步消元.
回代过程:
在最后的一方程中解出
n
x
,得:
)
1
(
)
1
(
/
n
nn
n
n
n
a
b
x
再将
n
x
的值代入倒数第二个方程,
解出
1
n
x
,
依次往上反推,
即
可求出方程组的
其通项为
3,...1
-
n
2,
-
n
k
/
)
(
)
1
(
1
)
1
(
)
1
(
k
kk
n
k
j
j
k
kj
k
k
k
a
x
a
b
x
五、实验内容:
function
maintest2
clc
clear
all
A=[1 3 4;2 4 5;1 4 6];
%
系数矩阵
b=[1 7 6]'
%
常数项
num=length(b)
for
k=1:num-1
for
i=k+1:num
if
A(k,k)~=0
l=A(i,k)/A(k,k);
A(i,:)=A(i,:)-A(k,:).*l;
b(i)=b(i)-b(k)*l;
end
end
end
A
b
%
回代求
x
x(num)=b(num)/A(num,num);
for
i=num-1:-1:1
sum=0;
for
j=i+1:num
sum=sum+A(i,j)*x(j);
end
x(i)=(b(i)-sum)/A(i,i);
end
x
End
六、实验结果:
A =
1.0000
3.0000
4.0000
0
-2.0000
-3.0000
0
0
0.5000
b =
1.0000
5.0000
7.5000
x =
16
-25
15
matlab高斯消元法列主元,matlab用列主元高斯消去法解方程组,相关推荐
- MATLAB 五对角矩阵 Jacobi迭代法 SOR迭代法 解方程组
% ----------------------------------------------------------------------------------------------- % ...
- 高斯消去法解线性方程组的fortran程序实现
高斯消去法解方程组的fortran程序实现 许多实际问题的解决,常常要化为求解线性代数方程组:例如,用最小二乘法处理测量结果和用差分法求解偏微分方程时,都会得到线性方程组:同时,很多物理学的问题最后也 ...
- matlab 高斯消去法程序,高斯消去法解线性方程的Matlab程序
1151091杨晨辉 高斯消去法解线性方程的Matlab程序 方法一: function x = gauss(A,b) n = length(b); for k = 1 : n-1 if A(k,k) ...
- matlab解方程组方法,第二章解线性方程组的直接方法matlab用法
第二章解线性方程组的直接方法matlab用法 第二章 解线性方程组的直接方法的 MATLAB 程序24. 在这章中我们要学习线性方程组的直接法,特别是适 合用数学软件在计算机上求解的方法. 2.1 方 ...
- 高斯消元法的MATLAB实现
这是一个基于最大主元的高斯消元法的matlab实现,代码中并未考虑对方程组是否有解以及解的唯一性的判断,具体原理可参考高等代数或<MATLAB数学建模>. function s=Gua ...
- matlab中 三种方法计算 Ax b,在MATLAB中,方程Ax=B的解可以用哪个命令求得? matlab 求助 解方程组...
matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=inv(A)*b - 采用求逆运算解方程组: (2)x=A\B - ...
- 【matlab】解方程组
在matlab中解方程组是很方便的 例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=inv(A)*b - 采用求逆运算解方程组: (2)x=A\b - ...
- MATLAB阶段性方程组,[转载]matlab 解方程组
1.解方程 最近有多人问如何用matlab解方程组的问题,其实在matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=i ...
- matlab水汽计算公式,[转载]matlab 解方程组
1.解方程 最近有多人问如何用matlab解方程组的问题,其实在matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=i ...
- MATLAB机器人机械臂运动学正逆解、动力学建模仿真与轨迹规划
MATLAB机器人机械臂运动学正逆解.动力学建模仿真与轨迹规划,雅克比矩阵求解.蒙特卡洛采样画出末端执行器工作空间 基于时间最优的改进粒子群优化算法机械臂轨迹规划设计 ID:4610679190520 ...
最新文章
- 双十一来这儿,华为昇腾的秘密都给你!
- 数据结构源码笔记(C语言):统计字符串中出现的字符及其次数
- 推荐系统的实践与思考(上篇)【转】
- scorecardpy库的使用简介
- 使用javassist动态注入代码
- 利用 Git OpenSSH 查看/生成 本机 ssh 公钥
- python如何输出两列数据_如何用python将一列数据分为两列?
- python3之后版本读取网页的内容
- 简述php语言的特点是_PHP语言有哪些优势和特点(一)
- 思科三层交换机开启ipv6路由功能_三层交换机实现路由功能配置示例与详解 (Cisco Packer Tracer 模拟器)...
- 程序员代码大全c语言,程序员大神教你,新手零基础学C语言编程代码训练
- GridView边框样式简单美化
- 一句批处理快速把优酷KUX转MP4
- 作品交流:锁相环环路滤波器系数、NCO增益单位、鉴相器输出之间的关系
- protal认证的稳定性测试脚本
- PHP之支付宝APP支付
- 概率论知识回顾(二):古典概型,几何概型
- 基本概念学习(9013)---通用寄存器、机器字长、数据通路
- 在Windows10环境下安装RabbitMAQ、Erlang的坑
- 64位win7系统下安装USB下载器驱动FriendlyArm-usb-dnw-driver-的解决方案
热门文章
- python大神年薪_我程序员年薪 80 万被亲戚鄙视不如在二本教书的博士生?
- windows设置定时任务(win10任务计划程序)
- 【领英如何一键批量有效地加好友?】
- 由对称性知定点一定在x轴上_圆锥曲线中的定点定值问题的四种模型.doc
- 《居里夫人自传》读后感
- 3 万亿美元,苹果问鼎全球市值最高公司宝座
- 计算机随机抽样的方法有,卫生统计学四种随机抽样方法
- 体验+营销+云,Adobe越来越“不务正业”了?
- 常用0x000000类型颜色代码表
- python爬楼梯递归_爬楼梯(Python3)