资源描述：
《MATLAB在信号处理中的应用中》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第一讲信号及其运算的MATLAB表示1.1连续信号1、指数信号形式：y=A*exp（a*t）2、正弦信号形式：y=A*cos（w0*t+f）或：y=A*sin（w0*t+f）13、抽样函数用sinc函数表示，其定义为：sinc(t)=sin(πt)/(πt)MATLAB语言调用形式:y=sinc(t)t=-3*pi:pi/100:3*pi;y=sinc(t/pi);plot(t,y),gridon24、矩形脉冲信号MATLAB语言调用形式:y=rectpuls(t,width)产生一个幅值为1,宽度为width、相对于t=0点左右对称的矩形波信号t=0:0.0
2、01:4;T=1;y=rectpuls(t-2*T,2*T);plot(t,y),gridonaxis([0,4,-0.5,1.5])以t-2*T=0,即t=2*T为对称中心3周期性矩形脉冲信号MATLAB语言调用形式:y=square(t,duty)产生一个周期为2π,幅值为±1的周期性方波,duty表示占空比(dutycycle)t=-0.0625:0.0001:0.0625;y=square(2*pi*30*t,75);plot(t,y),gridonaxis([-0.1,0.1,-1.5,1.5])占空比75%45、三角脉冲信号MATLAB语言调用形式
3、:y=tripuls(t,width,skew)产生一个最大幅值为1,宽度为width、斜度为skew的三角波信号.t=-3:0.0001:3;y=tripuls(t,4,0.5);plot(t,y),gridonaxis([-3,3,-1.5,1.5])5周期性三角脉冲信号MATLAB语言调用形式:y=sawtooth(t,width)产生一个周期为2π,最大幅值为1,最小幅值为-1周期性三角波信号.t=-5*pi:0.0001:5*pi;y=sawtooth(t,0.5);plot(t,y),gridonaxis([-16,16,-1.5,1.5])66、
4、一般周期性脉冲信号MATLAB语言调用形式:y=pulstran(t,d,’func’)依据名为func的连续函数并以之为一个周期,从而产生一串周期性的连续函数。用于指定周期性的偏移量，这个func函数会被计算length(d)次,最后值相当于:y=func(t-d(1))+func(t-d(2))+...更一般的调用形式为:y=pulstran(t,d,’func’,p1,p2)。其中p1,p2为需要传送给func函数的额外输入参数值,最后值相当于:y=func(t-d(1),p1,p2,...)+func(t-d(2),p1,p2,...)+...7t=0
5、:0.0001:1;d=0:1/3:1;y=pulstran(t,d,'rectpuls',0.1);figure(1),plot(t,y),gridonaxis([0,1,-0.1,1.5])yy=pulstran(t,d,'tripuls',0.1,-1);figure(2),plot(t,yy),gridonaxis([0,1,-0.1,1.5])81.2离散信号1、指数序列a.^k2、正弦序列离散的指数序列与正弦序列用MATLAB的表示与连续信号类似,只不过用stem函数而不是用plot函数来画出序列的波形.93、单位冲激序列单位冲激序列k=-50:5
6、0;y=[zeros(1,50),1,zeros(1,50)];stem(k,y)function[y,k]=impseq(k0,k1,k2);k=k1:k2;y=[k-k0==0]写成MATLAB函数文件形式若δ[k-k0]在k1≤k≤k2范围内104、单位阶跃序列单位阶跃序列k=-50:50;y=[zeros(1,50),ones(1,51)];stem(k,y),axis([-60,60,-0.5,1.5])写成函数文件形式若u[k-k0]在k1≤k≤k2function[y,k]=stepseq(k0,k1,k2);k=k1:k2;y=[k-k0>=0
7、]111.3信号运算的MATLAB实现2、离散序列的差分与求和1、信号的尺度变换、反转和平移差分∇f(k)=f(k)-f(k-1)用diff函数求和用sum函数123、连续信号的微分与积分微分可用diff函数积分用quad函数,调用形式:quad(‘function_name’,a,b)例:对于三角波，画出df(t)/dt和t-∞f(t)dt的波形functiony=funtrip(t)y=tripuls(t,4,0.5);13h=0.001,t=-3:h:3;y1=diff(funtrip(t))*1/h;figure(1),plot(t(1:length
8、(t)-1),y1);title('d
