C语言实现matlab的butter函数（免complex.h）

#include 
#include 
#include 

#include 
double pi = 3.1415926535897932384626433832795L;
#if 0
#include “complex.h“
#else
typedef struct _C_double_complex
{
	double _Val[2];
} _C_double_complex;

_C_double_complex cexp（_C_double_complex x）
{
	_C_double_complex r;
	double p;
	p = exp（x._Val[0]）;
	r._Val[0] = p*cos（x._Val[1]）;
	r._Val[1] = p*sin（x._Val[1]）;
	return r;
}

_C_double_complex conj（_C_double_complex x）
{
	_C_double_complex r;
	r._Val[0] = x._Val[0];
	r._Val[1] = -x._Val[1];
	return r;
}

_C_double_complex _Cmulcr（_C_double_complex x double y）
{
	_C_double_complex r;
	r._Val[0] = x._Val[0] * y;
	r._Val[1] = x._Val[1] * y;
	return r;
}

_C_double_complex _Cmulcc（_C_double_complex x _C_double_complex y）