gibbs抽樣求二重積分代碼的簡單介紹

本文目錄一覽：

• 1、用c語言求二重積分~~
• 2、matlab中如何用複合辛普森公式求二重積分
• 3、MATLAB如何求二重積分

用c語言求二重積分~~

這個你用手工積分就行了，當然常數的代入計算還是可以用計算器的。給你的提示是分0ya/2和a/2ya兩次脫絕對值符號，最後都會演變為對sqrt(x^2+u^2)/(x^2+u^2+v^2)^2的積分形式，或者，你只關心數值結果的話，直接採用矩形法或者辛普森積分法，劃分的細一點就不會有太大誤差了

matlab中如何用複合辛普森公式求二重積分

代碼如下：

function q=DblSimpson(f,a,A,b,B,m,n)

if(m==1 n==1) %辛普森公式

q=((B-b)*(A-a)/9)*(subs(sym(f),findsym(sym(f)),{a,b})+…

subs(sym(f),findsym(sym(f)),{a,B})+…

subs(sym(f),findsym(sym(f)),{A,b})+…

subs(sym(f),findsym(sym(f)),{A,B})+…

4*subs(sym(f),findsym(sym(f)),{(A-a)/2,b})+…

4*subs(sym(f),findsym(sym(f)),{(A-a)/2,B})+…

4*subs(sym(f),findsym(sym(f)),{a,(B-b)/2})+…

4*subs(sym(f),findsym(sym(f)),{A,(B-b)/2})+…

16*subs(sym(f),findsym(sym(f)),{(A-a)/2,(B-b)/2}));

else %複合辛普森公式

q=0;

for i=0:n-1

for j=0:m-1

x=a+2*i*(A-a)/2/n;

y=b+2*j*(B-b)/2/m;

x1=a+(2*i+1)*(A-a)/2/n;

y1=b+(2*j+1)*(B-b)/2/m;

x2=a+2*(i+1)*(A-a)/2/n;

y2=b+2*(j+1)*(B-b)/2/m;

q=q+subs(sym(f),findsym(sym(f)),{x,y})+…

subs(sym(f),findsym(sym(f)),{x,y2})+…

subs(sym(f),findsym(sym(f)),{x2,y})+…

subs(sym(f),findsym(sym(f)),{x2,y2})+…

4*subs(sym(f),findsym(sym(f)),{x,y1})+…

4*subs(sym(f),findsym(sym(f)),{x2,y1})+…

4*subs(sym(f),findsym(sym(f)),{x1,y})+…

4*subs(sym(f),findsym(sym(f)),{x1,y2})+…

16*subs(sym(f),findsym(sym(f)),{x1,y1});

end

end

end

q=((B-b)*(A-a)/36/m/n)*q;

叫我雷鋒！

MATLAB如何求二重積分

在MATLAB軟件中輸入二重積分的代碼即可求二重積分，具體操作請參照以下步驟，演示軟件版本為MATLAB 2014版。

1、將要使用MATLAB計算下圖中的二重積分，首先在電腦上打開MATLAB軟件。

2、新建腳本（Ctrl+N），輸入圖中框住的代碼內容。其中Q1=dblquad(f,0,2*pi,-pi,pi,1.0e-3)採用默認方法quad計算二重積分，絕對計算精度設為1.0e-3。plot3(x,y,ff,’r’,’LineWidth’,3)是繪製被積分函數ff=x.*sin(y)-cos(x)+y-3的圖像。

3、保存和運行上述腳本，在命令行窗口（Command Window）得到如下結果：Q1 =-118.4351。也就是說，該二重數值積分的結果為-118.4351。

4、同時得到被積分函數x.*sin(y)-cos(x)+y-3的圖像。

5、也可採用quadl法計算二重積分，在Q1命令後面再添加一行命令Q2=dblquad(f,0,2*pi,-pi,pi,1.0e-3,’quadl’)即可。