fseries.m

通过实现文件ft1调用函数fseries

%1．fseries.m
% [A,B,y]=fseries(f,x,n,a,b)
% f为待展开函数，
% x为自变量
% n为展开项数
% [a,b]为展开区间，省略则为 
% A,B为傅里叶系数
% y为返回的展开式
% 程序如下：
function [A,B,F]=fseries(f,x,n,a,b)
if nargin==3, a=-pi;  b=pi; end  %若输入为3个参数表示f(x)在[-pi,pi]内展开
L=(b-a)/2; 
A=int(f,x,-L,L)/L; %f(x)在[-L,L]内积分
B=[];  
F=A/2;
for i=1:n
   an=int(f*cos(i*pi*x/L),x,-L,L)/L;   %积分求系数an
   bn=int(f*sin(i*pi*x/L),x,-L,L)/L;   %积分求系数bn
   A=[A,an];   %记录所有an
   B=[B,bn];  %记录所有bn
   F=F+an*cos(i*pi*x/L)+bn*sin(i*pi*x/L);  %得到n阶傅里叶展开式
end