| 查看: 3771 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
lsqcurvefit最小二乘提示错误，望高手指点 已有1人参与
|
||
|
x =[770,850,860,960,990,1050,100,1060,1090,1090]; y =[0,0.3,0.7,1.0,1.3,1.7,2.0,2.3,0.7,3.0]; f=@(a, x) a(1)./a(2).*(a(3)./(1700+850.*a(1).*(1-a(4))) .*log((1700.*a(4)+a(1).*(1-a(4)).*(x-770))./(a(4).*a(3).*850-a(3).*(x-770)))+1./(a(1).*(1-a(4))) .*log(1700+(a(1).*(1-a(4)).*(x-770))./a(4))-a(3)./(1700+850.*a(1).*(1-a(4)) ).* log(1700/(850.*a(3) ))-1./(a(1)*(1-a(4))) .*log(1700)) [a,res]=lsqcurvefit(f,[1,1,1,1], x,y); y1=f(a, x); plot(x,y, x,y1,’x’) f = @(a,x)a(1)./a(2).*(a(3)./(1700+850.*a(1).*(1-a(4))).*log((1700.*a(4)+a(1).*(1-a(4)).*(x-770))./(a(4).*a(3).*850-a(3).*(x-770)))+1./(a(1).*(1-a(4))).*log(1700+(a(1).*(1-a(4)).*(x-770))./a(4))-a(3)./(1700+850.*a(1).*(1-a(4))).*log(1700/(850.*a(3)))-1./(a(1)*(1-a(4))).*log(1700)) 错误使用 snls (line 47) Objective function is returning undefined values at initial point. lsqcurvefit cannot continue. 出错 lsqncommon (line 149) [xC,FVAL,LAMBDA,JACOB,EXITFLAG,OUTPUT,msgData]=... 出错 lsqcurvefit (line 252) [xCurrent,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ... 本人菜鸟初学，求高手指点怎么修改啊？多谢啊 |
回复此楼» 猜你喜欢
filecode
已经有14人回复
-
关于Filecode分析方法
已经有4人回复
-
FileCode能看出啥？
已经有20人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有3人回复
-
不应该看fileCode
已经有4人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
4楼2015-12-29 18:01:28
文刀星
木虫 (小有名气)
- 应助: 13 (小学生)
- 金币: 3651.2
- 散金: 500
- 红花: 5
- 帖子: 60
- 在线: 100.9小时
- 虫号: 3693767
- 注册: 2015-02-22
- 性别: GG
- 专业: 导航、制导与传感技术
2楼2015-12-25 20:10:49
文刀星
木虫 (小有名气)
- 应助: 13 (小学生)
- 金币: 3651.2
- 散金: 500
- 红花: 5
- 帖子: 60
- 在线: 100.9小时
- 虫号: 3693767
- 注册: 2015-02-22
- 性别: GG
- 专业: 导航、制导与传感技术
3楼2015-12-25 20:28:21
