资源描述:
《信号与信号处理实验参考答案_新0_20070611》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、实验一熟悉MATLAB环境2、(2)粗略描绘下列各函数的波形说明:MATLAB中有函数④f(t)=sint/tt=-3*pi:0.01*pi:3*pi;t1=t/pi;y=sinc(t1);plot(t,y);%输出一个波形,以t为纵坐标,y函数为横坐标。holdon%不用的话前面的图会消失的,要在第一个波形后面用。plot(t,0)%输出一个参考坐标轴。⑤在自己的工作目录work下创建Heaviside函数的M文件,该文件如下:functionf=Heaviside(t)f=(t>0)%t>0时f为1,否则f为0在命令窗口输入如下语句,就能绘出u(t)的波形。t=-1:0.01:3
2、;f=Heaviside(t);plot(t,f)%绘图的加不加分号都行。axis([-13–0.21.2])%只显示前面两个是横坐标范围,后面两个纵坐标范围内的图形。⑥t=-1:0.01:2;g=Heaviside(t)-Heaviside(t-1);plot(t,g);axis([-12-0.21.2])holdonplot(t,0)4、分别用for和while循环结构编写程序,求出s==1+2+22+23+…+262+263并考虑一种避免循环的简洁方法来进行求和。程序如下:s=1;fork=1:63%么有逗号。s=s+2^k;ends运行结果是:s=1.8447e+019(2)
3、s=1;k=1;whilek<=63%没有分号。s=s+2^k;k=k+1;ends运行结果:s=1.8447e+019(3)k=0:63;s=sum(2.^k)实验二信号的卷积与系统的响应和阶跃响应1.n=0:20;hn=0.9.^n;xn=[0,0ones(1,8),0,0];%离散的所以是用矩阵赋值,连续的用HEAVISIDE。yn=conv(hn,xn);%求卷积,系统在前,输入在后。stem(yn)%离散输出波形。3.利用MATLAB绘制下列信号的卷积积分f1(t)*f2(t)的时域波形。13(1)f1(t)=2[u(t+1)-u(t-1)],f2(t)=u(t+2)-u(
4、t-2)(2)f1(t)=tu(t),f2(t)=u(t)(3)f1(t)=u(t)-u(t-4),f2(t)=sin(лt)u(t);(4)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:取样
5、时间间隔f=conv(f1,f2);%计算序列f1与f2的卷积ff=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)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
6、)plot(k,f);h=get(gca,'position');h(3)=2.5*h(3);set(gca,'position',h)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=k1k2=0:p:10;f2=ones(1,length(k2));[f,k]=sconv(f1,f2,k1,
7、k2,p)第(2)题图上,第(3)题图下,5.已知描述某连续系统的微分方程为:y’’(t)+5y’(t)+8y(t)=3f’’(t)+2f(t)13绘出系统的冲激响应波形,求出t=0.5s,1s,1.5s,2s系统冲激响应的数值解。6、step()函数的调用格式与impulse()函数类似。已知描述某连续系统的微分方程为:2y’’(t)+y’(t)+8y(t)=f(t)试用MATLAB绘出该系统的冲激响应与阶跃响应的波形。5.程序如下:a=[158];b=