W11-12 << Previous Next >> w15
w13

#include <stdio.h>
 
int main() {
// 開啟文件以寫入位移和速度數據
FILE *outputFile = fopen("motion_data.txt", "w");
if (!outputFile) {
fprintf(stderr, "無法創建數據文件。\n");
return 1;
}
 
// 模擬運動 10 秒並計算位移和速度，同時將數據寫入文件
double x = 0.2; // 初始位移
double v = 0.0; // 初始速度
double dt = 0.01; // 時間步長
double t = 0.0; // 時間
 
while (t <= 10.0) {
double acceleration = (-10.0 * x - 0.5 * v) / 1.0; // 這裡修改了系統的參數
v += acceleration * dt;
x += v * dt;
 
fprintf(outputFile, "%lf %lf %lf\n", t, x, v);
 
t += dt;
}
 
// 關閉數據文件
fclose(outputFile);
 
// 使用popen啟動Gnuplot進程
FILE *gnuplotPipe = popen("gnuplot -persistent", "w");
if (!gnuplotPipe) {
fprintf(stderr, "無法啟動Gnuplot。\n");
return 1;
}
 
// 使用Gnuplot繪圖指令，指定字型文件和輸出PNG
fprintf(gnuplotPipe, "set terminal png font 'default,12' size 800,400\n");
fprintf(gnuplotPipe, "set output './../images/motion_plot.png'\n");
fprintf(gnuplotPipe, "set title 'displacement and velocity vs. time'\n");
fprintf(gnuplotPipe, "set xlabel 'time (s)'\n");
fprintf(gnuplotPipe, "set ylabel 'displacement (m)'\n");
fprintf(gnuplotPipe, "plot 'motion_data.txt' using 1:2 with lines lw 2 title 'displacement', \
'motion_data.txt' using 1:3 with lines lw 2 title 'velocity'\n");
 
// 關閉Gnuplot進程
fprintf(gnuplotPipe, "exit\n");
pclose(gnuplotPipe);
 
return 0;
}
W11-12 << Previous Next >> w15