1樓：匿名使用者
樓上兩位的回答都很精彩，但也都存在一些問題：
1、chinasunsunsun的回答主要有以下問題（查了好半天才把問題查全）：
（1）樓主給的公式有兩項積分，但您只計算了一項，而且是把第一項的被積函式與第二項的積分限混在一起使用；
（2）樓主的n取值是0~5，您只取了0~4；
（3）畫圖的方式：很明顯θ作為角度量，畫圖應該按照柱面座標，但您是按照直角座標畫的（相應的，θ取值範圍也不恰當）。
2、aliksai回答的內容非常豐富，但樓主想要的是f對r、θ的三維影象而已，您前面的大量討論都不能反映樓主的真正需求，最後的那張圖似乎有那麼點意思了，但終究還是難見全貌，正確性也可疑，而且過程過於繁瑣。
以上點評純粹出於技術**，希望不要介意。
我寫的**如下，如果有問題歡迎批評：
f1=@(t,r,t1,r1)(-1./((r*cos(t)-r1.*cos(t1)).
^2+(r*sin(t)-r1.*sin(t1)).^2+0.
005^2).^(3/2)+...
1./((r*cos(t)-r1.*cos(t1)).
^2+(r*sin(t)-r1.*sin(t1)).^2+0.
025^2).^(3/2)).*(r*cos(t)-...
r1.*cos(t1)).*r1;
f2=@(t,r,t1,r1)(1./((r*cos(t)-r1.*cos(t1)).
^2+(r*sin(t)-r1.*sin(t1)).^2+0.
005^2).^(3/2)-...
1./((r*cos(t)-r1.*cos(t1)).
^2+(r*sin(t)-r1.*sin(t1)).^2+0.
025^2).^(3/2)).*(r*cos(t)-...
r1.*cos(t1)).*r1;
an = @(n,t,r) 8.2e4*dblquad(@(r1,t1)f1(t,r,t1,r1),0.09,0.
14,(n/3-1/12)*pi,(n/3+1/12)*pi) ...
+ 8.2e4*dblquad(@(r1,t1)f2(t,r,t1,r1),0.09,0.14,(n/3+1/12)*pi,(n/3+1/4)*pi);
f = @(t,r)sum(arrayfun(@(n)an(n,t,r),0:5));
r = linspace(0,.1,20);
t = linspace(0,2*pi,20);
[t,r]=meshgrid(t,r);
z = arrayfun(f,t,r);
[x,y,z]=pol2cart(t,r,z);
surf(x,y,z)
2樓：匿名使用者
^^rr= [0:0.001:0.01];
tt= [0:0.001:0.02];
f=@(r,t,r0,t0) (-1./((r0*cos(t0)-r.*cos(t)).
^2+(r0*sin(t0)-r.*sin(t)).^2+.
005^2).^1.5+1.
/((r0*cos(t0)-r.*cos(t)).^2+(r0*sin(t0)-r.
*sin(t)).^2+.025^2).
^1.5).*(r0*cos(t0)-r.
*cos(t)).*r;
y=zeros(length(rr),length(tt));
for i=1:length(rr)
for j=1:length(tt)
for k=0:4
y(i,j)=y(i,j)+dblquad(@(r,t)f(r,t,rr(i),tt(j)),0.09,0.14,pi/12+k*pi/3,pi/4+k*pi/3);
endend
endy=y*8.2*10^4;
surf(tt,rr,y)
計算二重積分，二重積分怎麼計算？
把積分割槽域分為三個x型區域，剩下的就是簡單的定積分的計算了，你把公式代進去算就行了，望採納。根據對稱性可知，積分項中的3x 與2x積分結果為零，所以積分項可以簡化為 x y 2y x y 1 1 再結合右圖分割槽域積分。二重積分怎麼計算？化為二次積分。x y dxdy 0 1 dx 1 2 x y...
二重積分運算，計算二重積分？
1,2 dx 1,x xydy 1,2 xdx 1,x ydy 1,2 xdx y 2 2 1,x 1 2 1,2 x x 2 1 dx 1 2 1,2 x 3 x dx 1 2 x 4 4 x 2 2 1,2 1 2 4 2 1 4 1 2 9 8 1 2 dx x2 2 x2 1 x x2 y ...
二重積分，d是什麼意思,二重積分是什麼？
1 在積分中，無論是定積分，還是不定積分 無論是一重積 分 二重積分 還是多重積分 d 都表示微分的概念，d differentiation derivative。2 d的用法 有時表示積分割槽域，d domian d integral area region 有時表示求導符號 dy y 這是尤拉 ...