1 实验四 循环结构程序设计 1.根据n3212222211116,求 的近似值。当n 分别取100、1000、10000 时,记过是多少? n=100 程序设计: n=100; i=1:n; f=1./i.^2; x=sum(f); y=sqrt(6*x) 运行结果: x = 3.1321 n=1000 程序设计: n=1000; i=1:n; f=1./i.^2; x=sum(f); y=sqrt(6*x) 运行结果: x = 3.1406 n=10000 程序设计: n=1000; i=1:n; f=1./i.^2; x=sum(f); y=sqrt(6*x) 运行结果: x = 3.1415 2.根据y=1+1-n215131,求: (1)y<3 时的最大 n 值。 (2)与(1)的n 值对应的y值。 程序设计: y=1;n=1; while(y<3) n=n+2; y=y+1/n; end 2 y=y-1/n n=(n+1)/2 运行结果: y = 2.9944 n = 57 3. 考虑以下迭代公式: xxnnba1 其中a、b 为正的常数。 (1)编写程序求迭代的结果,迭代的终止条件为10 51 xxnn,迭代初值x0=1.0,迭代次数不超过 500 次。 程序设计: x=1.0; a=2; b=3; x=a/(b+x); n=1; while(abs(x-(a/x-b))> 10e-5&n<=500) x=a/(b+x); n=n+1; end x n 运行结果: x = 0.5616 n = 6 (2)如果迭代过程收敛于 r,那么 r 的准确值是242abb ,当(a,b)的值取(1,1)、(8,3)、(10,0.1)时,分别对迭代结果的准确值进行比较。 (a,b)取(1,1) 程序设计: x=1.0; a=1; b=1; x=a/(b+x); 3 n=1; while(abs(x-(a/x-b))>10e-5&n<=500) x=a/(b+x); n=n+1; end x r1=(-b+sqrt(b^2+4*a))/2 r2=(-b-sqrt(b^2+4*a))/2 运行结果: x = 0.6181 r1 = 0.6180 r2 = -1.6180 (a,b)取(8,3) 程序设计: x=1.0; a=8; b=3; x=a/(b+x); n=1; while(abs(x-(a/x-b))>10e-5&n<=500) x=a/(b+x); n=n+1; end x r1=(-b+sqrt(b^2+4*a))/2 r2=(-b-sqrt(b^2+4*a))/2 运行结果: x = 1.7016 r1 = 1.7016 r2 = -4.7016 (a,b)取(10,0.1) 程序设计: x=1.0; a=10; b=0.1; x=a/(b+x); 4 n=1; while(abs(x-(a/x-b))>10e-5&n<=500) x=a/(b+x); n=n+1; end x r1=(-b+sqrt(b^2+4*a))/2 r2=(-b-sqrt(b^2+4*a))/2 运行结果: x = 3.1127 r1 = 3.1127 r2 = -3.2127 4.已知3,23,12,01,1321321nnnnfffffffnnnn 求ff1001~中: (1)最大值、最小值、各数之和。 程序设计: f(1)=1; f(2)=0;...