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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 最小二乘法解方程组系数
51/1返回列表
查看: 2241  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

aimforao

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-10-19 19:41:02
已阅   关注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的回帖
查看全部 11 个回答

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的回帖
查看全部 11 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[材料工程] 材料调剂推荐 +9 蛋糕x2 2026-04-07 9/450 2026-04-13 09:07 by lhj2009
[考研] 2本，初试303，0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[硕博家园] 新一代电子信息294求调剂 不挑学校 +7 Ytyt11 2026-04-09 8/400 2026-04-12 16:57 by ajpv风雷
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 347求调剂 +4 mhyqyy 2026-04-06 4/200 2026-04-12 02:27 by 秋豆菜芽
[考研] 267求调剂 +8 再忙也要吃饭啊 2026-04-09 8/400 2026-04-11 21:42 by cfdbai
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[考研] 085402通信工程调剂，有4项学科竞赛国奖（电赛国二），硕士研究生调剂自荐信。 +5 m永o不v言o弃m 2026-04-09 5/250 2026-04-11 09:33 by zhq0425
[考研] 337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 08:57 by zhq0425
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 282，电气工程专业，求调剂，不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[考研] 一志愿中科大070300化学，314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 264求调剂 +11 麦小叮当 2026-04-07 11/550 2026-04-08 16:05 by 一只好果子?
[考研] 307求调剂 +14 超级伊昂大王 2026-04-06 14/700 2026-04-08 07:03 by 无际的草原
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 085602调剂 初试总分335 +3 19123253302 2026-04-06 3/150 2026-04-07 18:00 by jp9609
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 材料调剂 +17 小刘同学吖吖 2026-04-06 18/900 2026-04-07 11:41 by 诗与自由
[考研] 307求调剂 +3 所念及所望 2026-04-06 3/150 2026-04-06 17:30 by 土木硕士招生
信息提示
关闭
请填处理意见
关闭