连续信号卷积
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
连续信号卷积
————————————————————————————————作者: ————————————————————————————————日期:
ﻩ
实验二连续信号卷积
一、实验目的
卷积积分可理解为某线性时不变系统在给定激励下的零状态响应。理解和掌握卷积运算对于线性系统分析来说至为关键。本实验的主要目的就是学习在
MATLAB环境中如何计算和分析连续时间信号的卷积
二、实验内容
在Matlab中,连续信号f1(t)与f2(t)的卷积可按下述过程求解:
1)构造两离散序列,f1(k)和f2(k),对应的时间向量k1, k2;
2)对两连续信号进行等间隔取样,得到离散序列f1(k)和f2(k);
3)调用Matlab提供的conv()函数计算两序列的卷积和f(k);
4)构造离散序列f(k)对应的时间向量k
三、学生实验内容:
1)已知两连续信号如下图所示,求它们的卷积近似,并记录波形
p=0.1;
k1=0:p:2;
f1=1/2*(k1);
k2=k1;
f2=1/2*(k2);
[f,k] = s co nv(f1,f2,k1,k2,p);[f,k] = sconv(f1,f2,k1,k2,p );
0.5
1 1.52
0.5
1
f1(t)
t
f 1(t )
00.5
1 1.52
0.5
1
f2(t)
t
f 2(t )
00.51 1.5
2
2.53
3.54
0.20.40.60.8f(t)=f1(t)*f2(t)
t
f (t )
2) 计算输入分别为ε(t)和t 时, 冲激响应为ε(t)的线性时不变系统的输出,验证该
系统为积分器 (将积分器的理论输出与sc onv()函数提供的卷积近似绘制在一张图上)
3)选择任意连续信号,验证f(t-τ1)*h(t-τ2)= f(t) *h(t-τ1-τ2)=
f(t-τ1-τ2) * h(t)
4)
5)计算并比较有限长序列的线性卷积和循环卷积。计算循环卷积可调用Matlab提
供的cconv(A,B,N)函数,其中A,B为参与卷积的序列,N为拓延周期。N的缺省值为length(A)+ length(B) -1,即缺省情况下无混叠发生。改变N的取值,以观察拓延周期(或者说频域抽样间隔)对循环卷积的影响。
6)
7)试用fft()和ifft()函数(快速离散傅里叶变换和反变换)计算两序列的循环
卷积
Circonvt.m
function y=circonvt(x1,x2,N)
x1=[x1,zeros(1,N-length(x1))];
x2=[x2,zeros(1,N-length(x2))];
X1=fft(x1,N);
X2=fft(x2,N);
X=X1.*X2;
y=ifft(X,N);
y=real(y);
A=[1 -7 2-38 -1 -3 0];
B=[002-3 5-1 2 -4];
C=circonvt(xn,yn,8);
stem(C);
四、实验要求
学习和掌握在Matlab环境中如何定义一个新函数。尝试换用不同信号测试编写的连续信号卷积程序sconv(),依据运算结果进一步理解函数卷积的物理意义;认真完成实验报告
五、实验小结
⑴通过实验我们学习和掌握了在Matlab环境中如何定义一个新函数。
⑵在实验中,我们应用Matlab实现对一些较难画出的图形信号函数进行画图。
⑶通过这次仿真实验,是我对进一步的理解了函数卷积的物理意义,使得我们对函数积学起来更容易些。