信号与系统实验报告
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
信号与系统实验报告
一、 实验内容
1. 利用MATLAB 产生下列连续信号并作图。
(1)51),1(2)(<<---=t t u t x
程序:
t=-1:0.02:5;
x=((t-1)>0);
y=-2*x;
plot(t,y);
axis([-1,5,-3,3]);
title('111000401阶跃函数');
(2)300),3
2sin()(3.0<<=-t t e t x t
t=0:0.001:30;
x=exp(-0.3*t)*sin(2*t/3);
plot(t,x);
=t
-
t
t
x
t
<
+
3000
(<
1.0
cos
,
100
cos
)
t=-0.1:0.001:0.1;
x=cos(100*t)+cos(3000*t);
plot(t,x);
title('111000401x=cos(100*t)+cos(3000*t)')
实验二周期信号的傅里叶级数展开
一、实验内容
1、计算以上所列的四个基本周期信号的傅里叶级数展开。
2、在计算机屏幕上画出周期信号的时域波形。
3、分别计算出5、10、15次谐波迭加的值。
4、把各次谐波迭加的波形重叠画在周期信号的时域波形上。
5、观察逼近情况。
(1)矩形脉冲
function squ_fourier_seq(n)
E=1;
T1=0.5;
Fn=2000;
t0=T1/2;
T=2;
t = -T/2:1/Fn:T/2;
n=5;
squ = E*(1+square(2*pi/T1*t))/2;
squ=cat(2,squ(T1/4*Fn:end),squ(1:T1/4*Fn-1));
a0=E*t0/T1;
N=1:n;
COSIN=cos(2*pi/T1*N'*t);
AN=(2*E)*sin(N*pi*t0/T1)./(N*pi);
AN_1=meshgrid(AN,ones(size(t)))';
COSIN_1=AN_1.*COSIN;
ft=a0+sum(COSIN_1);
plot(t,squ,'r');
plot(t,ft,'b');
grid on;
legend('原始脉冲信号','傅里叶级数逼近'); hold off
title('矩形脉冲5次谐波111000401白明辉')
function squ_fourier_seq(n)
E=1;
T1=0.5;
Fn=2000;
t0=T1/2;
T=2;
t = -T/2:1/Fn:T/2;
squ = E*(1+square(2*pi/T1*t))/2;
squ=cat(2,squ(T1/4*Fn:end),squ(1:T1/4*Fn-1));
a0=E*t0/T1;
N=1:n;
COSIN=cos(2*pi/T1*N'*t);
AN=(2*E)*sin(N*pi*t0/T1)./(N*pi);
AN_1=meshgrid(AN,ones(size(t)))';
COSIN_1=AN_1.*COSIN;
ft=a0+sum(COSIN_1);
plot(t,squ,'r');
hold on
plot(t,ft,'b');
grid on;
legend('原始脉冲信号','傅里叶级数逼近');
hold off
title('矩形脉冲10次谐波111000401白明辉')
function squ_fourier_seq(n)
E=1;
T1=0.5;
Fn=2000;
t0=T1/2;
T=2;
t = -T/2:1/Fn:T/2;
n=15;
squ = E*(1+square(2*pi/T1*t))/2;
squ=cat(2,squ(T1/4*Fn:end),squ(1:T1/4*Fn-1)); a0=E*t0/T1;
N=1:n;
COSIN=cos(2*pi/T1*N'*t);
AN=(2*E)*sin(N*pi*t0/T1)./(N*pi);
AN_1=meshgrid(AN,ones(size(t)))';
COSIN_1=AN_1.*COSIN;
ft=a0+sum(COSIN_1);
plot(t,squ,'r');
hold on
plot(t,ft,'b');
grid on;
legend('原始脉冲信号','傅里叶级数逼近'); hold off
(2)锯齿波
function saw_fourier_seq(n)
E=1;
T1=0.5;
Fn=2000;
T=2;
t = -T/2:1/Fn:T/2;
n=5;
saw= E*(sawtooth(2*pi/T1*t))/2; saw=fliplr(saw);