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

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 使用matlab实现批量线性回归并出图（需要得到每条拟合直线斜率、截距以及R）
【调剂】北京石油化工学院2024年16个专业接受调剂
141/2返回列表12下一页
查看: 5288  |  回复: 13
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

hxyonline

木虫 (小有名气)

[求助] 使用matlab实现批量线性回归并出图（需要得到每条拟合直线斜率、截距以及R）

使用matlab实现批量线性回归并出图（需要得到每条拟合直线斜率、截距以及R）

我只会使用origin对多组数据进行线性回归，不会用matlab。现在数据比较多，大约100多组，我想请大家帮忙，能不能帮我编段代码，用matlab实现一下。

以下是19组x数据
第一列是x1，第二列是x2，第三列是x3，第四列是x4。
0.002022556        0.001988032        0.001969762        0.001964628
0.001996335        0.001964189        0.001943093        0.001935219
0.001981611        0.001949346        0.001927129        0.001917188
0.001970649        0.00193815        0.001915423        0.001904113
0.001961535        0.001928879        0.001905915        0.001893661
0.001953477        0.001920729        0.00189764        0.001884743
0.001946023        0.001913234        0.001890102        0.001876756
0.001938926        0.001906103        0.001882979        0.001869328
0.001931971        0.001899125        0.001876041        0.001862156
0.001924985        0.001892187        0.001869134        0.001855023
0.001917828        0.001885019        0.001861987        0.001847673
0.001910269        0.00187738        0.001854392        0.001839875
0.001902071        0.001869042        0.001846104        0.001831356
0.001892955        0.001859669        0.001836815        0.001821706
0.001882585        0.001848774        0.00182603        0.001810434
0.001870366        0.001835692        0.001813        0.001796783
0.001855466        0.001819289        0.001796539        0.001779488
0.001836159        0.001797598        0.001774305        0.00175647
0.001807239        0.00176581        0.001741439        0.00172295

以下是19组y数据
第一列是y1，第二列是y2，第三列是y3，第四列是y4。
-3.67690886        -2.82263223        -2.54038545        -2.3374172
-2.88747038        -2.2224454        -1.92551946        -1.76636473
-2.53216404        -1.89237371        -1.56284418        -1.40323708
-2.31213051        -1.67715533        -1.3273451        -1.14904952
-2.18366325        -1.52771779        -1.17044401        -0.98115865
-2.09629845        -1.41768361        -1.06197465        -0.8556566
-2.0261498        -1.37633075        -1.00242748        -0.77572852
-2.00359223        -1.34177795        -0.96967448        -0.72675601
-2.00739724        -1.34464356        -0.96082092        -0.72327488
-2.03091777        -1.36084567        -0.97708542        -0.73964904
-2.07702435        -1.41028593        -1.03703842        -0.78959161
-2.14686422        -1.49418987        -1.11865095        -0.87128426
-2.25176376        -1.60933943        -1.22500753        -0.99028035
-2.37892058        -1.76497984        -1.37653368        -1.14793724
-2.53168303        -1.93794198        -1.5470687        -1.32228376
-2.72859257        -2.13983856        -1.77128415        -1.55105839
-2.96375081        -2.42463548        -2.05080296        -1.82164519
-3.28871856        -2.76514356        -2.41175178        -2.16648122
-3.83274982        -3.23537955        -2.93338179        -2.65588527

(x1,y1), (x2,y2),(x3,y3),(x4,y4)分别是坐标上2个点，这4个点可以拟合出一条直线，我需要得到拟合出直线的斜率，截距和相关系数。
第一行做完之后再对第二行进行直线拟合。。。。。然后直到把19行数据都就完。

