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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 信息科学 » 程序算法 » 利用FT和DFT计算卷积
71/1返回列表
查看: 1311  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 北科281学硕材料求调剂 +6 tcxiaoxx 2026-03-20 6/300 2026-03-22 20:23 by edmund7
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 354求调剂 +7 Tyoumou 2026-03-18 10/500 2026-03-22 11:11 by 人来盛
[考研] 资源与环境 调剂申请(333分) +5 holy J 2026-03-21 5/250 2026-03-21 22:42 by Catalysis25
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 初始318分求调剂（有工作经验） +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 295材料求调剂，一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
信息提示
关闭
请填处理意见
关闭