4.1 傅里叶级数的MATLAB计算
设周期信号x(t)的基本周期为T1，且满足狄里克利条件，则其傅里叶级数的系数可由式2.4计算得到。式2.4重写如下：
2/2/1110)(1TTtjkkdtetxTa
基本频率为：
102T
对周期信号进行分析时，我们往往只需对其在一个周期内进行分析即可，通常选择主周期（Principle period）。假定x1(t)是x(t)中的主周期，则
2/2/11110)(1TTtjkkdtetxTa
计算机不能计算无穷多个系数，所以我们假设需要计算的谐波次数为N，则总的系数个数为2N+1个。在确定了时间范围和时间变化的步长即T1和dt之后，对某一个系数，上述系数的积分公式可以近似为：
ntjknTTtjkkTdtetxdtetxTa12/2/11/)()(10110
121/],,[)](),(),([02010TdteeetxtxtxMtjktjktjkM
对于全部需要的2N+1个系数，上面的计算可以按照矩阵运算实现。MATLAB实现系数计算的程序如下：
dt = 0.01;
T = 2; t = -T/2:dt:T/2; w0 = 2*pi/T;
x1 = input(‘Type in the periodic signal x(t) over one period x1(t)=’);
N = input(‘Type in the number N=’);
k = -N:N; L = 2*N+1;
ak = x1*exp(-j*k*w0*t’)*dt/T;
需要强调的是，时间变量的变化步长dt的大小对傅里叶级数系数的计算精度的影响非常大，dt越小，精度越高，但是，计算机计算所花的时间越长。
例题2-1：给定一个周期为T1 = 2s的连续时间周期方波信号，如图所示，其一个周期内的数学表达式为：
短时傅里叶变换（Short-Time Fourier Transform, STFT）是一种在信号处理领域广泛应用的技术，特别在语音处理、音频处理、地震学、通信系统等领域有着重要的作用。STFT可以将信号从时间域转换到频率域，从而能够以时间和频率的双重视角来分析信号的特性。在MATLAB中，我们可以使用内置的STFT函数来实现信号的时频分析，以及一些功能强大的工具箱来进行更深入的信号处理和分析。
1. STFT的原理
STFT可以看作是对信号在一段时间内进行傅里叶变换的过程。在传统的傅里叶变换中，我们是对整段信号进行傅里叶变换，从而得到信号在整个时间范围内的频率特性。然而，STFT允许我们对信号进行局部的傅里叶变换，这样就可以观察到信号在不同时间段内的频率变化，从而更加全面地理解信号的特性。
2. MATLAB中的STFT函数
在MATLAB中，我们可以通过调用stft函数来实现对信号的短时傅里叶变换。该函数可以指定窗口长度、重叠长度以及窗口函数等参数，从而灵活地调整STFT的分辨率和精度。通过这些参数的设置，我们可以得到不同粒度和分辨率的时频分析结果，从而更好地理解信号的时频特性。
3. STFT的应用
在实际的工程和科研中，STFT有着广泛的应用。在语音信号处理中，可以利用STFT来进行语音的特征提取和分析，从而实现语音识别、语音合成等功能。在音频处理领域，STFT可以用于音乐信号的谱分析和乐器识别。STFT还可以应用于地震学领域的地震信号处理，通信系统中的信号调制解调等多个领域。
4. MATLAB工具箱的应用
除了内置的stft函数外，MATLAB还提供了丰富的工具箱来支持STFT相关的功能。Signal Processing Toolbox提供了丰富的时频分析工具函数，可以对信号进行更加深入的分析和处理。另外，Wavelet Toolbox也可以用于时频分析，提供了小波变换等功能，能够更好地适应不同频率分量的信号分析。
——————————————————————————————————————————————— MATLAB傅里叶级数
河南工程学院
理
学 院
MATLAB
姓名：杨标
班级：信息1241
学号：201211002126
（傅立叶级数的模拟）
实验目的： （1）用MATLAB模拟傅立叶级数
F(x)=a0/2+?(ancosnx?bnsinnx)
n?1?
求a0,an,bn;
（2）对矩形函数，三角函数求傅立叶级数；（作图）
（3）绘制函数及其傅立叶逼近
实验分析（1）：
将F(x)=a0/2+?(ancosnx?bnsinnx) 展开为傅立叶级数，既是求
n?1?
ai,bi 根据三角函数的正交性，可得：
a0=1
? ——————————————————————————————————————————————— an=1
???f(x)dx; ????f(x)cosnxdx; ??
bn=1
???f(x)sinnxdx应用积分求解公式int ??
源程序：
clear
syms x n
“函数方程”） f=（
a0=int(f,x,-pi,pi)/pi
an=int(f*cos(n*x),x,-pi,pi)/pi bn=int(f*sin(n*x),x,-pi,pi)/pi 例：求函数f(x)=x^3+5 clear
syms x n
f=x^3+5
a0=int(f,x,-pi,pi)/pi
an=int(f*cos(n*x),x,-pi,pi)/pi bn=int(f*sin(n*x),x,-pi,pi)/pi运行后得：
线性部分：
两边同时对
非线性部分 :U :z 二
inU 分步傅里叶算法的MATLAB程序实现举例
模型:
2
.仁 U
I 一
；z 2 ：：x
2
.2 d n n = U (x, z j +d (z )—2 ex
其中:
d z i=exp z 或 d z = exp -、z
dU = I C2U
；z 2 ；x2
x变量作傅里叶变换
两边积分
x h 2 1 ~ x h 2 I 2
x =dU = % - r dz U 2
最后有
U (x + h/2)=U (x)exp I；(险 f *
再对x变量作作傅里叶逆变换
U (x + h/2)=FFT-1 ]FFT[U(x)]expgz )
两边积分
x」h 1 x "h x jdU = x Indz 折射率部分 |2
n = U (x, z) + d(z)马 ex
n =FFT u2 d z * i
n=FFT-1 FFT _U
2 1 2d z lnU2L±.inh U x
最后有
U x h = U x exp inh
两边同时对x变量作傅里叶变换
FFT U 2
1 2d z
再对x变量作作傅里叶逆变换
MATLAB程序实现：
clear all
delta=1;
x0=1;
%% ------------------
n=2048;
hx=0.06;
x=(- n/2:n /2-1)*hx;
hw=2*pi/( n*hx);
w=fftshift((- n/2: n/2-1)*hw);
%% ------------------
q=exp(-1*(x-x0).A2/2)+ exp(-1*(x+x0).A2/2); % q=sech(x);
u1(:,1)=(abs(q)A2)：
% ------------------
L=500;
nm=L*100;
h=L/ nm;
% ------------------
