资源描述：
《234 葛雷+数学实验周15》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、淮阴工学院数学软件专业实践周班级：金融1131姓名：葛雷学号：1134104134教师：王小才 胡平教师王小才胡平答辩分总分得分数理学院17数学实验周题目一、代数与矩阵综合题生成3阶矩阵，其中，系数向量利用下述方法给定：；m为你的学号最后2位。(1)求的秩，行列式，的伴随矩阵。17程序代码：clccleara=fliplr(-(randi(8,3,1)-4)')A=[0,1,0;0,0,34;a]save('filename1','a','A')disp('A的秩为:')disp(rank(A))disp('A的行列式')disp(det(A))disp('A的伴随
2、矩阵')disp(inv(A)*det(A))运行结果：A=01000342-2-2A的秩为:3A的行列式68A的伴随矩阵68234680002017(2)向量由生成，求方程的解。程序代码：17b=randi(20,3,1)-4x=Ab运行结果：x=7.1765-1.00000.1765(3)求的特征多项式，计算，并对结果进行必要的讨论；17PA=poly(A)fx=poly2str(PA,'x')FA=A^3+2*A^2+68*A-68运行结果：PA=1.00002.000068.0000-68.0000fx=x^3+2x^2+68x–68FA=0-68-68-
3、680-6817-68-68017(4)确定一个区间，使该区间包含特征方程所有的根，并画出曲线的图像（显示的效果尽可能地好）；程序代码：17c=0,d=2x=-10:0.05:10y=x.^3+2*x.^2+68*x-68plot(x,y,'color','r','Linewidth',3)17(5)使用函数roots或solve求方程最小的实根或唯一的实根；程序代码：symsxx1=solve(x^3+2*x^2+68*x-68==0)x=vpa(x1)结果：x=0.95989-1.4799+8.2856*i-1.4799-8.2856*i结论：有唯一实根x=0.
4、95989(6)求的特征值和化成Jordan标准型的变换矩阵；运行代码D=1/2*(A+A')disp('D的特征值')disp(poly(D))[V,J]=jordan(D)Q=V结果：D的特征值1.00002.0000-257.2500-16.5000Q=-0.0312-0.0000i0.1014+0.0000i-32.2741-0.0000i-0.9381+0.0000i1.0627-0.0000i2.1381-0.0000i1.00001.00001.000017(7)根据计算的结果，指出矩阵具有的性质，并说明是否与你所了解的结论一致；程序：K=Q'*Q结果
5、：K=1.0e+03*0.0019-0.0000+0.0000i0.0000-0.0000i-0.0000-0.0000i0.00210.0000+0.0000i0.0000+0.0000i0.0000-0.0000i1.0472K是一个对角矩阵，一致的(8)求二元函数的具体表达式，说明该二元函数是那种类型的二次曲面，编写函数型m文件定义该二元函数；程序代码：symsxyzz=[xy1]*D*[x;y;1]结果：z=x+16*y+x*(y/2+1)+y*(x/2+16)–2类似于马鞍面的二次曲面functionz=mfun1(x,y)z=x+16*y+x.*(y./
6、2+1)+y.*(x./2+16)-2end(9)选取适当的圆形区域画出该二次曲面，并使显示效果尽可能地好。clcclf[x,y]=meshgrid(-5:0.01:5);z=mfun1(x,y);k=find(x.^2+y.^2>25);x(k)=NaN;mesh(x,y,z)17二、绘图基础题要求对图像作必要的修饰。1.用subplot命令把图像分成4块，分别绘出下面4个函数的图像①②作出极坐标方程为的对数螺线的图形.③④自己画一个分段函数的图像其中是由randint（或randi）命令产生的1到20之间的正整数。randi(20,9,1)ans=1938122
7、6111020可设为：c1=19;c2=3;c3=8;c4=12;n=2;a=6;b=12;c=10;d=20;17程序代码clcclearc1=19;c2=3;c3=8;c4=12;n=2;a=6;b=12;c=10;d=20;subplot(2,2,1)symsxytx=cos(c1*t).*cos(c2*t);y=sin(c3*t).*cos(c4*t);ezplot(x,y)subplot(2,2,2)theta=linspace(-2*pi,2*pi);h=polar(theta,exp(theta/5));set(h,'Linewidth',2,'c
