數(shù)學(xué)模型習(xí)題解答

1、上機(jī)練習(xí)題一班級： 姓名： 學(xué)號： 1.建立起始值=3,增量值=5.5,終止值=44的一維數(shù)組x 答案： x=(3:5.5:44)2.寫出計(jì)算 Sin(30o)的程序語句. 答案： sin(pi*30/180) 或 sin(pi/6) 3.矩陣,矩陣;分別求出及A與B中對應(yīng)元素之間的乘積. 答案：A = 3,2,3; 4,2,6; 7,8,1 B = 1,1,1; 2,2,2; 3,3,3 A*B；A.*B4計(jì)算行列式的值。答案：det(A)5對矩陣 進(jìn)行下述操作。(1)求秩。答案：rank(A)(2)求轉(zhuǎn)置。答案：A(3) 對矩陣求逆，求偽逆。答案：inv(A) ，pinv(A)(4) 左右

2、反轉(zhuǎn),上下反轉(zhuǎn)。答案：fliplr(A)，flipud(A)(5) 求矩陣的特征值. 答案：u,v=eig(A)(6) 取出上三角和下三角. 答案：triu(A) tril(A)(7)以A為分塊作一個(gè)3行2列的分塊矩陣。答案：repmat(a)6 計(jì)算矩陣與之和。 a=5 3 5;3 7 4;7 9 8; b=2 4 2;6 7 9;8 3 6; a+b7 計(jì)算與的數(shù)組乘積。 a=6 9 3;2 7 5; b=2 4 1;4 6 8; a.*bans = 12 36 3 8 42 408 已知：，分別計(jì)算a的數(shù)組平方和矩陣平方，并觀察其結(jié)果。 a=1 2 3;4 5 6;7 8 9; a.2a

3、ns = 1 4 9 16 25 36 49 64 81 a2ans = 30 36 42 66 81 96 102 126 150上機(jī)練習(xí)題二班級： 姓名： 學(xué)號： 1 對于，如果，求解X。 A=4 9 2;7 6 4;3 5 7; B=37 26 28; X=ABX = -0.5118 4.0427 1.33182 角度，求x的正弦、余弦、正切和余切。 x=30 45 60; x1=x/180*pi; sin(x1)ans = 0.5000 0.7071 0.8660 cos(x1)ans = 0.8660 0.7071 0.5000 tan(x1)ans = 0.5774 1.0000

4、1.7321 cot(x1)ans = 1.7321 1.0000 0.57743 將矩陣、和組合成兩個(gè)新矩陣：（1）組合成一個(gè)43的矩陣，第一列為按列順序排列的a矩陣元素，第二列為按列順序排列的b矩陣元素，第三列為按列順序排列的c矩陣元素，即 （2）按照a、b、c的列順序組合成一個(gè)行矢量，即 答案：a=4 2;5 7; b=7 1;8 3; c=5 9;6 2;% (1) d=a(:) b(:) c(:) d = 4 7 5 5 8 6 2 1 9 7 3 2% (2) e=a(:);b(:);c(:) e = 4 5 2 7 7 8 1 3 5 6 9 2 或利用（1）中產(chǎn)生的d e=re

5、shape(d,1,12) ans = 4 5 2 7 7 8 1 3 5 6 9 24求解在x=8時(shí)多項(xiàng)式(x-1)(x-2) (x-3)(x-4)的值。 p=poly(1 2 3 4); polyvalm(p,8)% polyval(p,8) ans = 8405求方程的全部根。p=3,7,9,0,-23; %建立多項(xiàng)式系數(shù)向量x=roots(p) %求根上機(jī)練習(xí)題三班級： 姓名： 學(xué)號： 1、 設(shè)x是數(shù)組，求均值和方差解：函數(shù)文件如下： function xx,s=func1(x)n=length(x);xx=sum(x)/n;s=sqrt(sum(x.2)-n*xx2)/(n-1);命

