%Samp3_15
clf;
fs=1000;N=1024;   %采樣頻率和數據點數
n=0:N-1;t=n/fs;   %時間序列
f1=50;f2=120;   %信號的頻率成分
x=sin(2*pi*f1*t)+2*sin(2*pi*f2*t);   %無噪聲信號
x=x+2*randn(1,length(t));  %含有噪聲的信號
y=fft(x,N);            %對含有噪聲信號進行Fourier變換
mag=abs(y);     %得到振幅值
f=n*fs/N;   %真實頻率序列
subplot(2,1,1),plot(f(1:N/2),mag(1:N/2)*2/N); %繪出Nyquist頻率之前的振幅
xlabel('頻率/Hz');
ylabel('振幅');title('N=1024');grid on;
x=sin(2*pi*f1*t)+2*sin(2*pi*f2*t);  %無噪聲信號
y=fft(x,N);     %對無噪聲信號進行快速Fourier變換
mag=abs(y);    %求取振幅
f=n*fs/N;   %真實頻率序列
subplot(2,1,2),plot(f(1:N/2),mag(1:N/2)*2/N); %繪出Nyquist頻率之前的振幅
xlabel('頻率/Hz');
ylabel('振幅');title('N=1024（無噪聲）');grid on;