河北工業(yè)大學(xué)材料的模擬與計(jì)算上機(jī)題

1、第一次：產(chǎn)生10個(gè)隨機(jī)數(shù)，計(jì)算其標(biāo)準(zhǔn)差和平均數(shù)r=rand(1,100)a=mean(r)b=var(r)c=std(r)第二次：判斷立方晶系晶體結(jié)構(gòu)。已知X射線衍射峰數(shù)據(jù)(20度)以下： 2theta=38.52,44.76,65.14,78.26,82.47,99.11,112.03,116.6,137.47,若 wavelength=0.154056nm，求晶格 常數(shù)。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)%經(jīng)過(guò)驗(yàn)證與面心立方結(jié)構(gòu)相符合f=

2、0.154056a=sqrt(3)*f/2/d(1)第三次：生成等差向量X，y=cos(x)，利用直線擬合，求出擬合方程，作圖，并計(jì)算最大誤差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)第四次：已知某溫度的擴(kuò)散系數(shù)，求擴(kuò)散激活能。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的球，四分之一的圓面積，求兀產(chǎn)生X，算出y= (1-x2) 0.5,取平均數(shù)為面積 x=rand(1,10000)y=sqrt(

4、1-x.2)a=mean(y)p=a*4產(chǎn)生(x，y)，看落在圓內(nèi)的幾率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 (攝氏度)時(shí)Zn在Cu中的擴(kuò)散系數(shù)為5*10(-13)m2/s，將純銅和黃銅(Cu-30Zn)對(duì)焊在 一起在1000攝氏度加熱擴(kuò)散，計(jì)算擴(kuò)散時(shí)界面附近的Zn濃度分布，可認(rèn)為距離界面2mm以為的濃度不受擴(kuò)散 影響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,+)熔點(diǎn)（攝氏度）10831453熔點(diǎn)（攝氏度）10831453Cu 13.26Ni 17.47求出不同溫度下的平衡相濃度，進(jìn)而繪制相圖。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)