matlab程序設(shè)計例題及答案

1、matlab 程序設(shè)計例題及答案1. 編寫程序：計算 1/3+2/5+3/7+ +10/21 法一： s=0;for i=1:10 s=s+i/(2*i+1); end ss = 法二：sum(1:10)./(3:2:21) ans =2. 編寫程序：計算 1100 中即能被 3 整除，又能被 7 整 除的所有數(shù)之和。s=0;for i=1:100if mod(i,3)=0&mod(i,7)=0 s=s+i; end,end ss =2103. 畫出 y=n！的圖 , 階乘的函數(shù)自己編寫，禁用 MATLAB 自帶的階乘函數(shù)。x=1:10; for i=1:10try y(i)=y(i-1)*i

2、;catch y(i)=1; end,endplot(x,y)106123456789104. 一個數(shù)恰好等于它的因子之和，這個數(shù)就稱為完數(shù)。例如， 6 的因子為 1， 2,3 ，而 6=1+2+3，因此 6 就是一個完 數(shù)。編程找出 20XX 以內(nèi)的所有完數(shù)。g=;for n=2:20XX s=0;for r=1:n-1if mod(n,r)=0 s=s+r; end endif s=ng=g n;end end gg =628 4965. 編寫一個函數(shù)，模擬 numel 函數(shù)的功能，函數(shù)中調(diào)用 size 函數(shù)。function y=numelnumel(x) m=size(x); y=m(

3、1)*m(2); numelnumel(1 2 3;4 5 6)ans =66. 編寫一個函數(shù)，模擬 length 函數(shù)的功能，函數(shù)中調(diào) 用 size 函數(shù)。function y=lengthlength(x) m=size(x); y=max(m(1),m(2);lengthlength(1 2 3;4 5 6)ans =37. 求矩陣 rand 的所有元素和及各行平均值，各列平均值。 s=rand(5);sum=sum(sum(s) mean2=mean(s,2) mean1=mean(s) sum = mean2 = mean1 =8. 編程判斷 1001,1003,1007 ，1009

4、， 1011 為素數(shù)，若 不是，輸出其約數(shù)。function sushupd(x) s=fix(sqrt(x); p=1;for i=2:sif mod(x,i)=0 p=0; break end endif p=1disp( 該數(shù)是素數(shù) ) else g=; for i=1:xif mod(x,i)=0 g=g i; end enddisp( 該數(shù)不是素數(shù) ) g end 是否sushupd(1001) 該數(shù)不是素數(shù) g =1 7 11 13 77 91 143 sushupd(1003) 該數(shù)不是素數(shù) g =100117591003 sushupd(1007)該數(shù)不是素數(shù) g =1 19

5、53 1007 sushupd(1009) sushupd(1011)1 3 337該數(shù)是素數(shù) 該數(shù)不是素數(shù) g = 10119. 將1563 35 4123 65 121233 356 57; 1223543 43 23132 345 45 31 從小到大排序，排成一列，禁用 sort 函數(shù) s=1563 35 4123 65 121233 356 57;1223543 43 23 132345 45 31; s=s(1,:) s(2,:); for i=1:length(s)-1 fora=s(i);i=1:length(s)-1 if s(i)s(i+1)s(i)=s(i+1)s(i+1

6、)=a; end end,end s s =Columns 1 through 823 31 3543455765132Columns 9 through 14345 356 15634123121233122354310. 兩個整數(shù) a 和 b，如果 a 的因子和等于b，b 的因子和也等于 a，且 a 不等于 b，則稱 a，b 為一對親密數(shù)。請尋找并輸出 20XX 以內(nèi)的所有親密數(shù)對g=;for a=1:20XXs1=0;s2=0;for r=1:a-1if mod(a,r)=0 s1=s1+r; end endif a=s1 continue endif find(g=a)=0 conti

7、nue endfor r=1:s1-1if mod(s1,r)=0 s2=s2+r;end endif s2=as=s1 s2;g=g;s;end end g g =284 220 1210 118411. 用二分法求函數(shù) y=exp(x)-10 的 0 點，要求誤差小 于x1=0; x2=10;while x2-x1if (exp(x1)-10)*(exp(x1+x2)/2)-10)0 x1=(x1+x2)/2; else x2=(x1+x2)/2;end endx=x1 x2 x =12. 現(xiàn)有一組數(shù)如下： x y 1 12 2 27 3 34 4 48 5 676 79 7 89 8 1

8、09 但 y 的數(shù)據(jù)不太準(zhǔn)確，已知其函數(shù)關(guān)系可 能是 y=12*x ， y=12*x+1 ， y=12*x-1 ，y=*x, 請用 matlab 判 斷這四個函數(shù)，哪個與數(shù)據(jù)的契合度最高？x=1:8;y(1,:)=12*x; y(2,:)=12*x+1; y(3,:)=12*x-1; y(4,:)=*x;y(5,:)=12 27 34 48 67 79 89 109; t=zeros(1,4);for i=1:4for j=1:8t(i)=t(i)+abs(y(i,j)-y(5,j);end endx,i=min(t) i =4 即為 y=*x13. 已知 rand 函數(shù)能生成 0 到 1 之

9、間的隨機(jī)數(shù)。而古人 用拋石頭的方法，根據(jù)落在在一個圓內(nèi)的石子數(shù)與在園外正 方形內(nèi)的石子數(shù)之比來求 ，請用 matlab 軟件模擬此過程， 并求出 的值，然后，分析所求 的值的誤差產(chǎn)生于什么地 方。s=rand(2,1000); pai=0;for i=1:1000if s(1,i).2+s(2,i).2 sushupd(1003)該數(shù)不是素數(shù) g =1 17 59 1003 sushupd(1007)該數(shù)不是素數(shù) g =1 19 53 1007 sushupd(1009) 該數(shù)是素數(shù) sushupd(1011) 該數(shù)不是素數(shù) g =1 3 337 10119. 將1563 35 4123 65 121233 356 57; 1223543 43 23132 345 45 31 從小到大排序，排成一列，禁用 so