MATLAB实验报告

数字信号处理及MATLAB

实验报告

班级：

学号：

姓名：

4.7.2 例4,2

设x(n)是由两个正弦信号及白噪声的叠加，试用FFT文件对其作频谱分析。程序清单

%产生两个正弦加白噪声

N=256;

f1=.1;f2=.2;fs=1;

a1=5;a2=3;

w=2*pi/fs;

x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N);

%应用FFT求频谱

subplot(2,2,1);

plot(x(1:N/4));

title('原始信号');

f=-0.5:1/N:0.5-1/N;

x=fft(x);

y=ifft(x);

subplot(2,2,2);

plot(f,fftshift(abs(x)));

title('频域信号');

subplot(2,2,3);

plot(real(x(1:N/4)));

title('时域信号');

例4.3

设x（n）为长度N=6的矩形序列，用MATLAB程序分析FFT取不同长度时x（n）频谱的变化。N=8,32,64,时x（n）的FFT MATLAB实现程序如下。

x=[1,1,1,1,1,1];

N=8;

y1=fft(x,N);

n=0:N-1;

subplot(3,1,1);stem(n,abs(y1),'.k');axis([0,9,0,6]);

N=32;

y2=fft(x,N);

n=0:N-1;

subplot(3,1,2);stem(n,abs(y2),'.k');axis([0,40,0,6]);

N=64;

y3=fft(x,N);

subplot(3,1,3);stem(n,abs(y3),'.k');axis([0,80,0,6]);