| 查看: 3372 | 回复: 11 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
matlab多元线性回归拟合求助 已有2人参与
|
||||
大家好，最近在尝试利用matlab进行多元线性回归拟合，由于是零基础感觉困难重重，请各位高手在百忙之中给予指点 ，该如何编程呢 T(温度)=[353.15 353.15 353.15 353.15 363.75 363.75 363.75 363.75 373.75 373.75 373.75 373.75 393.75 393.75 393.75 393.75 413.75 413.75 413.75 413.75]; d（密度）=[593.89 636.74 671.27 699.95 533.17 580.38 619.21 651.66 480.53 528.94 570.19 605.34 401.15 445.86 486.70 523.39 347.58 386.61 423.72 458.44]; S=[4.39 4.65 5.24 5.64 4.40 4.67 5.68 5.72 4.85 5.59 5.94 6.92 5.35 6.53 7.23 7.95 6.76 7.28 8.21 9.96]; 需要拟合的方程是：ln S=k ln d+A/T+B 如何通过拟合求出k，A，B？ 请大神们指点，感谢感谢！ |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于均匀设计实验数据的回归模型
已经有7人回复
- 求助用matlab 多元非线性拟合
已经有4人回复
- 菜鸟求助matlab习题两道，跪求指导，不胜感激！！！！
已经有4人回复
4楼2017-08-28 17:35:10
tieer
木虫 (正式写手)
- 程序强帖: 4
- 应助: 48 (小学生)
- 金币: 1528.2
- 散金: 565
- 红花: 26
- 帖子: 551
- 在线: 739.7小时
- 虫号: 761880
- 注册: 2009-05-02
- 专业: 高分子材料结构与性能
2楼2017-08-28 15:22:24
3楼2017-08-28 15:33:44
pdl9527
专家顾问 (小有名气)
-
专家经验: +8 - 应助: 100 (初中生)
- 金币: 2110.3
- 红花: 23
- 帖子: 282
- 在线: 112小时
- 虫号: 1227333
- 注册: 2011-03-09
- 性别: GG
- 专业: 分离过程
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
料峭春寒5: 金币+20, ★★★★★最佳答案, 我先研究下，不懂的地方再多多向您请教哈，谢谢！ 2017-08-28 19:36:47
感谢参与，应助指数 +1
料峭春寒5: 金币+20, ★★★★★最佳答案, 我先研究下，不懂的地方再多多向您请教哈，谢谢！ 2017-08-28 19:36:47
|
[k,A,B]= [1.46 -2914.37 0.36]; 下面是程序，你可以根据经验更改初值x0，也许能获得更好的结果。 function question_46 clear;clc;close all T=[353.15 353.15 353.15 353.15 363.75 363.75 363.75 363.75 373.75 373.75 373.75 373.75 393.75 393.75 393.75 393.75 413.75 413.75 413.75 413.75]; d=[593.89 636.74 671.27 699.95 533.17 580.38 619.21 651.66 480.53 528.94 570.19 605.34 401.15 445.86 486.70 523.39 347.58 386.61 423.72 458.44]; S=[4.39 4.65 5.24 5.64 4.40 4.67 5.68 5.72 4.85 5.59 5.94 6.92 5.35 6.53 7.23 7.95 6.76 7.28 8.21 9.96]; fun=@(x) exp(x(1).*log(d)+x(2)./T+x(3))-S; x0=[1,1,1]; [x,resnorm,residual,exitflag] = lsqnonlin(fun,x0) plot3(T,d,S,'ro-') hold on plot3(T,d,fun(x)+S,'b*-') |
5楼2017-08-28 19:06:18
200