求助：函数拟合！

jnxh · 发表于 2006-8-22 10:55

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

main.m
tdata=[0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 14 15 16]
ydata=[30 68 75 82 82 77 68 68 58 51 50 41 38 35 28 25 18 15 12 10 7 7 4]
x0=[10,10,10]
[x,resnorm]=curvefit('Yt',x0,tdata,ydata)

Yt.m
function F = myfun(x,tdata)
F=((2*x(3)/(x(2)-x(1)))*(exp(-x(1)*tdata)-exp(-x(2)*tdata)))/48

这个程序运行总出错，哪位高手可以帮忙看一下啊

[ 本帖最后由 cdwxg 于 2006-8-22 15:19 编辑 ]

dingd · 发表于 2006-8-22 12:08

初值给的不对，公式中有：x(3)/(x(2)-x(1))，x(1)和x(2)都给初值10，会造成分母为0，自然无法计算。做拟合等优化问题，还是1stOpt最好，不用猜初值

结果：
x3          5005.10308270522
x2          2.00793803886629
x1          0.185501959221012

微微 · 发表于 2006-8-22 15:12

恩，高手阿

happy · 发表于 2006-8-23 20:13

dingd是不是1stOpt的业务员阿？怎么到处推销1stOpt

happy · 发表于 2006-8-23 20:27

不过话说回来，这个拟合用matlab是能做出来，不过效果不如1stOpt好

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 求助：函数拟合！

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。