三对角矩阵逆C

时间:2013-10-19 21:28:15

标签: c matrix inverse

您好我正在尝试为c中的三对角矩阵创建反函数,但是我的代码给了我一个段错误。代码是

/* Inverse of a n by n matrix */
#include<stdio.h>
#include<math.h>
#include<stdlib.h>

double J[3][3];

double phi(i, n)
{
double p;
if(i==n)
{
    return 1; 
}
else if(i==n-1)
{
    return J[i][i];
}
else
{
    p=J[i][i]*phi(i+1,n)-J[i][i+1]*J[i+1][i]*phi(i+2,n);
}
return p;
}

double theta(i,n)
{
double p;
if(i==-1)
{
    return 1;
}
else if(i==0)
{
    return J[i][i];
}
else
{
    p=J[i][i]*theta(i-1,n)-J[i-1][i]*J[i][i-1]*theta(i-2,n);
}
return p;
}

void main(void)
{
double T[3][3];
int n=3;
int i, j, a, b;
J[0][0]=1;
J[0][1]=1;
J[0][2]=0;
J[1][0]=2;
J[1][1]=2;
J[1][2]=2;
J[2][0]=0;
J[2][1]=2;
J[2][2]=2;
for(i=0; i<n;i++);
{
    for(j=0; j<n;j++)
    {
    b=j;
    a=i;
    if(i<j)
    {
        T[i][j]=pow(-1, i+j);
        for(a; a<b; a++)
        {
            T[i][j]*=J[a][a+1];
        }
        T[i][j]*=theta(i-1,n)*phi(j+1)/theta(n-1,n);
    }
    else if(i==j)
    {
        if(i==0)
        {
            T[i][j]=0;
        }
        else
        {
            T[i][j]=theta(i-1,n)*phi(j+1,n)/theta(n-1,n);
        }
    }
    else
    {
        T[i][j]=pow(-1, i+j);
        for(b; b<i; b++)
        {
            T[i][j]*=J[a][b];
        }
        T[i][j]*=theta(j-1, n)*phi(i+1,n)/theta(n-1, n);
    }
    }
}

for(i=0;i<n;i++)
{
    for(j=0;j<n;j++)
    {
        printf("%lf\n", T[i][j]);
    }
}
}   

那里的J矩阵只是一个示例矩阵,我将把这个代码实现到另一个调用此函数矩阵的代码中。

1 个答案:

答案 0 :(得分:0)

首先,您应该用double phi(i, n)替换您的函数签名double phi(int i, int n)。 theta函数也是如此。此外,在您的代码中,您已经写了这一行: T[i][j]*=theta(i-1,n)*phi(j+1)/theta(n-1,n); 你的phi函数应该有2个参数,所以也许你想要这个:

T[i][j]*=theta(i-1,n)*phi(j+1, n)/theta(n-1,n);

我还建议您使用T[i][j]=pow((double)(-1), (double)(i+j));。 您可以看到signatures of pow here