matlab教程第四章數(shù)值計(jì)算(編輯修改稿)
2024-09-26 13:33 本頁面
【文章內(nèi)容簡(jiǎn)介】 AMAX = AMIN = AMED = AMEAN = ASTD = 【例 】 cov和 corrcoef的使用示例。 rand(39。state39。,1) 8 X=rand(10,3)。Y=rand(10,3)。 mx=mean(X)。Xmx=Xones(size(X))*diag(mx)。 CCX=Xmx39。*Xmx/(size(Xmx,1)1) CX=cov(X),CY=cov(Y) Cxy=cov(X,Y) PX=corrcoef(X) Pxy=corrcoef(X,Y) CCX = CX = CY = Cxy = PX = Pxy = 差分和累計(jì)指令 【例 】用一個(gè)簡(jiǎn)單矩陣表現(xiàn) diff和 gradient指令計(jì)算方式。 F=[1,2,3。4,5,6。7,8,9] Dx=diff(F) Dx_2=diff(F,1,2) [FX,FY]=gradient(F) [FX_2,FY_2]=gradient(F,) F = 1 2 3 4 5 6 7 8 9 Dx = 3 3 3 3 3 3 Dx_2 = 1 1 1 1 1 1 FX = 1 1 1 1 1 1 1 1 1 FY = 3 3 3 3 3 3 3 3 3 9 FX_2 = 2 2 2 2 2 2 2 2 2 FY_2 = 6 6 6 6 6 6 6 6 6 【例 】函數(shù) 22),( yxyxz ?? 的梯度 ]22[),( yxyxz ?? 和 4),(2 ?? yxz ，用數(shù)值計(jì)算驗(yàn)證，并圖示。（見圖 ） clear, dd=。 x=1:dd:1。y=x。 [X,Y]=meshgrid(x,y)。 Z=(X.^2)+(Y.^2)。 [DZx,DZy]=gradient(Z,dd,dd)。 DZ2=4*del2(Z,dd,dd)。 DDZx0=DZx2*X。DDZy0=DZy2*Y。DDZ20=DZ24。 subplot(1,3,1),stem3(X,Y,DDZx0) subplot(1,3,2),stem3(X,Y,DDZy0) subplot(1,3,3),stem3(X,Y,DDZ20) axis([1,1,1,1,]) xlabel(39。x39。),ylabel(39。y39。) 101101010110101011010xy 圖 理論計(jì)算和數(shù)值計(jì)算的差別圖示 【例 】求積分 ?? x dttyxs0 )()(，其中 ttety )( ?? ， 100 ??x 。 dt=。t=(0:dt:10)39。 y=exp(*t.*abs(sin(t)))。 ss=dt*cumsum(y)。 ss10=dt*sum(y)。ssend=ss(end)。 st=cumtrapz(t,y)。 st10=trapz(t,y)。stend=st(end)。 disp([blanks(5),39。sum39。,blanks(6),39。cumsum39。,blanks(4),39。trapz39。,blanks(5),39。cumtrapz39。]) disp([ss10,ssend,st10,stend]) plot(t,y,39。b:39。,t,ss,39。r39。,t,st,39。r.39。) legend(39。y(x)39。,39。cunsum39。,39。cumtrapz39。,0) sum cumsum trapz cumtrapz 10 0 2 4 6 8 100123y(x) cunsum cumtrapz 圖 矩形法和梯形法求積比較 MATLAB 泛函指令 求函數(shù)零點(diǎn) 【例 】通過求 tbettf at ?? ?)( s in)( 2 的零點(diǎn)，綜合敘述相關(guān)指令的用法。 P1=。P2=。 y_C=39。sin(x).^2.*exp(P1*x)P2*abs(x)39。 x=10::10。 Y=eval(y_C)。 clf,plot(x,Y,39。r39。)。hold on,plot(x,zeros(size(x)),39。k39。)。 xlabel(39。t39。)。ylabel(39。y(t)39。),hold off 10 5 0 5 105432101ty(t) 圖 函數(shù)零點(diǎn)分布觀察圖 zoom on [tt,yy]=ginput(5)。zoom off 圖 局部放大和利用鼠標(biāo)取值圖 11 tt tt = [t4,y4,exitflag]=fzero(y_C,tt(4),[],P1,P2) Zero found in the interval: [, ]. t4 = y4 =
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1