北京科技大学-材料成型自动控制基础大作业.docx

《北京科技大学-材料成型自动控制基础大作业.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、材料成形自动控制基础作业一、回归方程系数a,b1.Matlab程序：clcdisp('函数为y=ax^b');disp('取对数lny=lna+blnx')x=[0.111.92.83.74.66.47.38.29.11010.911.812.713.614.515.416.317.2]y=[0.2.418.7164264.73169157.9271316.9832912.91390.4552015.6522812.8313804.45012.66460.078173.03110175.03124

2、90.7615145.291816421573.9]c=log(x);d=log(y);sumc=sum(c);sumd=sum(d);n=length(c);e=c.*d;sume=sum(e);f=c.*c;sumf=sum(f);a1=(sume-sumc*sumd/n)/(sumf-(sumc)*(sumc)/n);c0=(sumc)/n;d0=(sumd)/n;b=a1;a0=d0-a1*c0;a=exp(a0);x1=0.01:0.1:21;y1=a.*x1.^b;plot(x,y,'

3、bp',x1,y1,'g-')legend('数据点','y=ax^b拟合曲线')disp('拟合得到的a和b的值为：')ab2.拟合结果拟合得到的a和b的值为：a=2.4001b=3.2000一、计算轧制压力1.Matlab程序clcC=[0.04160.04410.03670.03670.0330.0290.0290.0330.0290.00170.00130.00190.00190.0019]Mn=[0.250.230.260.260.260.250.250.260.250.150.140.

4、140.140.14]Cr=[0.030.030.020.020.040.040.040.040.040.020.020.020.020.02]V=[0.020.030.020.020.020.010.010.020.010.020.020.030.030.03]T=[948.4950.18944.26958.31954.99947.67954.43960.78950.46967.32960.46967.41954.13956.74]S=[1.391.381.331.2641.3721.3711.

5、3721.3131.3581.4521.5381.4691.4741.52]r=[0.41830.42960.4280.42480.40720.41150.41140.40850.41230.3930.36390.39470.39550.3974]H=[43.30346.6546.36145.92445.90845.94346.11346.04546.11546.4145.96446.28446.3646.317]D=[787.9787.9787.9787.9787.9787.9787.9787.

6、9787.9787.9787.9787.9787.9787.9]Nb=0;[C]=C+Mn/6+(Cr+V+Nb)/12;h=H.*(1-r);R=D./2;dh=H-h;T_=(T+273.15)./1000;f_h=sqrt(r./(1-r));e=(2./sqrt(3)).*log(H./h);e_r=e*1000.*S./(sqrt(R.*dh));ifT_>=0.95.*(([C]+0.41)./([C]+0.32))sigma_0=2.75.*exp(5./T_-0.01./([C]+

7、0.05));m=(-0.019.*[C]+0.126).*T_+(0.075.*[C]-0.05);elsesigma_0=2.75.*(30.*([C]+0.9).*(T_-0.95.*([C]+0.49)/([C]+0.42)).^2+([C]+0.06)./([C]+0.09)).*exp(([C]+0.32)./(0.19.*([C]+0.41))-0.01./([C]+0.05));m=(0.081.*[C]-0.154).*T_-0.019.*[C]+0.207+0.027./([C

8、]+0.32);endsigma_f=sigma_0.*(1.3.*(5.*e).^(0.41-0.77.*[C])-1.5.*e).*(e_r/10).^(m);fai=tan(pi./8.*log(1-r)./sqrt(R./h)+0.5.*atan(f_h))./sqrt(R./h);Qp=0.5.*pi./f_h.*atan(f_h)-pi/4-sqrt(R.*h)./f_h.*(log(2.*R.*(1-cos(fai))./h+1)+0.5.*log(1-r));for