赵志斌2012021241第二次作业试题.doc

2013—2014学年第二学期2012级物理学 专业计算物理基础上机安排 第二次练习 姓名：赵志斌 学号：2012021241 绘制函数表达式的二维图形。 答案：ezplot(x^2-y^4) 绘制误差函数的二维图形。 答案： Syms x ezplot(erf(x)) grid 在极坐标下绘制函数表达式的二维图形。 答案：ezpolar(1+cos(r)) 根据表达式，和，绘制三维曲线。 答案：t=0:0.01:15; x= sin (t); y=cos(t); z=t;plot3(x,y,z,*) plot3(x,y,z,*) 练习课本44页的图1.15和1.16，并对图1.16进行涂色。 答案：subplot(1,3,1) sphere (40); subplot(1,3,2) [X,Y]=meshgrid(-8:0.5:8); surf(X,Y,Z) subplot(1,3,3) A=[0,0,0]; B=[1,0,0]; C=[0.5,0.5,0.5]; D=[0.5,0.5,0]; Xabc=[A(1),B(1),C(1),A(1)]; Yabc=[A(2),B(2),C(2),A(2)]; Zabc=[A(3),B(3),C(3),A(3)]; Xbcd=[B(1),C(1),D(1),B(1)]; Ybcd=[B(2),C(2),D(2),B(2)]; Zbcd=[B(3),C(3),D(3),B(3)]; Xacd=[A(1),C(1),D(1),A(1)]; Yacd=[A(2),C(2),D(2),A(2)]; Zacd=[A(3),C(3),D(3),A(3)]; view(-44,-22) hold on fill3(Xbcd, Ybcd, Zbcd,y) fill3(Xacd, Yacd, Zacd,r) fill3(Xabc, Yabc, Zabc,g) 练习课本47页图1.18。 答案：[x,y,z]=meshgrid(-2:.2:2); v=x.*exp(-x.^2-y.^2-z.^2); slice(v,[5 15],15,10) hold on colorbar(horiz) colorbar(vert) view([-25 65]) 练习课本48、49页复数函数的图形。 subplot(1,3,1) subplot(1,2,1) u=cplxgrid(20); cplxmap(u,u.^2) colorbar(vert) subplot(1,2,2) u=cplxgrid(20); cplxmap(u,u.^2) colorbar(vert) 参照课本56页的等值面表现三维标量场，做出不同颜色，不同表面透明度的图形。 答案：[x y z v] = flow; p=patch(isosurface(x,y,z,v,-3)); isonormals(x,y,z,v,p) set(p,FaceColor,red,EdgeColor,none); daspect([1 1 1]) daspect([1 1 1]) view(3) axis tight; grid on camlight;lighting phong alpha(.5);