專題2 簡(jiǎn)單的程序設(shè)計(jì)(數(shù)學(xué)實(shí)驗(yàn))

1、2022年7月3日數(shù)學(xué)建模實(shí)驗(yàn)報(bào)告專題2 MATLAB基本使用目錄一、書上命令練習(xí)：3二、 MatLab電子版作業(yè)32.1練習(xí)1 :基本運(yùn)算32.2練習(xí)數(shù)組及矩陣運(yùn)算62.3 程序設(shè)計(jì)192.4 練習(xí)4繪圖30三、<數(shù)學(xué)實(shí)驗(yàn)>-預(yù)備實(shí)驗(yàn)作業(yè)參考解答36說(shuō)明：專題 2作業(yè)：Matlab 基本使用：電子版作業(yè)（基本運(yùn)算+數(shù)組+程序設(shè)計(jì)+作圖）+數(shù)學(xué)實(shí)驗(yàn)P35-實(shí)驗(yàn)內(nèi)容專題2一、書上命令練習(xí)：Diary filename20+30*rand(5)diary off2、 MatLab電子版作業(yè)2.1練習(xí)1 :基本運(yùn)算 注意: (1),請(qǐng)按要求作題,寫出程序并把結(jié)果附上,并保存文件.1, (

2、5*2+1.3-0.8)*102/25 答案：(5*2+1.3-0.8)*102/25ans = 422，y = 答案：sin(10*pi)*exp(-0.3+42)+log(23)/log(4)ans = 2.26183，z = 10sin(/3) . cos(/3) 答案： z = 10*sin(pi/3)* . cos(pi/3) z = 4.33014，x = sin(223/3)， y = x2，z = y*10 ； 求x+2y-5z 答案：x = sin(223/3); y = x2;z = y*10 ;>> x+2*y-5*zans = -37.60255，,計(jì)算 答

3、案：>> z1=2+7*iz1 = 2.0000 + 7.0000i>> z2=2*i;z3=5*exp(2*pi*i)z3 = 5.0000 - 0.0000i>> z=z1*z2/(z2+z3)z = -2.1379 + 1.6552i6, 建立起始值=3,增量值=5.5,終止值=44的一維數(shù)組x (填寫程序語(yǔ)句)答案：3：5.5：447, 建立等差一維數(shù)組x：首項(xiàng)為0,末項(xiàng)為,項(xiàng)數(shù)為15 (填寫程序語(yǔ)句) 答案：linspace(0,pi,15)8, 計(jì)算行列式的值 答案：>> A=3 2 3;4 2 6;7 8 1A = 3 2 3 4

4、2 6 7 8 1>> det(A)ans = -89, 矩陣,矩陣;分別求出A x B及A與B中對(duì)應(yīng)元素之間的乘積的程序語(yǔ)句。 答案：A=3 2 3;4 2 6;7 8 1;B=1 1 1;2 2 2;3 3 3;>> A*B,A.*Bans = 16 16 16 26 26 26 26 26 26ans = 3 2 3 8 4 1221 24 310,寫出計(jì)算 Sin(36o)的程序語(yǔ)句. 答案：sin(36*pi/180)ans = 0.58782.2 練習(xí)數(shù)組及矩陣運(yùn)算1.       輸入矩陣 。答案

