北邮22级信通院DSP:用C++程序实现DFT连续求取任意两个离散实序列从n到m点循环卷积

本文主要是介绍北邮22级信通院DSP:用C++程序实现DFT连续求取任意两个离散实序列从n到m点循环卷积,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

北邮22信通一枚~

跟随课程进度更新北邮信通院DSP的笔记、代码和文章,欢迎关注~

获取更多文章,请访问专栏:

北邮22级信通院DSP_青山入墨雨如画的博客-CSDN博客

目录

1.定义

2.求解方法(通过实例讲解)

3.程序思维分析

3.1选择存储结构

3.2双向队列的优点

3.2.1.自然倒置

3.2.2.快速头插

3.3实现过程分析

3.3.1.补零操作

3.3.2DSP中的“倒置”和程序中的“倒置”矛盾与解决方法

3.3.3移位与计算

3.3.4汇总

4.代码部分

4.1代码部分

4.2运行结果


1.定义

 

2.求解方法(通过实例讲解)

求解DFT的循环卷积有三种方法,分别是同心圆法、波形图法和解析式法。

例题:

同心圆解法

 

 

​ 

3.程序思维分析

注:这里的“倒置”是“程序倒置”,而不是DSP中规定的序列倒置,详情请参考下文

3.3.2DSP中的“倒置”和程序中的“倒置”矛盾与解决方法

3.1选择存储结构

其中,序列的存储和计算可以使用队列结构,综合考虑需要进行倒置操作,最终选定双向队列实现。

3.2双向队列的优点

3.2.1.自然倒置

如果使用普通队列,就涉及队列倒置函数的书写,其中时间复杂度至少为n;而使用双向队列可以从不同的方向输入,如果从队列头部开始输入,则在输入过程中就已经自然地实现了队列的倒置操作。

3.2.2.快速头插

如果使用普通队列,在队列头部添加元素是件很繁琐的工程,而双向队列有相应的内置函数,简便地完成头插操作。

3.3实现过程分析

3.3.1.补零操作

在DSP中,如果两个序列的长度小于要计算的卷积点数,需要进行补零操作。

补零操作之后,第二个序列才会进行倒置。

所以,对第一个序列(不用进行倒置的序列),直接采用双向队列的尾插法补零即可。

对第二个序列(需要进行倒置的序列),倒置之后直接采用双向队列的头插法补零即可。

3.3.2DSP中的“倒置”和程序中的“倒置”矛盾与解决方法

按照DSP的观点,序列x(n)=  {1,2,3,4}的倒置应为x(-n)=  {1,4,3,2};

而程序中倒置{1,2,3,4}的结果为{4,3,2,1};

所以在用程序实现之前,需要先进行一步移位操作。

多举几个例子:

DSP中序列x(n)=  {1,2,3,4,0}的倒置应为x(-n)=  {1,0,4,3,2};

而程序中倒置{1,2,3,4,0}的结果为{0,4,3,2,1};

所以在用程序实现之前,需要先进行一步移位操作。

3.3.3移位与计算

        补零、倒置之后两个序列长度相等。对应项相乘在相加,得到结果序列中的第一个元素。第二个序列再次将头元素移动到队列末尾构成新序列,重复上面的操作,得到结果序列的第二个元素。

        重复上面过程,最终输出结果。

3.3.4汇总

第一个序列尾插法输入;

第二个序列头插法输入(实现了程序倒置,但还没有成为DSP倒置);//{0,4,3,2,1}

补零操作;

第二个序列移位一次(这里才真正意义上完成了第一次DSP倒置);//{1,0,4,3,2}

两个序列对应项相乘,结果相加,尾插入结果序列;

第二个序列移位一次;

两个序列对应项相乘,结果相加,尾插入结果序列;//{2,1,0,4,3}

……

输出结果序列。

3.3.5多点循环卷积的更多例子

4.代码部分

4.1代码部分

#include<iostream>
#include<deque>
#include<cmath>
using namespace std;
int deque_min;
int deque_max = 100;
void deque_print(deque<int>a)
{cout << "打印队列如下:" << endl;for (int i = 0; i < a.size(); i++){cout << a.front() << " ";a.push_back(a.front());a.pop_front();}cout << endl;
}
void destruction(deque<int>&a)
{while(!a.empty())a.pop_back();
}
int main()
{system("color 0A");deque<int>a, b, rtn;int input, cnt_a = 0, cnt_b = 0, temp_sum = 0;cout << "请问您想从几点卷积开始计算?" << endl;cin >> deque_min;cout << "请问您想终止于几点卷积?" << endl;cin >> deque_max;cout << "请输入周期卷积序列1:" << endl;while (cin >> input){a.push_back(input);cnt_a++;if (getchar() == '\n')break;}cout << "请输入周期卷积序列2:" << endl;while (cin >> input){b.push_front(input);cnt_b++;if ((getchar() == '\n'))break;}//judgeif (cnt_a > deque_min || cnt_b > deque_min || deque_min > deque_max){cout << "error" << endl;return 0;}//calculatecout << endl;for(int i=deque_min;i<=deque_max;i++){cout << i << "点卷积结果如下:" << endl;//add zero to each onewhile (cnt_a < i){a.push_back(0);cnt_a++;}while (cnt_b < i){b.push_front(0);cnt_b++;}//cyclic convolutionint cnt = i;while (cnt--){b.push_front(b.back());b.pop_back();for (int j = 0; j < i; j++){temp_sum += (a.front() * b.front());a.push_back(a.front()); b.push_back(b.front());a.pop_front(); b.pop_front();}rtn.push_back(temp_sum);temp_sum = 0;}//printdeque_print(a);deque_print(b);deque_print(rtn);//destruct the result deque for next period of storagedestruction(rtn);cout << endl;}return 0;
}

