信号与系统（3）

实验内容

function [t,x]=sconv(x1,x2,t1,t2,dt)
%x1,x2为两个可近似连续信号的离散序列;t1,t2为两个信号的时间范围；dt为时间间隔
x=conv(x1,x2)*dt;%计算序列x1,x2的卷积和
t0=t1(1)+t2(1);%获取序列x非零样值的起点位置
l=length(x1)+length(x2)-2;%计算卷积和x的非零样值的宽度
t=t0:dt:(t0+l*dt);%确定卷积和x的非零样值的时间向量
end
%调用函数
dt=0.1;t1=-5:dt:5;t2=-5:dt:5;   
f1=3*(heaviside(t1+1)-heaviside(t1-2));
f2=2*exp(-2*t2);
[t,f]=sconv(f1,f2,t1,t2,dt);
plot(t,f),title('f2*f1'),grid

fs1=10; %采样频率
t=0:1/fs1:1000/fs1;
y1=-1+2*sin(0.2*pi*t)-3*cos(pi*t);
figure(1);
plot(t,y1);xlabel('时间t');ylabel('幅度值');title('原始信号');
 
fs2=10;
Y=fftshift(fft(y1));  %傅里叶
L=length(y1);  %长度
Y2=abs(Y)/L;  %幅度
X2=(0:L-1)*fs2/L-fs2/2;%频率
figure(2);
plot(X2,Y2);xlabel('频率f/hz');ylabel('幅度');title('幅度谱');