追赶法

原创 2012年03月22日 13:39:31
//追赶法   
#include <iostream.h>   
#include <stdio.h>   
#define N 4 //定义一个4*4的矩阵,改变N可以改变矩阵的大小   
void TDMA(float a[N],float b[N],float c[N-1],float f[N],float x[N])
{
	float u[N],l[N],y[N];

    u[0]=b[0];   
    y[0]=f[0];	
    for(int i=1;i<=N-1;i++)   
    {   
        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]; 		
    }    
	//计算x:   
    x[N-1]=y[N-1]/u[N-1];      
    for(i=N-2;i>=0;i--)   
    {   
        x[i]=(y[i]-c[i]*x[i+1])/u[i];   
    }      
}


//主函数   
void main()   
{        
float a[]={0,-1,-2,-3};
float b[]={2,3,4,5};
float c[]={-1,-2,-3};
float f[]={6,1,-2,1};  
float x[N];
TDMA(a,b,c,f,x);
  //输出x:   
    cout<<"方程组解x:"<<endl;   
    for(int i=0;i<=N-1;i++)      
    {   
        cout<<x[i]<<"\t";   
    }      
    cout<<endl;     
   
}  
/*
方程组解x:
5       4       3       2
Press any key to continue
*/

追赶法

给ab范围让a/b接近一个数,用追赶法。 直接枚举,取a,b为1每次对a,b求商,如果a / b > x,则a增加1,否则b增加1,每次记录下差值最小时a,b的值。...
  • ACpac
  • ACpac
  • 2014年12月30日 10:15
  • 352

计算方法之用追赶法求线性方程组

/************************************* * 用追赶法求线性方程组 * * |- -| |- -| |- -| * | 2 -1...
  • wzhg0508
  • wzhg0508
  • 2013年06月11日 12:57
  • 1187

追赶法求解三对角方程组

在这篇文章里,我们介绍追赶法的基本原理,以及用追赶法求解三对角方程组的算法....
  • u010450214
  • u010450214
  • 2015年12月04日 09:37
  • 7487

追赶法求解三对角方程组

在这篇文章里,我们介绍追赶法的基本原理,以及用追赶法求解三对角方程组的算法....
  • u010450214
  • u010450214
  • 2015年12月04日 09:37
  • 7487

数值分析 追赶法求解三对角线性方程组 MATLAB实现

函数主体部分编程算法  参考 数值分析 第四版 颜庆津 P27 运行结果截图: %追赶法求解三对角线性方程组,Ax=b,A用一维数组a,c,d存储。 function [L,U,x]=cr...
  • jingmiaa
  • jingmiaa
  • 2015年11月05日 11:48
  • 4727

Guass列主元、平方根法、追赶法求解方程组的C++实现

一,要解决的问题选用合适的算法,求解三种线性方程组:一般线性方程组,对称正定方程组,三对角线性方程组。 方程略。二,数值方法1,使用Guass列主元消去法求解一般线性方程组。Guass列主元是为了防...
  • qq_24421591
  • qq_24421591
  • 2016年03月28日 22:05
  • 1434

计算方法之追赶法

追赶法是解决稀疏矩阵的 一种方法,但应该满足对角占优,否则不能进行计算,下面是程序,语言为java import java.util.Scanner; /** * Created by q on...
  • baidu_38370610
  • baidu_38370610
  • 2017年11月30日 19:19
  • 80

追赶法解线性方程组

C语言代码实现; # include # define N 5 /* *实验二:用追赶法解方程组 */ main(){ float A[N][N]={0}; float f...
  • yyywww666
  • yyywww666
  • 2015年01月17日 14:11
  • 724

追赶法解方程

算法 详情见追赶法公式 这个博客只是用来储存。。 package 计算方法大作业; public class zhuiganfa { public static void main(Stri...
  • CrackMe_
  • CrackMe_
  • 2017年04月09日 17:33
  • 67

追赶法求三角方程组

#include #include #include using namespace std; class Mat { public: Mat(int x, int y, double d[...
  • qq_38875908
  • qq_38875908
  • 2017年11月06日 19:05
  • 77
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:追赶法
举报原因:
原因补充:

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