计算方法与MATLAB

求解，谢谢，计算方法与MATLAB. abcdysjshdbxbsxx

参考链接

Matlab入门教程--基本运算与函数(一)_落叶_201306的博客-CSDN博客 Matlab入门教程--基本运算与函数(一)在MATLAB下进行基本数学运算，只需将运算式直接打入提示号（>>）之後，并按入Enter键即可。例如：>>(5*2+1.3-0.8)*10/25ans =4.2000MATLAB会将运算结果直接存入一变数ans，代表MATLAB运算後的答案（Answer），并显示其数值於萤幕上。（为简便起见，在下述各例中，我们不再印出MAT...

https://blog.csdn.net/qq_23747821/article/details/79675687?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522166729337116782417065473%2522%252C%2522scm%2522%253A%252220140713.130102334..%2522%257D&request_id=166729337116782417065473&biz_id=0&utm_medium=distribute.pc_search_result.none-task-blog-2~blog~baidu_landing_v2~default-2-79675687-null-null.nonecase&utm_term=%E8%AE%A1%E7%AE%97%E6%96%B9%E6%B3%95%E4%B8%8EMATLAB&spm=1018.2226.3001.4450

这个问题的回答你可以参考下: https://ask.csdn.net/questions/7640048

希望能帮到你
https://b23.tv/DEly15S

1、雅可比迭代法求解线性方程组

https://blog.csdn.net/fanjufei123456/article/details/127641829?csdn_share_tail=%7B%22type%22%3A%22blog%22%2C%22rType%22%3A%22article%22%2C%22rId%22%3A%22127641829%22%2C%22source%22%3A%22fanjufei123456%22%7D

function [y,n]=jacobi(A,b,x0,ep)
%%输出的参数 y指方程的解 n为迭代的次数
% 输入的参数分别是系数矩阵 右端列向量 迭代的初值 精度
D=diag(diag(A))%%求对角矩阵
L=-tril(A,-1);%%求下三角
U=-triu(A,1);%%求上三角
B=D\(L+U);
f=D\b;
y=B*x0+f;
n=1;
while norm(y-x0)>=ep %%用2范数去逼近
    x0=y;
    y=B*x0+f;
    n=n+1;
end

方程求解实现：

A = [1 5 -3; 5 -2 1; 2 1 -5];
b = [2; 4; -11];
x0 = [0; 0; 0];
eps = 1.0e-4;
[x, n] = jacobi(A,b,x0,eps);

2、第二种方法：高斯-赛德尔送代法求解线性方程组

function [y,n]=gauseidel(A,b,x0,ep)
D=diag(diag(A))%%求对角矩阵
L=-tril(A,-1);%%求下三角
U=-triu(A,1);%%求上三角
B=(D-L)\U;
f=(D-L)\b;
y=B*x0+f;
n=1;
while norm(y-x0)>=ep %%用2范数去逼近
    x0=y;
    y=B*x0+f;
    n=n+1;
end

求解：

A = [-8 1 1; 1 -5 1; 1 1 -4];
b = [1; 16; 7];
x0 = [0; 0; 0];
eps = 1.0e-3;
[x, n] = gauseidel(A,b,x0,eps);

结果展示：

3、全国大学生数学建模专栏
数学方法涉及优化、评价、预测、分类模型，全面解析数学建模比赛常用算法，原理+代码双重实现，全面助力数学建模竞赛，欢迎您的关注！！！（想获奖必关注）https://blog.csdn.net/fanjufei123456/category_11322258.html