24小时热门版块排行榜

>论坛更新日志 (3431)
>虫友互识 (1653)
>休闲灌水 (745)
>导师招生 (86)
>文献求助 (62)
>基金申请 (33)
>论文投稿 (33)
>考博 (30)
>论文道贺祈福 (29)
>微米和纳米 (27)
>硕博家园 (27)
>考研 (27)
>博后之家 (26)
>教师之家 (26)
>公派出国 (19)
>找工作 (16)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

料峭春寒5

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 761.8
帖子: 63
在线: 23.2小时
虫号: 1079812
注册: 2010-08-22
专业: 纤维

[求助] 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？

请大神们指点，感谢感谢！

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2017-08-28 13:09:43

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

pdl9527

专家顾问 (小有名气)

专家经验: +8
应助: 100 (初中生)
金币: 2109.3
红花: 23
帖子: 282
在线: 112小时
虫号: 1227333
注册: 2011-03-09
性别: GG
专业: 分离过程
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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*-')