菜鸟之MATLAB学习——FM0波形生成及FFT变换

首先声明：我是MATLAB初学者，只做笔记记录。

clc;close all;fdata=160*1000;% 数据速率160kT=1/fdata;%信号周期N_sample=10;%每个周期的采样点数dt=T/N_sample;d0_1=[ones(1,N_sample/2),-1*ones(1,N_sample/2)];d0_2=[-ones(1,N_sample/2),ones(1,N_sample/2)];d1_1=[ones(1,N_sample)];d1_2=-1*ones(1,N_sample);N_data=100;fm0_seq=randi([0,1],1,N_data);% 产生N_data个01交替的序列FM0_WAVE=[];iffm0_seq(1)==0%先判断第一个数据是0还是1；FM0_WAVE=[FM0_WAVE d0_1];elseFM0_WAVE=[FM0_WAVE d1_1];endfori=2:N_dataifFM0_WAVE(length(FM0_WAVE))==-1iffm0_seq(i)==0FM0_WAVE=[FM0_WAVE d0_1];elseFM0_WAVE=[FM0_WAVE d1_1];endelseiffm0_seq(i)==0FM0_WAVE=[FM0_WAVE d0_2];elseFM0_WAVE=[FM0_WAVE d1_2];end end endfs=1/dt;%采样频率N_all=N_data*N_sample;%总的采样点数df=fs/N_all;t=0:dt:(N_all-1)*dt;f=0:df:(N_all-1)*df;sf=fft(FM0_WAVE);sf=fftshift(sf);%fft三件套f=f-N_all*df/2;subplot(2,1,1),plot(t,FM0_WAVE);axis([0(N_all-1)*dt-22]);subplot(2,1,2),plot(f/1000,abs(sf));xlabel("频率/kHz")

波形及FFT变换

可以多运行几次，因为数据是随机的，每次的频谱也不一样。

FFT运算的理解

注意3个要素：采样频率、采样点数、频率分辨率（三者中只有两个变量是独立的）。

MATLAB是如何构建FFT算法的呢？MATLAB在计算时并不关注采样频率，也不关注采样间隔，只关注采样点数，即做多少个点的FFT，而我们想分析信号的频谱，就必须规定采样频率或者频率分辨率（二者取其一），并用plot函数（构成一 一映射关系）将其绘制即可。