MATLAB
文章平均质量分 70
密函一封
可以不断进步乃人生一大幸事!
展开
-
关于MATLAB函数不可用问题。??? Undefined function or method 'minfo' for input arguments of type 'double'.
使用MATLAB过程中,需要使用vrho函数求矩阵谱半径,但运行此函数出现错误提示:??? Undefined function or method 'minfo' for input arguments of type 'double'.程序版本导致函数缺失,解决办法,借助开源代码,自己建了vrho.m函数文件。原创 2015-11-05 17:01:58 · 12579 阅读 · 1 评论 -
数值分析 追赶法求解三对角线性方程组 MATLAB实现
函数主体部分编程算法 参考 数值分析 第四版 颜庆津 P27运行结果截图:%追赶法求解三对角线性方程组,Ax=b,A用一维数组a,c,d存储。function [L,U,x]=crout(a,c,d,b)%数组a存储三角矩阵A的主对角线元素,c、d存储主对角线上边下边带宽为1的元素 n=length(a); n1=length(c);原创 2015-11-05 11:48:53 · 41046 阅读 · 5 评论 -
数值分析 jacobi迭代法求解线性方程组 MATLAB程序实现
数值分析 jacobi迭代法MATLAB程序实现 函数内容(保存为jacobi.m文件,在MATLAB中运行):%function [G,d,x,N]=jacobi(A,b)%Jacobi迭代法求解线性方程组。迭代公式x(k+1)=-Gx(k)+d,(k=0,1,…).%返回迭代矩阵G,d和方程组的解x,以及迭代次数N原创 2015-11-06 19:01:28 · 35583 阅读 · 4 评论 -
数值分析 Gauss-Seidel迭代法求解线性方程组 MATLAB程序实现
Gauss-Seidel迭代法 参考数值分析第四版 颜庆津著 P39运行输入运行结果函数内容(保存为gauss.m文件,在MATLAB中运行)%Gauss-Seidel迭代法求解线性方程组。迭代公式x(k+1)=Gx(k)+d,(k=0,1,…).%返回迭代矩阵G,d和方程组的解x,以及迭代次数N原创 2015-11-06 19:36:45 · 20118 阅读 · 5 评论 -
数组分析 jacobi法求实对称矩阵全部特征值与特征向量 MATLAB实现
Jacobi迭代法求解实对称矩阵的特征值与特征向量。MATLAB实现输入:eps允许的误差限,A实对称矩阵;输出:lam是n个特征值,U的列向量对应特征向量原创 2015-11-19 19:48:31 · 8948 阅读 · 4 评论 -
数值分析 幂法求矩阵A按模最大的特征值和相应的特征向量
%function [x,y]=jingmi(A,u0,ep)%x特征值,y特征向量,u0初始向量,ep允许误差。%幂法求矩阵A按模最大的特征值和相应的特征向量%参考数值分析第四版 颜庆津 P49%2015.11.7 密密编写 (*^__^*) 嘻嘻……function [x,y]=jingmi(A,u0,ep)%x特征值,y特征向量%xx=1:10原创 2015-11-07 17:53:12 · 12683 阅读 · 4 评论 -
数值分析 反幂法求矩阵按模最小特征值 MATLAB实现
值分析第四版 颜庆津 计算实习题P238结果:%function [lam]=jingmi(ep)%反幂法求矩阵按模最小特征值%2015.11.8 密密编写 (*^__^*) 嘻嘻…… function [lam]=jingmi(ep)%lam是A的按模最小的特征值%初始化 n=501; a=1:n;%存储A的主对角线元素 f原创 2015-11-11 18:29:05 · 11127 阅读 · 1 评论 -
选主元doolittle分解法求解n元线性方程组 MATLAB实现
算法参考数值分析第四版 颜庆津 P23 (选主元的doolittle分解法算法)运行结果截图:输入矩阵A,b。输出QA=LU的分解矩阵和Ax=b的解。A=[1 8 2 3;-6 -3 8 1;2 4 4 2;10 5 -5 6]b=[12;40;-50;80]doolittle.m文件函数内容:function [L,U,x] = doolittle(A,b)%原创 2015-11-04 22:18:08 · 8253 阅读 · 5 评论