1 / 8
文档名称:

Matlab 离散傅里叶变换 实验报告.doc

格式:doc   大小:231KB   页数:8页
下载后只包含 1 个 DOC 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

Matlab 离散傅里叶变换 实验报告.doc

上传人:mh900965 2016/9/2 文件大小:231 KB

下载得到文件列表

Matlab 离散傅里叶变换 实验报告.doc

文档介绍

文档介绍:陕西科技大学实验报告第页班级信工 082 学号 16 姓名刘刚实验组别实验日期室温报告日期成绩报告内容: ( 目的和要求,原理,步骤,数据,计算,小结等) 实验三离散傅立叶变换( DFT ) 1. 离散傅立叶级数给定有限长序列[123 4],延拓为周期 N=6 的周期序列,并求其 DFS 。代码: N1=6;x1=[1 234]; N2=length(x2); n1=0:5*N2-1; x2=[x1,zeros(1,(6-length(x1)))];k=0:5*N2-1;x3=x2(mod(n1,N2)+1) Xk=x3*exp(-j*2*pi/N1).^(n1'*k); subplot(321),stem(x1,'.');title(' 原序列') subplot(322),stem(x3,'.');title(' 原序列周期延拓') subplot(312),stem(Xk,'.');title('DFS') subplot(325),stem(abs(Xk),'.');title('DFS 模值') subplot(326),stem(angle(Xk),'.');title('DFS 相位') 结果: 附页第页 2. 求以下序列的 N=16,32,64 点的快速傅立叶变换 x(n)=exp(-i*pi/8*n);x2(n)=cos(pi/8*n);x3(n)=sin(pi/8*n) (1)代码: Fs=20; N1=16;n1=0:N1-1;xn1=exp(-i*pi/8*n1);D1=2*pi*Fs/N1; N2=32;n2=0:N2-1;xn2=exp(-i*pi/8*n2);D2=2*pi*Fs/N2; N3=64;n3=0:N3-1;xn3=exp(-i*pi/8*n3);D3=2*pi*Fs/N3; k1=floor(-(N1-1)/2:(N1-1)/2);k2=floor(-(N2-1)/2:(N2-1)/2); k3=floor(-(N3-1)/2:(N3-1)/2); X1=fft(xn1,N1);X2=fft(xn2,N2);X3=fft(xn3,N3); subplot(3,4,1);stem(n1,xn1);title(' 原信号');ylabel('C=16') subplot(3,4,2);stem(X1,'.');title('FFT 结果 X'); subplot(3,4,3);stem(k1*D1,abs(X1));title('X 的模值'); subplot(3,4,4);plot(k1*D1,angle(X1));title('X 的相位'); subplot(3,4,5);stem(n2,xn2);ylabel('C=32') subplot(3,4,6);stem(X2,'.'); subplot(3,4,7);stem(k2*D2,abs(X2)); subplot(3,4,8);plot(k2*D2,angle(X2)); subplot(3,4,9);stem(n3,xn3);ylabel('C=64') subplot(3,4,10);stem(X3,'.'); subplot(3,4,11);stem(k3*D3,abs(X3)); subplot(3,4,12);plot(k3*D3,an