绘制符号傅里叶级数

Question

我已经为傅里叶级数编写了代码。这是我目前所拥有的：

function FS = FourierSeries(f,degree)

cosCoefficients = zeros(1,degree);
sinCoefficients = zeros(1,degree);

syms x;

a0 = double((1/pi)*int(f,-pi,pi));


for n = 1:degree
    cosCoefficients = cosCoefficients + (1/pi)*int(f*cos(n*x),-pi,pi);
    sinCoefficients = sinCoefficients + (1/pi)*int(f*sin(n*x),-pi,pi);
end

for n = 1:degree
    FS = 0.5*a0 + cosCoefficients.*cos(n*x) + sinCoefficients.*sin(n*x);
end

然后，我还创建了以下函数文件：

function y = func1(x)
syms x
y = (x^2);
end

问题是当我尝试绘制 func1 和 FS1 = FourierSeries(func1, 4) 时，我不断收到一条错误消息

Data must be numeric, datetime, duration or an array convertible to double.

如何绘制这个傅里叶级数？

Answer 1

whos FS1
  Name      Size            Bytes  Class    Attributes

  FS1       1x4                 8  sym

这告诉您 FS1 是一个 符号函数 ，在绘制它之前需要对其求值：

FS1 = FourierSeries(func1, 4);
xIDX = -10:0.1:10;
array = zeros(size(xIDX));

for ii = 1:numel(xIDX)
    x = xIDX(ii);
    array(ii) = sum(double(subs(FS1)));
end

figure
plot(array)

subs 将您的符号表达式转换为字符串，前提是符号变量 x 在这种情况下存在于工作区中。 double 然后将字符串转换为实际数字，由于 FS 中有 4 个项，我们需要对它们求和，毕竟这就是傅里叶级数。

更短，不用数值计算，而是利用MATLAB内置的符号绘图函数fplot

figure
hold on
fplot(func1,'b')
fplot(sum(FS1),'r') % sum over the Fourier components
legend ('func1', 'FS1')

绘制符号傅里叶级数

Plotting a symbolic Fourier series

matlab

plot

symbolic-math

continuous-fourier