matlab傅里叶变换程序

clc
clear
fs=10000; % 采样频率>2倍的信号频率
N=256000; % 采样点数目（=2的幂次）
n=0:N-1; % 构造采样点序列
t=n/fs; % 得到采样时间序列，t=nT=n/fs
A=1sin(2pi300t);%载波
B=0.1sin(2pi15t);%调制信号
x=A+B; % 产生时域信号的样本值，向量
% x=A.*(1+B); % 产生时域信号的样本值，向量
Y=fft(x,N); % N点的DFT计算
mag=abs(Y); % FFT的振幅
phase=unwrap(angle(Y)); % FFT的相位
figure
plot(t,x)
xlim([0,0.1])
figure
% 1. 以下绘制物理频谱图（即正频部分）
fn=(0:N/2)fs/N; % 频率轴上的离散频率点，起始于0频(对应直流成分)，终
%于Nyquist频率fs/2，共N/2+1个频率点
% subplot(2,2,1)% 将图形窗口分割为2×2的子窗口，并指定第1个子窗口为绘图区
plot(fn,mag(1:N/2+1),‘b’) % 取出前N/2+1个振幅作图，即正频率分量
xlim([0,400])
figure
plot(fn,20log10(mag(1:N/2+1)),‘b’)%对数谱，以显示幅值较小的频率
xlabel(‘频率/Hz’);ylabel(‘振 幅’);
title(‘图1： 物理(正频)幅频图’);grid on % 加网格线
xlim([0,400])
% % 2. 以下绘制全频率的幅频图
% fn1=(0:N-1)*fs/N;
% subplot(2,2,2) %指定第2个子窗口为绘图区
% plot(fn1,mag);
% xlabel(‘频率/Hz’);ylabel(‘振 幅’);
% title(‘图2： 全频率的幅频图’); grid on
% %3. 以下绘制正频部分的相频图
% subplot(2,2,3) %指定第3个子窗口为绘图区
% plot(fn,phase(1:N/2+1));
% xlabel(‘频率/Hz’);ylabel(‘相 位’);title(‘图3： 相 频 图’);grid
% %4. 以下移动零频点
% Y1=fftshift(Y); % fftshift移动频率零点，并将Y的左右两部分交换
% mag1=abs(Y1); % 重新计算振幅
% fn2=fn1-fs/2; % 零点移动到fs/2处，故需重新标记频率轴
% subplot(2,2,4); %指定第4个子窗口为绘图区，最终4幅图绘制在一张图上了
% plot(fn2,mag1);
% xlabel(‘频率/Hz’);ylabel(‘振 幅’);
% title(‘图4： fftshift后的幅频图’);grid