Question

我正在尝试在Linux上进行实时绘图。一切运行正常，只是当程序结束时，一段时间内会有大量的CPU使用。此外，因为我在usleep上使用很短的时间（如0.3秒）我不知道为什么bash shell由gnuplot接口占主导地位，即使对于其他bash windows我尝试打开，所以我不能使用ctrl + c来停止该计划。

我可以使用其他功能代替usleep吗？或者我可以纠正的东西，以避免窗户上的gnuplot优势？

提前致谢并抱歉我的英文不好

我让代码在这里

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
float s=10.;
float r=28.;
float b=8.0/3.0;
/* Definimos las funciones */
float f(float x,float y,float z){
    return s*(y-x);
}
float g(float x,float y,float z){
    return x*(r-z)-y;
}
float h(float x,float y,float z){
    return x*y-b*z;
}
FILE *output;
FILE *gp;

int main(){
    gp = popen("gnuplot -","w");
    output = fopen("lorenzgplot.dat","w");
    float t=0.; 
    float dt=0.01;
    float tf=30;
    float x=3.;
    float y=2.;
    float z=0.;
    float k1x,k1y,k1z, k2x,k2y,k2z,k3x,k3y,k3z,k4x,k4y,k4z;
    fprintf(output,"%f %f %f \n",x,y,z);
    fprintf(gp, "splot '/home/david/documents/lorenzgplot.dat' with lines \n");
/* Ahora Runge Kutta de orden 4 */  
    while(t<tf){
        /* RK4 paso 1 */
        k1x = f(x,y,z)*dt;
        k1y = g(x,y,z)*dt;
        k1z = h(x,y,z)*dt;
        /* RK4 paso 2 */
        k2x = f(x+0.5*k1x,y+0.5*k1y,z+0.5*k1z)*dt;
        k2y = g(x+0.5*k1x,y+0.5*k1y,z+0.5*k1z)*dt;
        k2z = h(x+0.5*k1x,y+0.5*k1y,z+0.5*k1z)*dt;
        /* RK4 paso 3 */
        k3x = f(x+0.5*k2x,y+0.5*k2y,z+0.5*k2z)*dt;
        k3y = g(x+0.5*k2x,y+0.5*k2y,z+0.5*k2z)*dt;
        k3z = h(x+0.5*k2x,y+0.5*k2y,z+0.5*k2z)*dt;
        /* RK4 paso 4 */
        k4x = f(x+k3x,y+k3y,z+k3z)*dt;
        k4y = g(x+k3x,y+k3y,z+k3z)*dt;
        k4z = h(x+k3x,y+k3y,z+k3z)*dt;
        /* Actualizamos las variables y el tiempo */
        x += (k1x/6.0 + k2x/3.0 + k3x/3.0 + k4x/6.0);
        y += (k1y/6.0 + k2y/3.0 + k3y/3.0 + k4y/6.0);
        z += (k1z/6.0 + k2z/3.0 + k3z/3.0 + k4z/6.0);
        /* finalmente escribimos sobre el archivo */

        fprintf(output,"%f %f %f \n",x,y,z);
        fflush(output); 
        usleep(100000);
        fprintf(gp, "replot \n");
        fflush(gp);
        t += dt;
    }
    fclose(gp);
    fclose(output);
    return 0;
}

Answer 1

考虑：

gp = popen("nice gnuplot -","w");

会阻止gnuplot在某种程度上吃掉你的系统，但你有没有想过在启动gnuplot之前运行你的runge-kutta代码？将输出写入文件，然后将文件提供给后台作业中的gnuplot。

例如，如果你愿意的话，这一切都可以在一个你可以使用和/或使用的脚本中使用。

runge-kutta > rkfile.dat
cat rkfile.dat | gnuplot - &
wait

例如：

at -k now <<!
 cd /path/to/rk
 nice ./my_rkscript.sh
!

所有这些都是保持工作不吃现有互动环境的一些方法。这样一来，“窗口”就不会被束缚，可以这么说。示例脚本需要重定向stderr以及您根据需要添加的其他一些改进。

gnuplot和usleep on C

1 个答案: