MATLAB計算二重數值積分(dblquad)
當被積分函數爲二元函數f(x,y),積分變量爲dxdy時稱爲二重積分。MATLAB提供了dblquad( )函數命令計算二重數值積分。dblquad的調用格式爲Q=dblquad(function,xmin,xmax,ymin,ymax,tol,’Method’),其中function爲被積分的二元函數,xmin,xmax,ymin,ymax分別爲積分變量x,y的上下限,tol爲絕對計算精度,method是計算一維積分的方法,一般有Simpson法(即quad,默認)和Lobatto法(即quadl),Q是返回計算的二重數值積分結果。
操作方法
(01)第一,使用MATLAB計算下圖中的二重積分。
(02)第二,啓動MATLAB,新建腳本(Ctrl+N),輸入如下代碼:close all; clear all; clcfigure(x27;Position',[50,50,600,500],'Name','double integration')f=@(x,y) x*sin(y)-cos(x)+y-3;Q1=dblquad(f,0,2*pi,-pi,pi,1.0e-3)x=linspace(0,2*pi);y=linspace(-pi,pi);ff=x.*sin(y)-cos(x)+y-3;plot3(x,y,ff,'r','LineWidth',3)grid on; view(35,35)xlabel('x','FontSize',15); ylabel('y','FontSize',15)其中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的圖像。
(03)第三,儲存和執行上述腳本,在命令行視窗(Command Window)得到如下結果:Q1 =-118.4351也就是說,該二重數值積分的結果爲-118.4351.
(04)第四,同時得到被積分函數x.*sin(y)-cos(x)+y-3的圖像。
(05)第五,在Q1=dblquad(f,0,2*pi,-pi,pi,1.0e-3)後面再添加一行命令,Q2=dblquad(f,0,2*pi,-pi,pi,1.0e-3,'quadl')。也就是在絕對計算精度都是1.0e-3的情況下,使用quadl方法計算二重數值積分。
(06)第六,儲存和執行上述腳本,在命令行視窗(Command Window)得到採用quad方法(默認)和quadl方法計算的二重數值積分結果。Q1 =-118.4351Q2 =-118.4355也就是說,採用採用quad方法(默認)計算的二重數值積分結果爲-118.4351,而採用quadl方法計算的二重數值積分結果爲-118.4355,兩者在小數點第四位略有差別。
特別提示
按照積分變量dxdy的順序,x稱爲內積分變量,y稱爲外積分變量。
-
Nvidia顯卡設定/Nvidia控制面板設定
下面簡單介紹一下如何進行Nvidia顯卡設定,Nvidia控制面板的設定。操作方法(01)在搜尋欄中輸入【控制面板】,回車。(02)開啟控制面板介面,點擊檢視方式出的【倒三角】,然後選擇【小圖標】。(03)找到【NVIDIA控制面板】,並開啟。(04)開啟後,在左邊一框的【3D設定】下列,點擊【通...
-
怎麼在Excel列A數據中篩選出列B數據?
已知excel列A數據和列B數據有部分內容相同,怎樣才能篩選出哪些混在列A中的列B數據呢?最笨的方法就是複製列B的數據查找。這裏,我給大家分享一種方法,既準確又快捷的方法。不信,您嘗試一下。操作方法(01)開啟需要編輯的Excel檔案,觀察列A和列B,列A列B有部分內容重複,現需...
-
如何利用ps去水印
當你在網絡上看到一張圖片的時候很想儲存下來的時候,卻無奈的發現這張圖片有水印,通常我們使用的方法就是截圖,用裁剪的方式去掉有水印的地方,但是如此就會破壞整張圖的佈局美感,下面就教你如何用ps的填充功能快速去除水印。操作方法(01)相信大家對於如何安裝ps已經很...
-
怎麼做音樂相冊?
把照片做成音樂相冊,不僅可以動感地展示照片,還可以配上好聽的音樂來營造氛圍。做一份婚禮、畢業、家人、兒童等等主題的音樂相冊,可以爲自己和他人留下一份珍貴美好的記憶~從朋友那裏學到了做音樂相冊的方法,今天就換我來教大家怎麼輕鬆做一份好看、好聽的音樂相冊...