电子科技大学信号与系统上机实验报告docWord文档格式.docx
《电子科技大学信号与系统上机实验报告docWord文档格式.docx》由会员分享，可在线阅读，更多相关《电子科技大学信号与系统上机实验报告docWord文档格式.docx（17页珍藏版）》请在冰豆网上搜索。
subplot（3,1,1）;
stem（n,x1）
title（'
figuresofsignals'
）
xlabel（'
n'
ylabel（'
x1[n]'
subplot（3,1,2）;
stem（n,x2）
x2[n]'
subplot（3,1,3）;
stem（n,x3）
x3[n]'
（二）实验图像：
三、实验结论
x1的信号周期是4，x2信号周期是4，x3周期是16
第二题
一、实验目的：
1.确定离散时间信号的性质
2.证明不满足线性性质
3.证明
不是因果的
（一）
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）
figureofy3'
y3'
subplot（2,1,2）;
stem（n,y4）
figureofy4'
y4'
2.实验图像
（二）
1.实验程序：
nx=[-5:
9];
x1=[zeros（1,5）1ones（1,9）];
x2=[zeros（1,4）1ones（1,10）];
y=x1+x2;
stem（nx,x1）;
）;
stem（nx,x2）;
subplot（3,1,3）;
stem（nx,y）;
y[n]'
2.实验图象：
1.系统不是线性的。
2.系统
不是因果的。
第三题
用conv函数计算y，并验证与解析法所求是否一致
5];
x=[1ones（1,5）];
h=n;
subplot（2,2,1）;
stem（n,x）
figures'
x'
subplot（2,2,2）;
stem（n,h）
h'
y=conv（x,h）;
ny=[0:
10];
stem（ny,y）
ny'
y'
3.解析法所求
函数所求结果与解析法所求一致
第四题
1.学习使用filter函数，并与上题比较结果是否一致
2.用freqz函数求差分方程的频率响应
3.用freqs函数画出微分方程的频率响应
x=[ones（1,6）zeros（1,5）];
h=[012345zeros（1,5）];
b=h;
y=filter（b,1,x）;
n=[1:
11];
subplot（1,1,1）;
stem（n,y）
b=[20-1];
a=[1-0.8];
[Homega]=freqz（b,a,4）
2.实验结果：
>
H=
5.0000
2.8200-1.3705i
1.8293-1.4634i
0.9258-0.9732i
omega=
0
0.7854
1.5708
2.3562
（三）
b=[01];
a=[13];
w=linspace（0,3*pi）;
H=freqs（b,a,w）;
plot（H）
频率响应'
三、实验结论：
1.用filter函数所得结果，与上题一致
2.用freqz函数求差分方程的频率响应如上
3.用freqs函数画出微分方程的频率响应图如上
第五题
1.掌握fft和ifft函数的用法。
2.观察虚部与实部相比能否忽略。
3.思考虚部成分是由什么引入的
1.程序：
N=32;
x=[11zeros（1,30）];
a=（1/N）*fft（x）;
stem（n,real（a））
Realpartofa'
real（a）'
stem（n,imag（a））
Imagepartofa'
image（a）'
x=N*ifft（a）;
stem（n,real（x））
Therealpartofx'
real（x）'
stem（n,imag（x））
Theimagpartofx'
imag（x）'
1.图象如上图所示
2.虚部与实部相比不能忽略
第六题
1.求频率响应，并用a和b分别表示其分母与分子。
2.用residue函数求出频率响应的部分分式展开式。
3.重新组合展开式，确认分式展开的正确性。
b=[1-2]
a=[13/21/2]
[r,p,k]=residue（b,a）
b=
1-2
a=
1.00001.50000.5000
r=
6
-5
p=
-1.0000
-0.5000
k=[]
1.展开分式是正确的
2.数据如上
第七题
1.用roots求每个函数的零极点，用plot画出零极点图。
2.用zplane函数画出零极点图。
（A）
b1=[15];
a1=[123];
zs1=roots（b1）;
ps1=roots（a1）;
plot（real（zs1）,imag（zs1）,'
o'
holdon
plot（real（ps1）,imag（ps1）,'
grid
axis（[-63-63]）
H1零极点图'
b2=[2512];
a2=[1210];
zs2=roots（b2）;
ps2=roots（a2）;
plot（real（zs2）,imag（zs2）,'
plot（real（ps2）,imag（ps2）,'
axis（[-22-55]）
H2零极点图'
b3=[2512];
a3=[141420];
zs3=roots（b3）;
ps3=roots（a3）;
plot（real（zs3）,imag（zs3）,'
plot（real（ps3）,imag（ps3）,'
axis（[-33-44]）
H3零极点图'
（B）
b=[125];
a=[1-30];
zplane（b,a）
如图所示
第八题
1.用解析法求出傅立叶级数系数表达示，并画出系数的图
2.画出出有限项近似的波形
k=[-80:
80];
a（k~=0）=sin（pi*k（k~=0）/2）./（k（k~=0）*pi）;
a（k==0）=0.5;
stem（k,a）
ak'
k'
t=[-3;
0.01;
3]
升级会员 