MyTDMA 数值分析P183

原创 2012年03月22日 13:06:47
function x=MyTDMA
%首先说明:追赶法是适用于三对角矩阵的线性方程组求解的方法,并不适用于其他类型矩阵。
%定义三对角矩阵A的各组成单元。方程为Ax=f
%    b为A的对角线元素(1~n),a为-1对角线元素(2~n),c为+1对角线元素(1~n-1)。
%    A=[2  -1  0    0
%       -1  3  -2   0
%       0  -2  4   -3
%       0  0   -3   5]
a=[0 -1 -2 -3];b=[2 3 4 5];c=[-1 -2 -3];f=[6 1 -2 1];
%消去过程
n=length(b);
u=zeros(n,1);
l=zeros(n,1);
y=zeros(n,1);

u(1)=b(1);
y(1)=f(1);
for i=2:n
    l(i)=a(i)/u(i-1);
    u(i)=b(i)-l(i)*c(i-1);
    y(i)=f(i)-l(i)*y(i-1);
end
%回带过程
x=zeros(n,1);
x(n)=y(n)/u(n);
for i=(n-1):-1:1
    x(i)=(y(i)-c(i)*x(i+1))/u(i);
end

相关文章推荐

SGU 183 Painting the balls (DP优化)

183. Painting the balls time limit per test: 0.25 sec.  memory limit per test: 4096 KB inpu...

Codeforces Round #183 (Div. 1) / 303A Lucky Permutation Triple (强大的数学&想法题)

http://codeforces.com/problemset/problem/303/A when n is odd, A[i] = B[i] = i when n is even...

LintCode 183 木材加工

本人电子系,只为一学生。心喜计算机,小编以怡情。第一次:用遍历的方式,然后超时,死在了[2147483644,2147483645,2147483646,2147483647], 4上public i...

codeforces 183 B Zoo--计算几何+枚举

/* 题意:有n个望远镜,被放置在X轴1~n的位置上,现有m只鸟,望远镜可以以任何角度观察,能看到在那条直线上的所有鸟, 问所有望远镜能看到的鸟的数量之和 计算几何+暴力 */ #include #...

LeetCode 数据库181,182,183,184,185整理

181. Employees Earning More Than Their ManagersThe Employee table holds all employees including thei...
  • WuLex
  • WuLex
  • 2017年07月02日 18:33
  • 135

NYOJ-183 赚钱啦【最短路】

题目链接:http://acm.nyist.net/JudgeOnline/problem.php?pid=183 解题思路: 最短路裸题。 最长路+判断正环 学会的东西是如何用SPF...

MATLAB·提取图像中多个目标

基于matlab工具箱提取图像中的多目标特征(代码如下): 代码前面部分为提取图像的边界信息,调用了后面的遍历函数Pixel_Search,函数实现方法见后~ %%ROI Testing c...

盈建科加密狗1.8.3 183 驱动

  • 2017年10月13日 17:16
  • 1.71MB
  • 下载

jquery183chm_jb51net.rar

  • 2014年12月07日 12:57
  • 4.26MB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:MyTDMA 数值分析P183
举报原因:
原因补充:

(最多只允许输入30个字)