灵鹊做题网matlab编程求解双重积分
来源：学生作业帮 编辑：灵鹊做题网作业帮 分类：综合作业 时间：2024/04/28 18:52:27
matlab编程求解双重积分
我想求解f=a*b*cos(x-y)/(a^2+b^2+d^2-2*a*b*cos(x-y)+2*l*(a*cos(x)-b*cos(y)))^(1/2);对x、y的双重积分,编程如下：
编写m文件：
function f=fun(x,y)
f=a*b*cos(x-y)/(a^2+b^2+d^2-2*a*b*cos(x-y)+2*l*(a*cos(x)-b*cos(y)))^(1/2);
命令窗口中：
a=50;
b=10;
d=2;
l=0;
Q=dblquad('fun',0,2*pi,0,2*pi)；
结果有错误,本人菜鸟求指点.
我想求解f=a*b*cos(x-y)/(a^2+b^2+d^2-2*a*b*cos(x-y)+2*l*(a*cos(x)-b*cos(y)))^(1/2);对x、y的双重积分,编程如下：
编写m文件：
function f=fun(x,y)
f=a*b*cos(x-y)/(a^2+b^2+d^2-2*a*b*cos(x-y)+2*l*(a*cos(x)-b*cos(y)))^(1/2);
命令窗口中：
a=50;
b=10;
d=2;
l=0;
Q=dblquad('fun',0,2*pi,0,2*pi)；
结果有错误,本人菜鸟求指点.
a=50;b=10;d=2;l=0;
f=inline('a*b*cos(x-y)./(a^2+b^2+d^2-2*a*b*cos(x-y)+2*l*(a*cos(x)-b*cos(y))).^(1/2)','x','y');
Q=dblquad(f,0,2*pi,0,2*pi)
Q=39.9838
再问： 数值积分，怎么解？
再答： function f=fun(x,y,a,b,d,l) f=a*b*cos(x-y)./(a^2+b^2+d^2-2*a*b*cos(x-y)+2*l*(a*cos(x)-b*cos(y))).^(1/2); end clear all clc a=50;b=10;d=2;l=0; f=@(x,y)fun(x,y,a,b,d,l); Q=dblquad(f,0,2*pi,0,2*pi)
f=inline('a*b*cos(x-y)./(a^2+b^2+d^2-2*a*b*cos(x-y)+2*l*(a*cos(x)-b*cos(y))).^(1/2)','x','y');
Q=dblquad(f,0,2*pi,0,2*pi)
Q=39.9838
再问： 数值积分，怎么解？
再答： function f=fun(x,y,a,b,d,l) f=a*b*cos(x-y)./(a^2+b^2+d^2-2*a*b*cos(x-y)+2*l*(a*cos(x)-b*cos(y))).^(1/2); end clear all clc a=50;b=10;d=2;l=0; f=@(x,y)fun(x,y,a,b,d,l); Q=dblquad(f,0,2*pi,0,2*pi)