实验5 Matlab程序设计2.docVIP

实验5 Matlab程序设计2 实验目的： 掌握建立和执行M文件的方法； 掌握实现选择结构的方法； 掌握实现循环结构的方法。 熟悉利用向量运算来代替循环操作的方法。 实验内容： 根据，求π的近似值。当n分别取100、1000、10000时，结果是多少？ 要求：分别用循环结构和向量运算（使用sum函数）来实现。 根据，求： y3时的最大n值。 与（1）的n值对应的y值。 先用函数的递归调用定义一个函数文件求，然后调用该函数文件求 已知 求中： 最大值、最小值、各数之和。 正数、零、负数的个数。 12.写出下列程序的输出结果。 （1）s = 0; a = [12,13,14;15,16,17;18,19,20;21,22,23]; for k = a for j = 1:4 if rem(k(j),2)~=0 s = s +k(j); end end end s (2) 命令文件exe.m global x x = 1:2:5;y = 2:2:6; exsub(y); x,y 函数文件sub.m function fun = sub(z) global x z = 3*x;x = x+z; 详细实验内容： 1根据，求π的近似值。当n分别取100、1000、10000时，结果是多少？ function [y]=sy5_1_1(n) y=0; % n=100; for i=1:n; y=y+1/i/i end y function [y]=sy5_1_2(n) % n=100; i=1:n; f=1./i.^2; y=sum(f); 输出结果： sy5_1_2(100) ans = 1.6350 sy5_1_2(1000) ans = 1.6439 sy5_1_2(10000) ans = 1.6448 要求：分别用循环结构和向量运算（使用sum函数）来实现。 2根据，求： （1）y3时的最大n值。 （2）与（1）的n值对应的y值。 function[]=sy5_2(n) y=0; i=1; while (y3) f=1/(2*i-1); y=y+f i=i+1; end y=y-f; i=i-1; y i 结果： y = 2.9944 i = 57 3先用函数的递归调用定义一个函数文件求，然后调用该函数文件求 function[y]=sy5_3(n,m) y=0; for i=1:n; a=i.^m; y=y+a; end y=sy5_3(100,1)+sy5_3(50,2)+sy5_3(10,-1) y = 4.7978e+004 4已知 求中： （1）最大值、最小值、各数之和。 （2）正数、零、负数的个数。 function y = fun(n) if n==1 y = 1; elseif n==2 y = 0; elseif n==3 y = 1; else a = zeros(1,n); a(1:3) = [1,0,1]; for i = 4:n a(i)=a(i-1)-2*a(i-2)+a(i-3); end y = a(n); % y = fun(n-1)-2*fun(n-2)+fun(n-3); end maxa = max(a); mina = min(a); suma = sum(a); zero=0; zhengshu=0; fushu=0; for i=1:n if a(i) == 0 zero=zero+1; elseif a(i) 0 zhengshu=zhengshu+1; else fushu=fushu+1; end end maxa mina suma zero zhengshu fushu 执行： sy5_5(10) 输出： maxa = 5 mina = -7 suma = 2 zero = 2 zhengshu = 5 fushu = 3 ans = -7 5写出下列程序的输出结果。 （1）输出结果： s = 108 (2) 输出结果： x = 4 12 20 y = 2 4 6 雄鹰 实验五 Matlab程序设计2 第1页 共5页