作者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)