C程序设计谭浩强第四版终极程序
发布网友
发布时间：2022-08-23 07:30
共3个回答
热心网友
时间：2024-11-26 09:02
#include "stdio.h"
#include "REG52.H"
float max=0,sum=0; //计算失真度时所需要的全局变量
struct compx
{
float real;
float imag;
} ;
/*****************************************************************/
//函数名称: struct compx EE(struct compx b1,struct compx b2)
//返回值：旋转因子 即：x(n)W-N-kn(N为下标，kn为上标)
//功能： 生成旋转因子
//参数1：struct compx b1 即：x(n)
//参数2：struct compx b2 即: W-N-kn
/*****************************************************************/
struct compx EE(struct compx b1,struct compx b2) //旋转因子
{
struct compx b3;
b3.real=b1.real*b2.real-b1.imag*b2.imag;
b3.imag=b1.real*b2.imag+b1.imag*b2.real;
return(b3);
}
/*****************************************************************/
//函数名称: void FFT(struct compx *xin,int N)
//返回值：无
//功能： 实现FFT变换
//参数1：struct compx *xin 即：需要FFT变换的数组，类型为compx，对
// 于ADC采集后的数组其实部为AD值；虚部为0
//参数2：N 即: 转换的点数
/*****************************************************************/
void FFT(struct compx *xin,int N) //*xin为需要FFT变换的数组，N为转换的点数
{
int f,m,LH,nm,i,k,j,L;
double p , ps ;
long int le,B,ip;
float pi;
struct compx w,t;
LH=N>>1;
f=N;
for(m=1;(f=f>>1)!=1;m++){;}
nm=N-2;
j=N>>1;
for(i=1;i<=nm;i++) // 变址运算
{
if(i<j){t=xin[j];xin[j]=xin[i];xin[i]=t;}
k=LH;
while(j>=k){j=j-k;k=k/2;}
j=j+k;
}
for(L=1;L<=m;L++) //实现蝶形运算
{
le=(long int)pow(2,L);
B=le>>1;
pi=3.14159;
for(j=0;j<=B-1;j++)
{
p=pow(2,m-L)*j;
ps=2*pi/N*p;
w.real=cos(ps);//生成Wnk
w.imag=-sin(ps);
for(i=j;i<=N-1;i=i+le)
{
ip=i+B;
t=EE(xin[ip],w);
xin[ip].real=xin[i].real-t.real;
xin[ip].imag=xin[i].imag-t.imag;
xin[i].real=xin[i].real+t.real;
xin[i].imag=xin[i].imag+t.imag;
}
}
}
}
/*****************************************************************/
//函数名称:float THDcount(unsigned int *result)
//返回值：失真度值
//功能：求出失真度
//参数1：unsigned int *result 即：需求失真度数组地址
/*****************************************************************/
float THDcount(unsigned int *result)
{
unsigned int i=0,n=2,m,t=0;
float THD=0;
max=result[1];
for(i=2;i<256;i++)
{
if(result[i]>max)
{
max=result[i];
t=i;
}
}
m=t*n;
while(m<256)
{
sum=result[m]*result[m]+sum;
n++;
m=t*n;
}
sum=sum/(max*max);
THD=100*sqrt(sum);
return THD;
}
/*****************************************************************/
//函数名称:void main()
//返回值：无
//功能： 测试{2,3,3,2} 理论FFT转换后的值{10,-1-j,0,-1+j}
//参数： 无
//在Targe中设置Memory Mode为：Large:Variables in XDATA否则无法编译
/*****************************************************************/
void main()
{
struct compx aa[4];
aa[0].real=2;aa[0].imag=0;
aa[1].real=3;aa[1].imag=0;
aa[2].real=3;aa[2].imag=0;
aa[3].real=2;aa[3].imag=0;
FFT(aa,4);
while(1);
}
热心网友
时间：2024-11-26 09:02
热心网友
时间：2024-11-26 09:02
