matlab求DFT 画出幅度谱,matlab作业1

Maltab数字信号处理上机作业

实验任务

XXXXXX

(1)阅读例子程序，观察输出波形，理解每条语句的含义。

(2)已知有限长序列x(n)=[7,6,5,4,3,2],求DFT和IDFT，要求：画出序列傅立叶变换对应的幅度谱和相位谱；画出原信号与傅立叶逆变换IDFT[X(k)]的图形进行比较。

(3)已知周期序列的主值x(n)=[7,6,5,4,3,2],求x(n)周期重复次数为3次时的DFS和IDFS。要求：画出原信号序列的主值和周期序列的图形；画出离散傅立叶变换对应的幅度谱和相位谱。

(4)求x(n)=[7,6,5,4,3,2], 0=

%%

%题(2)

close all;clear all;clc;

xn=[7,6,5,4,3,2];

N=length(xn);

n=0:N-1;k=0:N-1;

Xk=xn*exp(-j*2*pi/N).^(n'*k);

x=(Xk*exp(j*2*pi/N).^(n'*k))/N;

subplot(2,2,1),stem(n,xn,'k');

title('x(n)');axis([-1,N,0.8*min(xn),1.1*max(xn)]);grid; subplot(2,2,2),stem(n,x,'k');

title('IDFT|X(k)|');

axis([-1,N,0.8*min(x),1.1*max(x)]);grid;

subplot(2,2,3),stem(k,abs(Xk),'k');

title('|X(k)|');

axis([-1,N,min(abs(Xk)),1.1*max(abs(Xk))]);

subplot(2,2,4),stem(k,angle(Xk),'k');

title('arg|X(k)|');

axis([-1,N,min(angle(Xk)),1.1*max(angle(Xk))]);

%%

%题(3)

close all;clear all;clc;

xn=[7,6,5,4,3,2];

N=length(xn);

n=0:3*N-1;k=0:3*N-1;

xn1=xn(mod(n,N)+1);

Xk=xn1*exp(-j*2*pi/N).^(n'*k);

subplot(2,2,1),stem(xn,'k');

title('原主值信号x(n)');