sunyok 发表于 2008-3-12 10:49

内联函数中怎样处理字母系数

比如：对附件1中双重积分的处理——

function SS=double_int(fun,innlow,innhi,outlow,outhi)
y1=outlow;
y2=outhi;
x1=innlow;
x2=innhi;
f_p=fun;
SS=quad(@G_yi,y1,y2,[],[],x1,x2,f_p);


function f=G_yi(y,x1,x2,f_p)
y=y(:);
n=length(y);
if isnumeric(x1)==1;
    xx1=x1*ones(size(y));
else xx1=feval(x1,y);
end
if isnumeric(x2)==1;
    xx2=x2*ones(size(y));
else xx2=feval(x2,y);
end
for i=1:n
    f(i)=quad(f_p,xx1(i),xx2(i),[],[],y(i));
end
f=f(:);


function f=x_hi(y)
f=y;


ff=inline('x.^2+y.^2','x','y');      
SS=double_int(ff,3,@x_hi,3,2)

那对附件2中积分怎样处理呢？？？？？？

PS:我也不想重复发帖，只是刚才发第一次时有错误，抱歉。

[ 本帖最后由 sunyok 于 2008-3-12 10:59 编辑 ]

eight 发表于 2008-3-12 10:52

原帖由 sunyok 于 2008-3-12 10:49 发表 http://www.chinavib.com/forum/images/common/back.gif
比如：对附件1中双重积分的处理——

function SS=double_int(fun,innlow,innhi,outlow,outhi)
y1=outlow;
y2=outhi;
x1=innlow;
x2=innhi;
f_p=fun;
SS=quad(@G_yi,y1,y2,[],[],x1,x2,f_p);
...
发帖前，先看看置顶区所有帖子，特别是会员守则（对比一下你的标题，还有重复发帖）。如果不是见你的帖子内容较多、重写一次比较麻烦的话，我肯定删除。

[ 本帖最后由 eight 于 2008-3-12 10:55 编辑 ]

sigma665 发表于 2008-3-12 10:53

回复 楼主 的帖子

不清楚，不过可以加个循环啊
页: [1]
查看完整版本: 内联函数中怎样处理字母系数