Matlab在材料力学的应用

《Matlab在材料力学的应用》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、Matlab在材料力学中的应用在工业生产中，为研究某种材料应力与应变的关系，假设我们测得一组数据如下表：应力σ925112516252125262531253625应变ε0.110.160.350.480.610.710.85如果假定应力与应变有如下关系（σ为应力值，ε为应变值）：ε=a+blnσ试计算a、b的值。MATLAB的表达形式如下：x=[925,1125,1625,2125,2625,3125,3625];y=[0.11,0.16,0.35,0.48,0.61,0.71,0.85];plot(x,y,'o'

2、)[p,resid1]=polyfit(x,y,2)holdonxi=linspace(700,3700,3000);yi=polyval(p,xi);plot(xi,yi)x0=[0.1,0.1];fff=inline('a(1)+a(2)*log(x)','a','x');[a,resid2]=lsqcurvefit(fff,x0,x,y)plot(xi,fff(a,xi),'r')执行程序得到图3，图中蓝色曲线为利用polyfit()函数得到的曲线，红色曲线为利用lsqcurvefit()函数得到的曲线；其显示

3、的结果为：p=-0.00000.0004-0.2266resid1=R:[3x3double]df:4normr:0.0331a=-3.58100.5344resid2=0.0064其中a的值代表利用lsqcurvefit()函数得到的关系为：ε=-3.5810+0.5344+σresid1、resid2分别代表运用polyfit()函数、lsqcurvefit()函数得到的残差。可以看出利用lsqcurvefit()函数残差更小，即得到了更好的拟合效果。