电子科技大学信号与系统上机实验报告 3.docx
《电子科技大学信号与系统上机实验报告 3.docx》由会员分享，可在线阅读，更多相关《电子科技大学信号与系统上机实验报告 3.docx（23页珍藏版）》请在冰点文库上搜索。
电子科技大学信号与系统上机实验报告3
信号与系统
实
验
报
告
姓名：
学号：
学院：
指导教师：
实验时间：
1第一题
1.1实验目的
画出离散时间正弦信号并确定基波周期
1.2实验程序
n=[0:
31];
x1=sin（pi*n/4）.*cos（pi*n/4）;
x2=（cos（pi*n/4））.^2;
x3=sin（pi*n/4）.*cos（pi*n/8）;
subplot（3,1,1）;stem（n,x1）
title（'figuresofsignals'）
xlabel（'n'）
ylabel（'x1[n]'）
subplot（3,1,2）;stem（n,x2）
xlabel（'n'）
ylabel（'x2[n]'）
subplot（3,1,3）;stem（n,x3）
xlabel（'n'）
ylabel（'x3[n]'）
1.3实验图像
1.4实验结论
可得信号x1[n]与x2[n]的周期均为4，x3[n]的周期为16.
2第二题
2.1实验目的
确定离散时间系统的性质
证明不满足线性性质
证明
不是因果的
2.2实验程序
2.2.1系统1
x1=[1zeros（1,10）];
x2=[2zeros（1,10）];
y1=sin（（pi/2）.*x1）;
y2=sin（（pi/2）.*x2）;
x3=2*x1+6*x2;
y3=sin（（pi/2）.*x3）;
y4=2*y1+6*y2;
subplot（2,1,1）;stem（n,y3）
title（'figureofy3'）
xlabel（'n'）
ylabel（'y3'）
subplot（2,1,2）;stem（n,y4）
title（'figureofy4'）
xlabel（'n'）
ylabel（'y4'）
2.2.2系统2
nx=[-5:
9];
x1=[zeros（1,5）1ones（1,9）];
x2=[zeros（1,4）1ones（1,10）];
y=x1+x2;
subplot（3,1,1）;
stem（nx,x1）;
xlabel（'n'）;
ylabel（'x1[n]'）;
subplot（3,1,2）;
stem（nx,x2）;
xlabel（'n'）;
ylabel（'x2[n]'）;
subplot（3,1,3）;
stem（nx,y）;
xlabel（'n'）;
ylabel（'y[n]'）;
2.3实验图像
2.3.1系统1
2.3.2系统2
2.4实验结论
系统
不是线性的。
系统
不是因果的。
3第三题
3.1实验目的
计算卷积
3.2实验程序
[x,n1]=stepseq（0,0,5）;
h=[012345];
y=conv（x,h）;
3.3实验图像
3.4实验结论
计算得x[n]和h[n]的卷积如图所示。
4第四题
4.1实验目的
使用filter函数，并与上题比较结果是否一致
用freqz函数求差分方程的频率响应
用freqs函数画出微分方程的频率响应
4.2实验程序
4.2.1程序1
h=[012345];
x2=[11111100000];
y2=filter（h,1,x2）;
stem（y2）
4.2.2程序2
b=[20-1];a=[1-0.8];
[Homega]=freqz（b,a,4）
4.2.3程序3
a=[13];b=3;
figure;freqs（b,a）;
w=linspace（0,3*pi）;
h=freqs（b,a,w）;
figure;plot（w,abs（h））;
4.3实验图像或输出
4.3.1程序1
4.3.2程序2
>>
H=
5.0000
2.8200-1.3705i
1.8293-1.4634i
0.9258-0.9732i
omega=
0
0.7854
1.5708
2.3562
4.3.3程序3
4.4实验结论
用filter函数所得结果，与上题一致
用freqz函数求差分方程的频率响应如上
用freqs函数画出微分方程的频率响应图如上
5第五题
5.1实验目的
计算离散时间傅里叶级数。
5.2实验程序
n=[0:
31];
N=32;
x=[11zeros（1,30）];
a=（1/N）*fft（x）;
subplot（2,1,1）;stem（n,real（a））
title（'Realpartofa'）
xlabel（'n'）
ylabel（'real（a）'）
subplot（2,1,2）;stem（n,imag（a））
title（'Imagepartofa'）
xlabel（'n'）
ylabel（'image（a）'）
5.3实验图像
5.4实验结论
函数实部与虚部的图像如图所示。
6第六题
6.1实验目的
验证LTI系统的特征函数
6.2实验程序
n=-20:
100;
x1=exp（i*0.25*pi*n）;
x2=sin（pi*n/8+pi/16）;
x3=power（0.9,n）;
x4=n+1;
figure
subplot（2,1,1）
stem（real（x1）,'.'）
subplot（2,1,2）
stem（imag（x1）,'.'）
figure
stem（x2,'.'）
figure
stem（x3,'.'）
figure
n=0:
100;
x1=exp（i*0.25*pi*n）;
x2=sin（pi*n/8+pi/16）;
x3=power（0.9,n）;
x4=n+1;
a=[1,-0.25];
b=[1,0.9];
y1=filter（b,a,x1）;
y2=filter（b,a,x2）;
y3=filter（b,a,x3）;
y4=filter（b,a,x4）;
figure
subplot（2,1,1）
stem（real（y1）,'.'）
subplot（2,1,2）
stem（imag（y1）,'.'）
figure
stem（y2,'.'）
figure
stem（y3,'.'）
figure
stem（y4,'.'）
stem（x4,'.'）
h1=y1./x1;
h2=y2./x2;
h3=y3./x3;
h4=y4./x4;
figure
subplot（2,1,1）
stem（real（h1）,'.'）
title（'therealpartofh1'）
xlabel（'n'）
ylabel（'h1'）
subplot（2,1,2）
stem（imag（h1）,'.'）
title（'theimagpartofh1'）
xlabel（'n'）
ylabel（'h1'）
figure
stem（h2,'.'）
xlabel（'n'）
ylabel（'h2'）
figure
stem（h3,'.'）
xlabel（'n'）
ylabel（'h3'）
figure
stem（h4,'.'）
xlabel（'n'）
ylabel（'h4'）
6.3实验图像
x1
x2
x3
x4
y1
y2
y3
y4
h1
h2
h3
h4
6.4实验结论
x1,x2,x3,x4,y1,y2,y3,y4,h1,h2,h3,h4的图像分别如图所示。
7第七题
7.1实验目的
用离散时间傅里叶级数综合信号
7.2实验程序
x1=[ones（1,8）ones（1,8）ones（1,8）ones（1,8）ones（1,8）ones（1,8）ones（1,8）ones（1,8）];
x2=[ones（1,8）zeros（1,8）ones（1,8）zeros（1,8）ones（1,8）zeros（1,8）ones（1,8）zeros（1,8）ones（1,8）zeros（1,8）];
x3=[ones（1,8）zeros（1,24）ones（1,8）zeros（1,24）];
stem（x1）
title（'singalx1'）
xlabel（'n'）;
ylabel（'x1'）
stem（x2,'.'）
title（'singalx2'）
xlabel（'n'）;
ylabel（'x2'）
stem（x3,'.'）
title（'singalx3'）
xlabel（'n'）;
ylabel（'x3'）
x1=[ones（1,8）];
x2=[ones（1,8）zeros（1,8）];
x3=[ones（1,8）zeros（1,24）];
a1=1/8*fft（x1）;
a2=1/16*fft（x2）;
a3=1/64*fft（x3）;
stem（abs（a1））
title（'a1Amplitude'）
xlabel（'k'）;
ylabel（'abs（a1）'）;
stem（abs（a2））
title（'a2Amplitude'）
xlabel（'2'）;
ylabel（'abs（a2）'）;
stem（abs（a3））
title（'a3Amplitude'）
xlabel（'k'）;
ylabel（'abs（a3）'）;
7.3实验图像
升级会员 