连续信号积分 matlab程序,连续信号的微积分与卷积.ppt

实验15 连续信号的微积分与卷积 实验目的 学习Matlab的符号运算功能及编程方法。 掌握用Matlab计算微积分的数值方法和解析方法。并应用于求解信号的功率或能量。 掌握用Matlab进行卷积运算的数值方法和解析方法。加深对卷积积分的理解。 实验原理与说明 微分和积分的数值解连续信号的微分是用差分来近似的，当步长(时间间隔)越小时，用差分表示微分就越精确。见图15-1所示。 实验原理与说明 微分和积分的解析解Matlab除了数值计算以外，还有强大的符号运算功能。在数值计算过程中，参与运算的变量都是被赋了值的数值。而在符号运算的整个过程中，参与运算的是符号变量。在符号运算中所出现的数字都是当做符号来处理的。Matlab中，对符号表达式微分的函数是diff()。利用这个函数，可以求符号表达式的一阶导数、n阶导数。该函数有三种调用格式diff(f), diff(f,a),diff(f,n)或diff(f,a,n)其中，f为符号表达式，a说明对其求导，n表示求导次数。积分运算的函数是int()，它也有三种调用形式int(f),int(f,v),int(f,a,b)或int(f,v,a,b)其中，f为符号表达式，v说明对其求积分，a,b 表示积分的区间。显然前两个是求不定积分，后两个是求定积分。 实验原理与说明 卷积积分的计算利用系统的冲激响应和叠加原理来求系统对任意激励信号作用时的零状态响应，这就是卷积方法的原理。因此，在时域内，卷积是求解线性非时变系统零状态响应的重要方法。特别是激励信号为时限信号时尤其如此。卷积积分计算实际上可用信号的分段求和来实现，即若只求当(n为整数)时的值，则由上式可得：当时间间隔足够小时，就是的数值近似。Matlab的conv(x,h)函数可以用来计算卷积。用Matlab的符号计算方法int函数还可以计算卷积的解析式。 计算示例1 计算示例2 计算示例2 计算示例2 计算示例2 计算示例3 计算示例4 计算示例4 计算示例5 计算示例5 实验内容1 实验内容2 实验内容 3 实验内容 4 实验步骤与方法 实验报告要求 根据求出的数学模型所编写出的程序。以及绘出的各种波形图。 上机调试程序的方法。 根据实验归纳、总结出用Matlab计算积分和卷积的方法。 心得体会及其他。 * * 故求导数就是近似求差分与步长之比。 图15-1 Matlab中用diff函数来计算差分 其调用格式为y=diff(f) 连续信号的定积分可由Matlab中的quad函数和Quadl来实现，其调用格式为 quad(‘function_name’,a,b) 采用自适应Simpson算法quadl(‘function_name’,a,b) 采用自适应Lobatto算法 其中，function_name为被积函数名，a,b为指定的积分区间。 计算定积分 解： 用三种积分函数计算如下>> format long; >> y=inline('exp(-0.5*t).*sin(t+pi/6)'); >> s=quad(y,0,3*pi) s = 0.100646 >> s1=quadl(y,0,3*pi) s1 = 0.775646 >> sf=int('exp(-0.5*t)*sin(t+pi/6)',0,3*pi) sf = . 符号积分计算是精确值。 画出如图15-2所示信号的一阶导数和积分的波形图。 图15-2 解: (a)用Matlab的数值计算方法的程序如下 % 画微分和积分的程序(数值计算) exp15_1.m t0=-3;t1=3;dt=0.01; t=t0:dt:t1; f1=tripuls(t,4,0.5);% 定义三角波 max_f=max(f1);min_f=min(f1); df=diff(f1)/dt;% 求导 max_df=max(df);min_df=min(df); f=inline('tripuls(t,4,0.5)');% 定义在线函数对象 for x=1:length(t)intf(x)=quad(f,-3,t(x));% 求积分 end max_if=max(intf);min_if=min(intf); subplot(3,1,1),plot(t,f1,'linewidth',2);% 画f(t) grid;line([t0 t1],[0 0]);ylabel('f(t)') axis([t0,t1,min_f-0.2,max_f+0.2]) subplot(3,1,2),plot(t(1:length(t)-1),df,'linewidth',2); grid;line([t0 t1],[0 0]);ylabel('df(t)/dt') % 画f