matlab教程第四章數(shù)值計(jì)算-wenkub
2022-09-01 13:33:14 本頁(yè)面
【正文】 equ A = 9 11 21 63 252 70 69 141 421 1684 575 575 1149 3451 13801 3891 3891 7782 23345 93365 1024 1024 2048 6144 24572 DJ = 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 D = Columns 1 through 4 + 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 Column 5 0 0 0 0 c_eig = +010 * c_equ = +017 矩陣的譜分解和矩陣函數(shù) 【例 】數(shù)組乘方與矩陣乘方的比較。 A_Ap=A.^ A_Mp=A^ A_Ap = A_Mp = + + + + + 4 【例 】標(biāo)量的數(shù)組乘方和矩陣乘方的比較。) A_sinA = A_sinM = 多項(xiàng)式和卷積 多項(xiàng)式 一 多項(xiàng)式表 達(dá)方式的約定 二 多項(xiàng)式運(yùn)算函數(shù) 【例 】求 1 )1)(4)(2(32 ?? ??? ss sss 的“商”及“余”多項(xiàng)式。 cq=39。 disp([cq,poly2str(q,39。)]) 商多項(xiàng)式為 s + 5 余多項(xiàng)式為 5 s^2 + 4 s + 3 【例 】求 3階方陣 A的特征多項(xiàng)式。 PA=poly(A) PPA=poly2str(PA,39。 P=poly(R) PR=real(P) PPR=poly2str(PR,39。PP=poly2str(P,39。 b=[3,2,5,4,6]。 n=3。b39。MarkerSize39。 x0=0::1。cubic39。,x0,y0,39。,20),xlabel(39。（ B） 假設(shè) A 序列中最后 4 個(gè)非零值未知，而成為截尾序列，求卷積并圖示。 b=ones(1,8)。nc1=n1+n3。nn1=3。 nx=nn2+n4。N=length(kcc)。filled39。g39。,0) A=randn(1000,4)。,1) 8 X=rand(10,3)。 CCX=Xmx39。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)證，并圖示。 [X,Y]=meshgrid(x,y)。 DDZx0=DZx2*X。x39。 dt=。 ss10=dt*sum(y)。stend=st(end)。cumsum39。cumtrapz39。r39。y(x)39。cumtrapz39。 y_C=39。 clf,plot(x,Y,39。k39。)。zoom off 圖 局部放大和利用鼠標(biāo)取值圖 11 tt tt = [t4,y4,exitflag]=fzero(y_C,tt(4),[],P1,P2) Zero found in the interval: [, ]. t4 = y4 = 0 ex
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1