数字信号处理实验报告(自己的实验报告).doc

数字信号处理
实验报告
西南交通大学信息科学与技术学院
姓名:伍先春
学号:20092487
班级:自动化1班
指导老师:张翠芳
实验一序列的傅立叶变换
实验目的
进一步加深理解DFS,DFT算法的原理;
研究补零问题;
快速傅立叶变换(FFT)的应用。
实验步骤
复习DFS和DFT的定义,性质和应用;
熟悉MATLAB语言的命令窗口、编程窗口和图形窗口的使用;
利用提供的程序例子编写实验用程序;
按实验内容上机实验,并进行实验结果分析;
写出完整的实验报告,并将程序附在后面。
实验内容
周期方波序列的频谱
已知
试画出下面四种情况下的的幅度频谱,并分析补零后,对信号频谱的影响。
有限长序列x(n)的DFT
已知:
取x(n)(n=0:10)时,画出x(n)的频谱X(k) 的幅度;
将(1)中的x(n)以补零的方式,使x(n)加长到(n:0~100)时,画出x(n)的频谱X(k) 的幅度;
取x(n)(n:0~100)时,画出x(n)的频谱X(k) 的幅度。
利用FFT进行谱分析
已知:模拟信号
以t=(n=0:N-1)进行采样,求N点DFT的幅值谱。
请分别画出N=45; N=50;N=55;N=60时的幅值曲线。
数字信号处理实验一
1.(1)
L=5;N=20;
n=1:N;
xn=[ones(1,L),zeros(1,N-L)];
Xk=dfs(xn,N);
magXk=abs([Xk(N/2+1:N) Xk(1:N/2+1)]);
k=[-N/2:N/2];
figure(1)
subplot(2,1,1);
stem(n,xn);xlabel('n');ylabel('xtide(n)');
title('DFS of :L=5,N=20');
subplot(2,1,2);
stem(k,magXk);
axis([-N/2,N/2,0,16]);
xlabel('k');ylabel('Xtide(k)');
(2)
L=5;N=40;
n=1:N;
xn=[ones(1,L),zeros(1,N-L)];
Xk=dfs(xn,N);
magXk=abs([Xk(N/2+1:N) Xk(1:N/2+1)]);
k=[-N/2:N/2];
figure(2)
subplot(2,1,1);
stem(n,xn);xlabel('n');ylabel('xtide(n)');
title('DFS of :L=5,N=40');
subplot(2,1,2);
stem(k,magXk);
axis([-N/2,N/2,0,16]);
xlabel('k');ylabel('Xtide(k)');
(3)
L=5;N=60;
n=1:N;
xn=[ones(1,L),zeros(1,N-L)];
Xk=dfs(xn,N);
magXk=abs([Xk(N/2+1:N) Xk(1:N/2+1)]);
k=[-N/2:N/2];
figure(3)
subplot(2,1,1);
stem(n,xn);xlabel('n');ylabel('xtide(n)');
title('DFS of :L=5,N=60');
subplot(2,1,2);
stem(k,magXk);
axis([-N/2,N/2,0,16]);
xlabel('k');ylabel('Xtide(k)');
(4)
L=7;N=60;
n=1:N;
xn=[ones(1,L),zeros(1,N-L)];
Xk=dfs(xn,N);
magXk=abs([Xk(N/2+1:N) Xk(1:N/2+1)]);
k=[-N/2:N/2];
figure(4)
subplot(2,1,1);
stem(n,xn);xlabel('n');ylabel('xtide(n)');
title('DFS of :L=7,N=60');
subplot(2,1,2);
stem(k,magXk);
axis([-N/2,N/2,0,16]);
xlabel('k');ylabel('Xtide(k)');
2. (1)
M=10;N=10;n=1:M;
xn=cos(*pi*n)+cos(*pi*n);
n1=[0:1:N-1];y1=[xn(1:1:M),zeros(1,N-M)];
figu