%  生成三個正弦函數并計算其幅度譜
N=120;  %  采樣點數
fs=400; %  采樣頻率
n=1:N;  %  離散時間
%  三個正弦函數和
x1=sin(2*pi*50*n/fs);x2=sin(2*pi*100*n/fs);x3=sin(2*pi*150*n/fs);x=x1+x2+x3;
%  時域波形顯示
figure(1);plot(n/fs,x);xlabel('離散時間 (s)');ylabel('振幅');
%  頻域幅度譜顯示
figure(2);stem(linspace(-fs/2,fs/2,N),abs(fftshift(fft(x))));xlabel('離散頻率 (Hz)');ylabel('幅度譜');