【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 1066  |  回复: 6

daijiantao

新虫 (正式写手)

[求助] 请教大家,这个公式中怎么拟合出n的值呀?求用MATLAB帮编一下,感激不尽!!已有1人参与

请教大家,这个公式中怎么拟合出n的值呀?已知D(t),D0,和t,未知是k和n。其中D(t)和t是变化的,D0就相当于一个常数,一直是281.27.
数据如下:


t               d(t)

0                281.27      d0
100                309.41      d1
1000        544.85         d2
43285        801.98         d3



哪位大神能帮忙用matlab帮编写一下,或者其他程序能算也行,感激不尽!!

请教大家,这个公式中怎么拟合出n的值呀?求用MATLAB帮编一下,感激不尽!!
微信截图_20181217212136.png
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2018-12-18 15:18:59
daijiantao: 金币+100, ★★★★★最佳答案 2018-12-18 20:56:22
xdata=[0,100,1000,43285]
ydata=[281.27,309.41,544.85,801.98]
y=@(n,x)(n(1).*x+281.27^n(2)).^(1/n(2)) %n(1)=k1, n(2)=n
n0=[1,1]
for ii=1:10
n=lsqcurvefit(y,n0,xdata,ydata)
n0=n
end
plot(xdata,ydata,xdata,y(n,xdata),'o')
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
2楼2018-12-18 00:15:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daijiantao

新虫 (正式写手)

引用回帖:
2楼: Originally posted by hollya at 2018-12-18 00:15:44
xdata=
ydata=
y=@(n,x)(n(1).*x+281.27^n(2)).^(1/n(2)) %n(1)=k1, n(2)=n
n0=
for ii=1:10
n=lsqcurvefit(y,n0,xdata,ydata)
n0=n
end
plot(xdata,ydata,xdata,y(n,xdata),'o')...

谢谢您,您这个程序看起来没问题呀,只是我跑一遍出来的结果为啥不是一个完整抛物线平滑的曲线而是两段分段直线,如下图。
还有第三个点的位置显示有点问题,请教一下
请教大家,这个公式中怎么拟合出n的值呀?求用MATLAB帮编一下,感激不尽!!-1
微信图片_20181218132459.png

3楼2018-12-18 13:26:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

啊,那个点数据fit出来的,线是原数据,有点误导啊,写成这样你再看下。

发自小木虫IOS客户端
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
4楼2018-12-18 19:45:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

最后一行:plot(xdata,ydata,’o’,xdata,y(n,xdata))

发自小木虫IOS客户端
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
5楼2018-12-18 19:46:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daijiantao

新虫 (正式写手)

引用回帖:
5楼: Originally posted by hollya at 2018-12-18 19:46:48
最后一行:plot(xdata,ydata,’o’,xdata,y(n,xdata))

哦哦,我明白了,那个线是原始数据点是拟合的线。。缺失被误导了,谢谢您的指点
6楼2018-12-18 20:56:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

引用回帖:
6楼: Originally posted by daijiantao at 2018-12-18 20:56:09
哦哦,我明白了,那个线是原始数据点是拟合的线。。缺失被误导了,谢谢您的指点...

没错没错

发自小木虫IOS客户端
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
7楼2018-12-18 21:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 daijiantao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 24年 申博 化学/材料 一作6篇sci +3 wangyp123 2024-04-23 4/200 2024-04-23 23:38 by 莹亮的星空
[博后之家] 南京大学-广州大学联合招聘博士后 欢迎广大优秀人才!!! +4 黑魔变身啾 2024-04-20 8/400 2024-04-23 23:25 by 黑魔变身啾
[考博] 真的好想读博! +7 wangzhe_bs 2024-04-22 9/450 2024-04-23 22:20 by wangzhe_bs
[找工作] 工作问题求助 +3 向前?@ 2024-04-18 5/250 2024-04-23 21:25 by 向前?@
[论文投稿] 无under review,直接DIP,咋么办。 +10 lizhengke06 2024-04-19 16/800 2024-04-23 19:15 by 化学程序员
[有机交流] 紧急求助,有谁用过三甲基硅醇钾TMSOK吗? 1000+3 genius2008 2024-04-19 8/400 2024-04-23 16:57 by czyzsu
[考博] 重庆大学计算机学院2024博士招生 +4 y_y_y_yyy 2024-04-18 6/300 2024-04-23 16:31 by 让一切随风l
[基金申请] 前同事拿我以前未中的国自然标书申请今年的国自然,现在基金委发来调查函 +28 bjdxyxy 2024-04-19 44/2200 2024-04-23 16:29 by 淀粉搬运工
[考博] 山东大学化学与化工学院刘磊课题组博士研究生招生启事 +7 sdorganic 2024-04-17 9/450 2024-04-23 16:28 by 让一切随风l
[基金申请] 颜宁:基础研究应顶天立地,立地即从需求端作为问题的出发点,立项背景就是面向真实需求 +10 zju2000 2024-04-18 10/500 2024-04-23 11:46 by 稻花香
[论文投稿] 编辑返稿让改格式,这个时候能大修内容吗? +6 双倍好运锦鲤 2024-04-17 9/450 2024-04-23 09:49 by fslaser
[考博] 研二光催化6月底4篇2区 +7 wjtab 2024-04-22 11/550 2024-04-23 06:59 by byron2012
[教师之家] 河海大学常州校区新进的讲师待遇如何?考核要求如何? +3 德国小球迷 2024-04-19 18/900 2024-04-22 19:23 by ZHONGWU_U
[考博] 青岛科技大学顾玉超教授课题组招收博士生 +3 yanqingwang 2024-04-19 3/150 2024-04-22 15:21 by 理科生yes
[访问学者] CSC的访问学者申请,没有个评审意见,也不知道怎么改,还有必要申请吗 +3 flyingship 2024-04-20 3/150 2024-04-22 13:52 by louistom_521
[考研] 问题已经解决。 +13 lekinna 2024-04-19 30/1500 2024-04-22 13:09 by 矛dei到
[考研] 博士生招生!! +3 2458827884 2024-04-19 3/150 2024-04-22 12:37 by 钰波love
[考研] 考研英语带学 +6 choice666 2024-04-18 6/300 2024-04-19 17:46 by charles-c
[有机交流] 紫外光谱 50+3 54胡 2024-04-17 3/150 2024-04-19 10:48 by Nanamiwww
[论文投稿] 投稿求助 5+3 我是洲洲啊 2024-04-17 5/250 2024-04-18 17:13 by topedit
信息提示
请填处理意见