应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 信息科学 » 程序算法 » 利用FT和DFT计算卷积
71/1返回列表
查看: 1740  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 yjjyjjhh 的 10 个金币

yjjyjjhh

银虫 (小有名气)

[交流] 利用FT和DFT计算卷积

我在将cos变换转成卷积形式的积分，因为：（1）采用对数坐标采样适合我后续的工作；（2）转换后余弦项的傅里叶变换有解析表达式。
但在一些区间上，计算结果和matlab命令conv的计算结果不太一致。求大佬帮我指点指点。

由于涉及公式，具体思路和结果放到图片里了。

感谢！金币有限，谢意无限！



图片无法上传。。。。。只能用网盘了

链接: https://pan.baidu.com/s/1zE0XMWYLrsBsP8upQHrA-Q 提取码: arav

链接: https://pan.baidu.com/s/1meCLiNqfNvWf6td0TKtXVQ 提取码: ygem
回复此楼

» 猜你喜欢
1楼 2021-12-17 11:15:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qqinly

专家顾问 (著名写手)

至人无己，神人无功，圣人无名
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
conv(a,b)=ifft(fft(a1).*fft (b1))，a和b要先补0到a和b的两者长度之和-1的长度，即补0的结果为a1和b1，再计算fft。如果想用2的整数次幂计算fft，还要多补0，计算结果丢弃掉多余的点

发自小木虫Android客户端
赞一下 回复此楼
学而时习之，不亦说乎？有朋自远方来，不亦乐乎？人不知而不愠，不亦君子乎？
2楼2021-12-18 09:31:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjjyjjhh

银虫 (小有名气)
引用回帖:
2楼: Originally posted by qqinly at 2021-12-18 09:31:12
conv(a,b)=ifft(fft(a1).*fft (b1))，a和b要先补0到a和b的两者长度之和-1的长度，即补0的结果为a1和b1，再计算fft。如果想用2的整数次幂计算fft，还要多补0，计算结果丢弃掉多余的点

...

你好。麻烦问一下，计算结果丢掉多余的点是指什么？比如a和b的原来数据点都为N0，补零到2*N0。需要丢到哪些点？另外，若a的傅里叶变换有解析解，有什么需要注意的地方？谢谢
赞一下 回复此楼
3楼2021-12-18 17:22:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qqinly

专家顾问 (著名写手)

至人无己，神人无功，圣人无名
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by yjjyjjhh at 2021-12-18 17:22:49
你好。麻烦问一下，计算结果丢掉多余的点是指什么？比如a和b的原来数据点都为N0，补零到2*N0。需要丢到哪些点？另外，若a的傅里叶变换有解析解，有什么需要注意的地方？谢谢...

我是说，如果你用2的整数次幂的点数计算fft和ifft，而你滑动计算出来的卷积结果小于这个点数的话，最后ifft计算出来的结果只保留前面一段和卷积长度一样的点。另外，卷积的结果的长度是两组数长度之和减1。
这个就是公式推导出来的卷积和fft的等价关系。
不清楚的就用matlab跑几个数验证一下。

发自小木虫Android客户端
赞一下 回复此楼
学而时习之，不亦说乎？有朋自远方来，不亦乐乎？人不知而不愠，不亦君子乎？
4楼2021-12-18 19:42:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjjyjjhh

银虫 (小有名气)
引用回帖:
4楼: Originally posted by qqinly at 2021-12-18 19:42:00
我是说，如果你用2的整数次幂的点数计算fft和ifft，而你滑动计算出来的卷积结果小于这个点数的话，最后ifft计算出来的结果只保留前面一段和卷积长度一样的点。另外，卷积的结果的长度是两组数长度之和减1。
这个就 ...

非常感谢耐心解答！我理解这个了。
我运算结果要移位后才能与conv计算结果在中间一定区间吻合，而在区间两端差异很大（链接2中的图），我以为您说的和这个有关。。这里面有两个问题我不知道为什么：（1）为什么要移位？（2）为什么只在一定区间吻合？
非常感谢！
若是方便，可以加我的QQ：610005697指导一下吗？若不嫌弃，可以有偿。非常感谢！！！
赞一下 回复此楼
5楼2021-12-19 14:36:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qqinly

专家顾问 (著名写手)

至人无己，神人无功，圣人无名
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
yjjyjjhh: 金币+10 2021-12-20 16:20:38
引用回帖:
5楼: Originally posted by yjjyjjhh at 2021-12-19 14:36:09
非常感谢耐心解答！我理解这个了。
我运算结果要移位后才能与conv计算结果在中间一定区间吻合，而在区间两端差异很大（链接2中的图），我以为您说的和这个有关。。这里面有两个问题我不知道为什么：（1）为什么要 ...

如果数据只有几千个点，直接一次卷积conv或用FFT/IFFT就搞定了，如果数据点数较多，是分段FFT/IFFT的话，分段取的数据，每段数据的前后要多取数据长度-1那么长的数据点再进行FFT/IFFT。
另外，看数据是否吻合，你还是再画两幅相减以后的图吧。
私聊给你联系方式
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

yjjyjjhh
学而时习之，不亦说乎？有朋自远方来，不亦乐乎？人不知而不愠，不亦君子乎？
6楼2021-12-19 22:19:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjjyjjhh

银虫 (小有名气)
送红花一朵
引用回帖:
6楼: Originally posted by qqinly at 2021-12-19 22:19:30
如果数据只有几千个点，直接一次卷积conv或用FFT/IFFT就搞定了，如果数据点数较多，是分段FFT/IFFT的话，分段取的数据，每段数据的前后要多取数据长度-1那么长的数据点再进行FFT/IFFT。
另外，看数据是否吻合，你 ...

您好，已私信。非常感谢！
赞一下 回复此楼
7楼2021-12-20 16:21:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yjjyjjhh 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-10-01 21:07 by Equinoxhua
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +7 IDs3scOF0tjC 2026-09-28 8/400 2026-10-01 19:08 by Equinoxhua
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:45 by Z9YWQ5EAO3qp
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
关闭
请填处理意见
关闭