如何用matlab進行多元函數(shù)偏導數(shù)計算
創(chuàng)新互聯(lián)是專業(yè)的江西網(wǎng)站建設公司,江西接單;提供網(wǎng)站建設、網(wǎng)站制作,網(wǎng)頁設計,網(wǎng)站設計,建網(wǎng)站,PHP網(wǎng)站建設等專業(yè)做網(wǎng)站服務;采用PHP框架,可快速的進行江西網(wǎng)站開發(fā)網(wǎng)頁制作和功能擴展;專業(yè)做搜索引擎喜愛的網(wǎng)站,專業(yè)的做網(wǎng)站團隊,希望更多企業(yè)前來合作!
可以調用 diff 函數(shù)求導。
舉例說明:
先定義符號 x、y 以及符號二元函數(shù)表達式 z,然后調用 diff 函數(shù)求偏導,代碼如下:
clc;clear;
syms x y
z=x^2+y^2+exp(x*y);
z_x=diff(z,x,1)
z_y=diff(z,y,1)
z_x2=diff(z,x,2)
z_y2=diff(z,y,2)
z_xy=diff(z,x,y)
結果如下:
z_x =
2*x + y*exp(x*y)
z_y =
2*y + x*exp(x*y)
z_x2 =
y^2*exp(x*y) + 2
z_y2 =
x^2*exp(x*y) + 2
z_xy =
exp(x*y) + x*y*exp(x*y)
把
for?i?=?2:100
a(i)=1;??%a2=a3=a4=...=a_num=1?
end
改成
for?i?=?2:100
eval(['a',num2str(i)?'=1;'])
end
你原來的寫法是把數(shù)組a的第2到第100個元素(即a(2)-a(100))賦值為1,而符號變量a2-a100仍然存在。
類型的異常在System.Data.dll中發(fā)生,但未在用戶代碼中進行處理,是設置錯誤造成的,解決方法如下:
1、首先,需要確保IIS是否啟用32位應用程序,因為OLEDB不支持64位。鍵盤上Win+R打開運行;輸入"inetmgr"打開IIS。
2、點擊IIS中的應用程序池,并選擇某一程序池后。
3、點擊右邊的“高級設置”選項。
4、在IIS高級設置對話框中找到“啟用32位應用程序”選項并設置為true。
5、設置完成后,點擊窗口下方的確定,問題就解決了。
matlab里一階導數(shù)用diff函數(shù)實現(xiàn),如:
syms t
diff(sin(t^2))
ans =
2*cos(t^2)*t
第二步在我們的電腦上打開matlab,在命令行窗口中輸入syms x,f(x)=sin(x)+x^2 ,diff(f(x)對f(x)函數(shù)進行求導
第三步按回車鍵,可以看見求導的結果是2*x+cos(x)
第四步在命令行窗口中輸入diff(f(x),3),按回車鍵求f(x)函數(shù)的3階導數(shù)
第五步也可以求偏導,輸入diff(f(x),x)求對
可以調用 diff 函數(shù)求導。
舉例說明:
先定義符號 x、y 以及符號二元函數(shù)表達式 z,然后調用 diff 函數(shù)求偏導,代碼如下:
clc;clear;
syms x y
z=x^2+y^2+exp(x*y);
z_x=diff(z,x,1)
z_y=diff(z,y,1)
z_x2=diff(z,x,2)
z_y2=diff(z,y,2)
z_xy=diff(z,x,y)
結果如下:
z_x =
2*x + y*exp(x*y)
z_y =
2*y + x*exp(x*y)
z_x2 =
y^2*exp(x*y) + 2
z_y2 =
x^2*exp(x*y) + 2
z_xy =
exp(x*y) + x*y*exp(x*y)
圖中調用了5次diff函數(shù),分別計算了 z 對 x 的一階偏導,z 對 y 的一階偏導,z 對 x 的二階偏導,z 對 y 的二階偏導,以及 z 對 x,y 都求一次偏導。
例子中所用的調用格式為: diff(f,var,n)
f 為符號表達式,也可以是符號函數(shù)(這個數(shù)據(jù)類型低版本的matlab沒有)。n 為求導次數(shù),缺省為1。var 為求導的符號變量,可以缺?。╩atlab會根據(jù)表達式自己選擇一個),但不建議缺省,除非表達式只含有一個符號變量。
此外,matlab還允許 diff(f,n,var)與 diff(S,v1,v2,...,vn)的調用形式。
diff(f,v1,v2,...,vn) 會把表達式 f 對變量 v1,v2 等 n 個變量都求一次偏導,得到 f 的 n 階偏導。
題名:科學與工程數(shù)值算法-Java版
作者:丁軍, 楊麗麗
出版社:清華大學出版社,2003
這本書挺好的,是去年我學《數(shù)值分析》的法寶。上面有微分方程的數(shù)值算法代碼,推薦參考一下