计算方法
文章平均质量分 67
lf144090
这个作者很懒,什么都没留下…
展开
-
计算方法中方程的近似解法中二分法matlab实现
计算方法方程的近似解法 二分法fun.m中代码如下:function fun(a,b,e)%f是自定义的函数%a为隔根区间左端点,b为隔根区间右端点,e为绝对误差限if nargin==2 e=1.0e-6;elseif nargin input(变量输入错误!); return;endif a>=b input(隔根区间输入错误!); ret原创 2008-03-12 21:07:00 · 6159 阅读 · 0 评论 -
线性方程的Seidel迭代法
function X=Seidel(A,b,x0,e)% 线性方程的Seidel迭代法% A为系数矩阵n*n;b为右值n*1;x0为初值n*1,默认为全0;e为精度,默认为10^(-6)% BJ为Jacobi迭代矩阵,B1为BJ下三角部分,B2为BJ上三角部分% 其基本迭代公式为 X = inv(E-B1)*B2*X + inv(E-B1)*g%原创 2008-07-09 20:20:00 · 1564 阅读 · 0 评论 -
线性方程的Jacobi迭代法
function X=Jacobi(A,b,x0,e)% 线性方程的Jacobi迭代法% A为系数矩阵n*n;b为右值n*1;x0为初值n*1,默认为全0;e为精度,默认为10^(-6)% BJ为Jacobi迭代矩阵,其基本迭代公式为 X = B*X + g% 编程过程中要用到的几个迭代收敛的充要条件:% 1,BJ的谱半径小于1% 2,BJ的第一范式或无穷范式原创 2008-07-09 20:15:00 · 3349 阅读 · 1 评论 -
求拉格朗日插值多项式
function y=LagrangePolynomial(xi,yi)% 函数说明:求拉格朗日插值多项式% 参数说明:xi为观测点,yi为相应观测值,均为向量形式% 返回值说明:y为关于x的多项式%% 作者:野渡无人% 最后修改日期:2008.4.24%%>> xi=[0 1 2];%>> yi=[1 2 3];原创 2008-07-09 20:31:00 · 2616 阅读 · 0 评论 -
用乘幂法求矩阵的主特征值
function [lamda,uk]=CM(A,u0,e)% 用乘幂法求矩阵的主特征值% 参数说明:A为所求主特征值的矩阵,要求为方阵 n*n;u0为初始向量 n*1;e为误差限,缺省值为10^(-6).% 返回值说明:lamda为主特征值,uk为lamda对应的特征向量%% 作者:野渡无人% 最后修改日期:2008.4.24原创 2008-07-09 20:26:00 · 8813 阅读 · 2 评论 -
矩阵的三角分解法之LU分解之Doolittle分解
function [L,U]=Doolittle(A)% 矩阵的三角分解法之LU分解之Doolittle分解 A=LU% Doolittle分解:LU分解中L为单位下三角阵,U为上三角阵% 说明:n阶(n>=2)矩阵A有惟一杜里特尔分解的充要条件是:A的前n-1个顺序主子式不为0% 编程思想:使用待定系数法%% 作者:野渡无人% 最后修改日期:2008.4原创 2008-07-09 20:13:00 · 11808 阅读 · 0 评论 -
矩阵的三角分解法之LU分解之Crout分解
function [L,U]=Crout(A)% 矩阵的三角分解法之LU分解之Crout分解 A=LU% Crout分解:LU分解中L为下三角阵,U为单位上三角阵% 说明:n阶(n>=2)矩阵A有惟一Crout分解的充要条件是:A的前n-1个顺序主子式不为0% 编程思想:使用待定系数法%% 作者:野渡无人% 最后修改日期:2008.4.9%% >>原创 2008-07-09 19:47:00 · 12587 阅读 · 0 评论 -
计算方法总体选主元素消去法求线性方程组的解
function x=zongti(A,b)% 用总体选主元素消去法求线性方程组的解% A为系数矩阵(要求解惟一,则必须为满秩方阵 n*n),b为右值 n*1% x为所求的根 n*1% 作者:野渡无人% 最后修改日期:2008.3.27% 如:% A=[2 4 -2 -2;1 2 4 -3;-3 -3 8 -2;-1 1 6 -3],b=[-4 5 7 7]% x=zongti原创 2008-04-08 19:37:00 · 1543 阅读 · 0 评论 -
计算方法列主元素消去法求线性方程组的解
function x=leizhuyuansu(A,b)% 用列主元素消去法求线性方程组的解% A为系数矩阵(要求解惟一,则必须为满秩方阵 n*n),b为右值 n*1% x为所求的根 n*1% 作者:野渡无人% 最后修改日期:2008.3.27% 如:% A=[2 4 -2 -2;1 2 4 -3;-3 -3 8 -2;-1 1 6 -3],b=[-4 5 7 7]% 结果为:% x原创 2008-04-08 19:35:00 · 2735 阅读 · 1 评论 -
计算方法Gauss-Jordan消去法求线性方程组的解
function x=GaussJordan(A,b)% 用Gauss-Jordan消去法求线性方程组的解% A为系数矩阵(要求解惟一,则必须为满秩方阵 n*n),b为右值 n*1% x为所求的根%作者:野渡无人% 最后修改日期:2008.3.27% 例如教材45页例3.1:% A=[1 4 -5;1 3 -2;6 -1 18], b=[3 2 2]% x=GaussJordan(A,b原创 2008-04-08 19:33:00 · 4338 阅读 · 1 评论 -
计算方法双点截弦法matlab实现
function Shuang(fun,a,b,e)% 使用双点截弦法求方程 fun=0 在区间[x1,x2]处的根 e为误差限,默认为 eps% 对于任意的 x 属于 [a,b] , diff(fun,x) 与 diff(fun,x,2) 连续且不变号% 选择x0,x1均为动点% 其中迭代公式为 x_n=x_0-(x_1-x_0)/(fun(x_1)-fun(x_0))*fun(x_0)原创 2008-04-08 19:31:00 · 3454 阅读 · 0 评论 -
计算方法单点截弦matlab实现
function Dan(fun,a,b,e)% 使用单点截弦法求方程 fun=0 在区间[x1,x2]处的根 e为误差限,默认为 eps% 对于任意的 x 属于 [a,b] , diff(fun,x) 与 diff(fun,x,2) 连续且不变号% 选择x0为不动点,xn为动点, 要求 fun(x_0)*diff(fun,x_0,2) > 0% 其中迭代公式为 y=a-(x-a)/(f原创 2008-04-08 19:25:00 · 2259 阅读 · 0 评论 -
用matlab得到有根区间向量
function Xb=brackPlot(fun,xmin,xmax,nx)% BrackPlot find subintervals on x that contain sign changes of f(x)% Synopsis: Xb=brackPlot(fun,xmin,xmax,nx)% Xb=brackPlot(fun,xmin,xmax)% Input:原创 2008-03-25 20:54:00 · 1244 阅读 · 0 评论 -
用拉格朗日插值多项式求值
function yi=Lagrange(x,y,xi) % 函数说明:用拉格朗日插值多项式求值% 参数说明:x为已知值,y为x对应的函数值,x,y均为向量,xi为所求插值点% 返回值说明:yi为xi对应的函数值%% 作者:野渡无人% 最后修改日期:2008.4.24%%>> x=[1 2 3];%>> x原创 2008-07-09 20:28:00 · 1852 阅读 · 0 评论