24小时热门版块排行榜

返回列表

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

料峭春寒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？

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

回复此楼

» 猜你喜欢

广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有9人回复
闲聊已经有3人回复
江苏省自然基金什么时候出结果已经有4人回复
如何查询2026年获批项目名称？已经有3人回复
两块石头已经有10人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有19人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有15人回复
面上没中，邀请各位路过的虫友分析一下分数已经有13人回复
初秋的晨风已经有6人回复
Ei源刊怎么投已经有4人回复

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

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

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

pdl9527

专家顾问 (小有名气)

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

引用回帖:

9楼: Originally posted by pdl9527 at 2017-08-29 19:58:54
对于这个问题来说使用非线性最小二乘法来计算没有任何问题，你这里用的左除（mldivide），如果你打开帮助文件仔细阅读你会发现这样一行字，If A is a rectangular m-by-n matrix with m ~= n, and B is a matrix w ...

function question_47
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];
[x1,resnorm,residual,exitflag] = lsqnonlin(fun,x0);
x2=[log(d) 1./T ones(numel(T),1)]\log(S);
norm1=sum(fun(x1).^2)
norm2=sum(fun(x2).^2)