追赶法

原创 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
*/

相关文章推荐

用追赶法解决数值积分的问题

  • 2010年11月27日 00:07
  • 324B
  • 下载

追赶法求解三对角矩阵

  • 2013年06月18日 21:12
  • 431B
  • 下载

追赶法解线性方程组

C语言代码实现; # include # define N 5 /* *实验二:用追赶法解方程组 */ main(){ float A[N][N]={0}; float f...

追赶法解三对角方程组

  • 2010年07月23日 20:40
  • 2KB
  • 下载

追赶法解方程

算法 详情见追赶法公式 这个博客只是用来储存。。 package 计算方法大作业; public class zhuiganfa { public static void main(Stri...

楚列斯基和追赶法的Matlab程序

  • 2014年11月13日 20:36
  • 29KB
  • 下载

追赶法-ZOJ-1601-Integer Approximation

#include#include #include #include #include #include #include #include #include #include #include #i...

计算方法-追赶法-C++编程实现

  • 2009年12月05日 17:55
  • 8KB
  • 下载

解三对角矩阵的追赶法

  • 2014年01月09日 11:41
  • 635B
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:追赶法
举报原因:
原因补充:

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