欢迎来到天天文库
浏览记录
ID:43442247
大小:21.19 KB
页数:3页
时间:2019-10-03
《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()函数残差更小,即得到了更好的拟合效果。
此文档下载收益归作者所有