matlab实验四循环结构程序设计,实验四循环结构程序设计

实验四循环结构程序设计

实验四 循环结构程序设计 实验内容： 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) y<3时的最大n值。 试验程序：y=0; n=1; while y<3 y=y+1/(2*n-1); n=n+1; end n=n-1 实验结果：n = 57 (2) 与(1)的n值对应得y值。 试验程序：y=0; n=1; while yp(n) min=p(n); end if max0 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=( 请输入一个复数： ) [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，输入时以角度为单位。 要求：定义一个线性方程组AX=B的函数文件，然后在命令文件中调用该函数文件。 试验程序：函数文件：function[A,B,X]=tim(m1,m2,theta) A=[m1*cos(theta),-m1,-sin(theta),eps;m1*sin(theta),eps,cos(theta),eps;eps,m2,-sin(theta),eps;eps,eps,-cos(theta),1]; B=[eps,m1*9.8,eps,m2*9.8] ; X=inv(A)\B; 函数文件调用：m1=( m1= ); m2=( m2= ); a=( a= ); theta=(a*pi/180); [A,B,X]=tim(m1,m2,theta); A B X 实验结果：m1=2 m2=2 a=30 A = 1.8356 -2.0000 -0.3970 0.0000 0.7940 0.0000 0.9178 0.0000 0.0000 2.0000 -0.3970 0.0000 0.0000 0.0000 -0.9178 1.0000 B = 0.0000 19.6000 0.0000 19.6000 X = -39.2000 0.0000 39.2000 19.6000 实验内容： 4：设f(x)=+，编写一个MATLAB函数文件fx.m，使得调用f(x)时，x可用矩阵带入，得出的f(x)为同阶矩阵。 试验程