matlab两序列卷积运算

需求规格说明书(老系统)

实验一 离散时间信号分析

一、实验目的

1．掌握各种常用的序列，理解其数学表达式和波形表示。

2．掌握在计算机中生成及绘制数字信号波形的方法。

3．掌握序列的相加、相乘、移位、反褶等基本运算及计算机实现与作用。

4．掌握线性卷积软件实现的方法。

5．掌握计算机的使用方法和常用系统软件及应用软件的使用。

6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。

二、实验原理

1．序列的基本概念

离散时间信号在数学上可用时间序列{x(n)}来表示，其中x(n)代表序列的第n个数字，n代表时间的序列，n的取值范围为 n 的整数，n取其它值x(n)没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号xa(t)进行等间隔采样，采样间

隔为T，得到{xa(nT)}一个有序的数字序列就是离散时间信号，简称序列。

2．常用序列

常用序列有：单位脉冲序列（单位抽样） (n)、单位阶跃序列u(n)、矩形序列RN(n)、实指数序列、复指数序列、正弦型序列等。

3．序列的基本运算

序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。

4．序列的卷积运算

y(n) x(m)h(n m)

m x(n) h(n)

上式的运算关系称为卷积运算，

需求规格说明书(老系统)

实验一 离散时间信号分析

一、实验目的

1．掌握各种常用的序列，理解其数学表达式和波形表示。

2．掌握在计算机中生成及绘制数字信号波形的方法。

3．掌握序列的相加、相乘、移位、反褶等基本运算及计算机实现与作用。

4．掌握线性卷积软件实现的方法。

5．掌握计算机的使用方法和常用系统软件及应用软件的使用。

6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。

二、实验原理

1．序列的基本概念

离散时间信号在数学上可用时间序列{x(n)}来表示，其中x(n)代表序列的第n个数字，n代表时间的序列，n的取值范围为 n 的整数，n取其它值x(n)没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号xa(t)进行等间隔采样，采样间

隔为T，得到{xa(nT)}一个有序的数字序列就是离散时间信号，简称序列。

2．常用序列

常用序列有：单位脉冲序列（单位抽样） (n)、单位阶跃序列u(n)、矩形序列RN(n)、实指数序列、复指数序列、正弦型序列等。

3．序列的基本运算

序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。

4．序列的卷积运算

y(n) x(m)h(n m)

m x(n) h(n)

上式的运算关系称为卷积运算，

DSP图像处理设计

实验1 序列的卷积运算 一、实验目的 二、实验原理 三、实验内容

四、实验思考题

1

北京交通大学 国家工科电工电子教学基地

DSP应用技术

DSP图像处理设计

一、实验目的1.掌握利用CCS建立工程、编译与调试代码的基本过程

2.掌握序列线性卷积和相关的基本原理3.掌握利用DSP实现序列线性卷积运算

2

北京交通大学 国家工科电工电子教学基地

DSP应用技术

DSP图像处理设计

二、实验原理序列x[k]和序列h[k]线性卷积的定义为y[k ] x[k ] h[k ]

n

x[n]h[k n]

3

北京交通大学 国家工科电工电子教学基地

DSP应用技术

DSP图像处理设计

三、实验内容1.工程的建立

2.工程的调试

4

北京交通大学 国家工科电工电子教学基地

DSP应用技术

DSP图像处理设计

1. 工程的建立CCS以工程(*.pjt)的方式管理DSP程序 在计算序列线性卷积和相关的工程中，需要如下文件：(1) 建立工程名test.pjt;

