matlab程序设计例题及答案

1、1.编写程序：计算1/3+2/5+3/7+10/21法一：s=0;for i=1:10 s=s+i/(2*i+1);endss =4.4096法二：sum(1:10)./(3:2:21)ans = 4.40962.编写程序：计算1100中即能被3整除，又能被7整除的所有数之和。s=0;for i=1:100 if mod(i,3)=0&mod(i,7)=0 s=s+i; end,endss = 2103.画出y=n！的图（1=n sushupd(1003)该数不是素数g = 1 17 59 1003 sushupd(1007)该数不是素数g = 1 19 53 1007 sushupd(100

2、9)该数是素数 sushupd(1011)该数不是素数g = 1 3 337 10119.将1563 35 4123 65 356 57; 43 23 132 345 45 31从小到大排序，排成一列，禁用sort函数。s=1563 35 4123 65 356 57; 43 23 132 345 45 31;s=s(1,:) s(2,:);for i=1:length(s)-1for i=1:length(s)-1 if s(i)s(i+1) a=s(i); s(i)=s(i+1) s(i+1)=a; endend,endss = Columns 1 through 8 23 31 35 4

3、3 45 57 65 132 Columns 9 through 14 345 356 1563 4123 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

4、=s1 s2; g=g;s; endendgg = 284 220 1210 118411.用二分法求函数y=exp(x)-10的0点，要求误差小于0.001x1=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; endendx=x1 x2x = 2.3022 2.302912.现有一组数如下：x12345678y12273448677989109但y的数据不太准确，已知其函数关系可能是y=12*x， y=12*x+1，y=12*x-1，y=12.5*x,请

5、用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); endendx,i=min(t)i = 4即为y=12.5*x13.已知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; endendpaipai = 3.112014.编写一个函数文件，函数的输入项为三维坐标的两个点，函数没有输出项，函数运行结果为画出这两个点所能确定的最小球function qiu(x1,x2)o=(x1+x2)/2;r=sqrt(sum(x1-x2).2)/2;X