%2017/8/23日 星期三
%工程北629
%功能:傅里叶分析——方波的合成——吉伯斯现象
%思想:由于f(x)是奇函数,所以系数a(k)=0,只有b(k)起作用
clear all
close all
%输入信号
N = 200; %数据点数
dt = 4/N; % 采样间隔0.02s N*dt = 200*0.02 = 4秒
for n = 1:N %得到方波 n = 1、2...198、199、200 n*dt = 0.02、0.04...3.96、3.98、4.00
if(n*dt >= 2)
x(n) = 0.8;
else
x(n) = -0.8;
end
end
%信号处理
a = zeros(1,N);
b = zeros(1,N);
nn = floor(N/2) + 1; % nn = 101
for k = 0:nn-1 % k = 0、1、2、3、4......98、99、100
a(k+1) = 0; % a(1)、a(2)、a(3)、a(4)......a(99)、a(100)、a(101)
b(k+1) = 0; % b(1)、b(2)、b(3)、b(4)......b(99)、b(100)、b(101)
for ii = 0:N-1 %0、1、2、3、4......198、199、200
a(k+1) = a(k+1) + 2/N*x(ii+1)*cos(2*pi*k*ii/N); %求解傅里叶系数ak
b(k+1) = b(k+1) + 2/N*x(ii+1)*sin(2*pi*k*ii/N); %求解傅里叶系数bk
end
c(k+1) = sqrt(
例3-2 用振幅为0.8的方波进行傅里叶分析,并用傅里叶分析得到的系数求解当K为不同值时的合成图。
最新推荐文章于 2021-07-26 23:28:54 发布