实验6 Matlab程序设计2
实验目的：
1、 掌握建立和执行M文件的方法； 2、 掌握实现选择结构的方法； 3、 掌握实现循环结构的方法。
4、 熟悉利用向量运算来代替循环操作的方法。 实验内容：
1. 根据
?26?1111?????，求π的近似值。当n分别取100、1000、10000时，122232n2结果是多少？
要求：分别用循环结构和向量运算（使用sum函数）来实现。 sum = 0;
n=input('input a number:'); for i=1:n
sum=sum+1/i^2; y=sqrt(6*sum); end y yi
input a number:100 y =
3.1321 >> yi
input a number:1000 y =
3.1406 >> yi
input a number:10000 y =
3.1415 2. 根据y?1?111????，求： 352n?1（1） y<3时的最大n值。 （2） 与（1）的n值对应的y值。
3. 先用函数的递归调用定义一个函数文件求
1005010?ii?1nm，然后调用该函数文件求
?k??k??2k?1k?11 k?1k
s=0;
>> s=factor(100,1)+factor(50,2)+factor(10,-1) s =
4.7978e+004 4. 已知
?f1?1,?????????????????????????????n?1?f?0,?????????????????????????????n?2?2 ?f?1,?????????????????????????????n?3?3??fn?fn?1?2fn?2?fn?3,?n?3求f1?f100中：
（1） 最大值、最小值、各数之和。
（2） 正数、零、负数的个数。 5. 写出下列程序的输出结果。 （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 s = 108
(2) 命令文件exe.m global x
x = 1:2:5;y = 2:2:6; sub(y); x,y
函数文件sub.m function fun = sub(z) global x z = 3*x;x = x+z; exe x =
4 12 20
y =
2 4 6