(2) 向工程中添加计算线性卷积和相关的C程序test.c(3) 向工程中添加连接命令文件test.cmd, (4) 向工程中添加库文件rts55x.lib (一般C语言程序需要C语言运行支持库，在C:\C

DSP图像处理设计

实验1 序列的卷积运算 一、实验目的 二、实验原理 三、实验内容

四、实验思考题

1

北京交通大学 国家工科电工电子教学基地

DSP应用技术

DSP图像处理设计

一、实验目的1.掌握利用CCS建立工程、编译与调试代码的基本过程

2.掌握序列线性卷积和相关的基本原理3.掌握利用DSP实现序列线性卷积运算

2

北京交通大学 国家工科电工电子教学基地

DSP应用技术

DSP图像处理设计

二、实验原理序列x[k]和序列h[k]线性卷积的定义为y[k ] x[k ] h[k ]

n

x[n]h[k n]

3

北京交通大学 国家工科电工电子教学基地

DSP应用技术

DSP图像处理设计

三、实验内容1.工程的建立

2.工程的调试

4

北京交通大学 国家工科电工电子教学基地

DSP应用技术

DSP图像处理设计

1. 工程的建立CCS以工程(*.pjt)的方式管理DSP程序 在计算序列线性卷积和相关的工程中，需要如下文件：(1) 建立工程名test.pjt;

(2) 向工程中添加计算线性卷积和相关的C程序test.c(3) 向工程中添加连接命令文件test.cmd, (4) 向工程中添加库文件rts55x.lib (一般C语言程序需要C语言运行支持库，在C:\C

实验一：离散时间序列卷积和MATLAB实现

（一）实验目的

1、掌握离散卷积计算方法；

2、学会线性常系数差分方程的迭代解法； 3、学会针对具体系统设计程序；

（二）实验内容

1、用MATLAB画出如下表达式的脉冲序列

f(n)?8.0?(n)?3.4?(n?1)?1.8?(n?2)?5.6?(n?3)?2.9?(n?4)?0.7?(n?5)

例如：x(1)=8.0; n=0:5;

fn=[8.0 3.4 1.8 5.6 2.9 0.7]; stem(n,fn); xlabel('n'); ylabel('fn');

876543210fn00.511.522.5n33.544.55

2、 计算卷积

用MATLAB计算序列{1 2 3 3 2 1 1}和序列{-1 2 0 1}的离散卷积。 首先用手工计算，然后用MATLAB函数conv(x,y)编程验证。 MATLAB部分参考程序如下： a=[1 2 3 3 2 1 1]; c=conv(a,b); M=length(c)-1; stem(n,c);

xlabel('n'); ylabel('幅度');

subplot(3,1,1); a=[1 2 3 3 2 1 1

信号与线性系统分析

实验四 使用matlab实现卷积的运算

一 实验目的

1、

2、

二 实验内容 学习MATLAB语言的编程方法及熟悉MATLAB指令； 深刻理解卷积运算，利用离散卷积实现连续卷积运算；

1、 完成f1(t)与f2(t)两函数的卷积运算

其中：f1(t) e 2tu(t),f2(t) u(t) u(t 4)在一个图形窗口中，画出f1(t)、f2(t)以及卷积结果。要求每个坐标系有标题、坐标轴名称。

p = 0.0001; %定义时间间隔 t= 0:p:10;

%定义时间向量

f1=exp(-2*t).*u(t); %将f（t）表示出来 f2=u(t)-u(t-4);

f=conv(f1,f2);

subplot(1,2,1);

plot(t,f1,t,f2); title('f1=e^-2t*u(t)'' / ''f2=u(t)-u(t-4)');

xlabel('t(sec)'); % 这行代码是给出x坐标的标签 ylabel('f(t)');

grid on；

subplot(1,2

. . . .. .. 连续时间信号卷积运算的MATLAB 实现

一、实验目的

（1） 理解掌握卷积的概念及物理意义。

（2） 理解单位冲击响应的概念及物理意义。

二、实验原理

根据前述知识，连续信号卷积运算定义为

1212()()()()()f t f t f t f f t d τττ∞-∞=*=-?

卷积计算可以通过信号分段求和来实现，即 1212120()()()()()lim ()()k f t f t f t f f t d f k f t k τττ∞∞

-∞?→=-∞=*=-=??-???∑?

如果只求当t n =?（n 为整数）时()f t 的值()f n ?，则由上式可得

1212()()()()[()]k k f n f k f n k f k f n k ∞∞=-∞=-∞?=??

???-?=????-?∑∑

上式中的

12()[()]k f k

1、卷积码编码

function [output]=cnv_encd(input)

%output=cnv_encd(g,k0,input) 卷积码编码函数 %g 生成矩阵 %k0 输入码长

%input 输入信源序列 %output 输出卷积编码序列 g=[1 1 1;1 0 1];编码矩阵 k0=1;

input=[1 1 0 1];

if rem(length(input),k0)>0

input=[input,zeros(size(1:k0-rem(length(input),k0)))]; end

n=length(input)/k0;

if rem(size(g,2),k0)>0

error('Error,g is not of the right size.') end

li=size(g,2)/k0; n0=size(g,1);

u=[zeros(size(1:(li-1)*k0)),input,zeros(size(1:(li-1)*k0))];

u1=u(li*k0:-1:1); for i=1:n+li-2

u1=[u1,u((i+li)*k0:-1:i*k0+1)]; end

uu=res

1、卷积码编码

function [output]=cnv_encd(input)

%output=cnv_encd(g,k0,input) 卷积码编码函数 %g 生成矩阵 %k0 输入码长

%input 输入信源序列 %output 输出卷积编码序列 g=[1 1 1;1 0 1];编码矩阵 k0=1;

input=[1 1 0 1];

if rem(length(input),k0)>0

input=[input,zeros(size(1:k0-rem(length(input),k0)))]; end

n=length(input)/k0;

if rem(size(g,2),k0)>0

error('Error,g is not of the right size.') end

li=size(g,2)/k0; n0=size(g,1);

u=[zeros(size(1:(li-1)*k0)),input,zeros(size(1:(li-1)*k0))];

u1=u(li*k0:-1:1); for i=1:n+li-2

u1=[u1,u((i+li)*k0:-1:i*k0+1)]; end

uu=res

龙源期刊网 http://www.qikan.com.cn

基于Ｍａｔｌａｂ的离散卷积

作者：刘国良

来源：《现代电子技术》2009年第05期

摘 要：卷积运算广泛用于通讯、电子、自动化等领域的线性系统的仿真、分析及数字信号处理等方面。在Matlab中可以使用线性卷积、圆周卷积和快速傅里叶运算实现离散卷积。线性卷积是工程应用的基础，但圆周卷积和快速傅里叶运算实现线性离散卷积具有速度快等优势，圆周卷积采用循环移位，在Matlab中没有专用函数，需要根据圆周卷积的运算过程编制程序代码；快速傅里叶运算(FFT)是DSP的核心算法，在序列比较长时FFT是一种最合适的方法，运算速度快、程序简单，序列越长其优势越明显。以同一个例子介绍了进行离散卷积仿真运算的两种方法与特点。

关键词：Matlab；数字信号处理；离散卷积运算；快速傅里叶运算 中图分类号：TP311文献标识码：B 文章编号：1004-373X(2009)05-125-02 Discreted Convolution Based on Matlab LIU Guoliang