用matlab求卷积不是稀罕的资料
工厂]墅的坛子皇里
axis([t(1)，t (l )，y_min—dy，y_max dy] )； grid:set(GCA，’y tick’，[y_min，yO，y_max]，’fontsize’，8 )调用这些函数，可以方便地计算并绘制任意有限长度波形和卷积波形。 用以下例子说明。
例1用Matlab描绘图中信号的卷积波形。
l，l‘f )
函数
y=cscoNVS(f，h，t_s，t_e，a，b )；
使用%计算卷积积分的解析解%f激励信号包括阶梯函数heaviside(t )；
%h不包含脉冲响应、卷积时混叠信号、阶跃函数，默认开始点为0；
%t_s，t_e是系统的零状态响应y=f*h波形的起点和终点，
%
a、b是卷积积分的基础上的下限。 例如：计算exp (一
t%h=sym(exp )一)2.t 2 ) ) :
工厂(2)。 r () ~
(a ) ) )。
-，l【"
2
-2
2t ) u(t 1 )与u ) t-3 )卷积。
%%%%
syms
f=sym (heav iside (t-2 ) ) ):CSCONVS(f ) f、h、1、5、2、t ) ) ) ) ) ) ) ) ) ) ) ) 652
(() ) )。
邻2
syms
t
tao
工厂]
一l
o
~
_
卅~
2
FTao=Subs(f )。 t tao； ht—Tao=Subs(h，t，t-tao )；
y=simple(int(Ftao*ht—tao，Tao，a，b ) )；
t=t—S:0.02:t—e；
LL
(b ) ) )。
图1
yt=subs(y ) y_max--max ) yt ) y_min=min ) yt ) dy=) y_ma
x—y_min(/I0 :
plot(t，yt，’line width’，2 )，axis ) [t_s11]；
xlabel(t ) sea ) )； 网格；
国际机场
y_min>；=O
t_e
Matlab程序如下。
fl=’2 * rect puls (t 1.5 )2* rect puls (t-1.5 ) ):F2=’2 * rect puls (t-2 )。 4 )’； 图形(1)。
sconv(fl，一)二，二，F2，0，4 ) fl=’2 * rect puls ) t，2 ) ):F2=’t.* rect puls ) t。 4 )’； 图形(2)。
csconv(fl，一1，l，f2，一2，2 ) )。
运行程序后显示的图表如图2所示。
y_min—dy
y_max dy]； %axis([一二] )
30
yo=(y_max—y—min )/2；
else
yO=O:end
set(GCA，’y tick’，[y_min，yO，y_max]，’fontsize’，8 ) titie’(卷积的波形)； disp’零状态响应’)，y活动实例说明其计算方法。
例2线性非时变系统的输入信号f(t )和脉冲响应h ) t )由以下的各式给出，系统的零状态响应)，)。
(a )、(f )=e“”I； o (一() (r一) )】、^ ()、) )=e一()“”
(b )解
、(f )=e-2、(O 3 )、^ ) r )=e-2t((1)1)。
(a )系统的零状态响应如下
(，)、)、) f ) ^ ) f )=j-ie -“f”【(o (一60-2 )】e-((t-.t ) e )卜t )出=j (根据E4“E1”。 L(O )碓(F(T ) dr(J ) E「E…1 ) e ) t ) t ) dr”
(a ) ) )。
图2
3
(b ) ) )。
=【e】【e】“出”e(f )-(p )出) e ) f(1)2) ) ) ) ) ) ) ) e ) f(1)2) )
=2e4"一E-T(S )，(一二e )，一e-O-I ) ) e ) T-2 ) "
(b )系统的零状态响应如下
卷积积分的符号计算
Matlab不仅有数值计算功能，还有强大的符号计算功能
可以。 也就是说，是具有推理功能的町得解析式。 卷积积分可以使用Matlab的积分int0函数计算不定积分和定积分。 调用格式为
int(f )，int(f，a，b ) )。
其中，f是被积函数，a、b是积分的上下限。 因此，计算卷积积分解析解的一般函数如下。
、 (f )=，(f )方向) f )=12
e一“eo 3 ) p一2‘’1 () (f—T—1 ) dr )”
=【P一2‘It，- I (打F2 ()的() r 2 )口一2 ) (()，2 ) ) ) ) ) ) ) ) )。
为了以CSCONVS ()函数计算，利用卷积的时间迁移特性，将被积函数变换如下。
万方数据
41