信號分析試驗

本文格式為Word版，下載可任意編輯——信號分析試驗>x1=ones(1,5);>>x2=0;>>fori=1:5

x2=x2+x1(i)*exp(j*(i-3)*a);end

>>plot(a,abs(x2));>>plot(a,angle(x2));>>

裝訂線

（2）程序：

a=linspace(0,2*pi,60);

>>n=0:10;>>x1=2.^n；>>x2=0;>>fori=1:5

x2=x2+x1(i)*exp(j*(i-1)*a);end

>>plot(a,abs(x2));>>plot(a,angle(x2));

2.

DFT：>>n=8;a=0:1:7;b=a;

x1=[8,7,9,5,1,7,9,5];x2=zeros(1,n);m=exp(-2*pi*j/n);forb=1:nfori=1:n

x2(b)=x2(b)+x1(i)*(m.^((i-1)*(b-1)));endend

stem(a,abs(x2))>>stem(a,angle(x2))

>>

FFT：

n=0:7;

xn=[8,7,9,5,1,7,9,5];Xk=fft(xn,8);stem(n,abs(Xk));stem(n,angle(Xk));結(jié)果與DFT一致如上圖3.FFT：

t=linspace(0,10,100);n=0:99;

xn=3*cos(8*pi*t);Xk=fft(xn,100);

subplot(1,2,1),stem(n,abs(Xk));subplot(1,2,2),stem(n,angle(Xk));

DFT：

a=linspace(0,10,100);N=100;

X1=3*cos(8*pi*a);n=0:1:(N-1);k=n;

X2=zeros(1,N);w=exp(-2*pi*j/N);fork=1:Nfori=1:N

X2(k)=X2(k)+x1(i)*(w.^((i-1)*(k-1)));endend

subplot(1,2,1),stem(n,abs(X2))

subplot(1,2,2),stem(n,angle(X2))結(jié)果同F(xiàn)FT。與X（w）基本一致。4.FFT：

a=linspace(0,10,100);

n=0:99;

y=randn(1,100);x1=3*cos(8*pi*a)+y;X2=fft(x1,100);

subplot(1,2,1),stem(n,abs(X2));subplot(1,2,2),stem(n,angle(X2));

>>DFT

基本與上題一致，在xn處加上randn(1,100)即可。結(jié)果與FFT一致。由結(jié)果可知，噪聲對于信號頻譜有著很大的影響。5.（1）程序N=6;n=0:5;

xn=[4,3,2,1,0,0];Xk=fft(xn,6);fork=1:6

Yk(k)=Xk(k)*exp(j*4*2*pi*(k-1)/N);end

yn=ifft(Yk,6)結(jié)果：yn=

Columns1through3

-0.0000+0.0000i-0.0000+0.0000i4.0000-0.0000i

Columns4through6