该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

%SOR迭代法

function [y,k,err,w]=sor(eps,a,n)

h=1/n;

A=zeros(n-1,n-1); %定义系数矩阵A

for i=1:n-1

A(i,i)=-(2*eps+h);

end

for i=1:n-1

for j=1:n-1

if i==j+1

A(i,j)=eps;

end

if i==j-1

A(i,j)=eps+h;

end

end

end

b=zeros(n-1,1); %定义常数项b

for i=1:n-2

b(i,1)=a*h^2;

end

b(n-1,1)=a*h^2-eps-h; %考虑边值问题

D=diag(diag(A));

L=triu(A)-A;

U=tril(A)-A;

B=D\(L+U);

x=eig(B); %B的特征值

p=abs(max(x)); %B谱半径

w=2/(1+sqrt(1-p^2)); %最佳松弛因子

if w<0||w>2

disp('迭代不收敛!');

return;

end;

L=(D-w*L)\((1-w)*D+w*U);

f=w*(D-w*L)\b;

q=1-h*h;

delta=(1.0e-4)*(1-q)/q;

y=zeros(n-1,1);

z=zeros(n-1,1);

k=0;

while 1

z=L*y+f;

if norm(z-y,inf)

break;

end

y=z;

k=k+1;

end

x=[h:h:(n-1)*h];

true=(1-a)/(1-exp(-1/eps))*(1-exp(-x/eps))+x*a;

t=true';

err=norm(t-y,inf);

sor迭代法matlab实例,SOR迭代法求线性方程组相关推荐

  1. 不动点迭代法matlab视频,不动点迭代法matlab

    不动点迭代法的 MATLAB 程序代码如下: Function [root,n]=StablePoint(f,x0,eps) %用不动点迭代法求函数的一个零点 %初始迭代向量:x0 %根的精度:eps ...

  2. sor迭代法matlab实例,Jacobi G-S SOR迭代法在matlab中例子.pdf

    [数值分析报告] [数学基地班 赵晨晓 2011301000007] 关于Jacobi G -S SOR 方法的收敛速度比较 摘要:本论文主要通过判断比较Jacobi G -S SOR 三种 计算方法 ...

  3. 二元牛顿迭代法matlab,2-8牛顿迭代法matlab

    4 实验七 牛顿迭代法 [实验目的] 1.了解牛顿迭代法的基本概念. 2.了解牛顿迭代法的收敛性和收敛速度. 3.学习掌握MATLAB 软件有关的命令. [实验内容] 用牛顿迭代法求方程0123=-+ ...

  4. picard迭代法matlab,从Picard迭代法来看欧拉等式

    上周,我们在<方程本质上是一个条件等式>中提到了迭代,本文是该文的扩展. 温习前文,点击链接: 方程本质上是一个条件等式 算式,通常是恒等式,有唯一解. 方程,本质是条件等式.条件不同,方 ...

  5. jacobi matlab程序,jacobi迭代法实验MATLAB程序数值分析

    jacobi迭代法实验MATLAB程序数值分析 例1. 求线性方程组 得近似解.精确解为x*=[3,2,1]'. 解:对方程进行移项就得 记为Ax=b,或写为x=B0 x+f,其中 取初始值,代入原方 ...

  6. 牛顿迭代法-matlab实现

    牛顿迭代法-matlab实现 牛顿迭代法简介: 牛顿迭代法又称为切线法,简单来说就是不断求切线与x轴的交点,来逐渐接近解的迭代过程.方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的 ...

  7. 牛顿迭代法例题 matlab,牛顿迭代法-matlab程序(解线性方程组)

    牛顿迭代法 matlab程序(解线性方程组) 作者:佚名来源:转载发布时间:2009-3-7 16:55:53 减小字体 增大字体 1.功能 本程序采用牛顿法,求实系数高次代数方程 f(x)=a0xn ...

  8. 牛顿迭代法解线性方程matlab程序,牛顿迭代法matlab程序(解线性方程组)

    <牛顿迭代法matlab程序(解线性方程组)>由会员分享,可在线阅读,更多相关<牛顿迭代法matlab程序(解线性方程组)(4页珍藏版)>请在金锄头文库上搜索. 1.牛顿迭代法 ...

  9. c语言编程雅可比解方程,求雅可比迭代法解方程组的C\C++程序,急需一个运用雅可比迭代法求线性方程组的C/C++程序!...

    问题标题 求雅可比迭代法解方程组的C\C++程序,急需一个运用雅可比迭代法求线性方程组的C/C++程序! 2019-7-8来自ip:14.137.150.56的网友咨询 浏览量:218 手机版 问题补 ...

最新文章

  1. Linux那些事儿之我是Sysfs(6)文件系统
  2. 2005年2月9日(星期三) 晚上,熱 - 大年初一,首次創作。
  3. Javascript编码规范,好的代码从书写规范开始,增强代码的可读性,可维护性,这是相当重要的!...
  4. 不同场景下MySQL的迁移方案
  5. Ubuntu LXC
  6. 删除右键文件下拉菜单(非打开方式)中的wine选项
  7. i标签content属性输入空白
  8. flutter持久化_开源项目【WanAndroid_Flutter】玩安卓 Flutter 版
  9. VMware:警惕 vSphere Web Client中的新漏洞
  10. iframe 滚动条不显示_BUG赏金 | 当我发现iFrame注入时的利用
  11. 解决wineQQ不能输入中文
  12. 设置Node.js脚本开机自启动
  13. 电力通信网与运营商5G网络融合方案
  14. JS边框图片基础知识
  15. Padded优化LinkedTransferQue并发性能是错误方向
  16. 2021 年,中台的“瓜”你还吃吗?
  17. 生活没那么复杂,“不要脸”就对了
  18. 如火如荼的人工智能现状
  19. setup factory打包工具的使用
  20. RoboCup中的机器人自主定位——从理论到实践(三):粒子滤波器

热门文章

  1. sounds音标软件Android,LearnEnglish Sounds Right
  2. 高级转录组分析和R数据可视化
  3. Perl 旁站查询(站长工具提取)
  4. Silverlight 3.0中文教程
  5. 【附源码】计算机毕业设计SSM网上拍卖系统
  6. iPad能装Android系统,我错了,原来iPad真的能装Windows和macOS系统!
  7. 如何在IPad上优雅地看移动硬盘中的视频
  8. Android - 框架使用
  9. 苏小红C语言程序设计第三版--设银行定期存款的年利率rate为2.25%,已知存款期为n年,存款本金为capital元,试编程计算并输出n年后的本利之和deposit
  10. ocaml-----hello world