计算物理(彭芳麟)课后部分答案
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
P96
8、t=[0:0.1:2*pi]; A1=5; A2=3; w1=2; w2=4;
x1=A1*sin(w1*t+pi/3); x2=A2*sin(w2*t+pi/4); plot(t,x1,'-r',t,x2,'-b');
1
2
3
4
5
6
7
-5-4-3-2-101234
5
P98 21 x=-5:0.1:5; y=0:0.1:10;
[X,Y]=meshgrid(x,y); z=X.^2.*Y + sqrt(Y)./X;
mesh(X,Y,z);
-5
5
5
10
-1000
100
200
300
-2
2
-2
2
-2
02
4
-2
2
-20
2
00.20.40.60.8
1
P97 20
subplot(1,2,1);
[X0,Y0,Z0]=sphere(20); X=2*X0;Y=3*Y0; Z=4*Z0+1; surf(X,Y,Z); axis equal subplot(1,2,2) t=-1:0.1:1;
[X,Y,Z]=cylinder(1+t.^2,20);%形成旋转曲面 surf(X,Y,Z); P195 1
x=[-1.0 -0.75 -0.50 -0.25 0 0.25 0.50 0.75 1.00];
y=[-0.2209 0.3295 0.8826 1.4392 2.0003 2.5645 3.1334 3.7061 4.2836]; a=polyfit(x,y,1);
a= 2.2516 2.0131
y=2.2516x+2.0131
P151 3
syms v
vp=1578;
f=int(4/pi^(1/2)*v^2/vp^3*exp(-v^2/vp^2),v,0,Inf);
VF=vpa(f);