ljh2011 发表于 2011-8-30 15:42

带有二个变参数的二重积分问题,以两变参数数值积分后绘图!

二元函数f(x,y) =sin(k.*x).*cos(m.*y).*x.^2,y从1积分到5,x从1积分到2;
代码:

kk=linspace(0,5);
mm=linspace(0,5);
y=zeros(size(kk));
ff=@(k,m) ['sin(',num2str(k),'*x).*cos(',num2str(m),'*y).*x.^2'];
f=@(k,m) quadl(@(x) arrayfun(@(xx) quadl(ff(k,m),1,5),x),1,2);
for ii=1:length(kk)
for jj=1:length(mm)
y(ii,jj)=f(kk(ii),mm(jj));
end
end
plot3(kk,mm,y)

不知错在哪,请各位指点!!

qibbxxt 发表于 2011-8-31 09:50

ezsurf(@(k,m)quad2d(@(x,y)sin(k.*x).*cos(m.*y).*x.^2,1,2,1,5))
页: [1]
查看完整版本: 带有二个变参数的二重积分问题,以两变参数数值积分后绘图!