angle.cpp

求取a+b*i的相位角的程序

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#define pi 3.1415926

void main()
{
	FILE *fp;
	int i;
	double w;
	double q[1001];
	double a;
	double b;
	for(i=0;i<=1000;i++)
	{
		w=-pi+pi*i/500;
		a=0.5+cos(w);
		b=-sin(w);
		if(a<0 && b>0)
		{
			q[i]=pi+atan(b/a);
		}
		else if(a<0 && b<0)
			{
				q[i]=atan(b/a)-pi;
			}
		else
		{
				q[i]=atan(b/a);
		}
	}
	if((fp=fopen("angle.txt","w"))!=NULL)
	{
			for(i=0;i<=1000;i++)
			
					fprintf(fp,"%f\n",q[i]);
					
			fclose(fp);
	}
}