C语言结合gnuplot绘制连续时间系统的频率响应曲线

C程序

#include<stdio.h>
#include<complex.h>
#define pi 3.1416
int main()
{ 
	double complex H; 
	double w;
	for(w = -6*pi; w <= 6*pi; w = w+0.01)
	{
		H = (I*w) / (0.2*cpow(I*w,3) + 0.2*cpow(I*w,2) + I*w);
		printf("%f\t%f\t%f\n", w, cabs(H), carg(H));
	}
}

gnuplot绘图

• 幅频特性曲线
  
• 相频特性曲线
  

matlab代码

w = -6*pi:0.01:6*pi;
b = [1 0];
a = [0.2 0.2 1 0];
H = freqs(b,a,w);
figure(1);
plot(w, abs(H));
figure(2);
plot(w, angle(H));

matlab绘图

• 幅频特性曲线
  
• 相频特性曲线