计算物理基础第一章习题(共6页).docx
精选优质文档-倾情为你奉上17.用fplot绘制曲线cos(tanx)。答案：fplot('cos(tan(pi*x)',0.01 10.1,1e-3) 。18用ezplot绘制曲线ey+cosxx+y=0。答案：ezplot('exp(y)+cos(x)/x+y')。19.用ezsurf绘制曲面。x=e-scosty=e-ssintz=t (0s8,0t5)Answer: ezsurf('exp(-s)*cos(t)','exp(-s)*sin(t)','t',0,8,0,5*pi)21.画曲面和等高线z=x2y+yx (-5x5, 0y10)的图像答案: x=linspace(-5,5,100);y=linspace(0,10,100);X,Y=meshgrid(x,y); %1.构造数据网格Z=X.2.*Y+sqrt(X)./X; %2.建函数 surf(X,Y,Z) CONTOUR(X,Y,Z)CONTOURF(X,Y,Z);26. 画出牛顿环的干涉图样。f.m在光学上，牛顿环是一个薄膜干涉现象。光的一种干涉图样，是一些明暗相间的同心圆环。在加工光学元件时，广泛采用牛顿环的原理来检查平面或曲面的面型准确度。%program 画出牛顿环的干涉图样。f.mlambda=589.3e-9;R=0.8551;x=-3:0.01:3;y=x; % 设定两轴的范围及间隔X,Y=meshgrid(x, y);r2=X.2+Y.2;% 产生n×n个网格点的坐标B=2*cos(2*pi*r2*1.e-6/R/lambda); % 光强的计算colormap(gray(255)W=B*255/5;image(x,y,W)title('牛顿环干涉图样');xlabel('x(mm)');ylabel('y(mm)');27. 画出等容、等压、等温曲线。f.m%program 画出等容、等压、等温曲线。f.m%由理想气体状态方程PV=nRT，令nR=1%等温图，令T=1figuresubplot(1,3,1)P=0.001:0.01:10;V=1./P;plot(P,V);xlabel('P');ylabel('V');set(gca,'xlim',0 10,'ylim',0 10)title('等温图(T=1)','fontsize',15)%等压图，令P=1subplot(1,3,2)V=0:10;T=V;plot(V,T);xlabel('V');ylabel('T');set(gca,'xlim',0 10,'ylim',0 10)title('等压图(P=1)','fontsize',15)%等容图，令V=1subplot(1,3,3)P=0:10;T=P;plot(P,T);xlabel('P');ylabel('T');set(gca,'xlim',0 10,'ylim',0 10)title('等容图(V=1)','fontsize',15)28. 画出表面水波的干涉图样。f.m%program 画出表面水波的干涉图样 f.m%缺陷：理论应标注a=4;w=4;A0=1;A1=1;k=1;x,y=meshgrid(-4*pi:pi/20:4*pi,-4*pi:pi/20:4*pi);r1=sqrt(x-a).2+y.2)+5;c1=1./r1;r2=sqrt(x+a).2+y.2)+5;c2=1./r2;t=0:0.1:10figurefor j=1:91 z=A0*c1.*cos(w*(t(j)-k*r1)+A1*c2.*cos(w*(t(j)-k*r2); surf(x,y,z); colormap(gray); axis equal; shading interp; view(50,36); axis off; title('剩余传播时间(秒)：',num2str(fix(91-j)*0.2),. 'fontsize',20,'Color','r') pause(0.2)end专心-专注-专业
