【正文】 g 10。 g 5。 g 5。 hessian。V。 r。 通過運(yùn)行Su rface E vo l ve r 軟件, 可以用有限元方法求解包含表面勢能 、 重力勢能以及外力作用勢能等作用下的系統(tǒng)在能量最小時(shí)的靜平衡狀態(tài), 得出表面形態(tài)的估計(jì)結(jié)果, 并能在運(yùn)行中有效地跟蹤和控制迭代過程 。 主要應(yīng)用于計(jì)算在表面勢能和其他形式的能量作用下的液態(tài)靜平衡形態(tài) 。 但建立這樣的函數(shù)方程是非常困難的 。wS39。 centroid = STATS(i).Centroid。和39。 % 設(shè)置圖像標(biāo)題%step 5 尋找邊界線，參數(shù)‘noholes’將加快圖像的處理速度，它可以防止bwboundaries搜索內(nèi)部邊界 [B,L] = bwboundaries(BW, 39。 % 閾值BW = im2bw(GRAY, threshold)。 % 讀入圖像figure, % 新建窗口imshow(RGB), % 顯示原始圖像title(39。練習(xí)縮放、旋轉(zhuǎn)和平移結(jié)構(gòu)。這樣使選定的重復(fù)單元以球棍模型顯示。在三維窗口中右擊鼠標(biāo)，從快捷菜單中選擇DisplayStyle。選擇Advanced（高級）選項(xiàng)頁，設(shè)置Torsion（轉(zhuǎn)矩）為60。此例中你將使用Materials Studio中創(chuàng)建高聚物的工具建立20個(gè)單元的全同立構(gòu)的PMMA，后繼部分將模擬研究它的性質(zhì)。(28MB) or*v1k1=int(x*v1v1,2)f1=int(x*v139。,[0::4],[,]) 。,39。) y =1/(sin(x) 1) syms x 》subs(y,x,[::]) ans = 求數(shù)值解第一種方法：先定義m函數(shù)文件，如下：function li=lichun(x,y)li=y^2*cos(x)。 x0 = 1 1 1 》x0=G*x0+f （重復(fù)多次）x0 = 實(shí)驗(yàn)04 基于matlab的微分方程數(shù)值解例43 用龐格庫塔公式，求解下列初值問題。 x0 = 1 1 1 》x0=B*x0+f (重復(fù)多次)x0 = 緊湊迭代計(jì)算機(jī)算法緊湊迭代公式：例32 用簡單迭代格式解下列方程組》clc 》clear all 》A=[2 1 1。))*subs([f1 f2],[x y],xy039。diff(f2,39。syms x y。%clear all。方法1：》fsolve(39。fc=c^*c^2+*。) x = 方法2：新建一個(gè)script文件，輸入如下語句。a39。a(1)*x(1,:)+a(2)*x(2,:)+a(3)*x(3,:)39。a(1)*exp(a(2)*x)39。*A)*(A39。) f3 = Inline function: f3(c,x) = c(1)*x(1,:).^c(2).*x(2,:).^c(3) 》[a3,resnormal3,RESIDUAL3]=lsqcurvefit(f3,[ ],x3,y3) Local minimum found.Optimization pleted because the size of the gradient is less thanthe default value of the function tolerance.a3 = resnormal3 = RESIDUAL3 = * 例14利用解矛盾方程的方法，用二次多項(xiàng)式函數(shù)擬合以下數(shù)據(jù).序號XY序號XY134511222622313735400》A=[1 3 9。x39。a39。 9例22求代數(shù)方程x^32*x5=0,在x0=2附近的根。盡管在用語上不完全一樣，各國在研究領(lǐng)域和特色上也不盡相同，但基本含義是相同的?！队?jì)算機(jī)在材料科學(xué)中的應(yīng)用》上 機(jī) 教 程（2011級高分子專業(yè)用）湖北工業(yè)大學(xué)材料科學(xué)與工程學(xué)院2013年2月《計(jì)算機(jī)在材料科學(xué)中的應(yīng)用》上機(jī)教程 李純清前 言當(dāng)前，計(jì)算機(jī)在材料科學(xué)中的應(yīng)用(Compnter Applicationln Materials Science)已經(jīng)發(fā)展成為材料科學(xué)與工程的一個(gè)新的學(xué)科分支，涉及很多的學(xué)科和領(lǐng)域，需要把材料科學(xué)的專業(yè)知識與計(jì)算機(jī)技術(shù)、計(jì)算數(shù)學(xué)、信號與系統(tǒng)、物理、圖形圖像學(xué)、自動控制技術(shù)、機(jī)械等各方面的知識技術(shù)進(jìn)行綜合應(yīng)用。高技術(shù)新材料是現(xiàn)代知識經(jīng)濟(jì)的重要組成部分，其發(fā)展日新月異，這為材料的計(jì)算與設(shè)計(jì)提供了發(fā)展機(jī)遇和廣闊空間。 10例23 10例24用牛頓迭代法求方程在x0=1附近的零點(diǎn)。,39。) f2 = Inline function: f2(a,x) = a(1)+a(2)*x+a(3)*x.^2 》x2=3:3 x2 = 3 2 1 0 1 2 3 》y2=[4 2 3 0 1 2 5] y2 = 4 2 3 0 1 2 5 》[a2,resnormal2]=lsqcurvefit(f2,[, 1, 0],x2,y2) Local minimum found.Optimization pleted because the size of the gradient is less thanthe default value of the function tolerance.a2 = resnormal2 = 多變量擬合例13》x3=[100 200 300 500 100 700 800。1 2 4。*Y) ans = 例15》x5=[3 2 1 2 4] x5 = 3 2 1 2 4 》y5=[ ] y5 = 》f5=inline(39。,39。,39。,39。clc。while abs(fc)*10^(4) c=(a+b)/2。x^32*x539。syms x。f1=4x^2+y^2。x39。)39。1 5 1。本實(shí)驗(yàn)用matlab中的相關(guān)命令來求解。再在windows mand下執(zhí)行如下語句：》[x1,y1]=ode45(39。x39。 [t ,y] ans = 0 plot(t,y(:,1),t,y(:,2))。,2) dv0=diff(v0,x)dv1=diff(v1,x) dv0v0=dv039。without JavaMaterials Studio提供建立均聚物、嵌端共聚物、無規(guī)共聚物和dendrimer的選項(xiàng)。這樣你就設(shè)置了創(chuàng)建20個(gè)單元的全同PMMA。在Atom選項(xiàng)頁，單擊Line選項(xiàng)。保持重復(fù)單元選定著，在文檔空白處右擊鼠標(biāo)，在快捷菜單中選Label（標(biāo)識）命令，打開Label對話框。關(guān)閉保存文檔。Original Image39。 % 圖像黑白轉(zhuǎn)換figure, % 新建窗口imshow(BW), % 顯示處理后的圖像title(39。noholes39。Extent39。 switch W(i) case 1 plot(centroid(1),centroid(2),39。)。 而且在建立函數(shù)方程的過程中表面張力 、 表面勢能及穩(wěn)定性也應(yīng)當(dāng)考慮到方程中去, 這就進(jìn)一步增加了確定表面形狀的難度 。 目前已經(jīng)得到了廣泛的應(yīng)用, 如: 用于計(jì)算火箭燃料箱中低重力作用下的燃料形態(tài) 、 計(jì)算毛細(xì)作用表面等 。（或自行copy）。 g 5。 r。 }// Evolution with 90 contact and high gravitygogo3 := { angle := 90。 hessian。 r。 t .1。 re。 r。 r。 hessian。 V。 g 5。 可處理體積約束 、 邊界約束 、方 程約束 、 預(yù)設(shè)平均曲率等約束條件 。該軟件是受美國國家科學(xué)基金資助的幾何超級計(jì)算工程項(xiàng)目中的能量部分軟件 。 因此, 可以通過建立表面能量與表面形狀間的函數(shù)方程來確定表面形狀 。 case 3 plot(centroid(1),centroid(2),39。 W(i) = W(i) + 2 * uint8((STATS(i).Extent 1) == 0 )。BoundingBox39。)。 % 設(shè)置圖像標(biāo)題%step 3 對圖像進(jìn)行黑白閾值化，目的是為圖像邊界檢測作準(zhǔn)備threshold = graythresh(GRAY)。)。3. 研究結(jié)構(gòu) 在the Polymethyl_methacrylate三維窗口中空白處單擊，取消選擇?；氐斤@示選項(xiàng)對話框的Atom選項(xiàng)頁，單擊Ball and Stick按鈕，關(guān)閉對話框。首先改變結(jié)構(gòu)的顯示方式。在Chain length（鏈長）文本框中輸入20。它是典型的自由基聚合物，使用臭氧、偶氮化合物、熱或光引發(fā)反應(yīng)。32bit Java,1,)v