6、令窗口： x=1 2 3 4 5;xx,s=func1(x)2、求滿足的最小m值s=0;n=0;while(se) k=k+1; a=abs(fun(k)/fun(k-1)-(1+sqrt(5)/2);end命令行： k,a=funTest(10-8)k = 21a = 9.7719e-009或者M(jìn)文件如下： clear; F(1)=1;F(2)=1;k=2;x=0; e=1e-8; a=(1+sqrt(5)/2; while abs(x-a)e k=k+1; F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1); end a,x,k4、分別用for和while循環(huán)結(jié)構(gòu)編寫程序，

7、求出，并考慮一種避免循環(huán)語句的程序設(shè)計(jì)，比較各種算法的運(yùn)行時(shí)間。解：for循環(huán)結(jié)構(gòu)：M文件loop.mk=0;for i=1:106 k=k+sqrt(3)*2-i;endkwhile循環(huán)結(jié)構(gòu)：M文件loop1.mk=0;i=1;while itic;loop;toc 循環(huán)結(jié)構(gòu)的執(zhí)行時(shí)間k = 1.7321Elapsed time is 1.813000 seconds. tic;nonLoop;toc 非循環(huán)結(jié)構(gòu)的執(zhí)行時(shí)間k = 1.7321Elapsed time is 1.094000 seconds.上機(jī)練習(xí)題四班級： 姓名： 學(xué)號： 1、作圖描述氣溫變化 x=0:24; y=15,1

8、4,14,14,14,15,16,18,20,22,23,25,28,31,32,31,29,27,25,24,22,20,18,17,16; plot(x,y)2、作出下列函數(shù)圖形（1） （分別使用plot和fplot完成）解： fplot(x2*sin(x2-x-2),-2 2) fplot方法 x=-2:0.1:2;y=x.2.*sin(x.2-x-2);plot(x,y) plot方法 如圖（4.1）（2） （橢圓 提示：用參數(shù)方程）解： r=-pi:0.1:pi;x=2*cos(r);y=3*sin(r);plot(x,y) 如圖（4.2）解法二 x=-2:1/100:2; y1=3

9、*sqrt(1-x.2/4); y2=-3*sqrt(1-x.2/4); plot(x,y1,r-,x,y2,r-); axis equal tight; 圖（4.1） 圖（4.2）（3） （拋物面） 解：（錯誤） x=-3:0.1:3;y=-3:0.1:3;z=x.2+y.2; plot3(x,y,z) 如圖（4.31）（正確） xa=-3:0.1:3;ya=-3:0.1:3;x,y=meshgrid(xa,ya); 如圖（4.32） z=x.2+y.2;mesh(x,y,z); surf(x,y,z) 圖（4.31）error 圖（4.32）（4）曲面解： xa=linspace(-3,3

10、,100);ya=linspace(-3,13,100); x,y=meshgrid(xa,ya); z=x.4+3*x.2+y.2-2*x-2*y-2*x.2.*y+6; mesh(x,y,z) surf(x,y,z) （5）空間曲線解： t=linspace(0,2,50);x=sin(t);y=cos(t);z=cos(2*t); plot3(x,y,z)（6）半球面解： a=linspace(0,pi/2,50);b=linspace(0,2*pi,50); a,b=meshgrid(a,b); x=2*sin(a).*cos(b);y=2*sin(a).*sin(b);z=2*cos

11、(a); surf(x,y,z) （7）三條曲線合成圖解： x=linspace(0,pi,50);y1=sin(x); plot(x,y1);hold on; y2=sin(x).*sin(10*x); plot(x,y2); y3=-sin(x); plot(x,y3); hold off; 3、作下列分段函數(shù)圖x=-5:0.1:5;for i=1:length(x) if x(i)1.1 y(i)=1.1; elseif x(i)1.1)+x.*( x =-1.1)-1.1.*(x1); %第a1列b1行對應(yīng)的x+y1 (x對應(yīng)列；y對應(yīng)行) %第a1列對應(yīng)的x值是xa(a1)；第b1行

