fortran：计算卷积

最新推荐文章于 2024-09-21 16:42:28 发布

chder_白南

最新推荐文章于 2024-09-21 16:42:28 发布

阅读量1.4k

点赞数 2

分类专栏： FortranNote Fortran教程数字信号处理文章标签： fortran

本文链接：https://blog.csdn.net/chd_lkl/article/details/90579499

版权

FortranNote 同时被 3 个专栏收录

81 篇文章 46 订阅

订阅专栏

Fortran教程

68 篇文章 57 订阅

订阅专栏

数字信号处理

7 篇文章 2 订阅

订阅专栏

计算卷积的物理意义与公式网上与书上很多，这里不再赘述。
下面是计算两个向量u和v的fortran代码
program test_conv
        !// 计算向量u和v的卷积
        implicit none
        integer            :: i, k
        integer, parameter :: m = 3, n = 8  !// 向量u和向量v的长度
        real               :: u(0:m-1) = [1, 1, 1], v(0:n-1) = [1, 1, 0, 0, 0, 1, 1, 1]
        real               :: conv(0:m+n-2)
        real               :: s1, s2 

        conv = 0.
        do i = 0, m+n-2

                do k = 0, m+n-2
                        if (k >= 0 .and. k <= m-1) then
                                s1 = u(k)
                        else
                                s1 = 0.
                        end if

                        if (i-k >=0 .and. i-k <= n-1) then
                                s2 = v(i-k)
                        else
                                s2 = 0.
                        end if

                        conv(i) = conv(i) + s1*s2
                end do

        end do

        print*, conv

end program test_conv