信号与系统实验报告.docx

姓名:刘伟权班级:12应用师1班学号20********** 实验一:产生信号波形的仿真实验
实验过程与结果:(1)用MATLAB软件中的funtool符合计算方法(图示化函数计算器)来产生并表示信号。
学生实验内容为产生以下信号波形:3sin(x)、5exp(-x)、sin(x)/x、1-2abs(x)/a、sqrt(a*x).
波形结果如下:x=[-2*pi, 2*pi] a= 1/2
(2)用MATLAB软件的信号处理工具箱来产生并表示信号向量表法。
t=-10::10;
f=sin(t)./t;
plot(t,f);
title(‘f(t)=Sa(t)’);
xlabel(‘t’);
axis([-10,10,-,]);
(3)离散时间序列波形产生并绘制
•f1(k)=cos(K*PI/8)
代码:
k=-8:8;
f=cos(pi/8*k);
stem(k,f,'filled');
axis([-4,4,-,]);
波形:
•f2(k)=cos(2k)
代码如下:
k=-8:8;
f=cos(2*k);
stem(k,f,'filled');
axis([-4,4,-,]); 波形:
实验二连续时间信号卷积及MATLAB实现
实验过程与结果:(1)输入程序代码function [f,k]=sconv(f1,f2,k1,k2,p)
f=conv(f1,f2);
f=f*p;
k0=k1(1)+k2(1);
k3=length(f1)+length(f2)-2;
k=k0:p:k3*p;
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)
plot(k,f);
h=get(gca,'position');
h(3)=*h(3);
set(gca,'position',h)
title('f(t)=f1(t)*f2(t)')
xlabel('t')
ylabel('f(t)')另存为文件名为matlab的文件
再新建空白输入:
p=;
k1=0:p:2;
f1=*k1;
k2=k1;
f2=f1;
[f,k]=sconv(f1,f2,k1,k2,p)另存数据后执行
结果:
(2)将(1)的程序代码中的k0改为0,其余不变
再输入:p=;
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)
结果:
实验三系统时域特性的仿真实验
实验过程与结果:(1)连续系统的冲激响应、阶跃响应
学生实验内容:绘出连续系统微分方程在0~
输入程序a=[2 1 8]; a=[2 1 8];
b=[1]; b=[1];
impulse(b,a,0::30) step(b,a,0::30)
结果:
(2)连续系统的单位脉冲响应
绘出离散时间系统分别以默认方式、0~60取样点范围内、-10~40离散时间范围内的波形图和-5~10离散时间范围内的数值解。
输入程序:a=[1 -1 ]; a=[1 -1 ]; a=[1 -1 ];
b=[1] b=[1]; b=[1];
impz(b,a) impz(b,a,0:60) impz(b,a,-10:40)
波形如下
输入 a=[1 -1 ];
b=[1];
y=impz(b,a,-10:40)
结果 y =
0
0
0
0
0

-


-


-


-
(3)对如下连续时间系统,通过仿真观察其在周期矩形信号作用下的零状态响应波形。
结果