看书看晕了，也没有找到类似例子，希望大家能给我一个基本代码，或者告诉我计算实现方法，我自己再找找书看。谢谢大家了。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-07-09 15:43:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
hxyonline: 金币+20, ★★★★★最佳答案, 非常感谢 2013-07-10 12:38:52
datax=[];  %所有的X写进去，按照x1-x19
datay=[];  %所有的Y写进去，按照y1-y19
for i=1:19
p=polyfit(x(i),y(i),1)
z=polyval(p,x);
figure
plot(x,y,'ro',x,z,'b');
end
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2013-07-09 17:18:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
2楼: Originally posted by 月只蓝 at 2013-07-09 17:18:27
datax=[];  %所有的X写进去，按照x1-x19
datay=[];  %所有的Y写进去，按照y1-y19
for i=1:19
p=polyfit(x(i),y(i),1)
z=polyval(p,x);
figure
plot(x,y,'ro',x,z,'b');
end

改正一下
datax=[];  %所有的X写进去，按照x1-x19
datay=[];  %所有的Y写进去，按照y1-y19
for i=1:19
p=polyfit(datax(i),datay(i),1)
z=polyval(p,datax(i));
figure
plot(x,y,'ro',x,z,'b');
end
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2013-07-09 17:19:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
3楼: Originally posted by 月只蓝 at 2013-07-09 17:19:40
改正一下
datax=[];  %所有的X写进去，按照x1-x19
datay=[];  %所有的Y写进去，按照y1-y19
for i=1:19
p=polyfit(datax(i),datay(i),1)
z=polyval(p,datax(i));
figure
plot(x,y,'ro',x,z,'b');
end...

不好意思， "plot(x,y,'ro',x,z,'b');"改成“plot(datax(i),datay(i),'ro',datax(i),z,'b')” ,或者去掉也行
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

hxyonline
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2013-07-09 17:22:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hxyonline

木虫 (小有名气)
送红花一朵
引用回帖:
4楼: Originally posted by 月只蓝 at 2013-07-09 17:22:11
不好意思， "plot(x,y,'ro',x,z,'b');"改成“plot(datax(i),datay(i),'ro',datax(i),z,'b')” ,或者去掉也行...

谢谢您，不过我想知道数据x和y写入的格式，是每个数一个逗号，没4个一个分号吗？
赞一下 回复此楼
5楼2013-07-09 17:31:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
5楼: Originally posted by hxyonline at 2013-07-09 17:31:38
谢谢您，不过我想知道数据x和y写入的格式，是每个数一个逗号，没4个一个分号吗？...

比如导入x数据的时候，可以在excel中复制19列x，然后直接粘贴进中括号中即可，不需要逗号、分号
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
6楼2013-07-09 18:08:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hxyonline

木虫 (小有名气)
引用回帖:
6楼: Originally posted by 月只蓝 at 2013-07-09 18:08:06
比如导入x数据的时候，可以在excel中复制19列x，然后直接粘贴进中括号中即可，不需要逗号、分号...

谢谢您的程序，我昨天晚上把程序又重新调试了一下。代码如下：
format long
[filename, pathname]= uigetfile('*.xls');
file=[pathname filename];
x=xlsread(file);
for i=1:19
xd1=x(i,1);xd2=x(i,2);xd3=x(i,3);xd4=x(i,4);yd1=x(i,6);yd2=x(i,7);yd3=x(i,8);yd4=x(i,9);
m=[xd1 xd2 xd3 xd4];  
n=[yd1 yd2 yd3 yd4];  
p=polyfit(m,n,1);
z=polyval(p,m);
poly2sym(p,'x');
hold on;
plot(m,n,'o',m,z)
end

可是我现在的问题是我想得到每一条直线的斜率和截距，可是输出的结果P只有一行，不知道我应该怎样能输出所有的结果。

还有一个问题，请问线性相关度怎么导出来？我非常需要这个数据。

麻烦您了
使用matlab实现批量线性回归并出图（需要得到每条拟合直线斜率、截距以及R）
QQ图片20130710082147.jpg
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com
  • 附件 1 : shuju.xls
    • 2013-07-10 08:22:37, 25.5 K
7楼2013-07-10 08:25:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hxyonline

木虫 (小有名气)
引用回帖:
6楼: Originally posted by 月只蓝 at 2013-07-09 18:08:06
比如导入x数据的时候，可以在excel中复制19列x，然后直接粘贴进中括号中即可，不需要逗号、分号...

