信号与系统上机实验new.pptVIP

信号与系统第1次上机实验 信号基本运算的Matlab实现 1、尺度变换 2、翻转 3、时移 4、求导 5、积分 6、差分 7、求和 信号的尺度变换、翻转、平移运算，实际上是函数自变量的运算。 尺度变换中只要自变量乘以一个常数，可用算术符号“*”来表示。 在信号翻转中自变量乘以一个负号，可以在MATLAB中直接写出，也可用函数fliplr(f)函数实现。 在信号时移中函数自变量加减一个常数，在MATLAB中用算术运算符“＋”或“－”来实现。 例2.1对一三角波f(t)，画出f(2t) 和f(-2t+1)的波形。 分析：f(t)?f(-2t+1)的步骤 clc Clear,close all t=-3:0.01:3; ft1=tripuls(t,4,0.5); subplot(2,2,1),plot(t,ft1);title(ft1) ft2=tripuls(2*t,4,0.5); ft3=tripuls(-2*t,4,0.5); ft4=tripuls(-2*t+1,4,0.5); subplot(2,2,2),plot(t,ft2);title(ft2) subplot(2,2,3),plot(t,ft3);title(ft3) subplot(2,2,4),plot(t,ft4);title(ft4) 信号的导数和积分 1、连续信号的导数和积分 diff函数用来求信号的导数 格式：y=diff(f(x)) /h,h为步长。 例如： quad函数用来求连续信号的定积分。 格式： q = quad(fun,a,b)。 %求导数 clc clear; close all h=0.01; t=-3:h:3; f2t=tripuls(t,4,0.5); ft2=diff(f2t)/h; plot(t(1:length(t)-1),ft2) axis([-3 3 -1.5 0.5]) title(导数) %求积分 t=-3:0.1:3; for x=1:length(t) y2(x)=quad(f2,-3,t(x)); end subplot(2,2,2),plot(t,y2,r) axis([-3 3 -0.2 2.2]) title(‘积分’) function yt=f2(t) yt=tripuls(t,4,0.5) 2、离散序列的差分与求和 离散序列的差分?f[k]=f[k]-f[k-1],在Matlab中用diff函数来实现，格式为： y=diff(f) 离散序列的求和 练习 ?1、? （1）编写如图Exercise1.1所示波形的MATLAB函数。 试画出f(t),f(0.5t),f(1-2t)的波形。 clc clear,close all t=-5:0.01:5 y=rectpuls(t,2)+tripuls(t-1.5,1,-1); figure,plot(t,y,r),axis([-5 5.2 0 1.3]),title(原始图形) 2、画出如图exercise2.2所示序列f[2k]、f[k+2],f[k-2]的波形。并求f[k]的和 ％用Matlab表达x(n),并求x(n)的和 clc clear,close all x=zeros(1,21);y1=zeros(1,21);y2=zeros(1,21);y3=zeros(1,21);y5=zeros(1,11); for n=-10:10 switch n case -2 x(n+11)=4; case -1 x(n+11)=1 case 0 x(n+11)=-2; case 1 x(n+11)=-1; case 2 x(n+11)=2; case 3 x(n+11)=-3; otherwise x(n+11)=0; end end xsum=sum(x(1:21)) %求x(2n) for n=1:21 k=(n+1)/2; if mod(n,2)==1 y5(k)=x(n); end end ％绘图 n=-10:10; subplot(2,2,1);stem(n,x);title(x[n]);xlabel(n);ylabel(x);grid on subplot(2,2,2);stem(n,y1);title(x[n-2]);xlabel(n);ylabel(x[n-2]);grid