信号与信号处理实验参考答案.docVIP

实验一 熟悉MATLAB环境 2、（2）粗略描绘下列各函数的波形 说明：MATLAB中有函数 ④ f(t)=sint/t t=-3*pi:0.01*pi:3*pi; t1=t/pi; y=sinc(t1); plot(t,y); hold on plot(t,0) ⑤在自己的工作目录work下创建Heaviside函数的M文件，该文件如下： function f=Heaviside(t) f=(t0) %t0时f为1，否则f为0 在命令窗口输入如下语句，就能绘出u(t)的波形。 t=-1:0.01:3; f=Heaviside(t); plot(t,f) axis([-1 3 –0.2 1.2]) ⑥t=-1:0.01:2; g=Heaviside(t)-Heaviside(t-1); plot(t,g); axis([-1 2 -0.2 1.2]) hold on plot(t,0) 分别用for和while循环结构编写程序，求出 s==1+2+22+23+…+262+263 并考虑一种避免循环的简洁方法来进行求和。 程序如下： s=1; for k=1:63 s=s+2^k; end s 运行结果是： s = 8447e+019 （2）s=1; k=1; while k=63 s=s+2^k; k=k+1; end s 运行结果： s = 8447e+019 k=0:63; s=sum(2.^k) 实验二 信号的卷积与系统的响应和阶跃响应 n=0:20; hn=0.9.^n; xn=[0,0 ones(1,8),0,0]; yn=conv(hn,xn); stem(yn) 3. 利用MATLAB绘制下列信号的卷积积分f1(t)*f2(t)的时域波形。 f1(t)=2[u(t+1)-u(t-1)], f2(t)=u(t+2)-u(t-2) f1(t)=tu(t), f2(t)=u(t) f1(t)=u(t)-u(t-4), f2(t)=sin(лt)u(t); f1(t)=e-2tu(t), f2(t)=e-tu(t) (1) 先编写实现连续信号卷积的通用函数sconv(),程序如下： function[f,k]=sconv(f1,f2,k1,k2,p) %计算连续信号卷积积分f(t)=f1(t)*f2(t) %f:卷积积分f(t)对应的非零样值向量 %k:f(t)的对应时间向量 %f1:f1(t)非零样值向量 %f2:f2(t)非零样值向量 %k1:f1(t)的对应时间向量 %k2:序列f2(t)的对应时间向量 %p:取样时间间隔 f=conv(f1,f2); %计算序列f1与f2的卷积f f=f*p; k0=k1(1)+k2(1); %计算序列f的非零样值的起点位置 k3=length(f1)+length(f2)-2; %计算卷积和f的非零样值的宽度 k=k0:p:((k3-(0-k0)/p)*p); %确定卷积和f非零样值的时间向量 subplot(2,2,1) plot(k1,f1) %绘制f1(t) title(f1(t)) xlabel(t) ylabel(f1(t)) subplot(2,2,2) plot(k2,f2) title(f2(t)) xlabel(t) ylabel(f2(t)) subplot(2,2,3) plot(k,f); h=get(gca,position); h(3)=2.5*h(3); set(gca,position,h) %将第三个子图的横坐标范围扩为原来的2.5倍 title(f(t)=f1(t)*f2(t)) xlabel(t) ylabel(f(t)) p=0.01; k1=-1:p:1; f1=2*ones(1,length(k1)); k2=-2:p:2; f2=ones(1,length(k2)); [f,k]=sconv(f1,f2,k1,k2,p) （2）p=0.01; k1=0:p:10; f1=k1 k2=0:p:10; f2=ones(1,length(k2)); [f,k]=sconv(f1,f2,k1,k2,p) 第（2）题图上 实验二 信号的卷积与系统的响应 1． n=0:20;hn=0.9.^n; xn=stepseq(2,0,20)-stepseq(10,0,20); yn=conv(hn,xn); stem(yn) 2．（1） p=0.01;k1=-2:p:2; f1=2*(u(k1+1)-u(k1-1)); k2=k1; f2=u(k2+2)-u(k2-2); [f,k]=sconv(f1,f2,k1,k2,p) （2） p=0.01;k1=-1:p:1