matlab对分段函数积分

2025-03-01 04:10:32
推荐回答(1个)
回答1:

function F=FX(L,h,n)
F=zeros(n+3);
f=@(x)1/6*((x+2).^3).*(x>=-2 & x<-1)+1/6*((x+2).^3-4*(x+1).^3).*(x>=-1 & x<0)+...
1/6*((2-x).^3-4*(1-x).^3).*(x>=0 & x<1)+1/6*((2-x).^3).*(x>=1 & x<=2);
g=@(x)1/6*((x+2).^3).*(x>=-2 & x<-1)+1/6*((x+2).^3-4*(x+1).^3).*(x>=-1 & x<0)+...
1/6*((2-x).^3-4*(1-x).^3).*(x>=0 & x<1)+1/6*((2-x).^3).*(x>=1 & x<=2);
for i=1:1:n+3
for k=1:1:n+3
fg=@(x)f(x/h-(i-2)).*g(x/h-(k-2));
F(i,k)=quad(fg,0,L);
end
end