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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 最小二乘法解方程组系数
51/1返回列表
查看: 2160  |  回复: 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 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 接收调剂 +5 津萌津萌 2026-03-02 12/600 2026-03-02 17:04 by 津萌津萌
[考研] 0856化工专硕求调剂 +15 董boxing 2026-03-01 15/750 2026-03-02 15:06 by 晃晃不许晃
[考研] 化工专硕348，一志愿985求调剂 +6 弗格个 2026-02-28 9/450 2026-03-02 14:09 by liyongv
[考博] 诚招农业博士 +3 心欣向荣 2026-02-28 3/150 2026-03-02 13:33 by 时间不狗
[考研] 338求调剂 +3 18162027187 2026-03-02 3/150 2026-03-02 13:12 by houyaoxu
[考研] 291 求调剂 +3 化工2026届毕业� 2026-03-02 3/150 2026-03-02 12:55 by houyaoxu
[考研] 275求调剂 +4 明远求学 2026-03-01 4/200 2026-03-02 12:45 by 无际的草原
[基金申请] 成果系统访问量大，请15分钟后再尝试。由此给您造成的不便，敬请谅解。 +5 xhuama 2026-03-02 5/250 2026-03-02 12:34 by stidwellNK
[考研] 求调剂 +3 熬夜的猫头鹰 2026-03-02 3/150 2026-03-02 11:45 by 刘兵
[考研] 0856材料与化工，270求调剂 +8 YXCT 2026-03-01 9/450 2026-03-02 11:01 by 无际的草原
[考研] 264求调剂 +4 巴拉巴拉根556 2026-02-28 4/200 2026-03-02 10:48 by yuchj
[考研] 材料类求调剂 +11 wana_kiko 2026-02-28 14/700 2026-03-02 08:46 by 聪明的大松鼠
[考研] 材料复试调剂 +4 学材料的点 2026-03-01 5/250 2026-03-02 08:26 by houyaoxu
[考研] 0857调剂 +4 一ll半 2026-02-28 5/250 2026-03-02 02:33 by 908055542
[考研] 279求调剂 +3 dua1 2026-03-01 4/200 2026-03-02 00:23 by 大脸蛋子
[考研] 272求调剂 +6 田智友 2026-02-28 6/300 2026-03-01 21:40 by 公瑾逍遥
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考博] 26申博 +4 想申博！ 2026-02-26 6/300 2026-03-01 17:32 by 想申博！
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
信息提示
关闭
请填处理意见
关闭