unwrap.c

自己编的数值计算方法的c语言实现源程序。

/*求解卷绕后的相频响应
y为一维数组，长度为len；*/
void unwrap(y,len)
int len;
double y[];
{	int i,k,j;
	double pi2,temp;
	pi2=8.0*atan(1.0);
	for(i=len-1;i>0;i--)
		{	temp=y[i]-y[i-1];
			temp=abs(temp);
			if( temp>=pi2/2)
			{	for(k=len-1;k>i-1;k--)  
				{	if(y[i-1]>0)
					{	
						y[k]=y[k]+pi2;
					}
					else
					{	
						y[k]=y[k]-pi2;
					}
				}
			}
		}
}