sor迭代法matlab实例,SOR迭代法求线性方程组
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
%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迭代法求线性方程组相关推荐
- 不动点迭代法matlab视频,不动点迭代法matlab
不动点迭代法的 MATLAB 程序代码如下: Function [root,n]=StablePoint(f,x0,eps) %用不动点迭代法求函数的一个零点 %初始迭代向量:x0 %根的精度:eps ...
- sor迭代法matlab实例,Jacobi G-S SOR迭代法在matlab中例子.pdf
[数值分析报告] [数学基地班 赵晨晓 2011301000007] 关于Jacobi G -S SOR 方法的收敛速度比较 摘要:本论文主要通过判断比较Jacobi G -S SOR 三种 计算方法 ...
- 二元牛顿迭代法matlab,2-8牛顿迭代法matlab
4 实验七 牛顿迭代法 [实验目的] 1.了解牛顿迭代法的基本概念. 2.了解牛顿迭代法的收敛性和收敛速度. 3.学习掌握MATLAB 软件有关的命令. [实验内容] 用牛顿迭代法求方程0123=-+ ...
- picard迭代法matlab,从Picard迭代法来看欧拉等式
上周,我们在<方程本质上是一个条件等式>中提到了迭代,本文是该文的扩展. 温习前文,点击链接: 方程本质上是一个条件等式 算式,通常是恒等式,有唯一解. 方程,本质是条件等式.条件不同,方 ...
- jacobi matlab程序,jacobi迭代法实验MATLAB程序数值分析
jacobi迭代法实验MATLAB程序数值分析 例1. 求线性方程组 得近似解.精确解为x*=[3,2,1]'. 解:对方程进行移项就得 记为Ax=b,或写为x=B0 x+f,其中 取初始值,代入原方 ...
- 牛顿迭代法-matlab实现
牛顿迭代法-matlab实现 牛顿迭代法简介: 牛顿迭代法又称为切线法,简单来说就是不断求切线与x轴的交点,来逐渐接近解的迭代过程.方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的 ...
- 牛顿迭代法例题 matlab,牛顿迭代法-matlab程序(解线性方程组)
牛顿迭代法 matlab程序(解线性方程组) 作者:佚名来源:转载发布时间:2009-3-7 16:55:53 减小字体 增大字体 1.功能 本程序采用牛顿法,求实系数高次代数方程 f(x)=a0xn ...
- 牛顿迭代法解线性方程matlab程序,牛顿迭代法matlab程序(解线性方程组)
<牛顿迭代法matlab程序(解线性方程组)>由会员分享,可在线阅读,更多相关<牛顿迭代法matlab程序(解线性方程组)(4页珍藏版)>请在金锄头文库上搜索. 1.牛顿迭代法 ...
- c语言编程雅可比解方程,求雅可比迭代法解方程组的C\C++程序,急需一个运用雅可比迭代法求线性方程组的C/C++程序!...
问题标题 求雅可比迭代法解方程组的C\C++程序,急需一个运用雅可比迭代法求线性方程组的C/C++程序! 2019-7-8来自ip:14.137.150.56的网友咨询 浏览量:218 手机版 问题补 ...
最新文章
- Linux那些事儿之我是Sysfs(6)文件系统
- 2005年2月9日(星期三) 晚上,熱 - 大年初一,首次創作。
- Javascript编码规范,好的代码从书写规范开始,增强代码的可读性,可维护性,这是相当重要的!...
- 不同场景下MySQL的迁移方案
- Ubuntu LXC
- 删除右键文件下拉菜单(非打开方式)中的wine选项
- i标签content属性输入空白
- flutter持久化_开源项目【WanAndroid_Flutter】玩安卓 Flutter 版
- VMware:警惕 vSphere Web Client中的新漏洞
- iframe 滚动条不显示_BUG赏金 | 当我发现iFrame注入时的利用
- 解决wineQQ不能输入中文
- 设置Node.js脚本开机自启动
- 电力通信网与运营商5G网络融合方案
- JS边框图片基础知识
- Padded优化LinkedTransferQue并发性能是错误方向
- 2021 年,中台的“瓜”你还吃吗?
- 生活没那么复杂,“不要脸”就对了
- 如火如荼的人工智能现状
- setup factory打包工具的使用
- RoboCup中的机器人自主定位——从理论到实践(三):粒子滤波器
热门文章
- sounds音标软件Android,LearnEnglish Sounds Right
- 高级转录组分析和R数据可视化
- Perl 旁站查询(站长工具提取)
- Silverlight 3.0中文教程
- 【附源码】计算机毕业设计SSM网上拍卖系统
- iPad能装Android系统,我错了,原来iPad真的能装Windows和macOS系统!
- 如何在IPad上优雅地看移动硬盘中的视频
- Android - 框架使用
- 苏小红C语言程序设计第三版--设银行定期存款的年利率rate为2.25%,已知存款期为n年,存款本金为capital元,试编程计算并输出n年后的本利之和deposit
- ocaml-----hello world