| 查看: 1032 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
matlab M文件 给规定参数设个边界
|
|||
|
Sample Text t=[0.0833 0.5 1 2 4 8]; c=[2.2505 1.6489 0.2789 0.2412 0.1077 0.1035]; 我用已知的两个函数模型拟合时出现了这种情况 （一）t=[0.0833 0.5 1 2 4 8]; c=[2.2505 1.6489 0.2789 0.2412 0.1077 0.1035]; myfun=inline('A(1)*exp(-A(2)*x)','A','x'); [A iter sse]=nlinfit(t,c,myfun,[1 1]) 结果：A = 2.6400 1.4099 iter = 6 sse = 0.2891 （二）：t=[0.0833 0.5 1 2 4 8]; c=[2.2505 1.6489 0.2789 0.2412 0.1077 0.1035]; myfun=inline('A(1)*exp(-A(2)*x)+A(3)*exp(-A(4)*x)','A','x'); [A iter sse]=nlinfit(t,c,myfun,[2.4 1.5 0.022 0.2]) 结果：A = 2.6266 1.4510 0.0227 -0.2002 iter = 19 sse = 0.2734 但是我的这个函数有个要求，就是参数A的值要大于零，这样的话A(4)就不符合了 （三）还有一种情况就是，用二的方式，但我给另外的随意的初值， 结果：A = 0.8571 1.4100 1.7829 1.4098 A(1)+A(3)=2.64，A(2)基本上等于A(4) iter = 6 sse = 0.2891 而这样的话就相当于把（一）的公式拆分成了两个一样模型的函数，与（一）没什么区别，我在用（二）计算式，我希望能够拟合出的结果中A(2)与A(4)不同，但是有还得是正值，所以我想给A(2)和A(4)设边界大于零 该怎么弄，我是matlab的菜鸟用户，请高手帮帮忙，谢谢！ [ Last edited by cheng1378653 on 2013-4-1 at 12:52 ] |
回复此楼» 猜你喜欢
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有281人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有28人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有22人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
月只蓝
禁虫
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68128.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
7楼2013-04-05 14:05:26
2楼2013-04-01 16:19:45
lgycjpcqu
主管区长
- 应助: 35 (小学生)
- 金币: 566.8
- 散金: 411
- 帖子: 371
- 在线: 176.9小时
- 虫号: 1327328
- 注册: 2011-06-20
- 性别: GG
- 专业: 电化学
★
fegg7502: 金币+1, 鼓励交流 2013-04-02 08:58:45
fegg7502: 金币+1, 鼓励交流 2013-04-02 08:58:45
|
function E=fconstrain(a,t,c) C=a(1).*exp(-a(2).*t)+a(3).*exp(-a(4).*t); E=c-C; 命令窗口 clear a0=[2.4 1.5 0.022 0]; LB=[-inf 0 -inf 0]; UB=[]; t=[0.0833 0.5 1 2 4 8]; c=[2.2505 1.6489 0.2789 0.2412 0.1077 0.1035]; options=optimset('lsqnonlin'); [A norm res ef]=lsqnonlin(@fconstrain,a0,LB,UB,options,t,c); cfit=A(1).*exp(-A(2).*t)+A(3).*exp(-A(4).*t); plot(t,c,'*',t,cfit,'rh') legend('原始数据','拟合数据') 结果是 A = 2.587176647067854 1.500899543303364 0.069042182557558 0.000000000000041  约束.jpg |
3楼2013-04-01 19:30:57
4楼2013-04-02 20:58:53
10