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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 最小二乘法解方程组系数
111/2返回列表12下一页
查看: 2227  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

aimforao

禁虫 (正式写手)
本帖内容被屏蔽

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-10-19 19:41:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mmjjwanm

金虫 (正式写手)
★
jjdg(金币+1): 感谢参与 2011-10-20 01:29:24
什么和什么的相关系数？？
R^2么？

拟合优度？衡量方程的拟合优劣的指标？
赞一下(1人) 回复此楼
2楼2011-10-19 21:48:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

aimforao

禁虫 (正式写手)
本帖内容被屏蔽
3楼2011-10-19 22:00:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mmjjwanm

金虫 (正式写手)
★ ★ ★
jjdg(金币+2): 感谢参与 2011-10-20 01:29:35
xzhdty(金币+1): 欢迎讨论交流 2011-10-20 07:19:07
引用回帖:
3楼: Originally posted by aimforao at 2011-10-19 22:00:47:
就是把数据带进去之后求出的ie和给出的ie的误差。你说的都对，

纠正你一个概念上的错误：
1）拟合优度，R^2，是衡量模型拟合的契合度，取值在0到1之间。值越大说明用此模型拟合数据越适合。
2）“数据代进去以及求出的ie和给出的ie之间的误差”，这个应该是残差平方和，此值越小说明拟合越好。

用SAS的NLIN过程（非线性模型拟合）采用牛顿迭代的方式进行参数估计，发现参数b和e的值在迭代过程中保持不变（初始值设为0）。数据、程序和结果见附件

不应助，建议LZ仔细检查下数据。。。
赞一下(2人) 回复此楼
4楼2011-10-19 23:06:06
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

aimforao

禁虫 (正式写手)
本帖内容被屏蔽
5楼2011-10-20 10:02:13
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

【答案】应助回帖

aimforao(金币+100): 40 2011-10-21 08:56:00
引用回帖:
1楼: Originally posted by aimforao at 2011-10-19 19:41:02:
以下是数据和方程，希望可以把相关系数R也求出来，或是用别的方法求出系数来。系数最好不能有零。

