matlab教程第四章數(shù)值計(jì)算-免費(fèi)閱讀
【正文】 0 ]。 18 P2=minreal(G2*P1/(1+G2*P1*H2/G4))。G2=tf(1,[1 1])。 B=[0 0 4。Input39。 y=filter(B,A,x)。在此 )()10 02c os ()102s in()( tntttx ??????? ?? ，而 ),0(~)( Ntn 。 plot(t,y) 0 2 4 6 8 10201510505101520 圖 受噪聲污染的信號 Y=fft(y)。error39。fill39。 AF=fft(a,M)。) 0 5 10 15 20 25 3021012x(t) 圖 微分方程解 （ 4） plot(yy(:,1),yy(:,2)) 13 2 1 0 1 2 432101234 圖 相平面軌跡 信號處理 快速 Fourier 變換和逆變換 【例 】利用 fft 和 ifft 指令重新計(jì)算兩序列的卷積。 ydot=[y(2)。 x0=[,1]。它即是著名的 Rosenbrock39。)。sin(x).^2.*exp(P1*x)P2*abs(x)39。,39。]) disp([ss10,ssend,st10,stend]) plot(t,y,39。 disp([blanks(5),39。t=(0:dt:10)39。DDZy0=DZy2*Y。（見圖 ） clear, dd=。Y=rand(10,3)。),hold off 4 6 8 10 12 14 16 18 20 22012345678910 圖 “完整”序列卷積和“截尾”序列卷積 數(shù)據(jù)分析函數(shù) 隨機(jī)數(shù)發(fā)生器和 統(tǒng)計(jì)分析指令 【例 】基本統(tǒng)計(jì)示例。 stem(kcc,cc(1:N),39。nn2=8。n3=2。x39。)。,20),xlabel(39。 P=polyfit(x0,y0,n) xx=0::1。s39。s39。s39。 p1=conv([1,0,2],conv([1,4],[1,1]))。 clear,A=[1 2 3。 A=[1,3。 td=toc, erd=norm(xxd), red=norm(A*xdb)/norm(b) td = erd = red = 奇異值分解和矩陣結(jié)構(gòu) 線性二乘問題的解 【例 】對于超定方程 Axy? ，進(jìn)行三種解法比較。 A=gallery(39。randsvd39。其中 A 取 MATLAB庫中的特殊函數(shù)生成。2,2/3]。4 5 6。 p2=[1 0 1 1]。)]),disp([cr,poly2str(r,39。) PA = PPA = s^3 45 s^2 18 s 【例 】由給定根向量求多項(xiàng)式系數(shù)向量。) 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 = * 【例 】部分分式展開。 yy=polyval(P,xx)。x39。 plot(xi,yi,39。) 0 120246810x 圖 通過三次多項(xiàng)式插值所得的曲線 卷積 一 離散序列的數(shù)值卷積 二 MATLAB 的“卷積”指令 7 【例 】有序列??? ?? e lsennA 12,4,301)( ? 和 ??? ?? e lsennB 9,3,201)( ?。n4=9。 cc=conv(aa,b)。r39。 randn(39。 mx=mean(X)。 x=1:dd:1。DDZ20=DZ24。 y=exp(*t.*abs(sin(t)))。sum39。b:39。cunsum39。 x=10::10。 xlabel(39。s Banana 測試函數(shù)，它的理論極小值是 1,1 ?? yx 。[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。mu*(1y(1)^2)*y(2)y(1)]。所給已知序列為 ??? ?? e lsenna 12,4,301)( ? ， ??? ?? e lsennb 9,3,201)( ?。BF=fft(b,M)。),grid,axis([0,31,0,9]) xlabel(39。) 14 0 5 10 15 20 25 3002468nncc0 5 10 15 20 25 30101error 圖 變 換法和直接法求卷積結(jié)果比較 【例 】 fft 在信號分析中的應(yīng)用。 Ts=t(2)t(1) Ws=2*pi/Ts。 clear,randn(39。 plot(t,x,39。,39。1 0 1]39。 G3=tf([1 1],[1 4 4])。
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1