12、對應(yīng)的y值是ya(b1)z(a1-1)*ny+b1)=0.5457*exp(-0.75*ya(b1).2-3.75*xa(a1).2-1.5*xa(a1); a2,b2=find(x+y-1);z(a2-1)*ny+b2)=0.7575*exp(-ya(b2).2-6*xa(a2).2);a3,b3=find(x+y pxy運(yùn)行結(jié)果如右圖：或者M(jìn)文件如下：clear;close; xa=-2:0.1:2;ya=-2:0.1:2;x,y=meshgrid(xa,ya); z=zeros(size(x); k1=find(x+y1); z(k1)=0.5457*exp(-0.75*y(k1).2-

13、3.75*x(k1).2-1.5*x(k1); k2=find(x+y-1); z(k2)=0.7575*exp(-y(k2).2-6*x(k2).2); k3=find(x+y1)+.0.7575*exp(-y.2-6*x.2).*(x+y-1)+.0.5457*exp(-0.75*y.2-3.75*x.2+1.5*x).*(x+ydemo2、查詢trapz的功能、用法、目錄、程序結(jié)構(gòu)、相同目錄下其它文件解： help trapz 功能用法 type trapz程序結(jié)構(gòu)，源碼 which trapz所在目錄 help C:MATLAB6p5toolboxmatlabdatafun該目錄下其它

14、文件3在0,4pi畫sin(x),cos(x)(在同一個(gè)圖象中); 其中cos(x)圖象用紅色小圓圈畫.并在函數(shù)圖上標(biāo)注 “y=sin(x)”, “y=cos(x)” ,x軸,y軸,標(biāo)題為“正弦余弦函數(shù)圖象”. x=linspace(0,4*pi,100);y=sin(x);plot(x,y);gtext(y = sin(x); % 圖形注解,注意要用鼠標(biāo)定位hold on;y=cos(x);plot(x,y,ro);gtext (y = cos(x); % 圖形注解xlabel(x軸); % x軸注解ylabel(y軸); % y軸注解 title(正弦余弦函數(shù)圖象); % 圖形標(biāo)題4從鍵盤

15、輸入若干個(gè)數(shù)，當(dāng)輸入0時(shí)結(jié)束輸入，求這些數(shù)的平均值和它們之和。sum=0;cnt=0;val=input(Enter a number (end in 0):);while (val=0) sum=sum+val; cnt=cnt+1; val=input(Enter a number (end in 0):);endif (cnt 0) sum mean=sum/cntend5若一個(gè)數(shù)等于它的各個(gè)真因子之和，則稱該數(shù)為完數(shù)，如6=1+2+3，所以6是完數(shù)。求1,500之間的全部完數(shù)。for m=1:500s=0;for k=1:m/2if rem(m,k)=0s=s+k;endendif m

16、=s disp(m);endend上機(jī)練習(xí)題六班級： 姓名： 學(xué)號： 1 假定數(shù)據(jù)點(diǎn)來源為：，試根據(jù)生成的數(shù)據(jù)進(jìn)行插值處理，得出較平滑的曲線。x=0:.12:1;y=(x.2-3*x+5).*exp(-5*x).*sin(x);plot(x,y,x,y,o)x1=0:.02:1;y0=(x1.2-3*x1+5).*exp(-5*x1).*sin(x1);y1=interp1(x,y,x1); y2=interp1(x,y,x1,cubic); y3=interp1(x,y,x1,spline); y4=interp1(x,y,x1,nearest); plot(x1,y1;y2;y3;y4,:

17、,x,y,o,x1,y0)或者plot(x1,y1,y2,y3,y4,:,x,y,o,x1,y0)2 用不同插值的方法計(jì)算sin(x)在pi/2的值clearclcx=0:0.2:pi;%產(chǎn)生包含被插值點(diǎn)的采樣點(diǎn)(做成一個(gè)向量)y=sin(x);%求出各采樣點(diǎn)對應(yīng)的樣本值y1=interp1(x,y,pi/2);%用默認(rèn)的linear方法計(jì)算sin(pi/2)y2=interp1(x,y,pi/2,nearest);%用默認(rèn)的linear方法計(jì)算sin(pi/2)y3=interp1(x,y,pi/2,cubic);%用三次多項(xiàng)式插值方法計(jì)算sin(pi/2)y4=interp1(x,y,pi

18、/2,spline);%用三次樣條插值方法計(jì)算sin(pi/2)y5=spline(x,y,pi/2);%直接用spline方法計(jì)算sin(pi/2),功能與y4相同disp(各種方法的插值結(jié)果:)out=y1=, num2str(y1) y2=,num2str(y2); y3=,num2str(y3) out2=y4=,num2str(y4) y5=,num2str(y5)3求 的通解x,y=dsolve(Dx=4*x-2*y,Dy=2*x-y, t)%P149(3)的求解x,y=dsolve(Dx=a*y,Dy=-b*x , x(0)=x0, y(0)=y0, t)y=dsolve(Dy=

19、(b*x)/(a*y), y(0)=y0,x);4 方法一：%微分方程的M函數(shù)funt.m文件function y=funt(t,y)y=(y2-t-2)/4/(t+1);%以下是求解的腳本m函數(shù)，可自由取名，然后在command窗口調(diào)用該函數(shù)求解.ts=0,10;%自變量的求解區(qū)間y0=2;%初值條件t,y=ode23(funt,ts,y0);%用2,3階龍格庫塔方法求funt文件里的微分方程y1=sqrt(t+1)+1;%求精確解對應(yīng)點(diǎn)上的函數(shù)值t,y,y1%以三列的形式顯示結(jié)果，其中第一列為采樣點(diǎn)t對應(yīng)的值，第二列為t對應(yīng)的數(shù)值解，第二列為t對應(yīng)的精確解。方法二：%或者直接用內(nèi)建函數(shù)編寫

20、待求的微分方程，在一個(gè)M腳本文件里執(zhí)行求解。funt=inline(y2-t-2)/4/(t+1) %內(nèi)建函數(shù)編寫待求的微分方程ts=0,10;%自變量的求解區(qū)間y0=2;%初值條件t,y=ode23(funt,ts,y0);%用2,3階龍格庫塔方法求funt文件里的微分方程y1=sqrt(t+1)+1;%求精確解對應(yīng)點(diǎn)上的函數(shù)值t,y,y15第一種方法建立lorenz函數(shù)模型的狀態(tài)方程function xp=lorenz(t,x)%建立lorenz函數(shù)模型的狀態(tài)方程xp=-8/3,0,x(2);0,-10,10;-x(2),28,-1*x%表明x是一個(gè)三維向量，前面是一個(gè)線性方程組的系數(shù)矩陣

21、，xp%是一個(gè)三維的輸出表示x(1),x(2),x(3)的一階導(dǎo)數(shù)向量。第二種方法建立lorenz函數(shù)模型的狀態(tài)方程function xp=lorenz1(t,x)%建立lorenz函數(shù)模型的狀態(tài)方程xp=-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(2)*x(1)+28*x(2)-x(3);%求解微分方程clear;clc;x0=0,0,eps;%三個(gè)初值條件構(gòu)成的向量t,x=ode23(lorenz,0,100,x0);t,xplot(t,x),grid,pause%繪制解x(1),x(2),x(3)各自相對于變量t的圖象，按任意鍵后繼續(xù)下面的程序figu

22、re;plot3(x(:,1),x(:,2),x(:,3);%繪制解x(1),x(2),x(3)的關(guān)于系統(tǒng)的相平面圖象axis(10,40,-20,20,-20,20);上機(jī)練習(xí)題七班級： 姓名： 學(xué)號： 1有一組測量數(shù)據(jù)如下表所示，數(shù)據(jù)具有y=x2+1的變化趨勢，用最小二乘法求解y。x11.522.533.544.55y-1.42.735.98.412.216.618.826.2 x=1 1.5 2 2.5 3 3.5 4 4.5 5 y=-1.4 2.7 3 5.9 8.4 12.2 16.6 18.8 26.2 e=ones(size(x) x.2 c=ey x1=1:0.1:5; y1

23、=ones(size(x1),x1.2*c; plot(x,y,ro,x1,y1,k)2求下列線性方程組的解（1）解： a=4 1 -1;3 2 -6;1 -5 3;b=9; -2; 1;x=ab 唯一解x = 2.3830 1.48942.0213（2）解： a=4 -3 3;3 2 -6;1 -5 3;b=-1;-2;1;x=ab 唯一解x = -0.4706 -0.2941 0（3）解： a=4 1;3 2;1 -5;b=1;1;1;x=ab 最小二乘近似解x = 0.3311 -0.1219（4），求通解解： a=2 1 -1 1;1 2 1 -1; 1 1 2 1;b=1;2;3;

24、rank(a),rank(a,b) ans = 3ans = 3 說明有無窮多解 rref(a,b) 行最簡化ans = 1.0000 0 0 1.5000 1.0000 0 1.0000 0 -1.5000 0 0 0 1.0000 0.5000 1.0000通解為：，3、求下列矩陣的行列式、逆、特征值和特征向量（1）解： a=4 1 -1;3 2 -6;1 -5 3;det(a),inv(a),v,d=eig(a)ans = -94 行列式 ans = 0.2553 -0.0213 0.0426 0.1596 -0.1383 -0.2234 矩陣的逆 0.1809 -0.2234 -0.0

25、532v = 0.0185 -0.9009 -0.3066 -0.7693 -0.1240 -0.7248 特征向量 -0.6386 -0.4158 0.6170d = -3.0527 0 0 0 3.6760 0 特征值 0 0 8.3766（2）解： a=4 3 1;3 3 -5;1 -5 3;det(a),inv(a),v,d=eig(a)ans = -124 行列式ans = 0.1290 0.1129 0.1452 0.1129 -0.0887 -0.1855 矩陣的逆 0.1452 -0.1855 -0.0242v = 0.3757 -0.8583 0.3496 -0.6881 -

26、0.0056 0.7255 特征向量 -0.6208 -0.5131 -0.5927d = -3.1480 0 0 0 4.6176 0 特征值 0 0 8.5304(3) 解： a=5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10;det(a),inv(a),v,d=eig(a)ans = 1 行列式ans = 68.0000 -41.0000 -17.0000 10.0000 -41.0000 25.0000 10.0000 -6.0000 矩陣的逆 -17.0000 10.0000 5.0000 -3.0000 10.0000 -6.0000 -3.0000 2.00

27、00v = 0.8304 0.0933 0.3963 0.3803 -0.5016 -0.3017 0.6149 0.5286 特征向量 -0.2086 0.7603 -0.2716 0.5520 0.1237 -0.5676 -0.6254 0.5209d = 0.0102 0 0 0 0 0.8431 0 0 0 0 3.8581 0 特征值 0 0 0 30.2887(4)n=5;tria = 6*eye(n-1); %構(gòu)造上三角rown = zeros(1,n-1);coln = zeros(n,1);tria = tria;rown;tria = coln,triaeyeb = 5*

28、eye(n); %構(gòu)造對角陣tric = 1*eye(n-1); %構(gòu)造下三角tric = rown;tric;tric = tric,colnA=tria+eyeb+tric 4求多項(xiàng)式的根，并分析誤差大?。?）解： A=1,1,1;roots(A)ans = -0.5000 + 0.8660i -0.5000 - 0.8660i y=polyval(A,ans) 驗(yàn)證y = 1.0e-015 * 0.3331 0.3331（2）解： A=3 0 -4 0 2 -1;x=roots(A)x = -0.9479 + 0.3845i -0.9479 - 0.3845i 1.0000 0.4479

29、 + 0.3435i 0.4479 - 0.3435i y=polyval(A,x) 驗(yàn)證y = 1.0e-013 * -0.0144 - 0.1138i -0.0144 + 0.1138i -0.0888 0.0044 - 0.0083i 0.0044 + 0.0083i（3）解： A=5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6 8 0 0 0 -5 0 0;x=roots(A) y=polyval(A,x)(4) 解：p1=2 3;p2=conv(p1,p1); p=conv(p2,p1); %(2x+3)3n=length(p);for (i=1:n-1) q(

30、i)=p(i);endq(n)=p(n)-4; %(2x+3)3-4x=roots(q)結(jié)果：x = -1.8969 + 0.6874i -1.8969 - 0.6874i -0.7063驗(yàn)證： y=polyval(q,x)y = 1.0e-014 * -0.7105 - 0.6217i -0.7105 + 0.6217i 0 上機(jī)練習(xí)題八班級： 姓名： 學(xué)號： 1有一正弦衰減數(shù)據(jù)y=sin(x).*exp(-x/10)，其中x=0:pi/5:4*pi，用三次樣條法進(jìn)行插值。 x0=0:pi/5:4*pi; y0=sin(x0).*exp(-x0/10); x=0:pi/20:4*pi; y=

31、spline(x0,y0,x); plot(x0,y0,or,x,y,b)2 編制一個(gè)解數(shù)論問題的函數(shù)文件：取任意整數(shù)，若是偶數(shù)，則用2除，否則乘3加1，重復(fù)此過程，直到整數(shù)變?yōu)?。 function c=collatz(n) % collatz % Classic “3n+1” Ploblem from number theory c=n; while n1 if rem(n,2)=0 n=n/2; else n=3*n+1; end c=c n; end3有一組學(xué)生的考試成績（見表），根據(jù)規(guī)定，成績在100分時(shí)為滿分，成績在9099之間時(shí)為優(yōu)秀，成績在8089分之間時(shí)為良好，成績在6079

32、分之間為及格，成績在60分以下時(shí)為不及格，編制一個(gè)根據(jù)成績劃分等級的程序。學(xué)生姓名王張劉李陳楊于黃郭趙成 績728356941008896685465方法一：clearName=王,張,劉,李,陳,楊,于,黃,郭,趙; Marks=72,83,56,94,100,88,96,68,54,65; % 劃分區(qū)域：滿分(100)，優(yōu)秀(90-99)，良好(80-89)，及格(60-79)，不及格(=90&Marks(i)=80&Marks(i)=60&Marks(i)=79 Rank(i,:)= 及格 ; else Rank(i,:)=不及格 ; end end % 將學(xué)生姓名，得分，級等信息打印出

33、來。 disp( ) disp(學(xué)生姓名 , 得分 , 等級); disp(-) for i=1:10; disp( ,Name(i), ,num2str(Marks(i), ,Rank(i,:); end 方法2：Name=王,張,劉,李,陳,楊,于,黃,郭,趙; Marks=72,83,56,94,100,88,96,68,54,65; % 劃分區(qū)域：滿分(100)，優(yōu)秀(90-99)，良好(80-89)，及格(60-79)，不及格(60)。 n=length(Marks); for i=1:10 ai=89+i; bi=79+i; ci=69+i; di=59+i; end; c=d,c; % 根據(jù)學(xué)