matlab分析傅里叶级数

前提

傅里叶级数求各次谐波的叠加

t1.m

clear all;
clc;
t=[1,3,5,9,19,39,79,159];
for i=t
    func(i)
end


func.m

function fn=func(m)
t=-2:0.0001:2;
y=square(2*pi*t,50);
T=1;
w=(2*pi)/T;
fn=0;
for n=1:2:m
    fn=fn+(4/pi)*((1/n)*sin(n*w*t))
end
figure(m)
plot(t,y,t,fn)
grid;

求傅里叶系数

t2.m


clear all;
clc;
tao=[1 1 2];
T=[10 5 10];
for i=1:length(tao)
    func2(tao(i),T(i))
end


func2.m

function fn=func2(tao,T)
A=1;
w=(-8*pi):(2*pi/T):(8*pi);
fn=A*tao*sinc(w/pi*tao/2);
figure
stem(w,fn)
grid
axis([-25 25 -0.5 2]);