用matlab画出连续周期信号的傅里叶级数图像，应该怎么写程序？

下面是一个示例程序，可以画出连续周期信号的傅里叶级数图像：

% 定义连续周期信号的函数
function y = f(x)
    y = sin(x) + 0.5*cos(2*x) + 0.25*sin(3*x);
end

% 定义傅里叶级数的函数
function y = fourier_series(x, N)
    a0 = (1/pi)*quad(@f, 0, 2*pi); % 计算a0
    a = zeros(1, N); % 初始化a
    b = zeros(1, N); % 初始化b
    for n = 1:N
        a(n) = (1/pi)*quad(@(t) f(t).*cos(n*t), 0, 2*pi); % 计算a
        b(n) = (1/pi)*quad(@(t) f(t).*sin(n*t), 0, 2*pi); % 计算b
    end
    y = a0/2 + sum(a.*cos(x*(1:N)) + b.*sin(x*(1:N))); % 计算傅里叶级数
end

% 画出傅里叶级数图像
N = 10; % 傅里叶级数的项数
x = linspace(0, 2*pi, 1000); % x轴范围
y = fourier_series(x, N); % 计算y值
plot(x, y); % 画图
xlabel('x'); ylabel('y'); title('Fourier Series'); grid on;

在上述程序中，`f(x)`是连续周期信号的函数，`fourier_series(x, N)`是计算傅里叶级数的函数。`N`是傅里叶级数的项数，`x`是x轴的范围，`y`是计算出来的傅里叶级数的值。最后用`plot`函数画出图像即可。