5、：A=3 4 -1 1 -9 10;6 5 0 7 4 -16;1 -4 7 -1 6 -8;2 -4 5 -6 12 -8;-3 6 -7 8 -1 1;8 -4 9 1 3 0A = 3 4 -1 1 -9 10 6 5 0 7 4 -16 1 -4 7 -1 6 -8 2 -4 5 -6 12 -8 -3 6 -7 8 -1 1 8 -4 9 1 3 02.       輸入一個(gè)與A同階的隨機(jī)矩陣B。要求元素為整數(shù)。答案：B=fix(100*rand(size(A)（其他方式取整也可以）B= 84 68 30 15 49 34

6、52 37 18 69 89 28 20 83 19 37 82 34 67 50 68 86 64 53 83 70 30 85 81 72 1 42 54 59 66 303.       輸入一個(gè)數(shù)組b=（1 3 5 7 9 11 ）答案：b=1 3 5 7 9 11 或b=1:2:11b = 1 3 5 7 9 114.       輸入與A同階的單位矩陣E 答案：程序代碼：E=eye(size(A)結(jié)果：E = 1 0 0 0 0 0 0 1 0 0 0 0 0 0

7、 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 15.       輸入與B同階的零矩陣Z 答案： 程序代碼： Z=zeros(size(B)Z = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 06.       構(gòu)建一個(gè)大矩陣 答案：D=A Z;E,BD = 3 4 -1 1 -9 10 0 0 0 0 0 0 6 5 0 7 4 -16 0

8、 0 0 0 0 0 1 -4 7 -1 6 -8 0 0 0 0 0 0 2 -4 5 -6 12 -8 0 0 0 0 0 0 -3 6 -7 8 -1 1 0 0 0 0 0 0 8 -4 9 1 3 0 0 0 0 0 0 0 1 0 0 0 0 0 83 69 17 13 28 51 0 1 0 0 0 0 56 62 97 1 46 33 0 0 1 0 0 0 37 79 27 89 6 43 0 0 0 1 0 0 70 95 25 19 98 22 0 0 0 0 1 0 54 52 87 29 58 57 0 0 0 0 0 1 44 88 73 66 42 767.

9、60;     將 46-70這25個(gè)整數(shù)填入一個(gè)五行五列的矩陣數(shù)表G中，使其各行、各列以及主對(duì)角線和次對(duì)角線的和等相等。 答案： G=45+magic(5)G = 62 69 46 53 60 68 50 52 59 61 49 51 58 65 67 55 57 64 66 48 56 63 70 47 548.       分別求d1=|A|、d2=|B|、d3=|D|、d4=|G|答案：d1=det(A),d2=det(B),d3=det(D),d4=det(G)結(jié)果：d1 =

10、 245295d2 = 4.9274e+010d3 =1.2087e+016d4 = 226200009.       對(duì)以上四個(gè)矩陣A、B、D、G中的可逆矩陣求逆。答案：inv(A),inv(B),inv(D),inv(G)ans = -0.0737 0.0604 -0.2297 0.0067 -0.0804 0.1042 0.3142 0.0036 0.2408 0.1605 0.1259 -0.1436 0.2099 -0.0395 0.3155 0.0364 0.0834 -0.0663 -0.0827 -0.0123 0.0

11、088 -0.0777 0.0779 0.0878 0.0134 -0.0335 -0.0159 0.1129 0.1061 0.0337 0.0377 -0.0525 -0.0110 0.0469 0.0698 0.0411ans = 0.0114 -0.0014 0.0144 -0.0030 0.0209 -0.0300 -0.0005 0.0155 -0.0059 0.0046 -0.0350 0.0219 -0.0060 0.0117 0.0040 -0.0041 0.0037 -0.0049 -0.0064 -0.0094 0.0187 0.0005 0.0208 -0.0179 -

12、0.0069 -0.0121 -0.0050 0.0101 0.0143 -0.0009 0.0092 -0.0135 -0.0189 -0.0056 -0.0011 0.0259ans = -0.0737 0.0604 -0.2297 0.0067 -0.0804 0.1042 0 0 0 0 0 0 0.3142 0.0036 0.2408 0.1605 0.1259 -0.1436 0 0 0 0 0 0 0.2099 -0.0395 0.3155 0.0364 0.0834 -0.0663 0 0 0 0 0 0 -0.0827 -0.0123 0.0088 -0.0777 0.077

13、9 0.0878 0 0 0 0 0 0 0.0134 -0.0335 -0.0159 0.1129 0.1061 0.0337 0 0 0 0 0 0 0.0377 -0.0525 -0.0110 0.0469 0.0698 0.0411 0 0 0 0 0 0 -0.0011 -0.0010 -0.0016 -0.0016 -0.0000 0.0004 0.0114 -0.0014 0.0144 -0.0030 0.0209 -0.0300 -0.0036 -0.0002 -0.0023 0.0010 0.0003 0.0018 -0.0005 0.0155 -0.0059 0.0046

14、-0.0350 0.0219 -0.0052 0.0003 -0.0054 -0.0025 -0.0020 0.0030 -0.0060 0.0117 0.0040 -0.0041 0.0037 -0.0049 -0.0010 0.0009 -0.0050 -0.0006 -0.0019 0.0006 -0.0064 -0.0094 0.0187 0.0005 0.0208 -0.0179 0.0050 0.0008 0.0030 0.0014 -0.0009 -0.0027 -0.0069 -0.0121 -0.0050 0.0101 0.0143 -0.0009 0.0074 0.0000

15、 0.0116 0.0013 0.0028 -0.0047 0.0092 -0.0135 -0.0189 -0.0056 -0.0011 0.0259ans = -0.0073 0.0488 -0.0378 -0.0012 0.0010 0.0408 -0.0397 -0.0070 0.0103 -0.0009 -0.0326 0.0007 0.0007 0.0007 0.0340 0.0023 -0.0089 0.0084 0.0411 -0.0394 0.0004 0.0026 0.0392 -0.0474 0.008710.    (1) 計(jì)算d5=A+B

16、答案： d5=A+Bd5 = 86 73 16 14 19 61 62 67 97 8 50 17 38 75 34 88 12 35 72 91 30 13 110 14 51 58 80 37 57 58 52 84 82 67 45 76 (2) 計(jì)算d6=A*B答案：d6=A*Bd6 = 460 883 384 372 258 511 780 189 -58 -724 644 -363 20 -113 -269 259 -200 -68 3 -365 91 173 -326 129 378 408 528 -467 918 -61 1005 1266 277 1007 366 8561

17、1.    取A的1，3，5行做一個(gè)子陣a1. 答案： a1=A(1:2:5,:)a1 = 3 4 -1 1 -9 10 1 -4 7 -1 6 -8 -3 6 -7 8 -1 1 12.    取A的2，4，6行，4，5，6列上的交叉元素做一個(gè)子陣a2. 答案：a2=A(2:2:6,4:6)或：a2=A(2,4,6,4,5,6)a2 = 7 4 -16 -6 12 -8 1 3 013.    由a2產(chǎn)生列向量a3. 答案：a3=a2(:)a3 = 7 -6 1 4 12 3 -16 -8 014.&#

18、160;   取a3絕對(duì)值大于4的元素構(gòu)成列向量a4答案：a4=a3(find(abs(a3)>4)a4 = 7 -6 12 -16 -815.    改變A的元素A（4，5）為0 值。答案：A(4,5)=0A = 3 4 -1 1 -9 10 6 5 0 7 4 -16 1 -4 7 -1 6 -8 2 -4 5 -6 0 -8 -3 6 -7 8 -1 1 8 -4 9 1 3 016.    改變A的第6列為b答案：A(:,6)=b'A = 3 4 -1 1 -9 1 6 5 0 7 4 3 1

19、 -4 7 -1 6 5 2 -4 5 -6 0 7 -3 6 -7 8 -1 9 8 -4 9 1 3 1117.    建立三角函數(shù)表陣S。第一列是自變量X的取值，從0-2共取20個(gè)值。第二列是對(duì)應(yīng)Y=cos X的值。答案：X=linspace(0,2*pi,20);>> Y=cos(X);>> S=X;YS = 0 1.0000 0.3307 0.9458 0.6614 0.7891 0.9921 0.5469 1.3228 0.2455 1.6535 -0.0826 1.9842 -0.4017 2.3149 -0.6773 2.6

20、456 -0.8795 2.9762 -0.9864 3.3069 -0.9864 3.6376 -0.8795 3.9683 -0.6773 4.2990 -0.4017 4.6297 -0.0826 4.9604 0.2455 5.2911 0.5469 5.6218 0.7891 5.9525 0.9458 6.2832 1.000018.以紅色的+號(hào)   畫y的圖形 答案：plot(X,Y,'r+')19.    分別取a2的1、2行作成兩個(gè)向量A1、A2 答案：A1=a2(1,:),A2=a2(2,:)A1 = 7 4

21、 -16A2 = -6 12 -820.    作向量的點(diǎn)積 C1=A1·A2 答案：C1=dot(A1,A2)C1 = 134或者：>> C1=A1*A2'C1 = 13421.    作向量的叉積 C2=A1×A2 答案： C2=cross(A1,A2)C2 = 160 152 1082.3 程序設(shè)計(jì) 1, 編輯函數(shù)文件hanshu.m,該函數(shù)是，并用它來(lái)計(jì)算、答案：函數(shù)文件hanshu.m如下：function y=hanshu(x)y=x.3-3*x.2-1*x+3；計(jì)算程序及結(jié)果如下：

22、x=0.34 3 6.87x = 0.3400 3.0000 6.8700>> hanshu(x)ans =2.3525 0 178.7820>> hanshu(3)-hanshu(2)*(hanshu(-4)2ans = 330752，編程求一個(gè)四位數(shù)，使。答案：(填寫程序語(yǔ)句和結(jié)果)3，利用二分法編程求方程 在0,3內(nèi)的根. 答案： (填寫程序語(yǔ)句和結(jié)果) 先編輯函數(shù)文件hanshu.m如下：function y=hanshu(x)y=x.3-3*x.2-x+3；然后編寫m文件如下：cleara=0; b=3; while abs(a-b)>10(-6) if

23、 hanshu2(a)*hanshu2(b)>0 b=(a+b)/2;else a=(a+b)/2;endendgen123=a結(jié)果： gen123 = 3.0000*4, 利用二分法編程求方程 tan(x)-x-0.5=0 在0,1.5內(nèi)的根. 答案：(填寫程序語(yǔ)句和結(jié)果)先定義函數(shù)：fin.mfunction y=fin(x)y=tan(x)-x-0.5;然后編輯m文件qiugen.m: fplot('fin',0 1.5)a=0;b=1.5;while abs(a-b)>10(-6) if fin(a)*fin(a+b)/2)>0 a=(a+b)/2;

24、else b=(a+b)/2; endenda結(jié)果： 0.97505,列出所有的水仙花數(shù), 水仙花數(shù)是一個(gè)三位數(shù),其各位數(shù)字立方和等于該數(shù)本身. 例如: 答案：(填寫程序語(yǔ)句和結(jié)果) 程序1：for a=1:9 for b=0:9 for c=0:9 abc=100*a+10*b+c; if abc=a.3+b.3+c.3 disp(abc) end end endend程序2：k=1;for a=1:9 for b=0:9 for c=0:9 abc=100*a+10*b+c; if abc=a.3+b.3+c.3 shuixianhua(k)=abc;k=k+1; end end ende

25、ndshuixianhua結(jié)果：shuixianhua = 153 370 371 407程序3：for m=100:999 m1=fix(m/100);%求m的百位數(shù)字 m2=rem(fix(m/10),10);%求m的10位數(shù)字 m3=rem(m,10); if m=m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) endend 6.利用 公式求的近似值,直到誤差小于。 k=1;s=0;while 4/(2*k-1)>10(-6) s=s+(-1)(k+1)/(2*k-1); k=k+1;endpai=4*s 結(jié)果：pai = 3.14159215358972【

26、中學(xué)時(shí)我們就已經(jīng)知道： ，從而，如果我們應(yīng)用泰勒公式將 展開(kāi)，就可以得到 把 代入，就得到 這個(gè)簡(jiǎn)單的公式就可以用循環(huán)實(shí)現(xiàn)來(lái)實(shí)現(xiàn)圓周率的近似計(jì)算（雖然它并不真正適用于計(jì)算圓周率）】7，寫出小于5000的、立方的末四位是8888的所有自然數(shù)的程序答案： (填寫程序語(yǔ)句和結(jié)果)clear;clcn=5000;ii=1;for kj=1:5000 if rem(kj3,10000)=8888 shu(ii)=kj;ii=ii+1; endendshu結(jié)果：shu = 1942 4442*8，先思考如何用推理的方法（1）求Fibonacc數(shù)列中第2007個(gè)數(shù)除以6 的余數(shù)（2）若將Fibonacc數(shù)列

27、按如下方式分組：（1）、（1，2）、（3，5，8）、（13，21，34，55）、問(wèn)2007組數(shù)的和除以6的余數(shù)是幾？然后再用matlab編程驗(yàn)證你的結(jié)果。答案： (填寫你的思考過(guò)程和你的結(jié)果以及matlab程序語(yǔ)句)提示：（1）可以發(fā)現(xiàn)Fibonacc數(shù)列中的數(shù)除以6 的余數(shù)呈周期性的變化（2）第2007組數(shù)開(kāi)始是幾？結(jié)尾是幾？考慮一個(gè)周期內(nèi)的和除以6的余數(shù)。9. 應(yīng)用近似求之值，要求誤差小于.程序：k=1;s=1;while 1/prod(1:k)>10(-6) s=s+1/prod(1:k); k=k+1;ende=s結(jié)果：e = 2.7182815255731910，定義一個(gè)函數(shù)s

28、hjiech計(jì)算雙階乘，要求對(duì)輸入的一個(gè)整數(shù)N，該函數(shù)能立即計(jì)算出其雙階乘。并用它來(lái)計(jì)算200！答案： (填寫程序語(yǔ)句和結(jié)果) function jiech=shjiech(n)if rem(n,2)=0 as=2:2:n;else as=1:2:n;endjiech=prod(as);11，定義多元函數(shù)，并用它來(lái)計(jì)算、，并畫出它的圖形。答案： (填寫程序語(yǔ)句和結(jié)果)函數(shù)文件hanshu4.mfunction fxy=hanshu4(a,b)fxy=(a2+b2)*exp(-(a2+b2);計(jì)算： >>hanshu4(-2,3)ans =2.938428229075371e-005

29、>> hanshu4(3,4),hanshu4(0,0)ans = 3.471985966241005e-010ans = 0>>繪圖程序：x=-2:0.01:2;y=x;xx,yy=meshgrid(x,y);zz=(xx.2+yy.2).*exp(-(xx.2+yy.2);mesh(zz)結(jié)果： *12，編寫函數(shù)文件SQRT1.m用迭代法求之值，迭代公式是，誤差小于10-6。答案： (填寫程序語(yǔ)句和結(jié)果)function genhaoa=SQRT1(a) an=1; while 1/2*abs(a/an-an)>10(-6) an=(an+a/an)/2;en

30、dgenhaoa=an;比如：SQRT1(2)ans = 1.41421356237469SQRT1(2456)ans = 49.55804687847163補(bǔ)充1. 求100,200之間的素?cái)?shù)。I.for循環(huán)程序： for m=101:2:200 k=fix(sqrt(m); for i=2:k+1 if rem(m,i)=0 break; end end if i>=k+1%除遍了所有（2到K+1） disp(int2str(m) end end 或者：（以數(shù)組的形式顯示）u=0;for m=101:2:200 k=fix(sqrt(m); for i=2:k+1 if rem(m,

31、i)=0 break; end end if i>=k+1 u=u+1; a(u)=m; end end b=a'II.while循環(huán)程序： m=101; while m<200 k=fix(sqrt(m); i=2; while i<=k+1 if rem(m,i)=0 break; end i=i+1; end if i>=k+1 disp(int2str(m) end m=m+2;end 或：i=1;x=;for j=100:200 if isprime(j)=1%此命令可直接判斷j是否為素?cái)?shù) x(i)=j; end i=i+1;end或：N=input(

32、'請(qǐng)輸入數(shù)值：m=')x=100:N; for j=2:sqrt(N); i=find(rem(x,j)=0); x(i)=; endx2.若一個(gè)數(shù)等于它各個(gè)真因子之和，則稱該數(shù)為完數(shù)，如6=1+2+3，所以6就是完數(shù)，求1,500之間的全部完數(shù)。程序：for m=1:500 s=0; for k=1:m/2 if rem(m,k)=0 s=s+k; end end if m=s disp(m) endend3.求雞兔同籠問(wèn)題，頭共36，腳共100.求雞、兔各多少？程序：i=1;while 1 if rem(100-i*2,4)=0&(i+(100-i*2)/4)=36

33、 break; end i=i+1;enda1=ia2=(100-2*i)/4a1 =22a2 = 144.用起泡法對(duì)10個(gè)數(shù)由小到大排序，即相鄰兩個(gè)比較，將小的調(diào)到前頭。程序1：x=input('x=')n=length(x);for i=1:n%比較n輪 for j=1:n-i %兩個(gè)比較第一輪比較n-1次，。 if x(j)>x(j+1) a=x(j); x(j)=x(j+1); x(j+1)=a; end endendy=x程序2.function y=qp(x)m,n=size(x);for i=1:n%比較n輪 for j=1:n-i %兩個(gè)比較第一輪比較n-

34、1次，。 if x(j)>x(j+1) a=x(j); x(j)=x(j+1); x(j+1)=a; end endendy=x5.有一個(gè)4×5矩陣，編程求出最大值及其所處的位置。程序1：x=input('x=')c,t=max(x);y,i=max(c);t(i);disp('最大值為：',num2str(y)disp('位置為：',num2str(t(i),'行',num2str(i),'列')程序2：function max-1(x)c,t=max(x);y,i=max(c);t(i);dis

35、p('最大值為：',num2str(y)disp('位置為：',num2str(t(i),'行',num2str(i),'列')2.4 練習(xí)4繪圖1, 在0,4pi畫sin(x),cos(x)(在同一個(gè)圖象中); 其中cos(x)圖象用紅色小圓圈畫.并在函數(shù)圖上標(biāo)注 “y=sin(x)”, “y=cos(x)” ,x軸,y軸,標(biāo)題為“正弦余弦函數(shù)圖象”. 答案：函數(shù)文件如下：t=0:0.02:4*pi;plot(t,sin(t),'m*',t,cos(t),'rO')text(pi/3,sin(pi

36、/3),'y=sin(x)')text(pi/3,cos(pi/3),'y=cos(x)')title('正線函數(shù)和余弦函數(shù)')xlabel('x軸')ylabel(' y軸')結(jié)果顯示：2, 畫出橢圓。 答案：函數(shù)文件如下：%畫橢圓tt=0:0.02:2*pi;xx=4*cos(tt);yy=3*sin(tt);axis equalplot(xx,yy,'m*')title('橢圓x2/16+y2/9=1')答案2：函數(shù)文件：xx=-4:0.02:4;yy=3*sqrt(1-xx.

37、2/16);plot(xx,yy,'m*')axis equalhold onplot(xx,-yy,'r+')title('橢圓x2/16+y2/9=1')3,繪制螺旋線. 答案：u=0:0.01:30;x=(1+2*cos(u);y=1+2*sin(u);z=3*u;plot3(x,y,z,'r*')4,繪制四葉玫瑰線。 答案：tt=0:0.03:2*pi;pp=sin(tt*2);polar(tt,pp,'m+')title('四葉玫瑰線')5，用subplot分別在不同的坐標(biāo)系下作出下列圖形

38、，并為每幅圖形加上標(biāo)題。第一副：標(biāo)題：“概率曲線” 第二副：標(biāo)題“雙扭線” 第三副：標(biāo)題“葉形線” 第四副：標(biāo)題：“曳物線” 答案：x=-5:0.1:5;y=exp(-x.2);subplot(2,2,1)plot(x,y)title('概率曲線')tt=-pi/4:0.03:pi/4;pp=2*sqrt(cos(tt*2);subplot(2,2,2)polar(tt,pp)hold on polar(tt,-pp)title('雙扭線')subplot(2,2,3)t=0:0.01:20;x=3*t./(1+t.3);y=3*t.2./(1+t.3);plo

39、t(x,y)title('葉形線')subplot(2,2,4)yyy=0:0.01:1;xxx=log(1+sqrt(1-yyy.2)./(yyy+eps)-sqrt(1-yyy.2);plot(xxx,yyy)hold onxxx=log(1-sqrt(1-yyy.2)./(yyy+eps)+sqrt(1-yyy.2);plot(xxx,yyy)title('曳物線') 三、<數(shù)學(xué)實(shí)驗(yàn)>-預(yù)備實(shí)驗(yàn)作業(yè)參考解答1設(shè)有分塊矩陣A = ，其中E, R, 0, S分別為單位矩陣、隨機(jī)矩陣、零矩陣和對(duì)角矩陣，試通過(guò)數(shù)值計(jì)算驗(yàn)證A2 = 程序：E=eye(3

40、);R=rand(3,2);O=zeros(2,3);S=diag(1,2);A=E,R;O,S; .B=A*A,C=E R+R*S;O S*S,a=B=C結(jié)果：B = 1.0000 0 0 0.8205 1.0586 0 1.0000 0 1.7873 2.4395 0 0 1.0000 0.1158 0.0296 0 0 0 1.0000 0 0 0 0 0 4.0000C = 1.0000 0 0 0.8205 1.0586 0 1.0000 0 1.7873 2.4395 0 0 1.0000 0.1158 0.0296 0 0 0 1.0000 0 0 0 0 0 4.0000a =

41、 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 12 用命令magic(n) 生成幻方矩陣，通過(guò)計(jì)算研究它的性質(zhì)，如行和、列和、兩條對(duì)角線和等(可以利用命令diag, sum, fliplr, flipud，其用法可以查閱MATLAB幫助系統(tǒng))程序：A=magic(4),s1=sum(A),s2=sum(A'),.s3=sum(diag(A),s4=sum(diag(fliplr(A)結(jié)果：A = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1s1 = 34 34 34 34s2 = 34 34 34 34

42、s3 = 34s4 = 343 設(shè)y1 = 1/(1+x2)，y2 = exp(x2/2)，y3 = sin(2x)，y4 =，x在2, 2內(nèi)適當(dāng)離散化，計(jì)算y1 + y2，y1y2，y3/y2，(5y4 y1)/y22程序：x=-2:0.5:2;y1=-1./(1+x.2);y2=exp(-x.2/2);y3=sin(2*x);y4=sqrt(4-x.2);>> a=y1+y2,b=y1.*y2,c=y3./y1,d=(5*y4-y1)./y2.2結(jié)果：a = -0.0647 0.0170 0.1065 0.0825 0 0.0825 0.1065 0.0170 -0.0647b = -0.0271 -0.0999 -0.3033 -0.7060 -1.0000 -0.7060 -0.3033 -0.0999 -0.0271c = -3.7840 0.4586 1.8186 1.0518 0 -1.0518 -1.8186 -0.4586 3.7840d = 10.9196 65.6748 24.9002 13.4597 11.0000 13.4597 24.9002 65.6748