| 查看: 1882 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
quan_zhang银虫 (小有名气)
|
[交流]
【求助】matlab用lsqnonlin非线性拟合的问题【已解决】 已有1人参与
|
||
|
大家好啊，我有个一个利用matlab的最小二乘法来求解参数的问题，困扰了我好久了。 具体情况如下：我有一个函数f（x）=346/15975*sin(4*pi*x/3.195)*pi-3196/3195*sin(2*pi*x/3.195)*pi且x=[0,1],另外一个带参数的函数为  且x=[0,1], a1+a2+a3=1，且a1,a2,a3=(0,1),我想用lsqnonlin来拟合。 我的做法如下: 1.先建立函数文件，并保存为myfunlin.m function y=myfunlin(x,xdata) xdata=[0 0.3195 0.639 0.9585 1.278 1.43775 1.5975 1.75725 1.917 2.2365 2.556 2.8755 3.195]; ydata=[0 -1.782446806 -2.948772499 -3.028762147 -1.911872775 -1.011104195 0 1.011104195 1.911872775 3.028762147 2.948772499 1.782446806 1E-15]; y=4.3688*((x(1)*(xdata-x(2)))./((xdata-x(2)).^2+x(3).^2)+(x(4))*(xdata-x(5))./((xdata-x(5)).^2+x(6).^2)+(1-x(1)-x(4))*(xdata-x(7))./((xdata-x(7)).^2+x(8).^2)); 2.然后调用优化程序： x0=[0.12 3.32 0.16 0.12 -0.13 0.17 0.16 1.2]; x=lsqnonlin('myfunlin',x0) 得到结果 -1.4874 2.8738 0.9386 -1.7023 0.3343 0.9831 1.5173 2.0017此时ai的范围不能控制，那么应该如何控制ai的范围呢？ [ Last edited by nono2009 on 2010-10-31 at 07:24 ] |
回复此楼» 猜你喜欢
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
今天务委会开完了，明天出结果吗
已经有11人回复
-
有没有大神帮我看看基金代码
已经有9人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
放榜前的不淡定
已经有18人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- lsqnonlin函数拟合微分方程组参数拟合问题
已经有10人回复
- matlab非线性参数拟合问题
已经有7人回复
- 50金币求一MATLAB的拟合程序
已经有11人回复
- MATLAB编程，非线性拟合。
已经有4人回复
- 【求助】急求关于matlab做非线性拟合的函数
已经有8人回复
- 【求助】二元函数拟合
已经有12人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
lidaxue
木虫 (正式写手)
之乎者也
- 应助: 10 (幼儿园)
- 金币: 1952.1
- 散金: 612
- 红花: 2
- 帖子: 527
- 在线: 856.1小时
- 虫号: 259962
- 注册: 2006-06-17
- 性别: GG
- 专业: 几何学
8楼2010-10-28 13:12:05
quan_zhang
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 365.4
- 散金: 12
- 红花: 1
- 帖子: 134
- 在线: 41.6小时
- 虫号: 864139
- 注册: 2009-10-06
- 专业: 金属结构材料
3楼2010-10-26 15:39:08
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
quan_zhang(金币+1): 2010-10-27 09:15:46
quan_zhang(金币+1): 2010-10-27 09:16:08
quan_zhang(金币+1): 2010-10-27 09:16:08
|
X = LSQNONLIN(FUN,X0,LB,UB) defines a set of lower and upper bounds on the design variables, X, so that the solution is in the range LB <= X <= UB. Use empty matrices for LB and UB if no bounds exist. Set LB(i) = -Inf if X(i) is unbounded below; set UB(i) = Inf if X(i) is unbounded above. |
4楼2010-10-26 17:40:52
lidaxue
木虫 (正式写手)
之乎者也
- 应助: 10 (幼儿园)
- 金币: 1952.1
- 散金: 612
- 红花: 2
- 帖子: 527
- 在线: 856.1小时
- 虫号: 259962
- 注册: 2006-06-17
- 性别: GG
- 专业: 几何学
quan_zhang(金币+2): 2010-10-27 09:16:03
quan_zhang(金币+2): 2010-10-27 21:15:11
quan_zhang(金币+2): 2010-10-27 21:15:11
|
function seqcurvefit11 clear all; load seqdata; beta0=[0.005 0.001]; lb=[0 0];ub=[inf inf]; [beta,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@seqfun,beta0,lb,ub,[],t,c); function y = seqfun(beta,t,c) ca=exp(-beta(1)*t); cb=beta(1)/(beta(2)-beta(1))*(exp(-beta(1)*t)-exp(-beta(2)*t)); y=[ca-c(:,1) cb-c(:,2)]; 楼主 依葫芦画瓢吧~ |
5楼2010-10-26 18:48:56
28