文章目录Jacobi迭代法matlab程序(《数值分析原理》)1、Jacobi迭代格式2、Jacobi迭代法的例子Jacobi迭代法matlab程序(《数值分析原理》)1、Jacobi迭代格式Jacobi迭代法是常见的几种迭代法之一,迭代格式如下图所示:(图片来自CHD的ztl老师的PPT)(具体内容详见《数值分析原理》)2、Jacobi迭代法的例子该例子使用matlab的命令文件格式,命名为jacobi.m。举例:设有方程组取初始向量为x(0)=(-3,1,1)(T),用Jacobi方法求解,要求||x(k+1)-x(k)||小于等于10……(-3)。%题目信息%A为线性方程组的系数矩阵A=
Jacobi正交多项式一些基本概念Jacobi正交多项式的定义证明[-1,1]上的正交多项式是Jacobi正交多项式注:本文的内容主要根据文末中的参考文档[1,2,3]中的内容进行整理完成。一些基本概念设I=[−1,1]I=[-1,1]I=[−1,1]是实轴上的标准区间,定义在III上的正函数:ωα,β(x)=(1−x)α(1+x)β,α>−1,β>−1\omega_{\alpha,\beta}(x)=(1-x)^{\alpha}(1+x)^{\beta},\alpha>-1,\beta>-1ωα,β(x)=(1−x)α(1+x)β,α>−1,β>−1是权函数。赋权的Sobolev空间记为H
1.要求考虑线性方程组Hx=b,其中H为n阶Hilbert矩阵,即通过先给定解(例如取x的各个分量为1),再计算出右端向量b的办法给出一个精确解已知的问题.(1)分别编写DoolittleLU分解法、Jacobi 迭代、Gauss-Seidel 迭代的一般程序;(2)取阶数n=6,分别用LU分解法、Jacobi 迭代、Gauss-Seidel 迭代去求解上述的病态方程组Hx=b;分别报告它们的数值结果(包括数值解、迭代步数)以及它们在1-范数下的计算误差。迭代法的停止条件均取为2.Matlab实现(取迭代初值为0)2.1.1 LU分解函数function[L,U,y,x]=LU(A,b)%LU
雅可比方法该方法是求解对称矩阵全部特征值和特征向量的一种方法,它基于以下结论:①任何实对称矩阵A可以通过正交相似变换成对角型,即存在正交矩阵Q,使得QTAQ=diag(λ1,λ2,…,λn)Q^TAQ=diag(λ1,λ2,…,λn)QTAQ=diag(λ1,λ2,…,λn)其中λi(i=1,2,…,n)是A的特征值,Q中各列为相应的特征向量。②在正交相似变换下,矩阵元素的平方和不变。即设A=(aij)n×nA=(a_{ij})_{n×n}A=(aij)n×n,Q为正交矩阵,记B=QTAQ=(bij)n×nB=Q^TAQ=(b_{ij})_{n×n}B=QTAQ=(bij)n×n,则∑