西安交大数字信号处理第二次实验报告文档格式.docx
- 文档编号:22035117
- 上传时间:2023-02-02
- 格式:DOCX
- 页数:12
- 大小:771.73KB
西安交大数字信号处理第二次实验报告文档格式.docx
《西安交大数字信号处理第二次实验报告文档格式.docx》由会员分享,可在线阅读,更多相关《西安交大数字信号处理第二次实验报告文档格式.docx(12页珍藏版)》请在冰豆网上搜索。
为了使系数为实数,可以将谐振器的共轭复根合并,不失一般性,假设N为偶数,于是可以得到如图1所示的结构。
三、程序及实验结果
1.
源程序:
f0=50;
N=16;
fs=N*f0;
L=2*N;
T=1/fs;
t=0:
T:
(L-1)*T
s=cos(2*pi*0*f0*t+0)+2*cos(2*pi*1*f0*t+pi/2)+cos(2*pi*2*f0*t+pi);
i=0:
1:
L-1;
stem(i,s(i+1))
实验结果:
2.
t=N*T:
x=fft(s);
mag=abs(x);
w=angle(x);
N-1
subplot(2,1,1)
stem(i,mag(i+1));
subplot(2,1,2)
stem(i,w(i+1));
3.
N=16
H=[1,exp(-j*pi*(N-1)/N),0.5*exp(-j*2*pi*(N-1)/N),0.25*exp(-j*3*pi*(N-1)/N),0,0,0,0,0,0,0,0,0,-0.25*exp(-j*13*pi*(N-1)/N),-0.5*exp(-j*14*pi*(N-1)/N),-exp(-j*15*pi*(N-1)/N)];
h=ifft(H);
h
y=abs(H);
w=angle(H);
N-1;
stem(i,y(i+1));
4.
(L-1)*T;
r=0.999;
x=s;
y=CombFilter(x,N,r);
x=y;
A=zeros(1,48);
fori=0:
(N/2)
H=[1,exp(-j*pi*(N-1)/N),0.5*exp(-j*2*pi*(N-1)/N),0.25*exp(-j*3*pi*(N-1)/N),0,0,0,0,0,0,0,0,0,-0.25*exp(-j*13*pi*(N-1)/N),-0.5*exp(-j*14*pi*(N-1)/N),-exp(-j*15*pi*(N-1)/N)];
y=Resonator2(x,N,r,i,H(i+1))
A=A+y;
end
y=A/N;
fori=N:
y1(i-N+1)=y(i+1);
i=0:
subplot(3,1,1);
stem(i,y1(i+1));
y=fft(y1);
mag=abs(y);
w=angle(y);
subplot(3,1,2)
stem(i,mag(i+1));
title('
·
幅频'
);
subplot(3,1,3)
stem(i,w(i+1));
相频'
5.
s0=cos(2*pi*0*f0*t+0);
s1=2*cos(2*pi*1*f0*t+pi/2);
s2=cos(2*pi*2*f0*t+pi);
s=s0+s1+s2;
y=CombFilter(s,N,r);
i=N:
y0=Resonator2(x,N,r,0,H
(1));
subplot(2,2,1);
stem(i,y0(i+1));
y1=Resonator2(x,N,r,1,H
(2));
subplot(2,2,2);
stem(i,y1(i+1));
y2=Resonator2(x,N,r,2,H(3));
subplot(2,2,3);
stem(i,y2(i+1));
y3=Resonator2(x,N,r,3,H(4));
subplot(2,2,4);
stem(i,y3(i+1));
6.
y=zeros;
y1=zeros;
y=CombFilter([1,zeros(1,N-1),1,zeros(1,N-1)],N,r);
(N/2);
A=A+Resonator2(y,N,r,i,H(i+1));
y1=A/N;
y2=zeros;
y2(i-N+1)=y1(i+1);
yfft=fft(y2);
mag=abs(yfft);
w=angle(yfft);
holdon
subplot(2,1,1)
stem(i,mag);
gridon;
subplot(2,1,2)
stem(i,w);
title('
思考题
1.在第2步的幅频特性中,各次谐波的幅度与相应的时域信号幅度有什么关系?
解:
零频分量的幅度等于时域波形中直流的幅度;
直流、二次、三
次谐波的谱线的幅度为时域波形相应谐波对应幅度的一半。
2.实验中为什么要观察第二个周期,如果直接观察第一个周期会怎么样?
第二个周期信号输入不会发生失真,与系统的单位脉冲响应卷积后,输出信号变周期。
直接观察第一个周期,卷积得到的信号不够完整,输出的就不是一个完整的周期信号,输出会发生失真。
3.如果取r=0.95,观察会出现什么情况。
r取0.95结果如图所示:
4.6步与第3步的滤波器幅频特性进行对比,观察并分析二者的关系。
两者均为低通滤波器,第三问的传递函数由内插得到,第六问中用冲激串得到幅频特性。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 西安 交大 数字信号 处理 第二次 实验 报告