实验四 循环结构程序设计.doc

实验四 循环结构程序设计 实验内容： 1：根据=+++……+，求π的近似值.当n分别取100,1000,10000时,结果是多少? 要求:分别用循环结构和向量运算(使用sum函数)来实现. 循环结构： 试验程序：(1)y1=0;n=100; for i=1:n y1=y1+1/i/i; end y=sqrt(6*y1) 实验结果：y = 3.1321 试验程序：(2) y1=0;n=1000; for i=1:n y1=y1+1/i/i; end y=sqrt(6*y1) 实验结果：y = 3.1406 试验程序：y1=0;n=10000; for i=1:n y1=y1+1/i/i; end y=sqrt(6*y1) 实验结果：y = 3.1415 向量运算： 试验程序：n=100; i=1:n; f=1./i.^2; y1=sum(f); y=sqrt(6*y1) 实验结果：y = 3.1321 试验程序：n=1000; i=1:n; f=1./i.^2; y1=sum(f); y=sqrt(6*y1) 实验结果：y = 3.1406 试验程序：n=10000; i=1:n; f=1./i.^2; y1=sum(f); y=sqrt(6*y1) 实验结果：y = 3.1415 实验内容： 2: y=1+++……+,求： （1） y3时的最大n值。 试验程序：y=0; n=1; while y3 y=y+1/(2*n-1); n=n+1; end n=n-1 实验结果：n = 57 （2） 与（1）的n值对应得y值。 试验程序：y=0; n=1; while y3 y=y+1/(2*n-1); n=n+1; end y=y-1/(2*(n-1)-1) 实验结果：y = 2.9944 实验内容： 4： 已知： 求f1~f100中： (1) 最大值，最小值，各数之和。 试验程序：p(1)=1; p(2)=0; p(3)=1; for n=4:100 p(n)=p(n-1)-2*p(n-2)+p(n-3); end min=p(1); max=p(1); sum=0; for n=1:100 if minp(n) min=p(n); end if maxp(n) max=p(n); end sum=sum+p(n); end min max sum 实验结果：min = -8.9941e+011 max = 4.3776e+011 sum = -7.4275e+011 (2) 正数，零，负数的个数。 试验程序：p=[]; p(1)=1; p(2)=0; p(3)=1; for n=4:100 p(n)=p(n-1)-2*p(n-2)+p(n-3); end zero=0; negative=0; positive=0; for n=1:100 if p(n)0 positive=positive+1; elseif p(n)0 negative=negative+1; else zero=zero+1; end end positive negative zero 实验结果：positive = 49 negative = 49 zero = 2 实验五 函数文件 实验内容： 1：定义一个函数文件，求给定复数的指数，对数，正弦和余弦，并在命令文件中调用该函数文件。 试验程序：函数文件：function[a,b,c,d]=exch(x) %sin 正弦函数 %cos 余弦函数 %log10 对数函数 %exp 指数函数 a=exp(x); b=log10(x); c=sin(x); d=cos(x); 函数文件调用：x=input(请输入一个复数：) [a,b,c,d]=exch(x) 实验结果：请输入一个复数：1+2i x = 1.0000 + 2.0000i a = -1.1312 + 2.4717i b = 0.3495 + 0.4808i c = 3.1658 + 1.9596i d = 2.0327 - 3.0519i 实验内容： 2：一物理系统可用下列方程组来表示： = 从键盘输入m1，m2和的值，求a1，a2，N1和N2的值。其中g取9.8，输入时以角度为单位。 要求