1、習(xí)題 l1 執(zhí)行下列指令， 觀 察其運(yùn)算 結(jié)果， 理 解 其 意 義 ：(I)I2;34 + 10 - 2 i(2)I2;3 4 .* 0 . 1 0.2; 0.3 0.4(3) (1 2; 3 4.(20 10;9 2(4) l 2; 3 4.2(5) exp(l 2; 3 4)(6)1og(lIO 100)(7)prod(l 2;3 4)(8)a,b=min( IO 20;30 40)(9)abs(l 2;3 4-pi)(10) 1 2;3 4=4,3;2 1(ll)find(lO 20;30 40=40,30;20 10)(12) a,b=find(IO 20;30 40=40,30;2

2、0 10)(提示： a 為行號(hào)， b 為列號(hào)）(13) all(l 2;3 41)(14) any(l 2;3 41)(15) linspace(3,4,5)生成線性等距向蜇(16) A= l 2;3 4;A(:,2)%Page20,ex I(5)等于exp( I),e xp(2);exp( 3),e xp( 4) (7) 3= I*3,8 = 2* 4(8)a 為各列最小伯， b 為品小伯所在的行號(hào)(10) I =4, fa lse , 2 = 3,fa lse , 3 =2 , ture, 4=I,cure(11) 答案表明： 編 劃 笫2元 余滿足不等式(30=20)和編劉第 4 元系滿

3、足不等式(40=10)(12) 答案表明：編 址 第 2 行第 1 列兀索滿足不等式(30=20)和編址第 2 行第 2 列兒索滿足不等式(40=10)2 執(zhí)行下列指令，觀 察 其運(yùn) 絆 結(jié)果、變恁類型和字節(jié)數(shù)， 理解 其意義：(1) clear; a=l ,b=num2s tr(a),c= a0 , a= =b, a= =c, b= =c(2) clear; fun=abs(x) ,x= -2,eval(fun),do uble(fun)%Page20, ex2(1) a = 1b =1c = Ians =0ans =1ans = 0a, b, c 的們盡管都是 1,但數(shù)據(jù)類別分別為數(shù)仇，字

4、符，邏輯，注 意 a 與 c 相等，但 他們不 等 千 b。(2) double(fun)輸出的分別是字符 a,b,s,(,x,)的 ASCII 碼3 本金 K 以每年 n 次，每 次 p % 的增值 率(n 與 p 的乘積 為每年增值額的百分比）增加，當(dāng)增加到 rK 時(shí) 所花 費(fèi)的 時(shí)間 為ln rT=（單位： 年 ）nln(l + O.Ol p )%Page20 ,ex3 r=2;p=0.5;n= l 2;T =log( r)/n/lo g( 1+0.0 1* p ) T =11.5 8 13用 M 燈L AB表 達(dá) 式 寫(xiě) 出 該公 式 并用下列數(shù)據(jù)計(jì)算： r=2, p=0.5, n=1

5、24.已知函數(shù)f(x)=x4 - 2x在(-2, 2)內(nèi)有兩個(gè)根。取步長(zhǎng) h=0.05,通過(guò)計(jì)算函數(shù)值求得函數(shù)的錄小值點(diǎn)和兩個(gè)根的近似解。（提示： 求 近似根 等價(jià)于求函數(shù)絕對(duì)值的最小值點(diǎn)）%Page20 ,ex4 x=-2:0 .05:2;f =x.A4 - 2 . Ax; ff min,min_indexl=min(t) fmin =-1.3907min_i ndex = 54 x(min_index)ans =尿小伯品小伯點(diǎn)編址0.6500% 最小伯點(diǎn) fl ,x l _inde x=min(abs(t) )% 求近似根 絕 對(duì) 但 最 小 的 點(diǎn)f1=0 . 03 28xJ _ ind

6、ex = 24 x(x l_ index)ans =-0.8500 x(x l _index)=fl;f = x.A4 - 2 . Ax;% 刪去絕對(duì)值展小的點(diǎn)以求函數(shù)絕對(duì)伯次小的點(diǎn) r僅 ，x2_indexl=min(abs(t)% 求另一近似根函 數(shù) 絕 對(duì) 值 次小的點(diǎn)t2 =0.0630x2_index = 65 x(x2_inde x)ans =1.25006 5.(1)用 z = magic( I0)得到 10 階魔方矩陣；(2) 求z 的 各 列 元素之和；(3) 求z的 對(duì) 角 線 元 素 之 和（提示 先 用 d i ag (z) 提取 z 的對(duì)角線）；(4) 將 z 的第二

7、列除以 3;(5) 將 z 的 第 3 行 元 素 加 到 第 8 行。6 6先 不 用 M ATLAB 判斷下面語(yǔ)句將顯示什么結(jié)果? size (B)又得出什么結(jié)果？ B1=1:9; David Beckham ;B2= 180:- 10:100; 100,80,75,;77,60,92;67 28 90;100 89 78; B=Bl , B2;B1,2(8)D=cell2s truct(B,f fl ,t2,2); a,b=D.fl然后用 MATL AB 驗(yàn)證你的判斷。進(jìn)一步， 察看變量類 型和字 節(jié)數(shù)， 并用 W or kspace 工具欄顯示 B 和 D 的具體內(nèi)容。習(xí)題 2I. 設(shè)

8、 x 為一個(gè)長(zhǎng)度為 n 的數(shù)組， 編程求下列均值和標(biāo)準(zhǔn)差艾 =-I Ix;,n i=IIIIs = . /- I;x/- n滬 ,nln - Ji=I%Page 40 e x l先在 編 鉛 器窗 口與下列 M 函數(shù) ， 保存為 eg2_ l.m扣nctio n xbar,s=ex2_ I( x) n=length(x);xba1-=su m(x)/n;s=sq rt(sum(x. 2)-n* xbar2)/(n-l );例如x= l81 70 65 51 76 66 90 87 61 77J;xbar,s=ex2_ l (x) xbar =72.4000S=12.11242.求滿足 f1o(

9、l + n) l OO 的 最 小 m 值 。,r=O% Pa g e 40 ex2 s=lo g ( I );n =O; while sek=k+ I;F( k) =F ( k- l) +F( k-2) ; x=F ( k)/F ( k- I) ;end a,x,k計(jì)算至 k=21 可滿足精度4. 分別用 for 和 while 循環(huán)結(jié)構(gòu)編寫(xiě)程序， 求出 K =106322 o并考慮一種避免循環(huán)語(yǔ)句的程序設(shè)計(jì)，比 較不同 算法的運(yùn)行時(shí)間。%Page 40 ex4 clear;tic;s=O; for i=l:1000000s=s+sqrt(3)/2i;end s,toctic;s=O;i=

10、l;while i= IO OOOOOs =s+sqrt(3)/2 i;i=i+ I ;i =Ie nd s,toc tic;s=O;i=l:1000000;s=sqrt(3)木sum( l./ 2. i); s,toc。5. 假定某天的氣溫變化記錄如下表，試 作圖 描述這一天的氣溫變化規(guī)律。時(shí)刻 t(h)I23456789101112溫度OC(t)15141414141516182022232528時(shí)刻 t(h)131415161718192021222324溫度 OC( t)313231292725242220181716%Page 40 ex5 t=0:24;c=l5 14 14 14

11、14 15 16 18 20 22 23 25 28 .31 32 31 29 27 25 24 22 20 18 17 161;plot(t,c)6. 作出下列函數(shù)圖象(i) 曲線 y = x2sin (i - x - 2), -2 f 心2(要求分別使用 plot 或 fplot 憲成）%(1)x=-2:0.1:2;y=x.A2 .*si n(x.A2 - x - 2) ; p lot(x,y)y=inline(xA2*s i n(xA2 - x-2);fplot(y,L- 2 2)(ii)橢圓 x2 ! 4+ /19= 1%(2) 參數(shù)方法t=linspace(0,2*pi,100);

12、x=2*cos(t); y=3*s in(t); plot(x,y)(iii)拋物而 z = x2 + y2, lxl 3, I計(jì):%(3)X=- 3: 0. 1 :3;y = x;x,y=meshgrid(x,y); Z=X. /2+ y. A2;S U頃x,y,z)(iv)曲面 z=x4 +3x 2 + /- 2x - 2y - 2x ,%(4)x=-3:0. 1: 3 ;y= - 3:0 . l : 1 3;x,y= meshg rid (x,y) ;3+6 , lxl3, -3y l3z =x. 4+Y x.2+y.2-2*x-2*y-2*x.2袞y+6; surf( x,y,z)(v

13、) 空間曲線 x=si nt, y=cost, z=cos (2t), 0t 2 兀%(5)t=0:0 .0 1:2* pi;x=sin(t);y=cos( t);z =cos (2* t) ; plo t3 (x ,y,z)(vi) 半球面 x=2s in t/Jcos0 , y =2si n汾si n 0 , z =2 co s 份，0在 360, 0綽華 90% ( 6 )th e ta =l i n s pace(0,2*pi,50);fai=linspace(O,p i/2 ,2 0) ; t heta,fai= meshgrid( thet a,fai) ;x =2 *si n(f

14、ai).*cos(theta);y=2*sin( fai).*s in(theta);z=2*cos(fai); surf(x,y,z)(vii) 三條曲線合成圖 y1 =si應(yīng)），2=si nxsin(l0 x), y3= - s inx , Oxl.1y =ixl x ll.l- 1.1x l .I )+x. (x=-1.1)-1.1*(x Ip (x, y) =0 .757 5 exp(- y2-6x 2)- I1);p=p+b*exp(-y.2-6 x.2). (x+y-1).*(x+y=1);p=p+a*exp(-0.75* y.2-3.75*x. 2+ l .5*x).*(x+y

15、A=4 1 - I;3 2 -6; 1 -5 3;b= 9;-2; 1 );rank(A), rank ( ans=A,bl)%fA,bl為增、丿矩 陣ans=33%可見(jiàn)方程組唯解 x=AbX=2.38301.48942.02134- 33x1 ( - 1(2)I 32- 6X2 I=I-2J- 53入3.丿J% E xercise 2(2) A=4 -3 3;3 2 -6; I -5 3l;b=r- I; -2; 11; r 叨 k(A), r 沺 k(A,bl) ans =ans =33% 可見(jiàn)方籽糾唯解 x=Ab x =。-0.4706-0.2941廠)4l 1(lX2(3)I 32I=

16、11l- 5丿U,%Exercise 2(3) A=4 1;3 2;1 -5l;b=f l ; l ; I); rank(A), rank(A,b) ans =ans =23% 可見(jiàn)方籽組無(wú)解 x=AbX =0.3311-0.1219 % 橄小二乘近似斛2l- 11 IIX i I(lX(4) 1l2l-1 1121 =1 2ll2l 丿1 X3 IuX4%Exercise 2(4)a=2, l ,-1,1;1,2,1,-1;1,1,2,1;b=1 2 3);%注怠 b 的寫(xiě)法 rank(a),rank(a,b) ans =3ans =3%rank(a)=ran k(fa , bl) ab a

17、ns =。0%個(gè)特解3求 第 2 題 第 (4) 小題的 通解。%Exercise 3a =2 , 1,-1 ,1; 1,2,1, -1;1,1,2,l ;b=1,2,3; x=null(a),xO= ab x=-0 .625 50.6255-0.20850.4170xO =。通 解 k x+ xO4. ( 人口流動(dòng)趨勢(shì)）對(duì)城鄉(xiāng)人口流動(dòng)作年度調(diào)查， 發(fā)現(xiàn) 有一個(gè)穩(wěn) 定的朝向城鎮(zhèn)流動(dòng)的趨勢(shì)， 每年農(nóng)村居民的 5%移居城鎮(zhèn)而城鎮(zhèn)居民的 1%遷出， 現(xiàn)在總?cè)丝诘?20 %位千城鎮(zhèn)。假如城鄉(xiāng)總?cè)丝诒３植蛔儯?且人口流動(dòng)的這種趨勢(shì)繼續(xù)下去，那 么( 1 )一年以后住在城鎮(zhèn)人口所占比例是多少？?jī)赡暌院竽?/p>

18、？ 十年以后呢？( 2 )很多年以后呢？( 3 )如果現(xiàn)在總?cè)丝?70%位于城鎮(zhèn)， 很多年以后城鎮(zhèn)人口所占比例是多少？( 4 )計(jì)尊轉(zhuǎn)移矩陣的最大特征值及對(duì)應(yīng)的特征向噩， 與問(wèn)題 (2) (3) 有何 關(guān)系？% E xercise 4 xO= f0 .2 0 .8l ; a = f0 .9 9 0 . 0 5 ;0.01 0.951; x I = a* x, x2=aA2 * x, x1O= aA1 O* x x=xO;fo r i= l : 10 00 ,x = a x;end,x x=0.83330.1667x0= 0.8 0.2); x=xO;fori= I : I0 00 ,x=a*x

19、;e nd, xX =0.83330.1667Lv,e=eig(a) v =0.9806-0.70710. 196 10.7071e =1.0000000.9400v(:,1)./x ans =1.17671.1767% 成比例， 說(shuō)明 x 是錄大特征伯對(duì)應(yīng)的特征向蕪5. ( 經(jīng)濟(jì)預(yù)測(cè)）在某經(jīng)濟(jì)年度內(nèi)， 各經(jīng)濟(jì)部門(mén)的投入產(chǎn)出表如下表 3.5 ( 單位： 億 元 ）消耗部門(mén)最后需求總產(chǎn)值工 業(yè)農(nóng) 業(yè)第三產(chǎn)業(yè)生產(chǎn)部門(mén)工業(yè)62j1625衣業(yè)2.2510.21.555第三產(chǎn)業(yè)30.21.81520假設(shè)某 經(jīng)濟(jì)年度工業(yè)，農(nóng) 業(yè)及第三產(chǎn)業(yè)的最后需求均為 17 億元，預(yù)測(cè) 該經(jīng)濟(jì)年度工業(yè)，衣 業(yè)及第三產(chǎn)業(yè)的

20、產(chǎn)出（提示： 對(duì)千一個(gè)特定的經(jīng)濟(jì)系統(tǒng)而言，直接消耗矩陣和 Leontief 矩陣可視作不變）。%Exercise 5用 到公式(3.11)(3.12)B=L6 , 2 , 1 ; 2.25, l , 0.2; 3,0. 2 , l. 8 ; x = 255 20 ;C=B/diag(x) C=0.24000.40000.05000.09000.20000.01000.12000.04000.0900A=eye(3,3)-C A =0.7600-0.0900-0.40000.8000-0.0500-0.0 JOO-0 . I200-0.04000.9 lOOD = f l 7 17 17 l;x

21、=ADX =37.569625.786224.76906. 求下列矩陣的行列式、逆、特征值和特征向晝.1 5765 7108741-n(ll-1 C I )I 32-6I(2)I O2-1 1 (3)I! %Exercise 6(1)1-53 丿l, 68109-120, 579IOa=4 1 -1;3 2 -6; 1 -5 3;de t(a),inv( a),v,d=eig(a) ans =-940.2553-0.02130.04260.1596-0.13 83-0.22340 .1809-0.2234-0.0 5320.0 1 85- 0 .9009-0 . 3066- 0 .7693-

22、0.1240-0 .72 4 8- 0 . 6 3 8 6- 0 .4 1580 . 6 170ans =V =3 . 676 0。d=- 3 . 0527。8.3766%Exercise 6(2)a=Lll- 1;02- 1; - l2 O;de t(a),inv(a),Lv , d = e ig (a) ans=1ans =2.0000-2.00001.00001.0000-1.00001.00002.0000-3.00002.0000v =-0.5773-0.5773-0.57740.5774 + O.OOOOi0.5774 - O.OOOOi0.57740.57740.5773 - O

23、.OOOOi0.5773+ O.OOOOid =1.0000。1 .0 00 0 + O.OOOOi。1.0000 - O.OOOOi%Exercise 6(3)A=5 7 6 5;7 10 8 7; 6 8 10 9;5 7 9 10 A =5765710876810957910de t(A),inv( A), v,d=eig( A) ans=Ians =68.0000-41.0000-17.000010.0000-41.000025.000010.0000-6.0000-17.000010.00005.0000-3.000010.0000-6.0000-3.00002.0000v =0.8

24、3040.09330.39630.3803-0.501 6- 0 . 30 170.6 1490 . 5286- 0 . 208 60 .7 6 03-0 .27 160 . 55200.1237-0.5676-0.62540.52090 .84 31。d =0.0102。3 .8 58 1。30 .28 8756156(4)n階方陣I1s, n 分別為5, 50,和 5006 5%Exercise 6(4)(以 n=5 為例）關(guān)鍵是矩陣的定義方法一（三個(gè) fo r) n=5;for i= l :n, a(i,i) =5;endfor i= l :(n-1),a(i,i+1)=6;endfor

25、 i= l :(n-1),a(i+ l ,i)= l ;enda方法-:. ( - 個(gè)for) n=5;a=ze ros(n,n);a( l ,1:2)=5 6;for i=2:(n- l ),a(i, i-1,i,i+ 1)=1 5 6;end a(n,n-1 n)=l 5;a方法三（不用 for) n=5;a=diag(5*ones(n, l);b=小ag(6*onse(n- l , l);仁 工iiag(ones(n-1,1);a=a+fze ros(n-1, l), b;ze ros(l ,n)+fzeros(l ,n);c,ze ros( n-1,1) l下列 計(jì)算det(a) an

26、s=6650.3173-0.58651.0286-1.62411.9489-0.09770.4887-0.85711.3534-1.62410.0286-0.14290.5429-0.85711.0286-0.00750.0376-0.14290.4887-0.58650.0015-0.00750.0286-0.09770.3173 in v(a) ans =lv,d=eig(a) V=-0.78430.5546-0.26140.0924-0.0218d=0.7574- 0.78 43-0 .92 3 70 . 9860- 0 .92 37-0 .5 546-0.3771- 0 .00 000

27、.3771- 0.26 140.0000-0 .164 30.0000-0.09240.0628-0.0000-0.0628-0.02180.02570.02740.0257oooo0 0。9.2426007.4495000005.0000000002.55057 判斷第 6 題各小題是否可以相似對(duì)角化， 如果是， 求出對(duì)角矩陣和對(duì)應(yīng)的相似變換矩陣。8 判斷第 6 題各小題是否為正定矩陣。9 求下列向量組的秩和它的一個(gè)最大線性無(wú)關(guān)組， 并將其余向量用該最大無(wú)關(guān)組線性表示。a,= (4, -3, 1,3), a2 = (2, -1, 3, 5), a 3= (1, -1, -1, -1), a4

28、 = (3, -2, 3, 4), a 5= (7, -6, -7, 0)10.( 二次型標(biāo)準(zhǔn)化）用正交變換化下列二次型為標(biāo)準(zhǔn)形f (Xi, X2心 ）= x產(chǎn)- 4 x 1 心+ 4x 1x 3- 2 x /+8心 X 3- 2 X 3%Exercise l0a= I -2 2; -2 -2 4;2 4 -2v,d=eig( a) v=0 .33330 .9339-0.12930 .6667-0 .3304-0 .6 681-0.66670 .1365-0 .7327d =-7.00000002.00000002.0000 v*v ans =1.00000.00000.00000.00001

29、.000000.000001.0000 %v 確實(shí)是正交矩陣泊 （電路網(wǎng)）圖 3.1 是連接三個(gè)電壓已知終端的電路網(wǎng)，求 a, b, c 點(diǎn)的電壓。20 4n3n20VaII bIov占c3n |5Vsn30圖 3.1 電路圖1236 12 .( Ha milton-Carle y 定理）就矩陣A =4561驗(yàn)證下列性質(zhì)780(i) 設(shè)八 A,z, ， 凡，為n 階方陣A 的特 征值， 則ItrA-; = (-1) l41;i = I(ii) 設(shè)f (x)為 A 的特 征多項(xiàng)式，則f (A) = 0。兒， =La; (A 的跡），”“i = Ii=I習(xí)題 41 求下列多項(xiàng)式的所有根，并進(jìn)行驗(yàn)算

30、。(1) x2+x + l ; ro ots( I I l)(2) 3x5 -4 x3+2x- 1; roots ( l3 0 -4 0 2 - I J) (3) 5x23 - 6x1+8x6- 5 武p=ze ros( l ,24);p( J 17 18 2 2) =5 - 6 8 -5 ;roots(p)(4) (2x+3)3 - 4 ( 提示： 先用 conv 展開(kāi)）p I = r2 31;% 們?yōu)?3 次方p2=conv(pl , pl );o/oc o nv 返回的足 P l *P l p3=conv(pl , p2);p3(end)= p3(end)-4; % 原 p3 品后個(gè)分臥

31、 4 roots(p3)2 求方程 x ln(f 了言+ x) - 心了習(xí)- 0 .Sx = O 的正根。fun=inline(x* log(sqrt(xA2 - l ) +x) -s q rt ( x A2 - l ) - 0 . S* x) ; fzero(fun,2)3 用 MATLAB 指令求解第一章習(xí)題 4。4 (超越方程）超越方程的解有時(shí)是很復(fù)雜的，作 出f (x) = x sin (1/x)在 - 0.1, 0.1內(nèi)的圖，可見(jiàn) 在 x = O 附近f (x) = 0 有無(wú)窮多個(gè)解，并設(shè)法 求出它們 的近似 解，使計(jì)算結(jié)果誤差不超過(guò)O.Ql 。5 求解下列非線性方程組在原點(diǎn)附近的根

32、9x2 + 36 y2+4 z2= 3 6x2-2 y2-2 0 z= 0 16x-x3 -2/-16z2 = 06 求解下列方程組在區(qū)域 Oa, /31 內(nèi)的解t-0 .7 s in a + 0 . 2 c o s jl/3 =0.7c o s a-0 . 2 s i n /3fu n=( x)x( 1) -0 .7 卞sin( x( I) ) -0.2* cos(x(2),x(2)-0.7*cos(x( I) ) +0 .2 卞sin( x(2) ) ); a,b,c) =fso lve ( ft m,0 .5 0.5)7 (橢園的 交點(diǎn)）兩個(gè)橢圓可能具有0 4 個(gè)交點(diǎn)， 求下列兩個(gè) 橢園

33、的所 有交點(diǎn) 坐標(biāo)(x - 2)2 + (), - 3 + 2x)2 = 52 (x-3)2 + (y/3)2=48 作 出下列函數(shù)圖形，觀察所 有的 局 部極 大，局部極小和全局最大，全局錄小值點(diǎn)的粗略位置， 并用 M ATLAB 函數(shù) fminbnd 和 fminsearch 求各極值點(diǎn)的確切位置(1) f(x)=x2 s i n ( x2 - x - 2) , - 2 ,2 );cl e a r ;扣 n=inline(x.A2.* s in(x.A2-x-2);fplo t(fu n,-2 2);grid on; %作圖觀察x( l )=-2;x(3)=fminbnd(fu n,- l

34、 ,-0.5); x(5)=fminbnd(fun,1,2);fun2=i nline(-x.A2 . * s i n( x .A2 - x - 2 ) ) ; x (2) = f m inbnd(fun2,-2,- I) ;x (4) = f m inbnd(fun2,-0.5,0.5);x(6)=2feval(fun,x)答案 以_r. x( 1)(3) (5)足屁部極小，x( 2) (4 )(6 ) 是局 部極大 ，從報(bào)后一句知道 x( l )全局品小， x(2 )垃大。(2)f(x ) = 3x5-20x3 + 1 0 , ( 3)f ( x) = Ix3 - x2 - x - 21

35、- 3, 3);0, 3).9 考慮函數(shù).f(x,y) = /!9+3x+9x2+y2切+9(1)作出 凡，y) 在- 2勺1 , -7 y l 的圖，觀 察極 值點(diǎn)的 位 置(2)用 M AT L A B 函數(shù) fminsearch 求極值點(diǎn)和極值。close;x=-2:0.1: 1;y=-7:0.1:1;x, y=mes hgrid( x,y) ;z=y.A3/9+3*x.A2.* y + 9* x.A2+y.A2 + x .*y+9; mesh(x, y,z);grid on;%作圖觀察扣n=inline(x(2)A3/ 9 + 3 * x( 1 )A2* x( 2) + 9 沁x( l

36、 )A2 + x( 2 )A2+ x ( I t x( 2) + 9 );x = f m i nsearch(fun,0 0)%求極小伯fun2=inline(-(x(2)A3/ 9 +3 *x ( 1)A2 x(2 )+9* x(1) A2 + x( 2) A2 + x( 1)*x (2)+ 9) );x=fminsearch(fun2,LO -5)%求極大伯10 假定某天的氣溫變化記錄如第二章習(xí)題 5, 試用最小二乘方法找出這一天的氣溫變化規(guī)律。考慮下列類型函數(shù)，作圖比較效果， 并計(jì)算均方誤差。(1) 二次函數(shù)；( 2 ) 三次函數(shù)；(3) 鐘 形 函數(shù)j(x) = a eb (r-14

37、) 1;(4) 函數(shù) f (x) = r s in(冗12t + 0)% E xerciseI 0 cl e a r;t=0:24;c=f l 5 14 14 14 14 15 16 18 20 22 23 25 28 31 3 2 31 2 9 27 25 24 22 20 18 17 16 ;p2=polyfit( t,c,2) p3=polyfit(t,c,3)fu n=i nline(a( l )*ex p(a(2)*(t- l4).A2 ) , a , t) ; a=ls q c u rvefit(fun, lO 0 ,t,c) % 初 伯 可 以試 探f=feval(fun, a,

38、t)norm(f-c)%擬合效朱plo t(t,c ,t,t) % 作圖檢驗(yàn)fu n2=inline(b(1)*sin(pi/121 t + b(2)+20,b,t);% 丿京題修改f(x)+20 b=lsqcurvefit(fun2,0 0 ,t,c)figuref2=feval( fu n2, b,t)norm(f2-c)%擬合效果plo t( t,c ,t, t2) % 作圖檢驗(yàn)11 (化 學(xué)反應(yīng)平衡）一等克分子數(shù)一氧化碳(CO)和氧氣(0 2)的 混合物在 300K 和 5bru 壓 力 下達(dá)到平衡， 理論反應(yīng)方程式為C 0 + 0 .5 0 2 今C O 2實(shí)際反應(yīng)方程式為CO+N2

39、x CO+0.5 (l+x) 0 2 + (l- x) CO2剩余co 比值 x 滿足化學(xué)平衡方程式(l-x)-Jl0.52 + XK=Px.J百飛壓O x I這里 Kp = 3.06, p = 5 bar 求 X.12(月還款額）作為房產(chǎn)公司的代理人，你要迅速準(zhǔn)確回答客戶各方面的問(wèn)題。現(xiàn)在有個(gè)客戶看 中了你公司一套建筑面積為 180 平方米，每平方單價(jià)7500 元的房子。他計(jì)劃首付 30%,其余70%用 20 年按揭貸款（貸款年利率 5.04%)。請(qǐng)你提供下列信息： 房 屋 總 價(jià)格、首付款額、月付還款額。如果其中 10 萬(wàn)元為公積金貸款（貸款年利率 4.05%)呢？%Exercise 12

40、1-=5.04 /12/1OO;N=20*12;x=7500*180 %房屆總價(jià)格y=歡 0.3 %竹付款額x0=x-y%貸款總額a=( I+r)N*r雙0/( I +r)N-1)%月付還款額rl =4.05/12/1OO;x I = I O* 10 000; %公積金貸款a I=(I +r I) N*rl*x1/(J +rI) N-I )x2=x0-x l %商業(yè)貸款a2=(1+r)N*r*x2/( l +r)N-1) a=al +a213 (栓牛鼻的繩子）農(nóng)夫老李有一個(gè)半徑 10 米的圓形牛欄， 里面長(zhǎng)滿了草， 老 李要 將家里一頭牛栓在一根欄樁上，但只 讓牛吃到一半草，他想讓上大學(xué) 的兒

41、 子告訴他，栓 牛鼻的 繩子應(yīng) 為多長(zhǎng)？t,14 (弦截法）牛頓迭代法是一種速度很快的迭代方法，但 是它需 要預(yù)先求得導(dǎo)函數(shù)。若用差商代替導(dǎo)數(shù)，可 得下列弦截法xk+l= xk -xk -x k- 1f伈 ）-f(xH)f ( Xk)這一迭代法需要兩個(gè)初值 Xo, X1 ,編 寫(xiě) 一 個(gè) 通 用 的 弦 截 法計(jì)算機(jī)程序并用以解習(xí)題 2。（提示函數(shù)參數(shù)求值用 MATLAB函數(shù) feval)!.1 5 (線性 迭 代） 迭代過(guò)程的 收 斂 性主要條件是在根的附近滿足 l,gX k+l = g (x k)Cx)I 1。從理論上證明線性迭代X k+i = a x k + l只有兩種極限形態(tài)：不動(dòng) 點(diǎn)

42、 或 無(wú)窮大 C分 別就a=0.9, -0.9, 1.1, -1.1(取Xo =1, 迭代 20 步）用圖形顯示迭代過(guò)程的不同表現(xiàn)（提示： 用 s ubplot 將 4 個(gè)子圖放在一個(gè)圖形窗口比較）!.1 6 ( 通道中的細(xì)桿）要運(yùn)送一根細(xì)桿子通過(guò)由寬 5cm 和寬 10cm 的通道垂直交叉口，在運(yùn)送過(guò)程中必須保持桿子是水平的（如圖 4.6), 問(wèn)這根細(xì)桿至多可有多長(zhǎng)？ 又通道為園柱形的且細(xì)桿不必保持水平， 細(xì)桿至多可有多長(zhǎng)？6 17證明當(dāng)且僅當(dāng) 3al + 石 ，Logi stic 映射有穩(wěn)定的周期 2 軌叫;10cm圖 4.6.l OOO,在 (x k, Yk) 處 亮一點(diǎn)（注意 不要連線

43、）可得所謂 Henon 引力線圖習(xí)題 51. 某河床的橫斷面如圖 5.8 所示， 為了計(jì)算最大的排洪蜀， 需 要計(jì) 算 它的斷面積，試 根據(jù)圖示測(cè)量數(shù)據(jù)（單位： 米 ）用梯 形 法計(jì)算其斷面積。x=lO 4 10 12 15 22 28 34 40;y=O I 3 6 8 9 5 3 OJ;trapz ( x, y)2. 求圖 5.8 各測(cè)盤(pán)點(diǎn)的坡度。x=O4 10 12 15 22 28 34 40; y=O I 3 6 8 9 5 3 OJ;di ff(y)./diff(x)圖 5.85O4101 21 522283 44023. 作圖 表 示函數(shù) z = x e- - y,(-lxl,

44、0y2), 沿 x 軸方向的梯度。xa= -1:0 . 1: I ;ya =0 : 0. 1:2;fx,yl=mes hgrid(xa,ya); z=x.*exp(-x. 2 -y.3);fp x,pyl = gradie nt(z,xa,ya); pxx = l n co s td ydy4 已知參 數(shù)方程,Otl.5,試取 t 的步長(zhǎng)O.Ql,求 和的數(shù)值解。y = cos t -t sin tdxd xx =- 1t=0:0.01:1.5;x=log(cos(t);y=cos(t)-t.*sin(t); dydx=gradient(y,x)x_ l ,id=min(abs(x-(-1);%找昆核近 x=-1 的點(diǎn)dydx(id)5 求下列積分的數(shù)值解(1) 1 1e-立，2心墨r(2)芢c o s3 (x)dx,(3)X ln(x4) arcsin dxX( 4 ) s in(x)dx ,X(5) 仁 d x,(6) 叫 1 + r2 sin(0)dr (7)ff(l +x+y2) dydx ,D 為 x2 +l傘D%Exercise 5(2)扣n=inline(exp(2*x).*cos(x).A3); quadl