信号与系统实验作业.doc

信号与系统实验 信号的描述、运算、绘图 用MATLAB生成下列函数，连续信号用plot，离散信号用stem绘图 (1) t=linspace(-2,5,1001); y=stepfun(t,0)?; plot(t,y)?; t=linspace(-2,5,1001); y=stepfun(t,0).*sin(t); plot(t,y); t=linspace(-2,5,1001); y=stepfun(t,0).*exp(-t); plot(t,y); t=linspace(-2,5,1001); y=stepfun(t,0).*exp(-t).*cos(t); plot(t,y); t=linspace(-2,2,1001); y=rectpuls(t,2); plot(t,y) t=linspace(-2,5,1001); y=sinc(t/pi*3); plot(t,y); (2) k=[-10:1:10]; y=cos(pi*k/2); stem(k,y) k=[-10:1:10]; y=cos(pi*k/8); stem(k,y) k=[-10:1:10]; y=cos(pi*k/4); stem(k,y) k=[-10:1:10]; y=cos(pi*k/1); stem(k,y) k=[-10:1:10]; y=cos(3*pi*k/2); stem(k,y) k=[-10:1:10]; y=cos(7*pi*k/4); stem(k,y) k=[-10:1:10]; y=cos(15*pi*k/8); stem(k,y) k=[-10:1:10]; y=cos(16*pi*k/8); stem(k,y) k=[-10:1:10]; y=cos(7*pi*k/4+pi/3); stem(k,y) (3) k=[-10:1:10]; y1=cos(pi*k/3); y2=sin(pi*k/3)*j; y=y1+y2; stem(k,y) k=[-10:1:10]; y1=cos(pi*k/3); y2=sin(pi*k/3)*j; y3=(0.5).^k; y=(y3).*(y1+y2); stem(k,y) k=[-10:1:10]; y=exp(pi/3.*k.*j) stem(k,y) 熟悉ezplot、polar、mesh等指令 用ezplot绘的图。 t=linspace(-5,5,1001); ezplot(exp(-t)*cos(t).*stepfun(t,0)) 用mesh绘的三维曲面。的取值范围是 clear;x=-8:0.5:8; y=x; X=ones(size(y))*x; Y=y*ones(size(x)); R=sqrt(X.^2+Y.^2)+eps; Z=sin(R)./R; mesh(X,Y,Z); colormap(hot) xlabel(x),ylabel(y),zlabel(z) 计算数值积分： t=linspace(0,pi,1001); x=exp(-t).*sin(t); y1=sum(x)*t(2); y2=trapz(x)*t(2), f=inline(sin(t)), y3=quad(f,0,pi) y2 = 0.5216 f = Inline function: f(t) = sin(t) y3 = 2.0000 t=linspace(-10*pi,10*pi,1001); x=sinc(2*t/pi); y1=sum(x)*t(502); y2=trapz(x)*t(502), f=inline(sinc(2*t/pi)), y3=quad(f,0,pi) y2 = 1.5549 f = Inline function: f(t) = sinc(2*t/pi) y3 = 0.7091 (2) ，，求 dt=0.01; k1=0:dt:100; f1=cos(3*k1).*stepfun(k1,0);%生成信号f1 k2=k1; f2=stepfun(k2,0)-stepfun(k2,4);%生成信号f2 f=dt*conv(f1,f2);%计算卷积结果f k0=k1(1)+k2(1); %计算序列f 非零样值的起点位置 k3=length(f1)+length(f2)-2; %计算卷积和f 的非零样值的宽度 k=k0:dt:k0+k3*dt; %确定卷积和f 非零样值的时间向量 subplot(2,2,1); plot(k1,f1);title(f1(t));