资源描述：
《利用matlab实现拉普拉斯变换及其逆变换》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、信号与系统实践报告有MATLAB实现连续时间周期函数学院：通信与电子信息工程学院班级：电子042班姓名：李瑞改学号：2004023082指导教师：朱恒军秦月2006年10月19日7摘要：本例的CTFShchsym.m函数文件有一定的通用性，用户只需编写好子函数time_fun_即可，但要注意，该函数是用符号表达式写成的。若要画出时间函数图形，用户需要另外编写一个子函数y=time_fun_e(t)。因为在MATLAB中，只定义了单位阶跃信号Heavisid作为一个符号对象，而不能把Heaviside看作MATLAB的函数加以调用
2、。同理，在信号与系统中，另一个十分重要的函数——单位脉冲函数Dirac(t)。它的使用方法可参照Heaviside进行。最后给出的数值是由完全准确解取32位有效数字后的简洁表示。关键字：单位阶跃信号MATLAB连续信号脉冲宽度目的:1.函数文件CTFStpshsym.m编写源程序72.有MATLAB实现连续时间周期函数一：函数step()将绘出连续系统的阶跃响应函数step()将绘出由向量a和b表示的连续系统的阶跃响应g(t)在指定时间范围内的波形图，并能求出其数值解。和impulse()函数一样，step()函数也有如下四种调
3、用格式：（1）step(b,a)（2）step(b,a,t)（3）step(b,a,t1:p:t2)（4）y=step(b,a,t1:p:t2)上述调用格式的功能和impulse()函数完全相同，所不同的是命令绘制的是系统的阶跃响应g(t)的曲线而不冲激响应h(t)的曲线。对上例，若执行命令step(b,a)则绘制的系统阶跃响应时域波形如图所示。连续系统的冲激响应连续系统的阶跃响应7二：函数文件CTFStpshsym.m编写源程序编写函数文件CTFStpshsym.m[CTFStpshsym.m]function[A_sym,B
4、_sym]=CTFStpshsym%采用符号计算求[0,T]内时间函数的三角级数展开系数，并绘制其双边频谱。%函数的输出为数值量%Nn输出数据的准确位数%A_sym第1元素是直流项，其后元素依次是1,2,3...次谐波cos项展开系数%B_sym第2,3,4,...元素依次是1,2,3...次谐波sin项展开系数%TT=m*tao，信号周期%Nf谐波的阶数%Nn输出数据的准确位数%m(m=T/tao)周期与脉冲宽度之比,如m=4,8,16,100等%tao脉宽:tao=T/msymstnyifnargin<3;Nf=input(
5、'pleasInput所需展开的最高谐波次数:Nf=');endT=input('pleasInput信号的周期T=');ifnargin<5;Nn=32;endy=time_fun_s(t);A0=2*int(y,t,0,T)/T;As=int(2*y*cos(2*pi*n*t/T)/T,t,0,T);Bs=int(2*y*sin(2*pi*n*t/T)/T,t,0,T);A_sym(1)=double(vpa(A0,Nn));fork=1:NfA_sym(k+1)=double(vpa(subs(As,n,k),Nn));7
6、B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn));endifnargout==0S1=fliplr(A_sym)%对A_sym阵左右对称交换S1(1,k+1)=A_sym(1)%A_sym的1*k阵扩展为1*(k+1)阵S2=fliplr(1/2*S1)%对扩展后的S1阵左右对称交换回原位置S3=fliplr(1/2*B_sym)%对B_sym阵左右对称交换S3(1,k+1)=0%B_sym的1*k阵扩展为1*(k+1)阵S4=fliplr(S3)%对扩展后的S3阵左右对称交换回原位置S5=S2-i
7、*S4;%用三角函数展开系数A、B值合成付里叶指数系数S6=fliplr(S5);N=Nf*2*pi/T;k2=-N:2*pi/T:N;S7=[S6,S5(2:end)];%形成-N:N的付里叶复指数对称系数subplot(3,3,3)x=time_fun_e%调用连续时间函数-周期矩形脉冲subplot(3,1,3)stem(k2,abs(S7));%画出周期矩形脉冲的频谱（T=M*tao）title('连续时间函数周期矩形脉冲的双边幅度谱')axis([-80,80,0,0.12])line([-80,80],[0,0])l
8、ine([0,0],[0,0.12])end%-------------------------------------------functiony=time_fun_s(t)%该函数是CTFStpshsym.m的子函数。它由符号变量和表达式写成。symsaa
