matlab教程第四章數值計算-文庫吧資料
2024-08-29 13:33本頁面
【正文】 t 在信號分析中的應用。),axis([0,31,1,1]) ylabel(39。) subplot(2,1,2),stem(nn,error,39。),ylabel(39。),grid,axis([0,31,0,9]) xlabel(39。 subplot(2,1,1),stem(nn,c,39。 c(M)=0。 cc=real(ifft(CF))。BF=fft(b,M)。 M=32。b([1,2])=0。a([1,2,3])=0。所給已知序列為 ??? ?? e lsenna 12,4,301)( ? ， ??? ?? e lsennb 9,3,201)( ?。x(t)39。 [tt,yy]=ode45(DyDt,tspan,y0)。 y0=[1。mu*(1y(1)^2)*y(2)y(1)]。（見圖 ） （ 1） （ 2） [] function ydot=DyDt(t,y) mu=2。 （ 2） Hf=fun。 數值積分 【例 】求 dxeI x? ?? 10 2，其精確值為 ? 。[sx,sfval,sexit,soutput]=fminsearch(ff,x0) Optimization terminated successfully: the current x satisfies the termination criteria using of and F(X) satisfies the convergence criteria using of sx = sfval = sexit = 1 soutput = iterations: 85 funcCount: 159 algorithm: 39。)。,39。 ff=inline(39。s Banana 測試函數，它的理論極小值是 1,1 ?? yx 。zoom off 圖 局部放大和利用鼠標取值圖 11 tt tt = [t4,y4,exitflag]=fzero(y_C,tt(4),[],P1,P2) Zero found in the interval: [, ]. t4 = y4 = 0 exitflag = 1 求函數極值點 【例 】求 222 )1()(1 0 0),( xxyyxf ???? 的極小值點。y(t)39。)。 xlabel(39。k39。)。 clf,plot(x,Y,39。 x=10::10。 y_C=39。 P1=。cumtrapz39。cunsum39。y(x)39。r.39。r39。b:39。cumtrapz39。trapz39。cumsum39。sum39。stend=st(end)。 st=cumtrapz(t,y)。 ss10=dt*sum(y)。 y=exp(*t.*abs(sin(t)))。 dt=。y39。x39。DDZ20=DZ24。 DDZx0=DZx2*X。 [DZx,DZy]=gradient(Z,dd,dd)。 [X,Y]=meshgrid(x,y)。 x=1:dd:1。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 【例 】函數 22),( yxyxz ?? 的梯度 ]22[),( yxyxz ?? 和 4),(2 ?? yxz ，用數值計算驗證，并圖示。 F=[1,2,3。 CCX=Xmx39。 mx=mean(X)。,1) 8 X=rand(10,3)。 rand(39。,0) A=randn(1000,4)。 randn(39。g39。b39。filled39。r39。N=length(kcc)。 kx=(ncc2+1):nx。 nx=nn2+n4。 cc=conv(aa,b)。nn1=3。 kc=nc1:nc2。nc1=n1+n3。n4=9。 b=ones(1,8)。n1=3。（ B） 假設 A 序列中最后 4 個非零值未知，而成為截尾序列，求卷積并圖示。) 0 120246810x 圖 通過三次多項式插值所得的曲線 卷積 一 離散序列的數值卷積 二 MATLAB 的“卷積”指令 7 【例 】有序列??? ?? e lsennA 12,4,301)( ? 和 ??? ?? e lsennB 9,3,201)( ?。,20),xlabel(39。,39。,x0,y0,39。 plot(xi,yi,39。cubic39。 xi=0::1。 x0=0::1。x39。MarkerSize39。.r39。b39。 yy=polyval(P,xx)。 n=3。（見圖 ） x0=0::1。 b=[3,2,5,4,6]。) PA=polyval(P,S) PM=polyvalm(P,S) S = 1 1 1 1 1 2 3 4 1 3 6 10
環(huán)評公示相關推薦
鄂ICP備17016276號-1