实验四循环结构程序设计

试验四循环构造程序设计试验内容：n2 1 1 1 1n1：依据 =6 2

+

+

+……+

，求π的近似值当n分别取2100,1000,10000时,结果是多少?要求:分别用循环构造和向量运算(sum函数)来实现.循环构造：试验程序：(1)y1=0;n=100;fori=1:ny1=y1+1/i/i;endy=sqrt(6*y1)试验结果：y=3.1321试验程序：(2)y1=0;n=1000;fori=1:ny1=y1+1/i/i;endy=sqrt(6*y1)试验结果：y=3.1406试验程序：y1=0;n=10000;fori=1:ny1=y1+1/i/i;endy=sqrt(6*y1)试验结果：y=3.1415向量运算：试验程序：n=100;i=1:n;f=1./i.^2;y1=sum(f);y=sqrt(6*y1)试验结果：y=3.1321i=1:n;f=1./i.^2;y1=sum(f);y=sqrt(6*y1)试验结果：y=3.1406i=1:n;f=1./i.^2;y1=sum(f);y=sqrt(6*y1)试验结果：y=3.1415试验内容：11 12: y=1+3+5+……+2n1,求：y<3n值。y=0;n=1;whiley<3y=y+1/(2*n-1);n=n+1;endn=n-1试验结果：n=57与〔1〕ny值。试验程序：y=0;n=1;whiley<3y=y+1/(2*n-1);n=n+1;endy=y-1/(2*(n-1)-1)试验结果：y2.9944试验内容：

f11f20

n1n2ff31fnf(n)2f(n)f(n)f1~f100中：最大值，最小值，各数之和。试验程序：p(1)=1;p(2)=0;p(3)=1;forn=4:100p(n)=p(n-1)-2*p(n-2)+p(n-3);

n3n3endmin=p(1);max=p(1);sum=0;forn=1:100ifmin>p(n)min=p(n);endifmax<p(n)max=p(n);end

endsum=sum+p(n);minmaxsum试验结果：min=-8.9941e+011max=4.3776e+011sum=-7.4275e+011正数，零，负数的个数。试验程序：p=[];p(1)=1;p(2)=0;p(3)=1;forn=4:100p(n)=p(n-1)-2*p(n-2)+p(n-3);endzero=0;negative=0;positive=0;forn=1:100ifp(n)>0positive=positive+1;elseifp(n)<0negative=negative+1;end

elseend

zero=zero+1;positivenegativezero试验结果：positive=49negative=49zero=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+2ix=1.0000+2.0000ia=-1.1312+2.4717ib=0.3495+0.4808ic=3.1658+1.9596id=2.0327-3.0519i试验内容：2：一物理系统可用以下方程组来表示：m1cos

m1 si0 1sin 0 cos 0a2 g = 0 m20 0

sico

0a3 0 1a4 m2g m1，m2和a1，a2，N1N2g9.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=input(”m1=”);m2=input(”m2=”);a=input(”a=”);theta=(a*pi/180);[A,B,X]=tim(m1,m2,theta);ABX试验结果：m1=2m2=2a=30A=1.8356-2.0000-0.39700.00000.79400.00000.91780.00000.00002.0000-0.39700.00000.00000.0000-0.91781.0000B=0.000019.60000.000019.6000X=-39.20230.000039.202319.6000试验内容：x设f= 1 x( )20.1

1(x)2

MATLABfx.m，使得调用f〔x〕时，x可用矩阵带入，得出的f〔x〕为同阶矩阵。试验程序：函数文件：function[Y]=f(x)Y=1./((x-2).^2+0.1)+1./((x-3).^4+0.01