24小时热门版块排行榜    

查看: 2714  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lihongru

金虫 (小有名气)

[求助] matlab非线性参数拟合问题

现已知x和y的一系列实验值,且已知x和y具有如下关系:
y=kln(x-ky)+a,期中k和a是待拟合的参数,我一开始采用线性拟和的方法,先给括号中的k赋初值k0,对y和ln(x-k0)线性回归出新的k和a,然后对k进行迭代和线性拟合,直到k收敛,现在的问题是有的情况下能够收敛出k,有时不收敛,请问有没有更好的非线性回归方法,求得k和a参数呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yqx1985

木虫 (著名写手)

云中仙

【答案】应助回帖

★ ★
微尘、梦想(金币+2): 谢谢参与应助 2011-10-24 17:46:23
lihongru(金币+10): 谢谢,我试着运行一下。 2011-10-24 21:00:32
你的表达式,经过转化不就是下面的形式吗?
x=e^[(y-a)/k]+ky;
你可以在matlab中建立一个函数
function y=fun(x,p)
a=p(1);k=p(2);
y=e^((x-a)./k)+k*x;
end
进而你可以用nlinfit函数进行非参数拟合
options = statset('Robust','on');
[beta,r,J,Sigma,mse] = nlinfit(x,y,@fun,rand(1,2),options);
即可得到参数的估计值,当然,你也可以重新定义函数,采用最小二乘的方法进行参数估计。
人面不知何时去,桃花依旧笑春风
3楼2011-10-24 11:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

lihongru(金币+10): 谢谢,我对1stOpt不是特别熟悉,回头试一下。 2011-10-24 20:58:45
你的问题已变成隐函数拟合了,可以试下1stOpt,处理的很好。
2楼2011-10-24 09:36:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihongru

金虫 (小有名气)

我试着用另一种方法解了一下,就是我知道k的大概范围,然后程序如下:
k=6:0.1:15%另一种求k值的方法
g=length(k);
for i=1:g
d0=d-k(i)*s;
lnd0=log(d0);
X0=[lnd0,DT,cons];
p=X0\lns;
m=p(1)-k(i);
f=abs(m);
if f>=delt
    return
else
end
end
p
可是运行不出结果,也不知程序错在什么地方。
4楼2011-10-24 21:04:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihongru

金虫 (小有名气)

原始的要拟合的方程是lns=kln(d-k*s)+a*DT+c,s、d和DT是已知的三个列向量,要拟合k、a和c,若DT固定后就成了一个隐函数了。不知采用上边的方法可行不可行。
5楼2011-10-24 21:13:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿东华大学化学070300,求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 328求调剂,英语六级551,有科研经历 +6 生物工程调剂 2026-03-17 10/500 2026-03-22 20:22 by edmund7
[考研] 求调剂一志愿海大,0703化学学硕304分,有大创项目,四级已过 +6 幸运哩哩 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[考研] 306求调剂 +6 chuanzhu川烛 2026-03-18 6/300 2026-03-22 19:39 by 公瑾逍遥
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 生物学一志愿985,分数349求调剂 +4 zxts12 2026-03-21 7/350 2026-03-22 09:57 by zxts12
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
信息提示
请填处理意见