ie=[0.365 0.479 0.592 0.427 0.575 0.711 0.595 0.701 0.891];
x1=[-7.022 -7.022 -7.022 -7.096 -7.096 -7 ...

方程可以转化成一下的线性形式：




计算步骤：
ie=[0.365 0.479 0.592 0.427 0.575 0.711 0.595 0.701 0.891]';
x1=[-7.022 -7.022 -7.022 -7.096 -7.096 -7.096 -6.648 -6.648 -6.648]';
x2=[-2.391 -2.391 -2.391 -2.524 -2.524 -2.524 -2.392 -2.392 -2.392]';
x3=[0.7963 0.7963 0.7963 1.9082 1.9082 1.9082 1.9277 1.9277 1.9277]';
x4=[106.444 106.444 106.444 98.851 98.851 98.851 135.257 135.257 135.257]';
x5=[10e-4 10e-3 10e-2 10e-4 10e-3 10e-2 10e-4 10e-3 10e-2]';

AA1 = x1.*x5;
AA2 = x2.*x5;
AA3 = x3.*x5;
AA4 = x4.*x5;
AA5 = x5;
AA = [AA1,AA2,AA3,AA4,AA5];
BB = 1./(1./ie-1);
x = AA\BB

计算结果：
Warning: Rank deficient, rank = 3,  tol =   3.9865e-014.
> In ddddd at 15

x =

   20.1351
         0
   19.7714
    1.3244
         0


>> [S,V,D]=svd(AA)

S =

   -0.0053   -0.0032    0.0078   -0.9992   -0.0358   -0.0030    0.0033    0.0078    0.0071
   -0.0534   -0.0319    0.0777   -0.0265    0.7190    0.0307    0.0563    0.3172   -0.6064
   -0.5338   -0.3186    0.7769    0.0126   -0.0715   -0.0030   -0.0057   -0.0318    0.0606
   -0.0050   -0.0062   -0.0060    0.0048   -0.1155   -0.9639   -0.0131   -0.0759   -0.2269
   -0.0496   -0.0623   -0.0596   -0.0035   -0.2080    0.2283   -0.1128   -0.7088   -0.6161
   -0.4960   -0.6232   -0.5963    0.0003    0.0220   -0.0132    0.0114    0.0716    0.0639
   -0.0068    0.0071   -0.0018    0.0061   -0.1055    0.0196    0.9900   -0.0633   -0.0656
   -0.0678    0.0707   -0.0176    0.0240   -0.6358    0.1306   -0.0609    0.6143   -0.4327
   -0.6776    0.7072   -0.1755   -0.0025    0.0646   -0.0133   -0.0038   -0.0608    0.0439


V =

   19.9909         0         0         0         0
         0    0.2101         0         0         0
         0         0    0.0871         0         0
         0         0         0    0.0000         0
         0         0         0         0    0.0000
         0         0         0         0         0
         0         0         0         0         0
         0         0         0         0         0
         0         0         0         0         0


D =

    0.0595    0.9413   -0.0662    0.2992   -0.1282
    0.0210    0.3092    0.0781   -0.6744    0.6656
   -0.0135   -0.0383   -0.9947   -0.0690    0.0651
   -0.9979    0.0641    0.0111   -0.0012   -0.0009
   -0.0086   -0.1128    0.0059    0.6715    0.7323

>> format long
>> AA

AA =

  -0.007022000000000  -0.002391000000000   0.000796300000000   0.106444000000000   0.001000000000000
  -0.070220000000000  -0.023910000000000   0.007963000000000   1.064440000000000   0.010000000000000
  -0.702200000000000  -0.239100000000000   0.079630000000000  10.644400000000001   0.100000000000000
  -0.007096000000000  -0.002524000000000   0.001908200000000   0.098851000000000   0.001000000000000
  -0.070960000000000  -0.025240000000000   0.019082000000000   0.988510000000000   0.010000000000000
  -0.709600000000000  -0.252400000000000   0.190820000000000   9.885100000000001   0.100000000000000
  -0.006648000000000  -0.002392000000000   0.001927700000000   0.135257000000000   0.001000000000000
  -0.066480000000000  -0.023920000000000   0.019277000000000   1.352570000000000   0.010000000000000
  -0.664800000000000  -0.239200000000000   0.192770000000000  13.525700000000001   0.100000000000000


可见系数矩阵的秩为3.
赞一下 回复此楼
6楼2011-10-20 21:06:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
★
jjdg:编辑内容 2011-10-21 00:08
jjdg(金币+1): 感谢参与 2011-10-21 00:08:22
原因是你的以下数据是每三组是重复的

x1=[-7.022 -7.022 -7.022 -7.096 -7.096 -7.096 -6.648 -6.648 -6.648];
x2=[-2.391 -2.391 -2.391 -2.524 -2.524 -2.524 -2.392 -2.392 -2.392];
x3=[0.7963 0.7963 0.7963 1.9082 1.9082 1.9082 1.9277 1.9277 1.9277];
x4=[106.444 106.444 106.444 98.851 98.851 98.851 135.257 135.257 135.257];

[ Last edited by jjdg on 2011-10-21 at 00:08 ]
赞一下(1人) 回复此楼
7楼2011-10-20 21:09:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mmjjwanm

金虫 (正式写手)
引用回帖:
6楼: Originally posted by change0618 at 2011-10-20 21:06:36:
方程可以转化成一下的线性形式：




计算步骤：
ie=[0.365 0.479 0.592 0.427 0.575 0.711 0.595 0.701 0.891]';
x1=[-7.022 -7.022 -7.022 -7.096 -7.096 ...

数学建模高人！
hiahia,您老人家滴头像很香艳
开个玩笑，呵呵，以后多向你学习
赞一下 回复此楼
8楼2011-10-20 21:50:12
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

aimforao

禁虫 (正式写手)
本帖内容被屏蔽
9楼2011-10-21 08:58:27
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
你还没有明白意思啊？你看似是9组数据，其实你只有3组，5个未知数3组数据，必然出现这种情况
赞一下 回复此楼
10楼2011-10-21 09:11:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aimforao 的主题更新
111/2返回列表12下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +11 月@163.com 2026-04-07 13/650 2026-04-11 22:55 by BruceLiu320
[考研] 电子信息270求调剂 +14 terminal469 2026-04-07 14/700 2026-04-11 19:44 by laoshidan
[考研] 269求调剂 +11 啊啊我我 2026-04-07 11/550 2026-04-11 16:45 by vgtyfty
[考研] 求调剂，一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 307求调剂 +8 tzq94092 2026-04-10 8/400 2026-04-10 17:33 by 286640313
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 生物与医药273求调剂 +18 荔题南墙 2026-04-05 19/950 2026-04-10 08:14 by kangsm
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +16 路痴小琪 2026-04-05 16/800 2026-04-10 08:08 by kangsm
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 一志愿中国科学院上海有机所，有机化学356分找调剂 +11 Nadiums 2026-04-09 11/550 2026-04-09 18:04 by lijunpoly
[考研] 322求调剂，08工科 +3 今天是个小号 2026-04-08 3/150 2026-04-09 15:53 by wp06
[考研] 生物学308分求调剂（一志愿华东师大） +13 相信必会光芒万� 2026-04-06 16/800 2026-04-09 13:54 by 徐良白眉大侠
[考研] 招收有机化学、化工，药学，食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 085404，334分，求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 275 求调剂 +8 Lei812514 2026-04-07 8/400 2026-04-08 12:46 by chemisry
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 277求调剂 +5 考研调剂lxh 2026-04-05 5/250 2026-04-05 19:03 by chy09050039
信息提示
关闭
请填处理意见
关闭