应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 【求助】求助:最小二乘法曲线拟合参数
111/2返回列表12下一页
查看: 1247  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

天边5328

新虫 (小有名气)

[交流] 【求助】求助:最小二乘法曲线拟合参数 已有4人参与

请高手帮忙编个程序，求出参数k，非常感谢！！！

方程：y=y0+(ylimt-y0)/(2*c0)*(x+1/k+c0-((x+c0+1/k)^2-4*c0*x)^0.5)

其中y0=0.723;   ylimt=0.330;   c0=0.00002


xdata=[0.00005 0.00015 0.0002 0.00025 0.0003 0.00035 0.00045 0.00055 0.0006 0.00075 0.00085 0.00095 0.001 0.00115 0.00125 0.00145 0.00155 0.00175 0.00205 0.00245]

ydata=[0.706 0.665 0.652 0.624 0.599 0.574 0.545 0.526 0.492 0.472 0.450 0.434 0.426 0.402 0.387 0.369 0.364 0.345 0.338 0.330]

[ Last edited by 天边5328 on 2010-4-3 at 17:35 ]
回复此楼

» 猜你喜欢
1楼 2010-04-03 17:33:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxlion811

金虫 (正式写手)
★
wangen994(金币+1):活动期间额外奖励 2010-04-09 18:51
请问楼主需要什么程序的软件？
能做拟合的软件太多了。
赞一下(1人) 回复此楼
开心努力一辈子
2楼2010-04-09 18:38:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天边5328

新虫 (小有名气)
jjdg:请记得升级你的论坛级别 2010-04-10 00:45
想用matlab软件
回复此楼
3楼2010-04-09 20:35:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

he866837

至尊木虫 (著名写手)
楼主学会做了吗
回复此楼
4楼2010-04-09 20:45:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天边5328

新虫 (小有名气)
我以前从没有接触过这个软件，不知道编程。我编的程序，软件提示错误。
赞一下 回复此楼
5楼2010-04-10 08:22:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sgyin

铁杆木虫 (正式写手)
★ ★ ★ ★
jjdg(金币+2):辛苦了 2010-04-11 03:21
wangen994(金币+2):活动期间额外奖励 2010-04-11 20:19
ff=inline('0.723+(0.330-0.723)./(2.*0.00002).*(x+1./k(1)+0.00002-((x+0.00002+1./k(1)).^2-4.*0.00002.*x).^0.5)','k','x')

[xx,res]=lsqcurvefit(ff,[1e3],xdata,ydata);


这是固定y0=0.723;   ylimt=0.330;   c0=0.00002 的结果！拟合只能这样了，楼主把这些参数放开试试

[ Last edited by sgyin on 2010-4-10 at 23:49 ]
赞一下(2人) 回复此楼
<ahref="http://www.TickerFactory.com/"><imgborder="0"src=&am
6楼2010-04-10 23:45:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sgyin

铁杆木虫 (正式写手)
ff=inline('k(2)+(k(3)-k(2))./(2.*0.00002).*(x+1./k(1)+0.00002-((x+0.00002+1./k(1)).^2-4.*0.00002.*x).^0.5)','k','x')
[xx,res]=lsqcurvefit(ff,[1e3,1,1],xdata,ydata);
fplot(@(x)ff(xx,x),[0,0.003],'r'),hold on,plot(xdata,ydata,'o')

这是拟合结果k= [1295.69223047089         0.756467370036168         0.171649191261954]
赞一下 回复此楼
<ahref="http://www.TickerFactory.com/"><imgborder="0"src=&am
7楼2010-04-10 23:54:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sgyin

铁杆木虫 (正式写手)
★ ★ ★ ★ ★
天边5328(金币+50):非常感谢您 2010-04-11 07:42
wangen994(金币+5):活动期间额外奖励 2010-04-11 10:52
ff=inline('k(2)+(k(3)-k(2))./(2.*k(4)).*(x+1./k(1)+k(4)-((x+k(4)+1./k(1)).^2-4.*k(4).*x).^0.5)','k','x')
[xx,res]=lsqcurvefit(ff,[1e3,1,1,0.0001],xdata,ydata);
fplot(@(x)ff(xx,x),[0,0.003],'r'),hold on,plot(xdata,ydata,'o')

k= [  3599.30404451002          0.73350054683102         0.263275701253758      0.000684697129278193]
赞一下(1人) 回复此楼
<ahref="http://www.TickerFactory.com/"><imgborder="0"src=&am
8楼2010-04-10 23:59:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天边5328

新虫 (小有名气)
我想请教您，为什么我输入最后那一行后，软件提示错误呢？
赞一下 回复此楼
9楼2010-04-11 07:58:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sgyin

铁杆木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wangen994(金币+1):活动期间额外奖励 2010-04-11 20:20
先要给出xdata, ydata的值
引用回帖:
Originally posted by 天边5328 at 2010-04-11 07:58:56:
我想请教您，为什么我输入最后那一行后，软件提示错误呢？

赞一下(2人) 回复此楼
<ahref="http://www.TickerFactory.com/"><imgborder="0"src=&am
10楼2010-04-11 12:13:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 天边5328 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 重磅！青年科学基金项目（C类）资助增幅预计超过50% +6 水和泥不是水泥 2026-05-13 8/400 2026-05-16 19:50 by Equinoxhua
[有机交流] 求助2，4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[基金申请] 精华III:  评审感受-评审感受-评审感受 +15 ferrarichen 2026-05-11 19/950 2026-05-16 19:44 by Equinoxhua
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-16 18:45 by j6b2pdz07o
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 3/150 2026-05-16 17:55 by j6b2pdz07o
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 3/150 2026-05-16 17:50 by j6b2pdz07o
[基金申请] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 v9tggjlwd0 2026-05-15 3/150 2026-05-16 17:37 by j6b2pdz07o
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 v9tggjlwd0 2026-05-15 3/150 2026-05-16 17:30 by j6b2pdz07o
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 16:35 by x28q7dxf75
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 6/300 2026-05-16 16:17 by 0i5p09z61n
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 xx7gd5zq4e 2026-05-15 5/250 2026-05-16 12:24 by h3oerqvkv9
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 11:36 by h3oerqvkv9
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 cjf4bx70cj 2026-05-14 6/300 2026-05-16 11:16 by h3oerqvkv9
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 04:29 by k37jurhrau
[教师之家] 上海大学实验技术岗位非升即走 +5 嘻嘻哈哈乐呵呵 2026-05-15 5/250 2026-05-16 00:17 by caiyun
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[高分子] 本人最近太闲了，谁有问题可以提，每天会统一回复 +8 一切都是空工 2026-05-12 19/950 2026-05-14 20:03 by 一切都是空工
[基金申请] 请问大佬b0816评完了吗 +3 市民华南虎 2026-05-12 7/350 2026-05-14 07:41 by 市民华南虎
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
关闭
请填处理意见
关闭