ft1.m

来自「通过实现文件ft1调用函数fseries」· M 代码 · 共 13 行

% 2．实现函数fourier.m
syms x;
f=abs(x)/x; %给出待展开函数
xx=[-pi:pi/200:pi];  %给出所有要作图的点
xx=xx(xx~=0);  %去掉为0的点
yy=subs(f,x,xx);  %计算f(x)的值
plot(xx,yy); %作f(x)图
hold on  

n=30;   %给出展开的阶数
    [a,b,f1]=fseries(f,x,n);% 求n阶展开 
    y1=subs(f1,x,xx);  %计算n阶展开在各点的函数值
    plot(xx,y1);  %作图