4.2运行结果

这篇关于北邮22级信通院DSP:用C++程序实现DFT连续求取任意两个离散实序列从n到m点循环卷积的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/912473

相关文章

springboot循环依赖问题案例代码及解决办法

《springboot循环依赖问题案例代码及解决办法》在SpringBoot中,如果两个或多个Bean之间存在循环依赖(即BeanA依赖BeanB,而BeanB又依赖BeanA),会导致Spring的... 目录1. 什么是循环依赖?2. 循环依赖的场景案例3. 解决循环依赖的常见方法方法 1:使用 @La

C++ 中的 if-constexpr语法和作用

《C++中的if-constexpr语法和作用》if-constexpr语法是C++17引入的新语法特性,也被称为常量if表达式或静态if(staticif),:本文主要介绍C++中的if-c... 目录1 if-constexpr 语法1.1 基本语法1.2 扩展说明1.2.1 条件表达式1.2.2 fa

C语言实现两个变量值交换的三种方式

《C语言实现两个变量值交换的三种方式》两个变量值的交换是编程中最常见的问题之一,以下将介绍三种变量的交换方式,其中第一种方式是最常用也是最实用的,后两种方式一般只在特殊限制下使用,需要的朋友可以参考下... 目录1.使用临时变量(推荐)2.相加和相减的方式(值较大时可能丢失数据)3.按位异或运算1.使用临时

C++中::SHCreateDirectoryEx函数使用方法

《C++中::SHCreateDirectoryEx函数使用方法》::SHCreateDirectoryEx用于创建多级目录,类似于mkdir-p命令,本文主要介绍了C++中::SHCreateDir... 目录1. 函数原型与依赖项2. 基本使用示例示例 1:创建单层目录示例 2:创建多级目录3. 关键注

C++从序列容器中删除元素的四种方法

《C++从序列容器中删除元素的四种方法》删除元素的方法在序列容器和关联容器之间是非常不同的,在序列容器中,vector和string是最常用的,但这里也会介绍deque和list以供全面了解,尽管在一... 目录一、简介二、移除给定位置的元素三、移除与某个值相等的元素3.1、序列容器vector、deque

C++常见容器获取头元素的方法大全

《C++常见容器获取头元素的方法大全》在C++编程中,容器是存储和管理数据集合的重要工具,不同的容器提供了不同的接口来访问和操作其中的元素,获取容器的头元素(即第一个元素)是常见的操作之一,本文将详细... 目录一、std::vector二、std::list三、std::deque四、std::forwa

Python循环缓冲区的应用详解

《Python循环缓冲区的应用详解》循环缓冲区是一个线性缓冲区,逻辑上被视为一个循环的结构,本文主要为大家介绍了Python中循环缓冲区的相关应用,有兴趣的小伙伴可以了解一下... 目录什么是循环缓冲区循环缓冲区的结构python中的循环缓冲区实现运行循环缓冲区循环缓冲区的优势应用案例Python中的实现库

C++字符串提取和分割的多种方法

《C++字符串提取和分割的多种方法》在C++编程中,字符串处理是一个常见的任务,尤其是在需要从字符串中提取特定数据时,本文将详细探讨如何使用C++标准库中的工具来提取和分割字符串,并分析不同方法的适用... 目录1. 字符串提取的基本方法1.1 使用 std::istringstream 和 >> 操作符示

C++原地删除有序数组重复项的N种方法

《C++原地删除有序数组重复项的N种方法》给定一个排序数组,你需要在原地删除重复出现的元素,使得每个元素只出现一次,返回移除后数组的新长度,不要使用额外的数组空间,你必须在原地修改输入数组并在使用O(... 目录一、问题二、问题分析三、算法实现四、问题变体:最多保留两次五、分析和代码实现5.1、问题分析5.

C++ 各种map特点对比分析

《C++各种map特点对比分析》文章比较了C++中不同类型的map(如std::map,std::unordered_map,std::multimap,std::unordered_multima... 目录特点比较C++ 示例代码 ​​​​​​代码解释特点比较1. std::map底层实现:基于红黑