1

因此，我编写了一些代码来从文件中获取矩阵，将其放入数组并计算其行列式。然而，在运行时，数组似乎没有被填充，程序返回一个 0 矩阵，因此行列式 = 0。我对编程很陌生，所以很难确定确切的原因，但我认为这是与从文件过程中读取的整个过程有关。

该文件只是一个 .data 文件，矩阵元素存储在空格分隔的列中，即：

1.7 1.2 0.5
0.0 -3.2 1.4
3.0 4.0 5.0

这是（我认为是）相关的代码部分，但如果有帮助，我可以发布整个内容。

int main(int argc, char* argv[])
{
    FILE          *input;
    int           record, i, j;
    int           curr_col;
    const int     dim = DIMENSION;
    double        entries[dim][dim];
    double        tmp, determinant;
    const char    inp_fn[]="matrix.data";

    /* Open file */
    input = fopen(inp_fn, "r");

    /* Check the pointer to file are not NULL */
    if(input != (FILE*) NULL)
    {
        for(j=0; j<dim; j++, fopen(inp_fn, "r"))
        {
            record = 0; i = 0;
            /* Read in records one by one and check the return value of fscanf */
            while(fscanf(input,"%lf",&tmp) == 1)
            {
                curr_col = (record % dim);
                if(curr_col==j)
                {
                    /* Copy data points to the array */
                    entries[i][j] = (double)(tmp);
                    i++;
                }
                record++;
            }
            fclose(input);
        }
    }
    else
        printf("*** Could not open input or output file! ***\n");

    /* Calculates determinant */
    determinant = det(dim, entries);

    printf("\nA = \n");
    for(i=0; i<dim; i++)
    {
        for(j=0; j<dim; j++)
        {
            printf("%.1lf ", entries[i][j]);
        }
        printf("\n");
    }
    printf("\n");

    printf("det(A) = %.3lf\n", determinant);

}

我得到“无法打开输入或输出文件！” 运行程序时出现错误和空矩阵...帮助！？

4

2 回答 2

3

肮脏的修复

从您的代码中，我看到您每次阅读不同专栏时都打算打开文件。它效率低下且笨重。您可以通过更改使其工作（仅读取输入部分，我不知道您的其余代码）：

for(j=0; j<dim; j++, fopen(inp_fn, "r"))

到

for(j=0; j<dim; j++, input = fopen(inp_fn, "r"))

您当前的代码将打开文件并浪费资源，同时fclose会遇到错误，因为文件在input上一次迭代中已关闭。

我上面建议的代码将分配新FILE*的 from fopento input。

当然，上面的方式效率非常低，正如我在开始时指出的那样。

更好的方法

一个更好的方法，在 if 语句中if(input != (FILE*) NULL)（用 删除循环j）：

record = 0;
// Read the file at most (dim * dim) times to fill up the array
// The reading also stops when it cannot read any double number
while(fscanf(input,"%lf",&tmp) == 1 && record < dim * dim)
{
    // Use some math to calculate the cell to put the new entry
    entries[record / dim][record % dim] = tmp; // entries and tmp are double, no need for casting

    record++;
}

// Close the file after done reading
fclose(input);

请注意，fopen在进入条件之前只调用一次if，并且一次读取所有内容。

您可能还想在阅读后添加检查以确保record == dim * dim- 以防提供的数据不够。

于 2012-11-22T02:17:43.957 回答
0

尝试perror("fopen");在打开失败时使用打印错误而不是 printf

另外，请注意，您每次都在 for 循环中重新打开文件。

改变

for(j=0; j<dim; j++, fopen(inp_fn, "r"))

到

for(j=0; j<dim; j++)
于 2012-11-22T01:50:21.070 回答