DSP中的卷积算法

卷积算法设计总框图卷积算法设计总框图 开 始 初始化 DSP 产生输入信号 将卷积信号X m 和 H m 输 入在同一个坐标系内 翻转 将 H m 以 m 0 的垂 直轴为轴翻褶成 H m 3 23 2 卷积算法设计的原理卷积算法设计的原理 移位 将 H m 移位 n 即 得 H n m 相乘 再将 H n m 和 X m 的相同 m 值的对应点值 相乘 相加 把以上所有点的对应点 的乘积叠加起来 即得 Y n 值 取移位值 n N 值取遍整 个坐标轴 结 束 YES 1 1 卷积算法基础理论 卷积算法基础理论 卷积的基本理论和公式 卷积和 对离散系统 卷积和 也是求线性时不变系统输出响应 零状态 响应 的主要方法 卷积和的运算在图形表示上可分为四步 A 翻褶 现在亚变量坐标 M 上作出 X m 和 H m 将 m 0 的垂直轴为轴翻 褶成 H m B 移位 将 H m 移位 n 即得 H n m 当 n 为正整数时 右移 n 位 当 n 为负整数时 左移 n 位 C 相乘 再将 H n m 和 X m 的相同 m 值的对应点值相乘 D 相加 把以上所有点的对应点的乘积叠加起来 即得 Y n 值 依上法 取 n 2 1 0 1 2 3 各值 即可得全部 Y n 值 程序流程图程序流程图 While TRUE 执行 dataIO2 子程序 执行 dataIO 子程序 执行 processing4 子程序 结束 YES NO 执行 processing4 子程序 开始 打印 volume example tarted n 执行 processing1 子程序 执行 processing2 子程序 执行 processing3 子程序 程序的自编函数及其功能 1 processing int input2 int output2 调用形式 processing1 int input2 int output2 参数解释 input2 output 为两个整型指针数组 返回值解释 返回一个 TURE 让主函数的 While 循环保持连续 功能说明 对输入的 input2 buffer 波形进行截取 m 点 再以零点的 Y 轴为对称轴进行翻褶 把生成的波形上的各点的值存入以 OUTPUT2 指针 开始的一段地址空间中 2 processing2 int output2 int output3 调用形式 processing2 int output2 int output3 参数解释 output2 output3 为两个整形指针数组 返回值解释 返回了一个 TREN 让主函数的 While 循环保持连续 功能说明 对输出的 output2 buffer 波形进行作 n 点移位 然后把生成 的波形上的各点的值存入以 OUTPUT3 指针开始的一段空间中 3 processing3 int input1 int output2 int output4 调用形式 processing3 int input1 int output2 int output4 参数解释 output2 output4 input1 为三个整型指针数组 返回值解释 返回了一个 TRUE 让主函数的 While 循环保持继续 功能说明 对输入的 input2 buffer 波形和输入的 input1 buffer 做卷 积和运算 然后把生成的波形上的各点的值存入以 OUTPUT4 指针开始的 地址空间中 4 processing4 int input2 int output1 调用形式 processing4 int input2 int output1 参数解释 output1 input 为两个整型指针数组 返回值解释 返回了一个 TRUE 让主函数的循环保持继续 功能说明 对输入的 input2 buffer 波形截取 m 点 然后把生成的波形上的各 点的值存入以 OUTPUT1 指针开始的一段地址空间中 源程序 源程序 include DSP281x Device h include DSP281x Examples h include f2812a h include stdio h include volume h int inp1 buffer BUFSIZE int inp2 buffer BUFSIZE int out1 buffer BUFSIZE int out2 buffer BUFSIZE int out3 buffer BUFSIZE int out4 buffer BUFSIZE 2 int size BUFSIZE int ain MINGAIN int zhy 0 int sk 64 unsigned int processingLoad 1 static int processing1 int output1 int output2 static int processing2 int output2 int output3 static int processing3 int input1 int output2 int output4 static int processing4 int input2 int output1 static void dataIO1 void static void dataIO2 void int input1 int input2 int output1 int output2 int output3 int output4 void main void int jishu 0 int input1 int input2 int output1 int output2 int output3 int output4 puts volume example started n while TRUE dataIO1 dataIO2 processing4 input2 output1 processing1 output1 output2 processing2 output2 output3 processing3 input1 output2 output4 jishu 在此处加断点 static int processing4 int input2 int output1 int m sk for m 0 m output1 input2 ain for size m 0 m output1 m 0 static int processing1 int output1 int output2 int m sk 1 for m 0 m output2 output1 ain static int processing2 int output2 int output3 int n zhy m size 64 for n0 i y x 0 z 0 f y for i 0 i g input1 z output2 f x x g z f output4 x y m sk y sk 1 w m zhy 1 for m 0 m y i y z sk 1 x 0 f sk y for i 0 i z f g input1 z output2 f x x g out4 buffer w x w static void dataIO1 return static void dataIO2 do data I O return 仿真设置仿真设置 1 在程序中 dataIO1 上单击鼠标右键选择 Toggle software breakpoint 设置软件断点 再在同一行上单击鼠标右键 选择 software breakpoint edit 来设置断点 此时打开了一个新 的窗口 最后设置成如下图所示 2 同理 在程序中 dataIO2 上同样操作设置 最后设置如下图 3 打开窗口菜单 Viwe Graph Frequency 进行如下设置 仿真图仿真图 1 当输入波为SIN SIN 时 SIN inp1 SIN inp2 SIN out 2 当输入波为SIN11 SIN11 时 SIN11 inp1 SIN11 inp2 SIN11 out 3 当输入波为SIN22 SIN22 时 SIN22 inp1 SIN22 inp2 SIN22 out 5 当输入波为SIN33 SIN33 时 SIN33 inp1 SIN33 inp2 SIN33 out 6 当输入波为SIN44 SIN44 时 SIN44 inp