MATLAB第3次上机作业.docxVIP

MATLAB第三次上机作业 数学与统计学院 071211班 时间：2014-4-14 晚上6:00 ~ 10:00 205机房 ? 第一部分：基本绘图 1. 在同一坐标系下面画出y=sinx, y=cosx, y=exsin5x在区间上[-2π,2π]的曲线图。 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:)? 运行结果： 2. 绘制三维螺旋线x=2costy=2sintz=0.5t, t∈[0, π] t=(0:0.1:pi); x=2*cos(t); y=2*sin(t); z=0.5*t; plot3(x,y,z) 运行结果： 3. 画出曲面z=sin?(xy)xy的网线图。 ?t=-10:0.01:10; ?[x,y]=meshgrid(t,t); ?z=sin(x.*y)./(x.*y+eps); ??mesh(x,y,z) 运行结果： 4. 画出曲面z=xe-(x2+y2)的图形。 t=-10:0.01:10; [x,y]=meshgrid(t,t); z=x.*exp(-(x.^2+y.^2)); mesh(x,y,z) 运行结果： 5. 作出下列曲面的3维图形， z=sin?(πx2+y2) 环面x=1+cosucos?vy=1+cosusinvz=sinu,u∈[0,2π]v∈[0,2π] 。 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) 运行结果： 6. 产生一个4阶的随机矩阵，执行下面的操作： 求其行列式，检验其是否可逆；若可逆，求其逆矩阵。 计算该矩阵的特征值、特征向量。 将该矩阵化为行最简的阶梯形。 验证矩阵的特征值之和等于矩阵主对角元之和，特征值之积等于矩阵的行列式。 a=rand(4) %产生一个4届随机矩阵； hls=det(a) %(1)计算行列式； njz=inv(a) %(1)计算逆矩阵； [tzxl,tzz]=eig(a) %(2)计算特征向量和特征值； Jtx=rref(a) %(3)化为行最简单的阶梯形； zzs=[]; for i=1:4 zzs=[zzs,det(a(1:i,1:i))]; end zzs %(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.141886338627215 hls = -0.026140719644426 njz = -15.299695954719372 3.076060796126461 14.723454499731773 9.644517041770815 -0.208808422147492 -1.844235608051758 1.036552712190957 1.871065860735144 14.569368265585791 -1.933730207042300 -14.649728277400790 -9.041327053613701 -0.369008852162959 0.534535526584963 1.437778589382914 -0.400838976479567 tzx