C语言结合gnuplot绘制sinc函数波形

C程序

#include<stdio.h>
#include<math.h>
#define pi 3.14159265
main()
{
	int i;
	double sinc, t;
	for(i = -5000; i<= 5000; i++)	//因为t=i/500.0，所以t的取值范围为[-10,10]，但严格的说t的取值范围为[-10,0)U(0,10]
	{
		if(i == 0)
			continue;
		t = i / 500.0;	//一个单位内取500个样本点
		sinc = sin(pi*t) / (pi*t);
		printf("%e\t%e\n", t, sinc);
	}
}

gnuplot绘图

如果上述程序中不加continue语句，则产生的波形如下图所示：

#include<stdio.h>
#include<math.h>
#define pi 3.14159265
main()
{
	int i;
	double sinc, t;
	for(i = -5000; i<= 5000; i++)	//因为t=i/500.0，所以t的取值范围为[-10,10]
	{
		t = i / 500.0;	//一个单位内取500个样本点
		sinc = sin(pi*t) / (pi*t);
		printf("%e\t%e\n", t, sinc);
	}
}

matlab代码

x = -10 : 0.01 : 10;
y = sin(pi*x) ./ (pi*x);
plot(x, y);
axis([-12 12 -0.5 1.5]);

matlab绘图