抱歉，忘了说了，附件的excel表格里面是我的原始数据。  图片中上面是结果P，只有一行。
赞一下 回复此楼
8楼2013-07-10 08:27:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by hxyonline at 2013-07-10 08:25:57
谢谢您的程序，我昨天晚上把程序又重新调试了一下。代码如下：
format long
= uigetfile('*.xls');
file=;
x=xlsread(file);
for i=1:19
xd1=x(i,1);xd2=x(i,2);xd3=x(i,3);xd4=x(i,4);yd1=x(i,6);yd2=x(i ...

p=polyfit(m,n,1); 这行程序中的分号去掉即可。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
9楼2013-07-10 11:12:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by hxyonline at 2013-07-10 08:25:57
谢谢您的程序，我昨天晚上把程序又重新调试了一下。代码如下：
format long
= uigetfile('*.xls');
file=;
x=xlsread(file);
for i=1:19
xd1=x(i,1);xd2=x(i,2);xd3=x(i,3);xd4=x(i,4);yd1=x(i,6);yd2=x(i ...

相关性系数的平方R2=1-sum((yexp-y).^2)./sum((yexp-mean(y)).^2); 其中y是计算值，yexp是你的实验值，代入上述即可。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
10楼2013-07-10 11:19:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hxyonline 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 某系主任被同校教授举报包养情妇、长期嫖娼 +11 zju2000 2024-04-13 13/650 2024-04-18 22:20 by jjsfb
[考博] 申请24博士 材料/化工/环境 +3 满目_星辰 2024-04-17 3/150 2024-04-18 17:35 by 我属驴核动力驴
[基金申请] 下雨了 +13 zju2000 2024-04-16 18/900 2024-04-18 17:02 by feng6531
[基金申请] 基金和生小孩 +32 Ausy 2024-04-15 34/1700 2024-04-18 12:13 by wangzhenyft
[公派出国] CSC在职申请求问 +5 Amber0919 2024-04-16 5/250 2024-04-17 16:15 by liaojiayan
[考研] 322求调剂 +7 本己上岸 2024-04-16 7/350 2024-04-17 11:49 by duanxz
[考研] 广州大学光电信息工程专业调剂，招收物理学专业学生 +5 txhx4010 2024-04-14 7/350 2024-04-16 10:52 by domax
[考研] 296求调剂 +3 Cclocomotive 2024-04-16 4/200 2024-04-16 10:04 by 19862091
[考研] 334求调剂 +4 学药救人 2024-04-14 4/200 2024-04-15 15:05 by hunanzang
[材料综合] 请问哪里可以测试低压吸氢的PCT曲线 +3 yunshengcd 2024-04-14 5/250 2024-04-15 13:59 by zqdsb
[考研] 284生物与医药求调剂（考数学没过A区国家线） +5 胖男孩 2024-04-12 5/250 2024-04-13 18:38 by yezixiaop
[考研] 298求调剂 +6 新阶段+有 2024-04-12 8/400 2024-04-13 17:42 by 新阶段+有
[考研] 331求调剂 +12 廖喆虓 2024-04-12 12/600 2024-04-13 17:10 by lincunhui
[考研] 086000生医没考数学304分，一志愿西农求调剂 +3 王有田有学上 2024-04-13 5/250 2024-04-13 16:18 by coco1981
[考研] 314求调剂 +18 阿洛法 2024-04-12 20/1000 2024-04-13 13:40 by kxxmaqq
[考研] 280求调剂 +7 黑皮冰棒 2024-04-12 7/350 2024-04-13 09:47 by lincunhui
[硕博家园] 生物炭采购 +3 锦鲤附体@ 2024-04-12 3/150 2024-04-12 22:10 by Danny614
[考研] 266求调剂 +3 ～酥糖 2024-04-12 15/750 2024-04-12 21:22 by 赵燕高仙兰
[考研] 343求调剂 +3 好好爱蓝宝 2024-04-12 3/150 2024-04-12 16:44 by 阿杰鲁呀
[考研] 生态学293求调剂，孩子想要读书！ +3 啊Q不乖 2024-04-12 9/450 2024-04-12 11:11 by 啊Q不乖
信息提示
关闭
请填处理意见
关闭