作者thereocr (眼睛闭上就可以睡着)
看板MATLAB
标题[发问] 二重积分的运算问题
时间Tue Dec 18 16:58:41 2012
2π 2π a1*a2*cos(α-β)
∫ ∫ _____________________________________dαdβ
0 0 (a1^2+a2^2+d^2 - 2*a1*a2*cos(α-β))
以上是题目!!
我试了好久用Matlab记算这个方程式都没算出来
不知道问题出在哪里...
希望高人可以指点我ㄧ下
--------------以下程式-------
----------------设定function----
function out=integrn(X,Y)
D=0.8;
a1=55;
a2=30;
up=(a1*a2*cos(X-Y));
down=sqrt(a1^2+a2^2+D^2-2*a1*a2*(cos(X-Y)));
out=up/down;
end
------------------------------------
---------运算--------------
Xmin=0;
Xmax=2*pi;
Ymin=0;
Ymax=2*pi;
ans=dbquad(@integrn,Xmin,Xmax,Ymin,Ymax)
--------------------------------------
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 150.69.51.201
※ 编辑: thereocr 来自: 150.69.51.201 (12/18 16:59)
※ 编辑: thereocr 来自: 150.69.51.201 (12/19 09:59)
1F：推 webbber:把integrn写在ablquad(@(x,y)函数,Xmin,Xmax,Ymin,Ymax) 12/19 14:09
2F：→ webbber:虽然这样会变得很长但是不会有错 函数里除了XY不能有其他 12/19 14:11
3F：→ webbber:变数 ex: a1,a2,D 12/19 14:12
clear
Xmin=0;
Xmax=2*pi;
Ymin=pi;
Ymax=3*pi;
M12=dblquad(@(X,Y)
(55*30*cos(X-Y))/sqrt(55^2+30^2+0.8^2-2*55*30*cos(X-Y)),Xmin,Xmax,Ymin,Ymax);
这是我的程式,还是会出现错误!!
请问我应该怎麽改呢?
如果方便的话W大可以帮我试试看吗??
※ 编辑: thereocr 来自: 150.69.51.201 (12/19 16:10)
4F：→ webbber:./(sqrt(55^2+30^2+0.8^2-2*55*30*cos(X-Y))) 12/19 22:02
5F：→ webbber:除号前面要有点 sqrt()外面要有()包起来 12/19 22:03
W大谢谢 我弄出来了XDDD 真的超级感谢你的
※ 编辑: thereocr 来自: 150.69.51.201 (12/20 11:25)