文章详情

短信预约-IT技能 免费直播动态提醒

请输入下面的图形验证码

提交验证

短信预约提醒成功

为什么我的两个正弦波的离散傅里叶变换在零处有一个巨大的尖峰?

2024-04-05 00:26

关注

偷偷努力,悄无声息地变强,然后惊艳所有人!哈哈,小伙伴们又来学习啦~今天我将给大家介绍《为什么我的两个正弦波的离散傅里叶变换在零处有一个巨大的尖峰?》,这篇文章主要会讲到等等知识点,不知道大家对其都有多少了解,下面我们就一起来看一吧!当然,非常希望大家能多多评论,给出合理的建议,我们一起学习,一起进步!

问题内容

我能够对两个正弦波(例如频率 1 和 10)相加进行傅里叶变换。我使用根据持续时间和样本数量计算出的采样频率将其转换为频率幅度图。

正如预期的那样,我在 1 和 10 处得到了两个不同的光点,但在 0 处也出现了我无法解释的巨大尖峰。

生成的代码将正弦波添加在一起以创建信号:

const input_num_samples = 1 << 10
const input_duration = 10
input_freqs := []float64{1, 10}

fftinput := make([]complex128, input_num_samples)
for _, freq := range input_freqs {
    for i := range fftinput {
        inputpoint := float64(i) / float64(input_num_samples) * input_duration
        fftinput[i] += complex(inputpoint, math.cos(inputpoint*freq*math.pi*2))
    }
}

执行傅立叶变换的代码:

fftoutput := fourier.newcmplxfft(input_num_samples).coefficients(make([]complex128, len(fftinput)), fftinput)

生成频率幅度图的代码

R := float64(INPUT_NUM_SAMPLES) / float64(INPUT_DURATION)
fmt.Printf("Sampling frequency: %f\n", R)
freqMag := make([]complex128, len(fftOutput))
for index, value := range fftOutput {
    freqMag[index] = complex((float64(index)/INPUT_NUM_SAMPLES)*R, math.Sqrt(real(value)*real(value)+imag(value)*imag(value)))
}

频率-幅度图(省略对称的另一半): 信号图:

我尝试了一些方法,包括:

  1. 使用不同的 fft 库
  2. 反转输入的实部和虚部,奇怪的是,这似乎为零差异。 我预计在 1 和 10 处有两个光点,就是这样(不计算对称性,输出的另一侧已从图中省略)


正确答案


线的右侧

fftinput[i] += complex(inputpoint, math.cos(inputpoint*freq*math.pi*2))

用实部 inputpoint 和虚部 math.cos(inputpoint*freq*math.pi*2) 构造一个复数值,例如

x(t) = 2 t + j (cos(t π 2) + cos(t π 4)),  0 ≤ t < 10.

因此,输入信号是一个斜坡(由于实部)加上两个正弦曲线(来自虚部)的和。斜坡是非负的,这解释了频率 0 (dc) 处频谱中的大尖峰。此外,斜坡在间隔端点处是不连续的,这解释了其他重要的低频分量。

根据您的描述,我相信您的意图是输入信号是两个实值正弦曲线。将行更改为:

fftInput[i] += complex(math.Cos(inputPoint*freq*math.Pi*2), 0.0)

本篇关于《为什么我的两个正弦波的离散傅里叶变换在零处有一个巨大的尖峰?》的介绍就到此结束啦,但是学无止境,想要了解学习更多关于Golang的相关知识,请关注编程网公众号!

阅读原文内容投诉

免责声明:

① 本站未注明“稿件来源”的信息均来自网络整理。其文字、图片和音视频稿件的所属权归原作者所有。本站收集整理出于非商业性的教育和科研之目的,并不意味着本站赞同其观点或证实其内容的真实性。仅作为临时的测试数据,供内部测试之用。本站并未授权任何人以任何方式主动获取本站任何信息。

② 本站未注明“稿件来源”的临时测试数据将在测试完成后最终做删除处理。有问题或投稿请发送至: 邮箱/279061341@qq.com QQ/279061341

软考中级精品资料免费领

  • 历年真题答案解析
  • 备考技巧名师总结
  • 高频考点精准押题
  • 2024年上半年信息系统项目管理师第二批次真题及答案解析(完整版)

    难度     813人已做
    查看
  • 【考后总结】2024年5月26日信息系统项目管理师第2批次考情分析

    难度     354人已做
    查看
  • 【考后总结】2024年5月25日信息系统项目管理师第1批次考情分析

    难度     318人已做
    查看
  • 2024年上半年软考高项第一、二批次真题考点汇总(完整版)

    难度     435人已做
    查看
  • 2024年上半年系统架构设计师考试综合知识真题

    难度     224人已做
    查看

相关文章

发现更多好内容

猜你喜欢

AI推送时光机
位置:首页-资讯-后端开发
咦!没有更多了?去看看其它编程学习网 内容吧
首页课程
资料下载
问答资讯