(答案)专：程序设计一

1、文件的命名：如：2022040108温玲玲程序设计题目、试编程实现以下问题公式求 的近似值,直到误差小于10 61 1 11. 利用一1 -4357k=1;s=0;while 4/(2*k-1)>10A(-6) s=s+(-1F(k+1)/(2*k-1); k=k+1;endpai=4*s结果：pai =112. 应用e 1 1 近似求e之值，要求误差小于102!n!程序：k=1;s=1;while 1/prod(1:k)>10A(-6)s=s+1/prod(1:k);k=k+1;ende=s结果：e =3. 求100,200之间的素数I.for循环程序：for m=101:2:2

2、00k=fix(sqrt(m);for i=2:k+1if rem(m,i)=0break;endendif i>=k+1% 除遍了所有( 2 到 K+1 ) disp(int2str(m)endend或者：(以数组的形式显示) u=0;for m=101:2:200k=fix(sqrt(m);for i=2:k+1if rem(m,i)=0 break;endendif i>=k+1u=u+1;a(u)=m;endendb=a'II.while 循环程序：m=101;while m<200k=fix(sqrt(m);i=2;while i<=k+1if rem

3、(m,i)=0break;endi=i+1;endif i>=k+1disp(int2str(m)endm=m+2;end或：i=1;x=;for j=100:200if isprime(j)=1% 此命令可直接判断 j 是否为素数x(i)=j;endi=i+1;end4. 假设一个数等于它各个真因子之和，那么称该数为完数，如 6=1+2+3 ，所以 6 就是完数，求 1,500 之间的全部完数。程序： for m=1:500s=0;for k=1:m/2if rem(m,k)=0 s=s+k;endendif m=sdisp(m)endend5. 列出所有的水仙花数 , 水仙花数是一个

4、三位数 , 其各位数字立方和等于该数 本身 .例如: 153 13 53 33程序： 程序 1 ：for a=1:9for b=0:9for c=0:9abc=100*a+10*b+c;if abc=a.A3+b.A3+c.A3disp(abc)endendendend程序 2 ：k=1;for a=1:9for b=0:9for c=0:9abc=100*a+10*b+c;if abc=a.A3+b.A3+c.A3shuixianhua(k)=abc;k=k+1;endendendend shuixianhua 结果： shuixianhua =153 370 371 407 程序 3： f

5、or m=100:999m1=fix(m/100);% 求 m 的百位数字m2=rem(fix(m/10),10);% 求 m的 10 位数字 m3=rem(m,10);if m=m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m)endend6. 写出小于 5000 的、立方的末四位是 8888 的所有自然数的程序 程序： clear;clcn=5000;ii=1;for kj=1:5000if rem(kjA3,10000)=8888 shu(ii)=kj;ii=ii+1;endendshu结果：shu =1942 44427. 求鸡兔同笼问题，头共 36，脚共 100.求

6、鸡、兔各多少？ 程序：i=1;while 1if rem(100-i*2,4)=0&(i+(100-i*2)/4)=36break;endi=i+1;end a1=i a2=(100-2*i)/4a1 =22 a2 =14、请编写函数 M 文件解决以下问题221，定义多元函数 f(x,y) (x2 y2)e (x y ) ，并用它来计算 f( 2,3) 、 f (3, 4) 、 f (0,0) ，并画出它的图形。答案： ( 填写程序语句和结果 ) 函数文件 hanshu4.mfunction fxy=hanshu4(a,b)fxy=(aA2+bA2)*exp(-(aA2+bA2);计算

7、： >>hanshu4(-2,3) ans =2.938428229075371e-005>> hanshu4(3,4),hanshu4(0,0)ans =ans =0>>绘图程序：x=-2:0.01:2;y=x;xx,yy=meshgrid(x,y);zz=(xx.A2+yy.A2).*exp(-(xx.A2+yy.A2);mesh(zz)结果：0 02.用起泡法对10个数由小到大排序，即相邻两个比拟，将小的调到 前头。程序1 :x=i nput('x=')n=len gth(x);for i=1:n% 比拟 n 轮for j=1:n-i

8、%两个比拟第一轮比拟n-1次，。if x(j)>x(j+1)a=x(j);x(j)=x(j+1);x(j+1)=a;endendendy=x程序2.function y=qp(x)m,n=size(x);for i=1:n% 比拟 n 轮for j=1:n-i % 两个比拟第一轮比拟 n-1 次，。 if x(j)>x(j+1) a=x(j);x(j)=x(j+1); x(j+1)=a;endendendy=x3.有一个4X5矩阵，编程求出最大值及其所处的位置 程序 1 ：x=input('x=') c,t=max(x);y,i=max(c);t(i);disp(' 最大值为： ',num2str(y)disp('位置为：,num2str(t(i),'行',num2str(i),'列')程序 2 ：fun