[问答] 关于傅里叶变换的简单C语言程序

[复制链接]

技术员

发表于 2016-11-8 11:17:03   1060 查看 10 回复 显示全部楼层 倒序浏览
分享
想做一个简单的傅里叶变换,只能显示32个给的原始数据,但是无法显示处理后的数据,不知道程序上哪里出错了,请大神看看,怀疑是子程序调用上有问题,请明示~~~坐等~~~
QQ截图20161108111543.png
#include "stdio.h"#include "math.h"
void dft(x,y,a,b,n,sign);
void main()
{
int i,j,n;
double a1,a2,c,c1,c2,d1,d2,q1,q2,w,w1,w2;
double x[32],y[32],a[32],b[32];
n=32;
a1=0.9;
a2=0.3;
x[0]=1.0;
y[0]=0.0;
for(i=1;i<n;i++)
        {
        x=a1*x[i-1]-a2*y[i-1];
        y=a2*x[i-1]+a1*y[i-1];
        }
printf("\n原始数据\n");
for(i=0;i<n/2;i++)
{
        for(j=0;j<2;j++)
                printf(" %10.7f +J %10.7f",x[2*i+j],y[2*i+j]);
        printf("\n");
}
q1=x[n-1];
q2=y[n-1];
printf("\n傅里叶变换后数据\n");
dft(x,y,a,b,n,1);
for(i=0;i<2/n;i++)
{
        for(j=0;j<2;j++)
                printf(" %10.7f +J %10.7f",a[2*i+j],b[2*i+j]);
                printf("\n");
}

system("pause");
}

/*傅里叶函数*/
void dft(x,y,a,b,n,sign)
int n,sign;
double x[],y[],a[],b[];
{int i,k;
double c,d,q,w,s;
q=6.28318530718/n;
for(k=0;k<n;k++)
        {w=k*q;
     a[k]=b[k]=0.0;
         for(i=0;i<n;i++)
                {d=i*w;
                 c=cos(d);
                 s=sin(d)*sign;
                 a[k]+=c*x+s*y;
                 b[k]+=c*y-s*x;
                }
        }
if(sign==-1)
        {
        c=1.0/n;
        for(k=0;k<n;k++)
        {a[k]=c*a[k];
        b[k]=c*b[k];
        }
        }

}


标签:C语言 程序
发表于 2016-11-8 13:27:48  

回帖奖励 +1

PCB在线计价下单

板子大小:

cm
X
cm

层数:

2

板子数量:

10

厚度:

1.6
你把n附一个初值看看,,,要不然没办法知道你n的值
回复

点赞

总工程师

发表于 2016-11-8 20:40:04  

回帖奖励 +1

我想问一下 当初在学校学过这个傅里叶基数  但是就是不知道这个到底有什么用呀
回复

点赞

技术员

发表于 2016-11-11 14:21:21    楼主|
王栋春 发表于 2016-11-8 20:40
我想问一下 当初在学校学过这个傅里叶基数  但是就是不知道这个到底有什么用呀

用处挺多,我是信号处理,把信号从时域变到频域,便于分析它的频域特性
回复

点赞

总工程师

发表于 2016-11-11 20:19:57  
dantehan 发表于 2016-11-11 14:21
用处挺多,我是信号处理,把信号从时域变到频域,便于分析它的频域特性

了解了  哎可惜没有学会
回复

点赞

实习生

发表于 2016-11-21 08:36:27  

回帖奖励 +1

加油!!!!!!!!!!!!!!!
回复

点赞

技术员

发表于 2016-11-27 00:05:24  
加油。
回复

点赞

noe

技术员

发表于 2017-3-12 13:00:40  

回帖奖励 +1

加油
回复

点赞

技术员

发表于 2017-3-15 11:58:47  

回帖奖励 +1

顶一个                     
回复

点赞

实习生

发表于 2017-3-15 16:07:24  
回复

点赞

等待验证会员

发表于 2017-3-26 22:13:08  
c++写的控制台啊
回复

点赞

高级模式
您需要登录后才可以回帖 登录 | 注册

专家问答 查看更多>>
关闭

站长推荐 上一条 /7 下一条

小黑屋|手机版|Archiver| 电子发烧友 ( 粤ICP备14022951号-2 )     

GMT+8, 2017-5-30 03:10 , Processed in 0.108976 second(s), 17 queries , Memcache On.

微信扫描
快速回复 返回顶部 返回列表
-

推荐专区

技术干货集中营

专家问答

用户帮助┃咨询与建议┃版主议事

我的提问

工程师杂谈

工程师创意

工程师职场

论坛电子赛事

社区活动专版

发烧友活动

-

嵌入式论坛

单片机/MCU论坛

FPGA|CPLD|ASIC论坛

DSP论坛

嵌入式系统论坛

-

电源技术论坛

电源技术论坛

-

硬件设计论坛

电路设计论坛

电子元器件论坛

控制|传感

总线技术|接口技术

-

测试测量论坛

LabVIEW论坛

Matlab论坛

测试测量技术专区

仪器仪表技术专区

-

EDA设计论坛

multisim论坛

PCB设计论坛

proteus论坛|仿真论坛

EasyEDA-中国人自已的EDA工具

-

综合技术与应用

电机控制

智能电网

光电及显示

参考设计中心

汽车电子技术论坛

医疗电子论坛

-

开源硬件

树莓派论坛

智能硬件论坛

开发快智能硬件开发平台

Intel物联网开发者专区

Waveshare

乐美客SBC专区

Arduino论坛

BeagleBone论坛

机器人论坛

创客神器NanoPi

小钢炮CANNON

比派科技banana pi专区

-

无线通信论坛

无线通信技术专区

天线|RF射频|微波|雷达技术

-

IC设计论坛

芯片测试与失效分析

Mixed Signal/SOC[数模混合芯片设计]

Analog/RF IC设计

设计与制造封装测试

-

个人版区

阿东Verilog技术专版

直流马达驱动电路设计

LabVIEW英雄联盟

特权同学FPGA专区

-

厂商专区

灵动微电子 MM32

盈鹏飞嵌入式

TI论坛

TI Deyisupport社区

芯灵思嵌入式论坛

Tisan

米尔科技

庆科社区

WIZnet技术专区

Cypress技术论坛

飞凌嵌入式

Qualcomm技术论坛

英创嵌入式

机智云GoKit论坛

-

检测技术与质量

电磁兼容(EMC)设计与整改

安规知识论坛

检测与认证

-

消费电子论坛

手机技术论坛

平板电脑/mid论坛

音视/视频/机顶盒论坛

-

电子论坛综合区

聚丰众筹官方社区

新人报道区

聚丰供应链

-

论坛服务区

-

供求信息发布

供需广告

招聘┃求职发布区

电子展览展会专区