[發問] 二重積分的運算問題
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)
推
12/19 14:09, , 1F
12/19 14:09, 1F
→
12/19 14:11, , 2F
12/19 14:11, 2F
→
12/19 14:12, , 3F
12/19 14:12, 3F
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)
→
12/19 22:02, , 4F
12/19 22:02, 4F
→
12/19 22:03, , 5F
12/19 22:03, 5F
W大謝謝 我弄出來了XDDD 真的超級感謝你的
※ 編輯: thereocr 來自: 150.69.51.201 (12/20 11:25)