基本 C 程序中的循环

Question

我无法找到并解决此代码中的错误。
我希望 Parray 和 Tarray 在循环中填充方程的迭代。我尝试了 for 和 while 循环，但似乎只有数组的第一个成员被赋予了正确的值，并且该值在新循环开始时丢失。

#include <stdio.h>
#include <stdlib.h>

int main() {

    int i;
    int j;
    float dt,P_init,t_init,t_end;
    dt = 0.03125;
    P_init = 30;
    t_init = 0;
    t_end = 30;

    int n_steps = 0;
    n_steps = t_end/(float)dt;

    int Parray[n_steps+1];
    int Tarray[n_steps+1];

    for (i=0; i<n_steps+1; i++)
       {
           Parray[i]=0;
           Tarray[i]=0;
       }



    Parray[0] = P_init;
    Tarray[0] = t_init;

    int x;
    float P,t,dpdt;

    x = 1;
    while (x < n_steps+1)
       {

        float s,d;
        P = Parray[x-1];
        t = Tarray[x-1];
        dpdt = 0.7 * P * (1-(P/750)) - 20;
        s = P + (dt * dpdt);
        Parray[x] = s;
        d = t + dt;
        Tarray[x] = d;
        x++;
        printf("%f  %f  ",s,d);

               }



    return(0);
}

Answer 1

#include <stdio.h>
#include <stdlib.h>

int main(void) 
{
    int i, j, x;
    float dt, P_init, t_init, t_end, P, t, dpdt;

    dt = 0.03125;
    P_init = 30;
    t_init = 0;
    t_end = 30;

    int n_steps = 0;
    n_steps = t_end / (float)dt;

    // using float arrays
    float Parray[n_steps+1];
    float Tarray[n_steps+1];

    for (i = 0; i < n_steps+1; i++)
           Parray[i] = 0;

    for (j = 0; j < n_steps+1; j++)
           Tarray[j] = 0;

    Parray[0] = P_init;
    Tarray[0] = t_init;

    for (x = 1; x < n_steps + 1; x++) {
        float s , d;
        P = Parray[x-1]; 
        t = Tarray[x-1];
        dpdt = 0.7 * P * (1 - (P / 750)) - 20;
        s = P + (dt * dpdt);
        Parray[x] = s;
        d = t + dt;
        Tarray[x] = d;
        printf("%f  %f\n", s, d);
    }

    return(0);
}

我将 int 数组更改为浮点数组。 你在这里实际计算的是什么？

Answer 2

这一行

dpdt = 0.7 * P * (1-(P/750)) - 20;

它似乎有一个浮点值，因此数组必须是浮点类型，并且您放入数组中的所有输入都必须是浮点类型。 改变这个：

int Parray[n_steps+1];
int Tarray[n_steps+1];

到

float Parray[n_steps+1];
float Tarray[n_steps+1];

和

P_init = 30;
t_init = 0;
t_end = 30;

它们都应该是 30.0, 0.0, 30.0 不要忘记这一点：

 Parray[i]=0;
 Tarray[i]=0;

将其更改为 0.0

简而言之，C 不擅长区分类型，所以如果你声明一个浮点数组，所有输入都必须是 x.0，而不仅仅是 x，其中 x 是整数。