【文章內(nèi)容簡介】 :pi/100:pi。y1=x.^2。y2=cos(2*x)。y3=y1.*y2。plot(x,y1,39。b39。,x,y2,39。r:39。,x,y3,39。k39。)運(yùn)行結(jié)果：（2）M文件：clc。x=pi:pi/100:pi。y1=x.^2。y2=cos(2*x)。y3=y1.*y2。subplot(1,3,1)。plot(x,y1,39。b39。)。title(39。y1=x^239。)。subplot(1,3,2)。plot(x,y2,39。r:39。)。title(39。y2=cos(2x)39。)。subplot(1,3,3)。plot(x,y3,39。k39。)。title(39。y3=y1*y239。)。.運(yùn)行結(jié)果：（3）M文件：clc。x=pi:pi/100:pi。y1=x.^2。y2=cos(2*x)。y3=y1.*y2。subplot(2,2,1)。plot(x,y1,39。b39。,x,y2,39。r:39。,x,y3,39。k39。)。subplot(2,2,2)。bar(x,y1,39。b39。)。title(39。y1=x^239。)。subplot(2,2,3)。bar(x,y2,39。r39。)。 title(39。y2=cos(2x)39。)。subplot(2,2,4)。bar(x,y3,39。k39。)。title(39。y3=y1*y239。)。由上面的M文件，只要依次將“bar”改為“stairs”、“stem”、“fill”,再適當(dāng)更改區(qū)間取的點(diǎn)數(shù)，運(yùn)行程序即可，即有下面的結(jié)果：3. 已知在5≤x≤5區(qū)間繪制函數(shù)曲線。解：M文件：clc。x=5::5。y=(x+sqrt(pi))/(exp(2)).*(x=0)+*log(x+sqrt(1+x.^2)).*(x0)。plot(x,y)運(yùn)行結(jié)果：由圖可看出，函數(shù)在零點(diǎn)不連續(xù)。4. 繪制極坐標(biāo)曲線ρ=asin(b+nθ)，并分析參數(shù)a、b、n對曲線形狀的影響。解：M文件如下：clc。theta=0:pi/100:2*pi。a=input(39。輸入a=39。)。b=input(39。輸入b=39。)。n=input(39。輸入n=39。)。rho=a*sin(b+n*theta)。polar(theta,rho,39。m39。)采用控制變量法的辦法，固定兩個參數(shù)，變動第三個參數(shù)觀察輸出圖象的變化。分析結(jié)果：由這8個圖知道， 當(dāng)a,n固定時，圖形的形狀也就固定了，b只影響圖形的旋轉(zhuǎn)的角度。當(dāng)a,b固定時，n只影響圖形的扇形數(shù)，特別地，當(dāng)n是奇數(shù)時，扇葉數(shù)就是n,當(dāng)是偶數(shù)時，扇葉數(shù)則是2n個。當(dāng)b,n固定時，a影響的是圖形大小，特別地，當(dāng)a是整數(shù)時，圖形半徑大小就是a。5. 繪制函數(shù)的曲線圖和等高線。其中x的21個值均勻分布[5，5]范圍，y的31個值均勻分布在[0，10]，要求使用subplot(2,1,1)和subplot(2,1,2)將產(chǎn)生的曲面圖和等高線圖畫在同一個窗口上。解：M文件：clc。x=linspace(5,5,21)。y=linspace(0,10,31)。[x,y]=meshgrid(x,y)。z=cos(x).*cos(y).*exp(sqrt(x.^2+y.^2)/4)。subplot(2,1,1)。surf(x,y,z)。title(39。曲面圖39。)。subplot(2,1,2)。surfc(x,y,z)。title(39。等高線圖39。)。運(yùn)行結(jié)果：6. 繪制曲面圖形，并進(jìn)行插值著色處理。解：M文件：clc。s=0:pi/100:pi/2。t=0:pi/100:3*pi/2。[s,t]=meshgrid(s,t)。x=cos(s).*cos(t)。y=cos(s).*sin(t)。z=sin(s)。subplot(2,2,1)。mesh(x,y,z)。title(39。未著色的圖形39。)。subplot(2,2,2)。surf(x,y,z)。title(39。shading faceted（缺?。?9。)。subplot(2,2,3)。surf(x,y,z)。shading flat。title(39。shading flat39。)。subplot(2,2,4)。surf(x,y,z)。shading interp。title(39。shading interp39。)。 運(yùn)行結(jié)果有：實驗七 低層繪圖操作二、實驗內(nèi)容1. 建立一個圖形窗口，使之背景顏色為紅色，并在窗口上保留原有的菜單項，而且在按下鼠標(biāo)器的左鍵之后顯示出Left Button Pressed字樣。解：M文件如下：clc。hf=figure(39。color39。,[1 0 0],... 39。WindowButtonDownF39。,39。disp(39。39。Left Button Pressed.39。39。)39。)。 運(yùn)行結(jié)果：左擊鼠標(biāo)后：2. 先利用默認(rèn)屬性繪制曲線y=x2e2x，然后通過圖形句柄操作來改變曲線的顏色、線型和線寬，并利用文件對象給曲線添加文字標(biāo)注。解：M文件：clc。x=2::2。y=x.^2.*exp(2*x)。h=plot(x,y)。set(h,39。color39。,[,],39。linestyle39。,39。39。,... 39。linewidth39。,2)。text(,^2*exp(2*),39。\leftarrow x^2exp(2x)39。,39。fontsize39。,9)。運(yùn)行結(jié)果:3. 利用曲面對象繪制曲面v(x,t)=(+π)。解：M文件：clc。x=0::2*pi。[x,t]=meshgrid(x)。v=10*exp(*x).*sin(2000*pi**x+pi)。axes(39。view39。,[37,30])。hs=surface(x,t,v,39。facecolor39。,... [,],39。edgecolor39。,39。flat39。)。grid on。xlabel(39。xaxis39。)。 ylabel(39。yaxis39。)。zlabel(39。zaxis39。)。title(39。meshsurf39。)。pause %按任意鍵繼續(xù)set(hs,39。FaceColor39。,39。flat39。)。text(0,0,0,39。曲面39。)。運(yùn)行結(jié)果：按任意鍵繼續(xù)：4. 以任意位置子圖形式繪制出正弦、余弦、正切和余切函數(shù)曲線。5. 生成一個圓柱體，并進(jìn)行光照和材質(zhì)處理。 解：M文件：[x,y,z]=cylinder(3,500)。 %cylinder是生成柱體的函數(shù)surf(x,y,z)。title(39。圓柱體的光照和材料處理39。)。Xlabel(39。Xaxis39。)。Ylabel(39。Yaxis39。)。Zlabel(39。Zaxis39。)。axis([5,5,5,5,0,1])grid off。light(39。Color39。,39。r39。,39。Position39。,[4,0,0],39。style39。,39。infinite39。)。shading interp。material shiny。view(0,10)。lighting phong。axis off。 運(yùn)行結(jié)果：實驗八 數(shù)據(jù)處理與多項式計算一、實驗?zāi)康?. 掌握數(shù)據(jù)統(tǒng)計和分析的方法。2. 掌握數(shù)值插值與曲線擬合的方法及其應(yīng)用。3. 掌握多項式的常用運(yùn)算。二、實驗內(nèi)容1. 利用MATLAB提供的rand函數(shù)生成30000個符合均勻分布的隨機(jī)數(shù)，然后檢驗隨機(jī)數(shù)的性質(zhì)：(1) 均值和標(biāo)準(zhǔn)方差。(2) 最大元素和最小元素。(3) 。解：M文件:clc。x=rand(1,30000)。mu=mean(x) %求這30000個均勻分布隨機(jī)數(shù)的平均值sig=std(x) %求其標(biāo)準(zhǔn)差σ1y=length(find(x))。 %p=y/30000 %運(yùn)行結(jié)果：mu = sig = p = 2. 將100個學(xué)生5門功課的成績存入矩陣P中，進(jìn)行如下處理：(1) 分別求每門課的最高分、最低分及相應(yīng)學(xué)生序號。(2) 分別求每門課的平均分和標(biāo)準(zhǔn)方差。(3) 5門課總分的最高分、最低分及相應(yīng)學(xué)生序號。(4) 將5門課總分按從大到小順序存入zcj中，相應(yīng)學(xué)生序號存入xsxh。提示：上機(jī)調(diào)試時，為避免輸入學(xué)生成績的麻煩，可用取值范圍在[45,95]之間的隨機(jī)矩陣來表示學(xué)生成績。解：M文件：clc。t=45+50*rand(100,5)。P=fix(t)。 %生成100個學(xué)生5門功課成績[x,l]=max(P) %x為每門課最高分行向量,l為相應(yīng)學(xué)生序號[y,k]=min(P)%y為每門課最低分行向列,k為相應(yīng)學(xué)生序號mu=mean(P) %每門課的平均值行向量sig=std(P) %每門課的標(biāo)準(zhǔn)差行向量s=sum(P,2) %5門課總分的列向量[X,m]=max(s)%5門課總分的最高分X與相應(yīng)學(xué)生序號m[Y,n]=min(s)%5門課總分的最低分Y與相應(yīng)學(xué)生序號n[zcj,xsxh]=sort(s) %zcj為5門課總分從大到小排序，相應(yīng)學(xué)生序號xsxh 運(yùn)行結(jié)果：3. 某氣象觀測得某日6:00~18:00之間每隔2h的室內(nèi)外溫度（0C）如實驗表1所示。實驗表1 室內(nèi)外溫度觀測結(jié)果（0C）時間h 6 8 10 12 14 16 18室內(nèi)溫度t1 室外溫度t2 試用三次樣條插值分別求出該日室內(nèi)外6:30~18:30之間每隔2h各點(diǎn)的近似溫度（0C）。解：M文件：clc。h=6:2:18。t1=[ ]。t2=[ ]。T1=interp1(h,t1,39。spline39。)%室內(nèi)的3次樣條插值溫度T2=interp1(h,t2,39。spline39。)%室外的3次樣條插值溫度 運(yùn)行結(jié)果：T1 = Columns 1 through 3 Columns 4 through 6 Column 7 T2 = Columns 1 through 3 Columns 4 through 6 Column 7 4. 已知lgx在[1,101]區(qū)間10個整數(shù)采樣點(diǎn)的函數(shù)值如實驗表2所示。實驗表2 lgx在10個采樣點(diǎn)的函數(shù)值x 1 11 21 31 41 51 61 71 81 91 101lgx 0 試求lgx的5次擬合多項式p(x)，并繪制出lgx和p(x)在[1,101]區(qū)間的函數(shù)曲線。解：M文件：x=1:10:101。y=lg10(x)。P=polyfit(x,y,5)y1=polyval(P,x)。plot(x,y,39。:o39。,x,y1,39。*39。) 運(yùn)行結(jié)果：Warning: Polynomial is badly conditioned. Add points with distinct X values, reduce the degree of the polynomial, or try centering and scaling as described in HELP POLYFIT. In polyfit at 80P = （這里出現(xiàn)警告是提示不必用5價函數(shù)就已經(jīng)可以完美擬合了，是可以降價擬合。）在[1,101]的區(qū)間函數(shù)圖像5. 有3個多項式P1(x)=x4+2x3+4x2+5，P2(x)=x+2，