实验11 综合实验 2new

《实验11 综合实验 2new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验11综合实验班级：　　　110501　　　　姓名：刘荣林　学号：11050126　一、计算题，要求给出程序和运行结果（每小题5分，总分25分）1、计算表达式的值，其中。代码：a=1+i;b=-7.9;y=2*pi*a^2*((1-pi/4)*b-(0.8333-exp(2)/2)*a)结果：y=-35.9553+14.6508i2、求多项式的根。代码：clcclearcloseallsymsxy='x^4-2*x^3+x^2+2*x-2';m=solve(y)结果：m=-111-ii+13、已知，求、代码：clcclearcloseallx=[892965129

2、];f=x^3-x+3第7页共7页det(f)结果：f=2420225516562278212417327267361096ans=-2.3023e+0074、已知，求A的特征值及对应特征向量。代码：clcclearcloseallA=[561891558];[x,y]=eig(A)结果：x=-0.4116-0.7576-0.0700-0.60880.6495-0.1834-0.67820.06570.9805y=15.5226000-0.23050006.70795、求线性方程组的解。代码：clc第7页共7页clearcloseallA=[145321246];b

3、=[7;3;9];B=Ab结果：B=0.6250-0.12501.3750二、已知（）。求n=100时表达式的值。（15分）代码：functionmain()n=100;s=s(n)functions=s(n)a=-1;s=0;d=1/8;e=1;fori=2:nb=2*i-1;c=2*i+2;a=-a;d=a*d*(b/c);s=s+d;end第7页共7页结果：s=0.0183三、利用牛顿迭代公式解方程f(x)=cosx-x=0。取初始值x0=3.1415/4，求使f(xn)<时的xn（15分）代码：functionmainy=sbb(5);i=0;f=0.78

4、54;whilef>1*10^(-6)i=i+1;f=abs(cos(y(i))-y(i));endi,y(i)functiony=sbb(n)fori=1:nifi==1sbb(i)=3.1415/4;endsbb(i+1)=sbb(i)-(cos(sbb(i))-sbb(i))/(-sin(sbb(i))-1);y(i)=sbb(i);end结果：i=3ans=0.7391第7页共7页四、将窗口按两行两列方式分割成四个子窗口，在各窗口依次按要求画出如下图型：1）用plot在上以绿色的虚线绘制；2）用fplot在上将曲线绘在同一窗口；3）用ezplot画出曲线，；

5、4）用ezpolar绘制曲线，。（20分）代码：1）clcclearx=-pi:0.01:pi;y=cos(2*x)+sin(3*x);plot(x,y,'g--')2）clcclearsubplot(121)fplot('atan(x)',[-10,10])subplot(122)fplot('acos(2*cos(x))',[-10,10])3）ezplot('2*t^2','3*t^3',[0,7])第7页共7页4）ezplot('r=cos(2*t)*log(1+t)',[0,6*pi])五、先执行rand('seed',100);A=rand(100)>0

6、.5;将矩阵A按“田”字形均分为4块，统计各块中的非零元的个数。（15分）代码：clcclearrand('seed',100);A=rand(100)>0.5;A1=A(1:50,1:50);A2=A(51:100,1:50);A3=A(1:50,51:100);A4=A(51:100,51:100);a1=find(A1,'0');aa1=numel(a1)a2=find(A2,'0');aa2=numel(a2)a3=find(A3,'0');aa3=numel(a3)a4=find(A4,'0');aa4=numel(a4)结果：第7页共7页aa1=48aa

7、2=48aa3=48aa4=48六、在自己的电脑上安装Notebook，新建一M-book文档，运行如下程序：A=[1,2,3;1,4,9;1,8,27];b=[5,-2,6]';x=inv(A)*b将运行结果复制如下。（10分）代码：结果：第7页共7页