MATLAB第三次上机作业.docxVIP

MATLAB第三次上机作业数学与统计学院071211班时间：2014-4-14晚上6:00 ~ 10:00 205机房第一部分：基本绘图在同一坐标系下面画出在区间上的曲线图。x=(-2*pi:0.1:2*pi);y1=sin(x);?y2=cos(x);?y3=(exp(x)).*sin(5.*x);?plot(x,y1,g:);?hold?on;?plot(x,y2,r:);?plot(x,y3,b:)?运行结果：绘制三维螺旋线, t=(0:0.1:pi); x=2*cos(t); y=2*sin(t); z=0.5*t; plot3(x,y,z)运行结果：画出曲面的网线图。?t=-10:0.01:10;?[x,y]=meshgrid(t,t);?z=sin(x.*y)./(x.*y+eps);??mesh(x,y,z)运行结果：画出曲面的图形。t=-10:0.01:10; [x,y]=meshgrid(t,t); z=x.*exp(-(x.^2+y.^2)); mesh(x,y,z)运行结果：作出下列曲面的3维图形，环面, 。t=-10:0.01:10;?????[x,y]=meshgrid(t,t);?????z=sin(pi.*sqrt(x.^2+y.^2));?????mesh(x,y,z)?t=[0:0.1:2*pi];?[u,v]=meshgrid(t,t);?x=(1+cos(u)).*cos(v);?y=(1+cos(u)).*sin(v);?z=sin(u);?mesh(x,y,z)运行结果：产生一个4阶的随机矩阵，执行下面的操作：求其行列式，检验其是否可逆；若可逆，求其逆矩阵。计算该矩阵的特征值、特征向量。将该矩阵化为行最简的阶梯形。验证矩阵的特征值之和等于矩阵主对角元之和，特征值之积等于矩阵的行列式。a=rand(4) %产生一个4届随机矩阵；hls=det(a) %(1)计算行列式；njz=inv(a) %(1)计算逆矩阵；[tzxl,tzz]=eig(a) %(2)计算特征向量和特征值；Jtx=rref(a) %(3)化为行最简单的阶梯形；zzs=[];fori=1:4zzs=[zzs,det(a(1:i,1:i))];endzzs %(5)求顺序主子式tzzh=sum(eig(a)) %（5）求特征值之和；zdjh=trace(a) %(5)求主对角元之和；tzzj=prod(eig(a)) %(5)求特征值之积运行结果：a = 0.814723686393179 0.632359246225410 0.957506835434298 0.957166948242946 0.905791937075619 0.097540404999410 0.964888535199277 0.485375648722841 0.126986816293506 0.278498218867048 0.157613081677548 0.800280468888800 0.913375856139019 0.546881519204984 0.970592781760616 0.141886338627215hls = -0.026140719644426njz = -15.299695954719372 3.076060796126461 14.723454499731773 9.644517041770815 -0.208808422147492 -1.844235608051758 1.036552712190957 1.87106586073514414.569368265585791 -1.933730207042300 -14.649728277400790 -9.041327053613701 -0.369008852162959 0.534535526584963 1.437778589382914 -0.400838976479567tzxl = -0.662059835624771 -0.714947788081498 0.174462593093463 0.182078143338004 -0.481933427298575 -0.029153849852601 -0.629136472171865 -0.928766659682403 -0.276571160291259 0.697193988674131 0.599472951677392 0.317787492557903 -0.502916831096