【文章內(nèi)容簡介】 d off legend(39。x(t)39。,39。dx/dt39。) / 17 plot(t,x,39。LineWidth39。,5)xlabel(39。t39。) ，采用 diff和 gradient計(jì)算該函數(shù)在區(qū)間 [0, 中的近似導(dǎo)函數(shù)，比較二者區(qū)別。 ，其中 ，并以符號計(jì)算結(jié)果為參考解，通過誤差分析比較上述方法的精度。 fminseach 函數(shù)求著名的 Rosenbrock39。s “Banana” 測試函數(shù) 。 的極小值點(diǎn)。 陣 ， ，采用三種不同的編程求這兩個(gè)矩陣的乘積 。并用 tic、toc 計(jì)時(shí)方法，在重復(fù) 1000次運(yùn)算的情況下比較運(yùn)行時(shí)間。 的矩陣，并計(jì)算出其行列式、跡、秩、條件數(shù)、逆以及特征分解。 N=100, p=。 ，并觀察正態(tài)分布標(biāo)準(zhǔn)差的幾何表示。 mu=3。sigma=。 x=mu+sigma*[3:1,1:3]。yf=normcdf(x,mu,sigma)。P=[yf(4)yf(3),yf(5)yf(2),yf(6)yf(1)]。xd=1::5。yd=normpdf(xd,mu,sigma)。clf for k=1:3 % xx=x(4k):sigma/10:x(3+k)。yy=normpdf(xx,mu,sigma)。 %subplot(3,1,k),plot(xd,yd,39。b39。)。hold on fill([x(4k),xx,x(3+k)],[0,yy,0],39。g39。)。hold off if k / 17 text(,39。[{mu}{sigma},{mu}+{sigma}]39。)kk=int2str(k)。text(,[39。[{mu}39。,kk,39。{sigma},{mu}+39。,kk,39。{sigma}]39。])else end text(,num2str(P(k)))。shg end xlabel(39。x39。)。 ，實(shí)現(xiàn)給定數(shù)據(jù)組 x0 , y0，求擬合三階多項(xiàng)式，并圖示擬合情況。改變多項(xiàng)式階數(shù)再次圖示并觀察結(jié)果變化。 x0=0::1。y0=[,,,]。n=3。P=polyfit(x0,y0,n) xx=0::1。yy=polyval(P,xx)。plot(xx,yy,39。b39。,x0,y0,39。.r39。,39。MarkerSize39。,20)legend(39。擬合曲線39。,39。原始數(shù)據(jù) 39。,39。Location39。,39。SouthEast39。)xlabel(39。x39。) 【實(shí)驗(yàn)報(bào)告要求】 1．整理實(shí)驗(yàn)結(jié)果。 2．總結(jié)實(shí)驗(yàn)心得體會 / 17 實(shí)驗(yàn)五 MATLAB可視化 【實(shí)驗(yàn)?zāi)康摹? MATLAB連續(xù)函數(shù)（信號）的可視化 MATLAB二維離散作圖的基本格式 MTALAB三維作圖的基本格式和高級作圖修飾 MTALAB動畫的制作過程 【實(shí)驗(yàn)原理】 ，基于 plot 命令而展開的一系列作圖命令格式為學(xué)習(xí)本節(jié)的基礎(chǔ) 【實(shí)驗(yàn)儀器】 ； 。 【實(shí)驗(yàn)內(nèi)容及步驟】 。 t1=(0:11)/11*pi。t2=(0:400)/400*pi。t3=(0:50)/50*pi。y1=sin(t1).*sin(9*t1)。y2=sin(t2).*sin(9*t2)。y3=sin(t3).*sin(9*t3)。subplot(2,2,1),plot(t1,y1,39。r.39。) % axis([0,pi,1,1]),title(39。(1)點(diǎn)過少的離散圖形 39。)subplot(2,2,2),plot(t1,y1,t1,y1,39。r.39。)% axis([0,pi,1,1]),title(39。(2)點(diǎn)過少的連續(xù)圖形 39。)subplot(2,2,3),plot(t2,y2,39。r.39。) % axis([0,pi,1,1]),title(39。(3)點(diǎn)密集的離散圖形 39。)subplot(2,2,4),plot(t3,y3) % axis([0,pi,1,1]),title(39。(4)點(diǎn)足夠的連續(xù)圖形 39。) 。 clf t=(0:pi/50:2*pi)39。k=::1。Y=cos(t)*k。 subplot(1,2,1),plot(t,Y,39。LineWidth39。,) / 17 title(39。By plot(t,Y)39。),xlabel(39。t39。)subplot(1,2,2),plot(Y,39。LineWidth39。,)title(39。By plot(Y)39。),xlabel(39。row subscript of Y39。) 及其包絡(luò)線。（圖 ） t=(0:pi/100:pi)39。y1=sin(t)*[1,1]。t3=pi*(0:9)/9。 % % % % % y2=sin(t).*sin(9*t)。y3=sin(t3).*sin(9*t3)。hold on plot(t,y1,39。r:39。,t,y2,39。bo39。)% plot(t3,y3,39。s39。,39。MarkerSize39。,10,39。MarkerEdgeColor39。,[0,1,0],39。MarkerFaceColor39。,[1,0])axis([0,pi,1,1]) % % hold off % %以下指令供讀者比較用。使用時(shí)，指令前的 % 號要去除。 %屬性影響該指令中的所有線對象中的離散點(diǎn)。 %plot(t,y1,39。r:39。,t,y2,39。bo39。,t3,y3,39。s39。,39。MarkerSize39。,10,39。MarkerEdgeColor39。,[0,1,0],39。MarkerFaceColor39。,[1,0]) 。演示采用長軸為 ，短軸為 。 t=0:2*pi/99:2*pi。x=*cos(t)。y=*sin(t)。title(39。Normal and Grid on39。)subplot(2,3,2),plot(x,y),axis equal,grid on,title(39。Equal39。)subplot(2,3,3),plot(x,y),axis square,grid on,title(39。Square39。)subplot(2,3,4),plot(x,y),axis image,box off,title(39。Image and Box off39。)subplot(2,3,5),plot(x,y),axis image fill,box off title(39。Image and Fill39。)subplot(2,3,6),plot(x,y),axis tight,box off,title(39。Tight39。) subplot(2,3,1),plot(x,y),axis normal,grid on, ，繪制二階系統(tǒng)階躍響應(yīng)，綜合演示圖形標(biāo)識。 clf。t=6*pi*(0:100)/100。y=1exp(*t).*cos(*t)。plot(t,y,39。r39。,39。LineWidth39。,3)hold on tt=t(find(abs(y1)))。ts=max(tt)。plot(ts,39。bo39。,39。MarkerSize39。,10)hold off axis([inf,6*pi,inf]) set(gca,39。Xtick39。,[2*pi,4*pi,6*pi],39。Ytick39。,[,1,max(y)])% / 17 % % % set(gca,39。XtickLabel39。,{39。2*pi39。39。4*pi39。39。6*pi39。})grid on text(,39。fontsize{12}{alpha}=39。)text(,39。fontsize{12}{omega}=39。) cell_string{1}=39。fontsize{12}uparrow39。% % % set(gca,39。YtickLabel39。,{39。39。39。139。39。39。39。max(y)39。}) % % cell_string{2}=39。fontsize{16} fontname{隸書 }鎮(zhèn)定時(shí)間 39。cell_string{3}=39。fontsize{6} 39。 % cell_string{4}=[39。fontsize{14}rmt_{s} = 39。 num2str(ts)]。text(ts,cell_string,39。Color39。,39。b39。,39。HorizontalAlignment39。,39。Center39。) title(39。fontsize{14}it y = 1e^{alpha t}cos{omegat}39。)xlabel(39。fontsize{14} bft rightarrow39。)ylabel(39。fontsize{14} bfy rightarrow39。) hold繪制離散信號通過零階保持器后產(chǎn)生的波形，驗(yàn)證疊繪命令 t=2*pi*(0:20)/20。y=cos(t).*exp(*t)。stem(t,y,39。g39。,39。Color39。,39。k39。)。hold on stairs(t,y,39。:r39。,39。LineWidth39。,3)hold off legend(39。fontsize{14}it stem39。,39。fontsize{14}it stairs39。)box on ，畫出函數(shù) 和積分驗(yàn)證雙縱坐標(biāo)。 clf。dx=。x=0:dx:4。y=x.*sin(x)。s=cumtrapz(y)*dx。 % % % % % a=plotyy(x,y,x,s,39。stem39。,39。plot39。)。 在區(qū)間 [0,4]上的曲線，text(,39。fontsize{14}ity=xsinx39。)sint=39。{fontsize{16}int_{fontsize{8}0}^{ x}}39。text(,3.5,ss)% ss=[39。fontsize{14}its=39。,sint,39。fontsize{14}itxsinxdx39。]。 set(get(a(1),39。Ylabel39。),39。String39。,39。被積函數(shù) ity=xsinx39。)% / 17 set(get(a(2),39。Ylabel39。),39。String39。,ss)% xlabel(39。x39。) ，演示 subplot指令對圖形窗的分割。 clf。t=(pi*(0:1000)/1000)39。y1=sin(t)。y2=sin(10*t)。y12=sin(t).*sin(10*t)。subplot(2,2,1),plot(t,y1)。axis([0,pi,1,1])subplot(2,2,2),plot(t,y2)。axis([0,pi,1,1])subplot(39。position39。,[,])plot(t,y12,39。b39。,t,[y1,y1],39。r:39。)axis([0,pi,1,1]) ，實(shí)現(xiàn)三維曲線繪圖。演示：三維曲線的參數(shù)方程；線型、點(diǎn)形和圖例。（圖 ） t=(0::2)*pi。x=sin(t)。y=cos(t)。z=cos(2*t)。plot3(x,y,z,39。b39。,x,y,z,39。bd39。)view([82,58]),box on xlabel(39。x39。),ylabel(39。y39。),zlabel(39。z39。)legend(39。鏈 39。,39。寶石 39。) ，體會 不同燈光、照明、材質(zhì)指令所表現(xiàn)的圖形。 clf。[X,Y,Z]=sphere(40)。colormap(jet) % % % % % % % subplot(1,2,1),surf(X,Y,Z),axis equal off,shading interp light(39。position39。,[010 ],39。style39。,39。infinite39。)lighting phong material shiny subplot(1,2,2),surf(X,Y,Z,Z),axis equal off,shading flat light。lighting flat % light(39。position39。,[1,1,2],39。color39。,39。y39。)% light(39。position39。,[1,1],39。style39。,39。local39。,39。color39。,39。w39。)set(gcf,39。Color39。,39。w39。) 11．透視圖。 / 17