MATLAB程序设计例题及答案.pdf

1.编写程序：计算 1/3+2/5+3/7+10/21 法一： s=0; for i=1:10 s=s+i/(2*i+1); end s s = 4.4096 法二： sum(1:10)./(3:2:21) ans = 4.4096 2.编写程序：计算 1100 中即能被 3 整除，又能 被 7 整除的所有数之和。 s=0; for i=1:100 if mod(i,3)=0 end,end s s = 210 3.画出 y=n！的图（1 sushupd(1003) 该数不是素数 g = 117591003 sushupd(1007) 该数不是素数 g = 119531007 sushupd(1009) 该数是素数 sushupd(1011) 该数不是素数 g = 133371011 9.将1563 35 4123 65 121233 356 57; 1223543 43 23 132 345 45 31从小到大排序， 排成一列，禁用 sort 函数。 s=1563 35 4123 65 121233 356 57;1223543 43 23 132 345 45 31; s=s(1,:) s(2,:); for i=1:length(s)-1 for i=1:length(s)-1 if s(i)s(i+1) a=s(i); s(i)=s(i+1) s(i+1)=a; end end,end s s = Columns 1 through 8 2331354345 5765132 Columns 9 through 14 34535615634123121233 1223543 10. 两个整数 a 和 b，如果 a 的因子和等于 b，b 的因子和也等于 a，且 a 不等于 b，则称 a，b 为 一对亲密数。请寻找并输出 2000 以内的所有亲 密数对。 g=; for a=1:2000 s1=0; s2=0; for r=1:a-1 if mod(a,r)=0 s1=s1+r; end end if a=s1 continue end if find(g=a)=0 continue end for r=1:s1-1 if mod(s1,r)=0 s2=s2+r; end end if s2=a s=s1 s2; g=g;s; end end g g = 284220 12101184 11.用二分法求函数 y=exp(x)-10 的 0 点，要求 误差小于 0.001 x1=0; x2=10; while x2-x10.001 if (exp(x1)-10)*(exp(x1+x2)/2)-10)0 x1=(x1+x2)/2; else x2=(x1+x2)/2; end end x=x1 x2 x = 2.30222.3029 12.现有一组数如下： x12345678 y12273448677989109 但 y 的数据不太准确，已知其函数关系可能是 y=12*x， y=12*x+1，y=12*x-1，y=12.5*x,请用 matlab 判断这四个函数， 哪个与数据的契合度最 高？（提示：根据函数在各个点的值与原数据差 异度来判断契合度） x=1:8; y(1,:)=12*x; y(2,:)=12*x+1; y(3,:)=12*x-1; y(4,:)=12.5*x; y(5,:)=12 27 34 48 67 79 89 109; t=zeros(1,4); for i=1:4 for j=1:8 t(i)=t(i)+abs(y(i,j)-y(5,j); end end x,i=min(t) i = 4 即为y=12.5*x 13.已知 rand 函数能生成 0 到 1 之间的随机数。 而古人用抛石头的方法， 根据落在在一个圆内的 石子数与在园外正方形内的石子数之比来求， 请用 matlab 软件模拟此过程， 并求出的值， 然 后，分析所求的值的误差产生于什么地方。 s=rand(2,1000); pai=0; for i=1:1000 if s(1,i).2+s(2,i).21 pai=pai+4/1000; end end pai pai = 3.1120 14.编写一个函数文件， 函数的输入项为三维坐标 的两个点，函数没有输出项，函数运行结果为画 出这两个点所能确定的最小球 function qiu(x1,x2) o=(x1+x2)/2; r=sqrt(sum(x1-x2).2)/2; X,Y,Z = sphere(50); X=r*(X+o(1)