C语言中的代数操作

时间:2015-03-16 10:02:13

标签: c algebraic-data-types computer-algebra-systems

现在我正在编写一个可以执行基本线性代数的程序,其中我有2个常数:A和B,我可以简单地输入它们的系数并在其上添加其他值。

示例:

我想找到fibbonaci序列的术语

A
B

我想拿A + B并将其附加到档案。

A
B
A+B

现在我想添加第二和第三个词

A
B
A + B
A + 2B

等等。

我有一个程序可以对所有数值做得很好,但是,我希望看到代数上升,不使用binets公式。

我唯一的猜测是将A和B的系数存储到2个单独的文件中以便计算它们然后打印出来,但是,如果C中的线性代数可用则更容易。

编辑:没关系,我忘了这是一个适合计算机书呆子的地方,而不是数学书呆子

1 个答案:

答案 0 :(得分:1)

作为一种传统的命令式语言,C没有对代数表达式的内置支持。要计算Fibonacci序列的系数,可以编写如下程序:

#include <stdio.h>

int main() {
    int i;
    int coeff_a_1 = 1;
    int coeff_b_1 = 0;
    int coeff_a_0 = 0;
    int coeff_b_0 = 1;

    printf("A\n");
    printf("B\n");

    for (i = 0; i < 15; i++) {
        int coeff_a = coeff_a_0 + coeff_a_1;
        int coeff_b = coeff_b_0 + coeff_b_1;

        printf("%dA + %dB\n", coeff_a, coeff_b);

        coeff_a_1 = coeff_a_0;
        coeff_b_1 = coeff_b_0;
        coeff_a_0 = coeff_a;
        coeff_b_0 = coeff_b;
    }

    return 0;
}

这将打印:

A
B
1A + 1B
1A + 2B
2A + 3B
3A + 5B
5A + 8B
8A + 13B
13A + 21B
21A + 34B
34A + 55B
55A + 89B
89A + 144B
144A + 233B
233A + 377B
377A + 610B
610A + 987B

对于更复杂的问题,我建议使用提供C绑定的computer algebra system