河北工业大学材料的模拟与计算上机题

1、第一次：产生10个随机数，计算其标准差和平均数r=rand(1,100)a=mean(r)b=var(r)c=std(r)第二次：判断立方晶系晶体结构。已知X射线衍射峰数据(20度)以下： 2theta=38.52,44.76,65.14,78.26,82.47,99.11,112.03,116.6,137.47,若 wavelength=0.154056nm，求晶格 常数。l=38.52,44.76,65.14,78.26,82.47,99.11,112.03,116.6,137.47b=l*pi/360c=sin(b).”2d=sin(b)e=b/b(1)%经过验证与面心立方结构相符合f=

2、0.154056a=sqrt(3)*f/2/d(1)第三次：生成等差向量X，y=cos(x)，利用直线拟合，求出拟合方程，作图，并计算最大误差x=linspace(0,pi/4,100)y=sin(x)n=100 x1=sum(x)y1=sum(y)x11=sum(x.”2)y12=sum(x.*y)a=(y1*x11-x1*y12)/(n*x11-x2)b=(n*y12-y1*x1)/(n*x11-x2)plot(x,y)hold on;y2=a+b*xplot(x,y2)t=max(abs(y2-y)第四次：已知某温度的扩散系数，求扩散激活能。T(摄氏度)800870 8301027D(m

3、2/s)*10(-12)2.51 8.01 17.6 51.5t=800;870;930;1027T=t+273d=2.51;8.01;17.6;51.5D=d*10(-12)y=log(D)x=1./Tx1=sum(x)y1=sum(y)x11=sum(x.”2)y12=sum(x.*y)a=(y1*x11-x1*y12)/(4*x11-x2)b=(4*y12-y1*x1)/(4*x11-x2)y2=a+b*xplot(x,y2)Q=-b*8.314第五次:半径为1的球，四分之一的圆面积，求兀产生X，算出y= (1-x2) 0.5,取平均数为面积 x=rand(1,10000)y=sqrt(

4、1-x.2)a=mean(y)p=a*4产生(x，y)，看落在圆内的几率clears=0;for a=1:1:10000 x=rand(1);y=rand(1);r(a)=x2+y2;if r(a) 0;i=i+1a=G*M/(R+h)2;v(i+1)=a*t+v(i);h=h-v(i)*t;plot(i,h);hold onend第八次：已知在1000 (摄氏度)时Zn在Cu中的扩散系数为5*10(-13)m2/s，将纯铜和黄铜(Cu-30Zn)对焊在 一起在1000摄氏度加热扩散，计算扩散时界面附近的Zn浓度分布，可认为距离界面2mm以为的浓度不受扩散 影响x=0.001/1000;t=0

5、.1;f=5*10(-13)*t/x2;c(501:1000)=0.3;c(1:500)=0;a(1)=0;a(1000)=0.3;for n=1:100000for i=2:999a(i)=f*c(i+1) + (1-2*f)*c(i)+f*c(i-1);end c=a;endplot(c,+)熔点（摄氏度）10831453熔点（摄氏度）10831453Cu 13.26Ni 17.47求出不同温度下的平衡相浓度，进而绘制相图。A=13260/8.314B=17470/8.314T1=1083+273T2=1453+273for T=T1:T2m=A*(1/T-1/T1)n=B*(1/T-1/T2)xAs=exp(m)*(1-exp(n)/(exp(m)-exp(n)