matlab教程第四章數(shù)值計(jì)算-全文預(yù)覽
【正文】 【例 】求方框圖（圖 ）所描述的那多環(huán)系統(tǒng)的傳遞函數(shù)。1 0 1]39。 A=[ 。,39。,39。 plot(t,x,39。 x=sin(2*pi*10*t)+cos(2*pi*100*t)+*randn(size(t))。 clear,randn(39。 plot(w(ii),Ya(ii)) grid,xlabel(39。 Ts=t(2)t(1) Ws=2*pi/Ts。state39。) 14 0 5 10 15 20 25 3002468nncc0 5 10 15 20 25 30101error 圖 變 換法和直接法求卷積結(jié)果比較 【例 】 fft 在信號(hào)分析中的應(yīng)用。) subplot(2,1,2),stem(nn,error,39。),grid,axis([0,31,0,9]) xlabel(39。 c(M)=0。BF=fft(b,M)。b([1,2])=0。所給已知序列為 ??? ?? e lsenna 12,4,301)( ? ， ??? ?? e lsennb 9,3,201)( ?。 [tt,yy]=ode45(DyDt,tspan,y0)。mu*(1y(1)^2)*y(2)y(1)]。 （ 2） Hf=fun。[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。s Banana 測(cè)試函數(shù)，它的理論極小值是 1,1 ?? yx 。y(t)39。 xlabel(39。)。 x=10::10。 P1=。cunsum39。r.39。b:39。trapz39。sum39。 st=cumtrapz(t,y)。 y=exp(*t.*abs(sin(t)))。y39。DDZ20=DZ24。 [DZx,DZy]=gradient(Z,dd,dd)。 x=1:dd:1。 F=[1,2,3。 mx=mean(X)。 rand(39。 randn(39。b39。r39。 kx=(ncc2+1):nx。 cc=conv(aa,b)。 kc=nc1:nc2。n4=9。n1=3。) 0 120246810x 圖 通過(guò)三次多項(xiàng)式插值所得的曲線 卷積 一 離散序列的數(shù)值卷積 二 MATLAB 的“卷積”指令 7 【例 】有序列??? ?? e lsennA 12,4,301)( ? 和 ??? ?? e lsennB 9,3,201)( ?。,39。 plot(xi,yi,39。 xi=0::1。x39。.r39。 yy=polyval(P,xx)。（見(jiàn)圖 ） x0=0::1。) PA=polyval(P,S) PM=polyvalm(P,S) S = 1 1 1 1 1 2 3 4 1 3 6 10 1 4 10 20 PP = s^4 29 s^3 + 72 s^2 29 s + 1 PA = +004 * PM = * 【例 】部分分式展開(kāi)。) P = PR = PPR = x^3 + x^2 + x + 【例 】?jī)煞N多項(xiàng)式求值指令的差別。) PA = PPA = s^3 45 s^2 18 s 【例 】由給定根向量求多項(xiàng)式系數(shù)向量。14 15 16。)]),disp([cr,poly2str(r,39。 cr=39。 p2=[1 0 1 1]。（ A取自例 ） A_sinA=sin(A) A_sinM=funm(A,39。4 5 6。 [V,D,c_eig]=condeig(A)。2,2/3]。 x=inv(A39。其中 A 取 MATLAB庫(kù)中的特殊函數(shù)生成。 cond(A) ans = +013 （ 2） tic xi=inv(A)*b。randsvd39。 1 第四章 數(shù)值計(jì)算 線性方程組的解 LU 分解、行列式、逆和恰定方程的解 【例 】“求逆”法和“左除”法解恰定方程的性能對(duì)比 （ 1） randn(39。 A=gallery(39。 b=A*x。 td=toc, erd=norm(xxd), red=norm(A*xdb)/norm(b) td = erd = red = 奇異值分解和矩陣結(jié)構(gòu) 線性二乘問(wèn)題的解 【例 】對(duì)于超定方程 Axy? ，進(jìn)行三種解法比較。y=[ ]39。 A=[1,3。 A=gallery(5) [VJ,DJ]=jordan(A)。 clear,A=[1 2 3。（ A取自例 ） pA_A=().^A pA_M=()^A pA_A = pA_M = 【例 】 sin的數(shù)組運(yùn)算和矩陣運(yùn)算比較。 p1=conv([1,0,2]
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1