 |
最新回复
#include <stdio.h>
#include <stdlib.h>
#include <csl.h>
#include <csl_timer.h>
#include <math.h>
/* ======================================================================== */
/* Include the DSPLIB header file */
/* ======================================================================== */
#include \"dsp_fft32x32.h\"
/* ======================================================================== */
/* Macro definition */
/* ======================================================================== */
#define PI (3.141592654)
#define NN (512)
/* ======================================================================== */
/* Input and output arrays */
/* ======================================================================== */
#pragma DATA_ALIGN(x, 8)
#pragma DATA_ALIGN(xx, 8)
#pragma DATA_ALIGN(y, 8)
#pragma DATA_ALIGN(w, 8)
int x[2*NN];
int w[2*NN];
int xx[2*NN];
int y[2*NN];
unsigned int z[2*NN];
short ch[NN];
/* ======================================================================== */
/* D2I -- Truncate a \'double\' to a \'int\', with clamping. */
/* ======================================================================== */
static int d2i(double d)
{
if (d >= 2147483647.0) return (int)0x7FFFFFFF;
if (d <= -2147483648.0) return (int)0x80000000;
return (int)d;
}
int gen_twiddle_fft32x32(int *w, int n, double scale)
{
int i, j, k, s=0, t;
for (j = 1, k = 0; j < n >> 2; j = j << 2, s++)
{
for (i = t=0; i < n >> 2; i += j, t++)
{
w[k + 5] = d2i(scale * cos(6.0 * PI * i / n));
w[k + 4] = d2i(scale * sin(6.0 * PI * i / n));
w[k + 3] = d2i(scale * cos(4.0 * PI * i / n));
w[k + 2] = d2i(scale * sin(4.0 * PI * i / n));
w[k + 1] = d2i(scale * cos(2.0 * PI * i / n));
w[k + 0] = d2i(scale * sin(2.0 * PI * i / n));
k += 6;
}
}
return k;
}
void main(void)
{
FILE *file;
short iRead;
int i;
int tempI,tempR;
/* ==================================================================== */
/* Initialize Chip Support Library */
/* ==================================================================== */
CSL_init();
file = fopen(\"test1.PCM\", \"rb\"
iRead = fread(ch, 1, NN, file);
//iRead = fread(ch, 1, NN, file);
fclose(file);
/* ==================================================================== */
/* Generate Q.31 input data */
/* ==================================================================== */
for(i=0; i<NN; i++)
{
x[2*i] = xx[2*i] = (int)ch;
x[2*i+1] = xx[2*i+1] = 0;
}
/* ==================================================================== */
/* Generate twiddle factors */
/* ==================================================================== */
gen_twiddle_fft32x32(w, NN, 2147483647.);
DSP_fft32x32(w, NN, x, y);
for(i=0; i<NN; i++)
{ tempI=y[i*2];
tempR=y[i*2+1];
z = tempI * tempI + tempR * tempR;
}
printf(\"END\\n\" );
}
/* ======================================================================== */
/* End of file: fft32x32_main.c */
/* ======================================================================== */
test1.pcm是8K采样的16BIT语音信号文件。用cool edit pro软件产生。实际是一个440hz正弦信号。
问题描述(用的是CCS3.1环境）：
1、当信号幅值很小的时候，运行程序到末尾后用view/graph看z数组。发现会出现两个对称的峰值。第一个落在z[28]附件，另外一个落在z[480]附件。为什么这样？这是一个正弦波
。虽然z没有开平方，但是应该是单一峰值啊。
2、当增大波形的幅值，频率不变后，再看z。发现峰值的波形变宽变乱。为什么这个算法和输入的动态范围有关？如何让FFT变换和输入幅值没有关系？
3、把iRead = fread(ch, 1, NN, file);语句重复，也就是ch的值取文件的第二次读写结果进行变换，z的结果完全乱了。按照道理，波形文件足够长，又是一个正弦波形。
应该没有变化才对。我怀疑什么地方错了。但是我用view/graph看xx的FFT，没有改变。这么说变换前数据是一样的。这是为什么？
详情
回复
| |
|
|
|
提升卡
变色卡
千斤顶
活动 更多>>
- 测评上线：国产FPGA 高云GW5AT-LV60 FPGA图像开发板
- 有奖直播｜ADMT4000 掀起多圈编码器设计革命
- 有奖直播 | 艾迈斯欧司朗 OSP 开放协议，从氛围灯动态照明到传感器交互融合
- 是德科技有奖直播 | EDICON 2025 《人工智能与机器学习》技术论坛
- 安富利、Nordic、TDK有奖直播报名 | AIoT“算法+芯片+终端”边缘计算解决方案
- TI有奖直播：MSPM0 系列 MCU 再添新成员:高性能与高性价比的优秀组合 MSPM0G351x/MSPM0L111x
- 【干货上新】电源解决方案和技术第二趴 | DigiKey 应用探索站
- Follow me第三季第1期来袭，与Digikey得捷一起解锁高性能嵌入式开发板STM32F429I-DISC1超能力
开源项目 更多>>
- FRDM-HB2001-EVM: 面向HB2001的Freedom平台，搭载10 A H桥、可编程有刷直流电机驱动器
- 35W 砖型适配器参考板
- STM32103ZET6拓展板LCD12864平台冰墩墩
- DC2874A-B，使用 LTC3310S 5V、30A 同步降压 Silent Switcher 2 的演示板
- LM431C 限流、电流源稳压器的典型应用
- 使用 Analog Devices 的 LTC3374AEUHF 的参考设计
- STM32F411 system1
- 使用 Seeed Technology Co.,Ltd 的 SSD1327Z 的参考设计
- 5 1811571128 徐晓剑
- DC1813A-H，用于 LTC2364CMS-18、18 位、250 ksps、15.5MHz 低功率、低噪声模数转换器的演示板
随便看看
查找数据手册?
EEWorld Datasheet 技术支持
热门标签
相关文章 更多>>
新帖速递
- STM32和无源蜂鸣器播放声音的问题
- 车规级AECQ200介绍，混合铝电解电容器的选择
- 嵌入式教程_DSP技术_DSP实验箱操作教程：2-28 搭建轻量级WEB服务器实验
- OPA847IDBVR运放器国产替代
- AG32VF407测试UART
- 【得捷电子Follow Me第二期】第一章 收到货物的分享
- 请问这个红外接收头是什么型号？能用哪个型号代替？谢谢
- 出售全新未拆封ZYNQ 7Z020 FPGA核心板
- 用在锂电池供电的水表设置上的LORA模块，当有100块水表集中安装在一个楼道内时，节能
- 请问一下，当某个端口被设置为 RX0后，这个端口的输入输出方向还有必要设置吗
- 今年怎么这么难，比疫情时还难，三十了面临失业好迷茫
- 请教稳压管测试问题
- 【小华HC32F448测评】关于小华半导体的UART中断发送和PRINTF构造和重定向
- 【BIGTREETECH PI开发板】 HDMI输出测试
- 【BIGTREETECH PI开发板】+08.音频测试（zmj）
关闭
站长推荐
1/10 
1/10
- 【TI有奖直播报名中】
- 直播主题：MSPM0 系列 MCU 再添新成员:高性能与高性价比的优秀组合 MSPM0G351x / MSPM0L111x
直播时间：4月24日（周四）上午10:00
直播奖励：好礼等您拿！
电子工程世界版权所有
京B2-20211791
京ICP备10001474号-1
电信业务审批[2006]字第258号函
京公网安备 11010802033920号
Copyright © 2005-2025 EEWORLD.com.cn, Inc. All rights reserved
京公网安备 11010802033920号
Copyright © 2005-2025 EEWORLD.com.cn, Inc. All rights reserved
