有限长序列频谱DTFT地性质Word下载.docx
《有限长序列频谱DTFT地性质Word下载.docx》由会员分享，可在线阅读，更多相关《有限长序列频谱DTFT地性质Word下载.docx（12页珍藏版）》请在冰豆网上搜索。
例如，a=0.5,length=10
a=0.9,length=10
a=0.9,length=20
2-1-3复指数序列
例如，a=0.4,b=0.6,length=20
2-1-4从正弦信号x（t）=sin（2πft+delta）抽样得到的正弦序列x（n）=sin（2πfnT+delta）。
如，信号频率f=1Hz，初始相位delta=0，抽样间隔T=0.1秒，序列长length=10。
2-1-5含两个频率分量的复合函数序列x（n）=sin（2πf1nT）+delta×
sin（2πf2nT+phi）。
如，
频率f1
（Hz）
频率f2
相对振幅
delta
初相位phi
（度）
抽样间隔T
（秒）
序列长
length
1
3
10
90
180
2-2用MATLAB，对上述各个序列，重复如下过程。
2-2-1画出一个序列的实部、虚部、模、相角；
观察并记录实部、虚部、模、相角的特征。
2-2-2计算该序列的幅度谱、频谱实部、频谱虚部；
观察和并记录它们的特征，给予解释。
2-2-3观察同种序列取不同参数时的频谱，发现它们的差异，给予解释。
_______________________________某某：
______________学号：
__________________
2.3用MATLAB，对下面系统进展分析
因果系统
求
（1）求出
画出零极点示意图
（2）画出
和
三、主要仪器设备
MATLAB编程。
四、操作方法和实验步骤
〔参见“二、实验内容和步骤〞〕
五、实验数据记录和处理
列出MATLAB程序清单，加注释。
2-1
2-1-1
N=15;
x=ones（1,N）;
%生成行矩阵
xn=5:
N-1+5;
%横坐标
figure
（1）;
xlabel（'
x（n）'
）;
ylabel（'
n'
stem（xn,2*x）;
%纵坐标变为2倍
axison;
axis（[-11404]）;
set（gca,'
YTick'
[0:
1:
4]）;
XTick'
[-1:
14]）;
title（'
单位阶跃序列'
2-1-2
N=10;
%可修改参数
a=0.5;
xn=0:
N-1;
stem（xn,a.^（xn））;
%画指数函数
axis（[-3N+30a^0]）;
N+2]）;
a=0.5length=10'
2-1-3
N=20;
a=0.4;
b=0.6;
stem（xn,（a+j*b）.^（xn））;
%画复指数函数
axis（[-3N+3-0.5（a+j*b）^0]）;
a=0.4b=0.6length=20'
2-1-4
f=1;
delta=0;
T=0.1;
T:
T*（N-1）;
stem（xn,sin（2*pi*f*xn+delta））;
%画采样后的结果
holdon;
%停留
t=-0.5:
0.001:
1.5;
%画原始函数
plot（t,sin（2*pi*f*t））;
axis（[-0.51.5-1.21.2]）;
gridon;
Sample'
2-1-5
f1=1;
f2=3;
delta=0.5;
phi=0;
stem（xn,sin（2*pi*f1*xn）+delta*sin（2*pi*f2*xn+phi））;
plot（t,sin（2*pi*f1*t）+delta*sin（2*pi*f2*t+phi））;
2-1-5
（1）'
2-2-1
（5）
phi=180;
yn=sin（2*pi*f1*xn）+delta*sin（2*pi*f2*xn+phi）;
yn0=real（yn）;
%实部
yn1=imag（yn）;
%虚部
yn2=abs（yn）;
%模
yn3=angle（yn）;
%角
subplot（2,2,1）;
stem（xn,yn0,'
filled'
real'
subplot（2,2,2）;
stem（xn,yn1,'
imag'
subplot（2,2,3）;
stem（xn,yn2,'
abs'
subplot（2,2,4）;
stem（xn,yn3,'
angle'
六、实验结果与分析
观察实验结果〔数据与图形〕的特征，做必要的记录，做出解释。
包括：
6-1各种序列的图形〔时域〕和频谱〔频域〕各有何特征，给予解释。
其中f1=1,f2=3,delta=0.5,phi=180
实部、虚部、模、角
升级会员 