- 博客(8)
- 收藏
- 关注
原创 平移反幂法
平移反幂法 求解距离p最近的特征值λ与单位特征向量xclear%% 程序说明% 平移反幂法求解距p最近的特征值及其对应的特征向量%% 输入端A=input('请输入目标矩阵A: \n');p=input('请输入平移单位p: \n');n=length(A);u=ones(n,1);pI=p*eye(n);k=0;% 精度选择disp('选择默认精度,请输入0;自定义请输入1');chr=input("");if (chr ~= 0)err=input('请输入指定精度:')
2021-12-28 21:41:56
387
原创 迭代法 求解线性方程组
迭代法 求解线性方程组 (MATLAB)统筹了 李庆扬《数值分析》第五版中关于求解Ax=b的四种常用迭代法一码多用Jacobi、Gauss-Seidel、SOR、SSOR四种迭代法可以自行选择迭代方法,自定义精度,选择收敛判定方案交互式软件般的体验代码如下cyc clear;%% 程序说明% 迭代法求解Ax=b. 可选求解器模型% A为n阶系数矩阵;b为常数项,err指定精度% R为迭代矩阵B的谱半径,k为迭代步数,x为最终解%% 输入端while 1A=input('请输入
2021-11-10 22:15:28
1837
原创 SSOR迭代法
SOR迭代法 求解线性方程组 (MATLAB)(数值分析作业自编)代码如下:*将代码分块,很容易理解。左端是输出项,右端输入自定义项function [R,k,x] = SSOR(A,b,error,w)% SSOR迭代% R为迭代矩阵B的谱半径;k为迭代步数;x为解向量% A为系数矩阵;b为常数项;error为指定相邻两次迭代精度;w为松弛因子%% w的输入控制while (1) if w>0 & w<=2 % 用谱半径判断收敛性 不满足
2021-11-10 17:27:03
1374
原创 SOR迭代法
SOR迭代法 求解线性方程组 (MATLAB)(数值分析作业自编)代码如下:*将代码分块,很容易理解。左端是输出项,右端输入自定义项function [R,k,x] = SOR(A,b,error,w)% Gauss-Seidel迭代% R为迭代矩阵B的谱半径;k为迭代步数;x为解向量% A为系数矩阵;b为常数项;error为指定相邻两次迭代精度;w为松弛因子%% w的输入控制while (1) if w>0 & w<=2 % 用谱半径判断
2021-11-09 22:14:37
2425
原创 Gauss-Seidel迭代法
Gauss-Seidel迭代法 求解线性方程组 (MATLAB)(数值分析作业自编)代码如下:*将代码分块,很容易理解。左端是输出项,右端输入自定义项function [R,k,x] = Gauss_Seidel(A,b,error)% Gauss-Seidel迭代% R为迭代矩阵B的谱半径;k为迭代步数;x为解向量% A为系数矩阵;b为常数项;error为指定相邻两次迭代精度%% 初始化n=length(b);x=zeros(n,1);x0=zeros(n,1);k=0;%%
2021-11-09 15:51:11
2329
原创 Jacobi迭代法
这里写`在这里插入代码片`自定义目录标题欢迎使用Markdown编辑器新的改变功能快捷键合理的创建标题,有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能,丰富你的文章UML 图表FLowchart流程图导出与导入导出导入欢迎使用Markdown编辑器你好! 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你
2021-11-09 11:56:18
1448
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人