matlab高斯消元,高斯消元法的MATLAB和PYTHON实现 | 学步园

超、凢脫俗 2023-01-15 03:51 157阅读 0赞

![Img\_line.gif][Img_line.gif]

![Img\_line.gif][Img_line.gif]

什么是高斯消元法？请见维基百科的定义

MATLAB实现：

function x = gauss(a，b)

n = length(b);

for i = 1 : n-1

for j = (i+1) : n

if a(i，i)~=0

lam = a(j，i)/a(i，i);

a(j，(i+1):n) = a(j，(i+1):n) - lam\*a(i，(i+1):n);

b(j) = b(j) -lam\*b(i);

end

end

end

for i = n:-1:1

b(i) = (b(i) - sum(a(i，(i+1):n).\*b((i+1):n)))./a(i，i);

end

x = b;

PYTHON实现：

import numpy

numpy import dot

def gauss(a，b):

n = len(b)

for i in range(0，n-1):

for j in range(i+1，n):

＃print a\[j，i\]

if a\[j，i\] != 0.0:

lam = float(a\[j，i\])/a\[i，i\]

a\[j，(i+1):n\] = a\[j，(i+1):n\] - lam\*a\[i，(i+1):n\]

b\[j\] = b\[j\] - lam\*b\[i\]

for k in range(n-1，-1，-1):

b\[k\] = (b\[k\] - dot(a\[k，(k+1):n\]，b\[(k+1):n\]))/a\[k，k\]

result = b

return result

x = numpy.matrix(\[\[2，1，-1\]，\[-3，-1，2\]，\[-2，1，2\]\]，dtype = numpy.float)

aa = numpy.matrix(\[\[2，1，-1\]，\[-3，-1，2\]，\[-2，1，2\]\]， dtype = numpy.float)

y = numpy.array(\[8，-11，-3\]，dtype = numpy.float)

bb = numpy.array(\[8，-11，-3\]， dtype = numpy.float)

c = gauss(aa，bb)

这次练习过程中发现的一个问题是关于矩阵变量的传递的问题，gauss(aa,bb)函数调用后，aa和bb的值都将会被修改，利用Python的id语句可以查看函数外的参数aa，bb和函数内的参数a和b的id，发现他们是一样的，这个问题应该怎么解决，怎样像C一样进行数值传递呢？目前能够想到的办法只有在函数中复制一下matrix，但是矩阵复制貌似和列表复制有所区别。列表的赋值使用b = a\[:\]就可以了，但是matrix不行…

终于找到答案了！利用matrix内置的method copy！

>>> x

matrix(\[\[ 2.,  1., -1.\],

\[-3., -1.,  2.\],

\[-2.,  1.,  2.\]\])

>>> y = x.copy()

>>> y

matrix(\[\[ 2.,  1., -1.\],

\[-3., -1.,  2.\],

\[-2.,  1.,  2.\]\])

>>> id(x)

32482800

>>> id(y)

39972208

list还是使用原来的方法进行复制，因为list里面没有copy这个method，在matrix或者dict中都有copy这个method用来实现真正的复制！

[Img_line.gif]: http://www.byywee.com/images/Img_line.gif

发表评论取消回复

表情：

评论列表（有 0 条评论，157人围观）

还没有评论，来说两句吧...

相关阅读

相关高斯消元 AcWing 883. 高斯消元解线性方程组

高斯消元 AcWing 883. 高斯消元解线性方程组原题链接 [AcWing 883. 高斯消元解线性方程组][AcWing 883.] 算法标签线性

川长思鸟来/ 2023年09月28日 13:11/ 0 赞/ 143 阅读

相关 matlab高斯消元,高斯消元法的MATLAB和PYTHON实现 | 学步园

![Img\_line.gif][Img_line.gif] ![Img\_line.gif][Img_line.gif] 什么是高斯消元法？请见维基百科的定义 MATL

超、凢脫俗/ 2023年01月15日 03:51/ 0 赞/ 158 阅读

相关【算法 | Python】高斯消元法

-------------------- 程序来源：[Gaussian Elimination][] -------------------- Arithmetic A

末蓝、/ 2022年12月03日 05:22/ 0 赞/ 208 阅读

相关线性代数——高斯消元法

文章目录先提问题高斯消元法简化模型先提问题我们先提出这样一个问题，对于如下这样一组方程，应该如何求出它的 x x x, y y y

￡神魔★判官ぃ/ 2022年09月16日 06:06/ 0 赞/ 220 阅读

相关高斯消元法的MATLAB代码

高斯消元法的原理详见于线性代数课本，其MATLAB代码如下： function[x]=gauss_elim(A,b) % A为系数矩阵且必须为方阵

港控/mmm°/ 2022年07月26日 00:21/ 0 赞/ 587 阅读

相关高斯消元训练

EXTENDED LIGHTS OUT <table> <tbody> <tr> <td><strong>Time Limit:</strong>

深碍√TFBOYSˉ_/ 2022年07月18日 05:27/ 0 赞/ 218 阅读

相关高斯消元

高斯消元法是用来对线性方程组进行求解的方法，模板如下 const int MAXN = 55; const double _inf = 1e-9;

Dear 丶/ 2022年06月12日 05:37/ 0 赞/ 246 阅读

相关 MATLAB之高斯消元法

第一种方法：文件名gauss1.m function x=gauss1(A,b) %参量说明：A,系数矩阵；B，常数列向量；zg，增广矩阵 %将增广矩

野性酷女/ 2022年05月15日 19:09/ 0 赞/ 245 阅读

相关高斯消元

https://www.luogu.org/problemnew/show/P3389 高斯消元用来解决n元一次方程组，复杂度为n ^ 3 流程和初中学的那个加减消元法差不

短命女/ 2022年01月07日 09:21/ 0 赞/ 325 阅读

相关高斯消元模板

　　[题目传送门][Link 1] 　　高斯消元其实很好理解，就是我们平常解线性方程组用的消元法，不过在代码中实现需要用到矩阵，这里就只放代码，思路网上的大佬们已经讲的

朱雀/ 2021年09月30日 15:28/ 0 赞/ 327 阅读