|
楼主 |
发表于 2010-6-22 06:55
|
显示全部楼层
请教matlab多重积分问题
求解上面积分问题, 参考rocwoods给的方法
clear all; theta=1; T=10; H=0.78;
ff1=@(s2,u2)@( s1,u1)exp(-theta*(s2-u2)-theta*(s1-u1)).*(s2-s1).^(2*H-2).*(u2-u1).^(2*H-2);
ff2=@(s2,u2)@( s1,u1)exp(-theta*(s2-u2)-theta*(s1-u1)).*(s2-s1).^(2*H-2).*(u1-u2).^(2*H-2);
ff3=@(s2,u2)@( s1,u1)exp(-theta*(s2-u2)-theta*(s1-u1)).*(s1-s2).^(2*H-2).*(u2-u1).^(2*H-2);
ff4=@(s2,u2)@( s1,u1)exp(-theta*(s2-u2)-theta*(s1-u1)).*(s1-s2).^(2*H-2).*(u1-u2).^(2*H-2);
ff5=@(s2,u2)@( s1,u1)exp(-theta*(s2-u2)-theta*(u1-s1)).*(s2-s1).^(2*H-2).*(u2-u1).^(2*H-2);
ff6=@(s2,u2)@( s1,u1)exp(-theta*(s2-u2)-theta*(u1-s1)).*(s2-s1).^(2*H-2).*(u1-u2).^(2*H-2);
ff7=@(s2,u2)@( s1,u1)exp(-theta*(s2-u2)-theta*(u1-s1)).*(s2-s1).^(2*H-2).*(u2-u1).^(2*H-2);
ff8=@(s2,u2)@( s1,u1)exp(-theta*(s2-u2)-theta*(u1-s1)).*(s1-s2).^(2*H-2).*(u2-u1).^(2*H-2);
[email=ff9=@(s2,u2]ff9=@(s2,u2[/email])@( s1,u1)exp(-theta*(u2-s2)-theta*(s1-u1)).*(s2-s1).^(2*H-2).*(u2-u1).^(2*H-2);
ff10=@(s2,u2)@( s1,u1)exp(-theta*(u2-s2)-theta*(s1-u1)).*(s2-s1).^(2*H-2).*(u1-u2).^(2*H-2);
ff11=@(s2,u2)@( s1,u1)exp(-theta*(u2-s2)-theta*(s1-u1)).*(s1-s2).^(2*H-2).*(u2-u1).^(2*H-2);
ff12=@(s2,u2)@( s1,u1)exp(-theta*(u2-s2)-theta*(s1-u1)).*(s1-s2).^(2*H-2).*(u1-u2).^(2*H-2);
ff13=@(s2,u2)@( s1,u1)exp(-theta*(u2-s2)-theta*(u1-s1)).*(s2-s1).^(2*H-2).*(u2-u1).^(2*H-2);
ff14=@(s2,u2)@( s1,u1)exp(-theta*(u2-s2)-theta*(u1-s1)).*(s2-s1).^(2*H-2).*(u1-u2).^(2*H-2);
ff15=@(s2,u2)@( s1,u1)exp(-theta*(u2-s2)-theta*(u1-s1)).*(s2-s1).^(2*H-2).*(u2-u1).^(2*H-2);
ff16=@(s2,u2)@( s1,u1)exp(-theta*(u2-s2)-theta*(u1-s1)).*(s1-s2).^(2*H-2).*(u2-u1).^(2*H-2);
现在请问对应的 Int1……Int16 应该怎么写呢? 谢谢
[ 本帖最后由 ChaChing 于 2010-6-22 23:43 编辑 ] |
|