求matlab程序,还有,stem画图和subplot,plot 画图的时候有什么区别1.利用FFT分析有限长序列xn={2,3,3,1,0,5; n=0,1,2,3,4,5} 的频谱,并绘出其幅度谱与相位谱.2.利用FFT分析无限长非周期序列 xn=(1/3)的n次幂的

Question

求matlab程序,还有,stem画图和subplot,plot 画图的时候有什么区别
1.利用FFT分析有限长序列xn={2,3,3,1,0,5; n=0,1,2,3,4,5} 的频谱,并绘出其幅度谱与相位谱.
2.利用FFT分析无限长非周期序列 xn=(1/3)的n次幂的频谱,n>=0开始,（1） 确定序列的长度 及窗函数的类型（序列为无限长时,需要根据能量分布进行截短）.
（2） 确定作FFT的点数 ；根据频域取样定理,为使时域波形不产生混叠,必须 .
（3） 使用fft函数作 点FFT计算 ,并绘出其幅度谱与相位谱.
（4） 进行理论值与计算值比较,讨论信号频谱分析过程中误差原因及改善方法.
.
第二题给出一部分答案也行,程序最好注解一下,分值再加

Answer 1

%问题一
xn=[2 3 3 1 0 5];
Xn=fft(xn);
figure(1);
subplot(211)
plot(abs(Xn));
xlabel('序列k');
ylabel('幅度');
title('序列幅频图');
subplot(212)
plot(angle(xn));
xlabel('序列k');
ylabel('幅度');
title('序列相频图');
%问题二
n=0:1000;%选取前1000个点 长度足够 简单选取规则 未进行能量分布域的分析
x1=(1/3).^n;
z=(hamming(length(n)))';%加窗 选为汉明窗 长度1001
x2=x1.*z;%时域加窗
X2=fft(x2);
figure(2);
subplot(211)
plot(abs(X2));
axis([0 1000 -0.4 0.4]);
xlabel('序列k');
ylabel('幅度');
title('序列幅频图');
subplot(212)
plot(angle(X2));
axis([0 1000 -0.4 0.4]);
xlabel('序列k');
ylabel('幅度');
title('序列相频图');
截为有限长,必然会产生频谱的混叠,因为信源的频率并不是带限的.