原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
matlab应用与建模方法2
问题化为：由实验数据来识别参数 记： 659 650 639 626 610 590 565 535 499 454 Cb 1000 900 800 700 600 500 400 300 200 100 t function f=fun507(x,t) f=x(1)+x(2)*exp(-0.02*x(3)*t); t=100:100:1000; cb=10^(-5)* [454,499,535,565,590,610,626,639,650, 659]; x0=[0.05,-0.005,0.05]; x=lsqcurvefit(‘fun507’,x0,t,cb) x = 0.0070 -0.0030 0.1011 例5.8 在某次阻尼振荡实验中测得18组数据点， 试确定其振动方程。 t：0 0.4 1.2 2 2.8 3.6 4.4 5.2 6 y: 1 0.85 0.29 -0.27 -0.53 -0.4 -0.12 0.17 0.28 t：7.2 8 9.2 10.4 11.6 12.4 13.6 14.4 15 y: 0.15 -0.03 -0.15 -0.07 0.06 0.08 0.032 -0.015 -0.02 由物理背景知： 采用自定义函数拟合 clc,clear x=[0,0.4,1.2,2,2.8,3.6,4.4,5.2,6,7.2,8,9.2,10.4,11.6,12.4,13.6,14.4,15]; y=[1,0.85,0.29,-0.27,-0.53,-0.4,-0.12,0.17,0.28,0.15,-0.03,-0.15,-0.07,0.059,0.08,0.032,-0.015,-0.02]; f=fittype(a*cos(k*x)*exp(w*x),independent,x,coeff,{a,k,w}); cfun=fit(x,y,f) xi=0:0.1:20; yi=cfun(xi); plot(x,y,r*,xi’,yi’,b-) 运行结果： General model: cfun(x) = a*cos(k*x)*exp(w*x) Coefficients (with 95% confidence bounds): a = 0.9988 (0.9837, 1.014) k = 1.001 (0.9959, 1.006) w = -0.2067 (-0.2132, -0.2003) L与下垂度x的关系： clc,clear a=100;m=1200;g=9.78;L=125; w=m*g/a; x=fzero(@(x)a*(1+(8/3)*(x/a)^2-(32/5)*(x/a)^4+(256/7)*(x/a)^6)-L,50) T=w*a/2*sqrt(1+(a/x)^2) 二、 利用插值与拟合方法的数据建模 已知函数y=f(x)的一数据集(xi,yi), i=1,2,…,m 而函数表达式未知，要从某函数类（如多项式函数）中求得一个函数F(x)作为f(x)的近似，这类数值计算问题称为数据建模。 常用方法：（1）插值（2）拟合 （一）插值法 已知函数y=f(x)在区间[a,b]上的n+1个不同点 的函数值 若存在一个简单函数F(x), 使 称F(x)为f(x)在区间[a,b]上的插值函数。 求差值函数的方法称为插值法，若F(x)为多项式，相应的问题称为多项式差值或代数插值； 若F(x)为分段函数，称为分段插值； 若F(x)为三角函数，称为三角插值。 n次代数插值： F(x)=a0xn+a1xn-1+…+a1x+a0 常用的代数插值方法有：拉格朗日插值，牛顿插值。 当插值次数过高时，常常会出现龙格现象， 因此在数据建模中不建议使用过高次的代数插值。 在区间[-1,1]上的代数插值 (2) S(x)在每个小区间 上是三次多项式， 为了避免龙格现象，常常使用分段低次插值多项式进行插值。 一般的分段低次插值具有收敛性和稳定性，但光滑性较差，样条差值不仅有具有收敛性和稳定性还具有很好的光滑性。是较为理想的插值方法。 设函数y=f(x)
文档评论（0）