【正文】 he function tolerance.a5 = resnorm5 = 例16》x6=[1 2 3 4 5 6 7 8] x6 = 1 2 3 4 5 6 7 8 》y6=[ ] y6 = 》f6=inline(39。,39。1 1 2 5] x7 = 1 1 2 3 1 3 5 1 1 1 2 5 》y7=[2 1 1 2] y7 = 2 1 1 2 》f7=inline(39。,39。,39。) fun8 = Inline function: fun8(a,x) = exp(a(1)a(2)./(x+a(3))) 》[a8,resnormal8]=lsqcurvefit(fun8,[14,2000,15],x8,y8) Local minimum possible.lsqcurvefit stopped because the final change in the sum of squares relative to its initial value is less than the default value of the function tolerance.a8 = +03 * resnormal8 = 》y88=exp(a8(1)a8(2)./(x8+a8(3))) y88 = +03 * 》(y88y8)./y8 ans = 實驗02：非線性方程求解例21 用對分法求f(x)=x^*x^2+*【1，2】之間的根。x39。a=1。c=(a+b)/2。 if fb*fc0 b=c。 endendformat longfx=fc, x=c 》half2root fx = x = 例22求代數(shù)方程x^32*x5=0,在x0=2附近的根。文件內(nèi)容如下：function F=nonlineequation(x)F=[*x(1)^2x(1)+*x(2)^2+1。 err=norm(x_tmpx0,inf) x0=x_tmp endx0err》Newton_method_solve_equation f =x^3 (77*x^2)/10 + (96*x)/5 153/10x0 = 1err = 1err = x0 = err = x0 = err = x0 = err = x0 = err = x0 = err = x0 = x0 = err = 例25用割線法求方程的根（取x0=,x1=4）。while(err) x_tmp=x1subs(f,x,x1)*(x1x0)/(subs(f,x,x1)subs(f,x,x0))。 endx1err f =x^3 (77*x^2)/10 + (96*x)/5 153/10x0 = x1 = 3x1 = err = 例26：用牛頓迭代法解下列非線性方程組%clear all。Jac=[diff(f1,39。)。y39。while(err) xy_tmp=inv(subs(Jac,[x y],xy039。 xy0=xy0+xy_tmp。 b = 5 8 11 》D=diag(diag(A)) D = 2 0 0 0 5 0 0 0 10 》B=eye(3)inv(D)*A B = 0 0 0 》f=inv(D)*b f = 》x0=[1 1 1]39。 b = 5 8 11 》D=diag(diag(A)) D = 2 0 0 0 5 0 0 0 10 》L=tril(A,1) L = 0 0 0 1 0 0 1 1 0 》U=triu(A,1) U = 0 1 1 0 0 1 0 0 0 》G=(DL)^1*U G = 0 0 0 》f=(DL)^1*b f = 》x0=[1 1 1]39。 b = 5 8 11 》omiga= omiga = 》D=diag(diag(A)) D = 2 0 0 0 5 0 0 0 10 》L=tril(A,1) L = 0 0 0 1 0 0 1 1 0 》U=triu(A,1) U = 0 1 1 0 0 1 0 0 0 》G=(Domiga*L)^1*((1omiga)*D+omiga*U) G = 》f=omiga*(Domiga*L)^1*b f = 》x0=[1 1 1]39。Dyy^2*cos(x)39。x39。,[0, ],[1]) 》[x2,y2]=ode45(39。y^2*cos(x)39。y39。exam4_439。本微分方程的物理意義是表示：在不考慮流體本身熱傳導(dǎo)時的套管傳熱微分方程。,1,)v1=[q10 q11 q12]v1v1=v139。*dv1 dk0=int(x^2*dv0v0,1,)dk1=int(x^2*dv1v1,2) kk0=k0+dk0kk1=k1+dk1 kkk0=zeros(5)kkk0(1:3,1:3)=kk0kkk1=zeros(5)kkk1(3:5,3:5)=kk1 k=kkk0+kkk1fff0=zeros(5,1)fff0(1:3)=f0fff1=zeros(5,1)fff1(3:5)=f1 f=fff0+fff1 syms u1d u2d realk_sym=sym(k)f_sym=sym(f)k_sym(1,1)=k_sym(1,1)(2u1d)/2k_sym(5,5)=k_sym(5,5)+(44*u2d)/4f_sym(1)=f_sym(1)2f_sym(5)=f_sym(5)+4u01234=inv(k_sym)*f_sym [u1d u2d]=solve(u01234(1)2,u01234(5)4,u1d,u2d) u01234=eval(subs(u01234))中篇 分析應(yīng)用篇實驗07：泡沫材料泡孔（顆粒）大小及形狀分析實驗?zāi)康模毫私鈭D像處理一般過程方法技巧及原理熟悉常用圖像處理軟件的用法及功能了解顆粒大小及形狀分析實驗原理：粒徑分析中的圖像處理方法及算法實驗軟件：WindowsDownload ImageJ 32bit Java或用繪圖軟件自行繪制一張。它是典型的自由基聚合物，使用臭氧、偶氮化合物、熱或光引發(fā)反應(yīng)。在Polymerize（聚合）選項頁，單擊Lirary下拉菜單，選擇acrylates（丙烯酸酯）。在Chain length（鏈長）文本框中輸入20。，在一個三維顯示出它的結(jié)構(gòu)。首先改變結(jié)構(gòu)的顯示方式。在PMMA結(jié)構(gòu)中選定任一原子（你可能需要放大一下才行）?；氐斤@示選項對話框的Atom選項頁，單擊Ball and Stick按鈕，關(guān)閉對話框。在Properties（屬性）中選擇Name（名稱）。3. 研究結(jié)構(gòu) 在the Polymethyl_methacrylate三維窗口中空白處單擊，取消選擇。請你按照期刊論文的格式及排版布局，用office、origin等相關(guān)軟件仿摹成一篇word電子檔的文件。)。 % 設(shè)置圖像標(biāo)題%step 2 將RGB圖像轉(zhuǎn)為灰度圖像GRAY = rgb2gray(RGB)。 % 設(shè)置圖像標(biāo)題%step 3 對圖像進(jìn)行黑白閾值化，目的是為圖像邊界檢測作準(zhǔn)備threshold = graythresh(GRAY)。)。)。 % 尋找邊界%step 6 測定形狀屬性STATS = regionprops(L, 39。BoundingBox39。%step 7 根據(jù)形狀屬性對形狀進(jìn)行分類識別figure, % 新建窗口imshow(RGB), % 顯示原始圖像title(39。 W(i) = W(i) + 2 * uint8((STATS(i).Extent 1) == 0 )。)。 case 3 plot(centroi