哈工大 数学实验报告(完整版)

《哈工大 数学实验报告(完整版)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验一Matlab的使用1.输入A,B>>A=[-3508;1-82-1;0-593;-70-45]A=-35081-82-10-593-70-45>>B=[02-16]'B=02-16(a)方法>>ABans=-0.6386-0.4210-0.35290.0237(b)方法>>inv(A)*Bans=-0.6386-0.4210-0.35290.02372.分块矩阵>>E=eye(3);>>R=rand(3,2);>>O=zeros(2,3);>>S=diag([12]);>>A=[ER;OS]A=1.0000000.81470.913401.000000.90580.6324001

2、.00000.12700.09750001.0000000002.0000>>A2=[ER+R*S;OS^2];>>err=A^2-A2>>norm(err)ans=03.多项式>>r=[2-31+2*i1-2*i0-6];>>p=poly(r)p=15-9-172-1800计算0.8�-1.2处的值�>>y1=polyval(p,0.8)y1=-100.2179>>y1=polyval(p,-1.2)y1=293.29004.特征方程和友元矩阵>>p=[10-63-8];>>A=compan(p)A=06-38100001000010求A的特征值>>lamda=eig(A)lamda=

3、-2.83742.46920.1841+1.0526i0.1841-1.0526iroots(p)结果>>roots(p)ans=-2.83742.46920.1841+1.0526i0.1841-1.0526i结论�一维向量p友元矩阵的特征值即为其降幂多项式的根5.做函数�����������的图形>>x=linspace(0,1);>>X=[x.^2;x.^3;x.^4;x.^5];>>plot(x,X(1,:),'r',x,X(2,:),'g',x,X(3,:),'b',x,X(4,:),'k')>>legend('y=x^2','y=x^3','y=x^4','y=x^5')用s

4、ubplot画四个函数图形的指令>>subplot(2,2,1);plot(x,X(1,:));title('y=x^2');>>subplot(2,2,2);plot(x,X(2,:));title('y=x^3');>>subplot(2,2,3);plot(x,X(3,:));title('y=x^4');>>subplot(2,2,4);plot(x,X(4,:));title('y=x^5');6.fplot的使用>>fplot('funfplot',[-0.10.1],2e-4)7.用作图法求�����������的根的近似值>>x=linspace(-6,4);>>y=4*

5、sin(x)-x-2;>>plot(x,y)>>[y0,n]=find(abs(y)<0.05);>>plot(x,y,x(n),y(n),'ro')对x,y轴加标志并加网格线�>>xlabel('xaxis')>>ylabel('yaxis')>>gridon8.作曲面�������的三维图形>>[X,Y]=meshgrid(-5:0.1:5);>>Z=X.^2-Y.^2;>>meshc(Z)或使用surfc(Z)9.建立m函数(a)自然数阶乘functiony=Factorial(x)ifround(x)~=x

6、x<0error('x必须是非负整数');endifx==0y=1;el

7、sey=prod(1:x);end(b)n中取m的组合functiony=Combination(n,m)ifround(m)~=m

8、round(n)~=nerror('m,n必须为整数');elseifn<1

9、m<0

10、n

11、(x)y0(i)=feval(inline(y),x(i));y1(i)=feval(inline(dy),x(i));y2(i)=feval(inline(ddy),x(i));endplot(x,y0,x,y1,x,y2)legend('y','dy','ddy')11.绘制极坐标曲线functionPolarPlot(a,b,n)theta=linspace(0,2*pi);r=a*cos(b+n*theta);polar(t