产生采样率为8000次每秒的正弦信号,时长为1秒(8000点)。正弦信号的幅度为0.6,频率为500Hz,或50Hz,将500Hz和50Hz的两个正弦信号相加。分别用gnuplot画出相应的波形图。
编译器:TCC
第一步:c语言产生数据
(1)频率为500Hz的c程序
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define pi 3.1415926
int main()
{
int i;
double t,s;
for (i=0;i<8000;i++)
{
t=i/8000.0 ; //采集时间=采样点数/采样率
s=0.6*sin(2*pi*500*t);//0.6是振幅,500是频率
printf("%.16lf %.16lf\n",t,s);
}
system ("pause");
return 0;
}
(2)频率为50Hz的c程序
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define pi 3.1415926
int main()
{
int i;
double t,s;
for (i=0;i<8000;i++)
{
t=i/8000.0 ; //采集时间=采样点数/采样率
s=0.6*sin(2*pi*50*t);
printf("%.16lf %.16lf\n",t,s);
}
system ("pause");
return 0;
}
(3)频率为50Hz和500Hz相加的c程序
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define pi 3.1415926
int main()
{
int i;
double t,s;
for (i=0;i<8000;i++)
{
t=i/8000.0 ; //采集时间=采样点数/采样率
s=0.6*sin(2*pi*500*t)+0.6*sin(2*pi*50*t);
printf("%.16lf %.16lf\n",t,s);
}
system ("pause");
return 0;
}
第二步:用gnuplot画出数据波形
方法一:重定向
TCC编译后生成的.exe文件变成.dat 文件
在gnuplot命令行提示符下输入作图指令(注意自己tcc路径,必要是需要将.exe和.dat复制到用户下面)
[0:0.02]定义了横坐标的范围,[-2:2]是纵坐标的显示范围,u 1:2
表示使用数据文件的1,2列数据,w l表示使用折线作图。
方法二:直接将"500hezi.dat"写成"<500hezi.exe"(注意<符号)
gnuplot> plot [0:0.02] [-2:2] "<500hezi.exe" u 1:2 w l
第三步:图像显示
(1)500Hz图像
(2)50Hz图像
(3)500Hz图像+50Hz图像