中國(guó)海洋大學(xué)MATLAB期末復(fù)習(xí)2014_答案

1、一、計(jì)算題：1. 設(shè)，計(jì)算：(1) ；(2) ；(3) 參考答案： u=2; v=3;（1） 4*u*v/log(v)ans = 21.8457（2） (exp(u)+v)2/(v2-u)ans = 15.4189（3） sqrt(u-3*v)/(u*v)ans = 0 + 0.4410i2. 計(jì)算cos 60- cos(pi/3)-(9-2(1/2)(1/3)ans = -1.46493. 計(jì)算下列積分。(1) (2) (3) 參考答案：（1） f = (x)x+x.3+x.5; q = quad(f,-1,1)q = 2.2204e-016(2) f5 = (x,y)sin(y).*(x+

2、y)./(x.2+4); q = dblquad(f5,1,10,1,10)q = 5.5254(3) f = inline(sin(2*(x+pi/3)+cos(x+pi/6),x); Q = quad(f,0,2*pi)Q = 6.7851e-0094. 創(chuàng)建一個(gè)含10000個(gè)元素的數(shù)組，其值依次為1到10000之間的整數(shù)。用邏輯數(shù)組計(jì)算大于5000的元素的平方根。a = 1:10000; %Declare array ab = a 5000; %Create maska(b) = sqrt(a(b); %Take square root5. 對(duì)進(jìn)行因式分解參考答案： factor(sym

3、(3*a2*(x-y)3 - 4*b2*(y-x)2) ans = (x-y)2*(3*a2*x-4*b2-3*a2*y)6. 求解線性方程組 A=3 4 -7 -12; 5 -7 4 2; 1 0 8 -5; -6 5 -2 10A = 3 4 -7 -12 5 -7 4 2 1 0 8 -5 -6 5 -2 10 B=4;-3;9;-8B = 4 -3 9 -8 x=ABx = -1.4841 -0.6816 0.5337 -1.24297. 計(jì)算下列各式（1）（2），求（3），求，（4），參考答案：（1） limit(sym(tan(x) - sin(x)/(1-cos(2*x) ans

4、 = 0（2） y = sym(x3 - 2*x2 + sin(x); diff(y) ans = 3*x2-4*x+cos(x)（3） f = x*y*log(x+y); fx = diff(f,x) fx = y*log(x+y)+x*y/(x+y) fy = diff(f,y) fy = x*log(x+y)+x*y/(x+y) f2xy = diff(fx,y) f2xy = log(x+y)+y/(x+y)+x/(x+y)-x*y/(x+y)2（4） syms t y = log(1+t); int(y) ans = log(1+t)*(1+t)-t-1 int(y,0,27) an

5、s = 56*log(2)+28*log(7)-278. 已知變量：A=ilovematlab;B=matlab，請(qǐng)找出：（1）B在A中的位置；（2）把B放在A后面。解：Lb=strfind(A,B)Lb= 6Le=Lb+length(B)-1Le= 119. 煉鋼基本上是一個(gè)氧化脫碳的過程，鋼液中原含碳量多少直接影響到冶煉時(shí)間的長(zhǎng)短，下表是某平爐的熔鋼完畢碳（x）與精煉時(shí)間（y）的生產(chǎn)記錄：X（0.01%）134150180104190163200Y（min）135168200100215175220現(xiàn)希望從上表的數(shù)據(jù)中找出x與y變化規(guī)律的經(jīng)驗(yàn)公式，用多項(xiàng)式進(jìn)行曲線擬合，并給出相應(yīng)的曲線。x

6、=104 134 150 163 180 190 200; y=100 135 168 175 200 215 220; f=polyfit(x,y,2)f = -0.0016 1.7666 -68.3091 plot(x,y,o,x,polyval(f,x)x與y的關(guān)系式為：y=-0.0016*x2+1.7666*x-68.309110. 有如下數(shù)據(jù)：X11.11.21.31.4y1.000001.233681.552711.993722.61170利用本章介紹的幾種插值方法對(duì)其進(jìn)行插值，得到每隔0.05的結(jié)果。編寫腳本文件，文件內(nèi)容為：% Interpolation using the f

7、our methodsx=1 1.1 1.2 1.3 1.4;y=1.00000 1.23368 1.55271 1.99372 2.61170;length_of_x=length(x);scalar_x=x(1):0.05:x(length_of_x);length_of_sx=length(scalar_x);y_nearest = zeros(length(scalar_x),1);y_linear = zeros(length(scalar_x),1);y_spline = zeros(length(scalar_x),1);y_cubic = zeros(length(scalar

8、_x),1);for i=1:length_of_sx y_nearest(i)=interp1(x,y,scalar_x(i),nearest); y_linear(i) =interp1(x,y,scalar_x(i),linear); y_spline(i) =interp1(x,y,scalar_x(i),spline); y_cubic(i) =interp1(x,y,scalar_x(i),cubic);endsubplot(2,2,1),plot(x,y,*),hold on,plot(scalar_x,y_nearest),title(method=nearest);subpl

9、ot(2,2,2),plot(x,y,*),hold on,plot(scalar_x,y_linear),title(method=linear);subplot(2,2,3),plot(x,y,*),hold on,plot(scalar_x,y_spline),title(method=spline);subplot(2,2,4),plot(x,y,*),hold on,plot(scalar_x,y_cubic),title(method=cubic);得到結(jié)果為： 二、繪圖題：1. 繪制雙曲拋物面：，參考答案： X,Y = meshgrid(-16:0.4:16,-4:0.1:4);

10、 Z = X.2/16 - Y.2/4; surf(X,Y,Z)2. 繪制函數(shù)在，上的表面圖參考答案： syms x y z = 1/(2*pi)*exp(-(x2+y2); ezsurf(x,y,z,-3,3,-3,3);3. 求函數(shù)cos(x)的5次擬合多項(xiàng)式p(x),并繪制出cos(x)和p(x)在0,2pi區(qū)間的函數(shù)曲線。4. 編寫MATALAB程序，完成下列任務(wù)：（1）在區(qū)間0,4*pi上均勻地取100個(gè)點(diǎn)構(gòu)成向量；（2）分別計(jì)算函數(shù)y1=sin(t)與y2=2cos(2t)在向量t處于區(qū)間0,4*pi上的函數(shù)值；（3） 在同一圖形窗口繪制曲線y1=sin(t)與y2=2cos(2t

11、)，要求y1曲線為黑色點(diǎn)畫線，y2曲線為紅色虛線圓圈；并在圖中恰當(dāng)位置標(biāo)注兩條曲線的圖例；給圖形加上標(biāo)題“y1 and y2”。參考答案：(1) t=linspace(0,4*pi,100);(2) y1=sin(t); y2=2*cos(2*t);(3) plot(t,y1,k-.); text(1.2,sin(1.2)-0.1,y1rightarrow,Fontsize,12); hold on; plot(t,y2,r-o); text(4.8,1.5,y2rightarrow,Fontsize,12); title(y1 and y2);5. 編寫程序，實(shí)現(xiàn)功能為：創(chuàng)建圖形窗口，并且設(shè)

12、置其默認(rèn)背景為黃色，默認(rèn)線寬為 4 個(gè)像素，在該窗口中繪制橢圓 的圖像，其中的和任選。參考答案：figure(Color,y);set(gca,DefaultLineLineWidth,4);a = 4; b = 3;x = linspace(-a,a,100);y1 = sqrt(1-x.2/(a2)*b2);y2 = -sqrt(1-x.2/(a2)*b2);h1=line(x,y1);hold on;h2=line(x,y2);set(h1,LineWidth,default);set(h2,LineWidth,default);三、編程題： 1. 編寫一個(gè)腳本，判斷輸入字符串中每個(gè)單詞

13、的首字母是否為大寫，若不是則將其修改為大寫，其他字母為小寫。參考答案：腳本文件內(nèi)容為：str = this is the string to be converted; nlength = length(str);for k=1:nlength if (k=1 | isspace(str(k-1) & (str(k)=a) str(k) = char(double(str(k) - 32); endenddisp(str);2. 在MATLAB中使用一個(gè)循環(huán)確定：如果用戶最初在一個(gè)銀行帳戶中存儲(chǔ)$10000，并且在每年的年終再存儲(chǔ)$10000(銀行每年支付6%的利息)，那么賬戶上要積累$100

14、0000要需要多長(zhǎng)時(shí)間。參考答案： 33年。n=0;z=10000;y=1+0.01*6;while z=1000000z=z*y+10000;n=n+1;enddisp(n)3. 編寫matlab語(yǔ)句，計(jì)算的值已知t從-9到9每隔0.5取一次值。運(yùn)用循環(huán)和選擇語(yǔ)句進(jìn)行計(jì)算。4. 編制程序，從鍵盤輸入一系列的測(cè)量數(shù)，計(jì)算它們的平均值以及標(biāo)準(zhǔn)差。這些數(shù)可以是正數(shù)，負(fù)數(shù)或0。要求給出輸入值的個(gè)數(shù)，然后用for 循環(huán)讀取所有數(shù)值。% Initialize sumssum_x = 0; sum_x2 = 0;% Get the number of points to input.n = input(E

15、nter number of points: );% Check to see if we have enough input data.if n 2 % Insufficient datadisp (At least 2 values must be entered.);else % we will have enough data, so lets get it.% Loop to read input values.for ii = 1:n% Read in next valuex = input(Enter value: );% Accumulate sums.sum_x = sum_x + x;sum_x2 = sum_x2 + x2;end% Now calculate s