1.信号的MATLAB表示与信号的基本运算.pptVIP

信号的MATLAB表示 指数信号Aeat y = A*exp(a*t); 指数序列ak 幂运算a.^k实现 正弦型信号 内部函数cos( ) 和sin( ) 抽样函数Sa(t) sinc(t) 矩形脉冲信号 y = rectpuls(t,width) 三角波脉冲信号 y = tripuls(t, width,skew) %decaying exponential t=0:001:10; A=1; a=-0.4; ft=A*exp(a*t); plot(t,ft) % rectpuls t=0:0.001:4; T=1; ft=rectpuls(t-2*T,T); plot(t,ft) axis([0,4,-0.5,1.5]) % tripuls t=-3:0.001:3; ft=tripuls(t,4,0.5); plot(t,ft) % unit impuls sequence k=-50:50; delta=[zeros(1,50),1,zeros(1,50)]; stem(k,delta) % unit step sequence k=-50:50; uk=[zeros(1,50), ones(1,51)]; stem(k,uk) 信号基本运算的MATLAB实现 1．信号的尺度变换、翻转、时移（平移） t=-3:0.001:3; ft1=tripuls(2*t,4,0.5); subplot(2,1,1) plot(t,ft1) title(f(2t)) ft2=tripuls((2-2*t),4,0.5); subplot(2,1,2) plot(t,ft2) title(f(2-2t)) 已知三角波f(t)，用MATLAB画出的f(2t)和f(2-2t) 波形 信号的相加与相乘 相加用算术运算符“+”实现 相乘用数组运算符“.*”实现 例画信号Aeatcos(w0t+f)的波形 t=0:0.001:8; A=1; a=-0.4; w0=2*pi;phi=0; ft1=A*exp(a*t).*sin(w0*t+phi); plot(t,ft1) 离散序列的差分与求和 差分 y=diff(f); 求和 y=sum(f(k1:k2)); [例]已知三角波f(t)，画出其微分与积分的波形 function yt=f2_2(t) yt=tripuls(t,4,0.5); %differentiation h=0.001;t=-3:h:3; y1=diff(f2_2(t))*1/h; plot(t(1:length(t)-1),y1) title(df(t)/dt) %integration t=-3:0.1:3; for x=1:length(t) y2(x)=quad(f2_2,-3,t(x)); end plot(t,y2) title(integral of f(t)) 三角波f(t)微分与积分的波形 * * t=0:0.1:10; A=1; a=-0.4; ft=A*exp(a*t); stem(t,ft) ft=tripuls(t,4,1); function [f,k]=impseq(k0,k1,k2) %产生 f[k]=delta(k-k0)；k1=k=k2 k=[k1:k2];f=[(k-k0)==0]; k0=0;k1=-50;k2=50; [f,k]=impseq(k0,k1,k2); stem(k,f) function [f,k]=stepseq(k0,k1,k2) %产生 f[k]=u(k-k0);k1=k=k2 k=[k1:k2];f=[(k-k0)=0]; k0=0;k1=-50;k2=50; [f,k]=stepseq(k0,k1,k2); stem(k,f) 连续信号的微分与积分 微分 y=diff(f)/h; h为数值计算所取时间间隔 定积分 quad(‘function_name’,a,b); function_name为被积函数名，a和b指定积分区间。 * * * * *