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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 最小二乘法解方程组系数
51/1返回列表
查看: 2161  |  回复: 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 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290求调剂 +6 ErMiao1020 2026-03-02 6/300 2026-03-02 18:14 by lature00
[考研] 化学，材料，环境类求调剂 +5 考研版棒棒 2026-03-02 5/250 2026-03-02 18:07 by lature00
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 14/700 2026-03-02 17:44 by 好好好1233
[考研] 接收调剂 +5 津萌津萌 2026-03-02 12/600 2026-03-02 17:04 by 津萌津萌
[考研] 一志愿华南理工大学材料与化工326分，求调剂 +3 wujinrui1 2026-02-28 3/150 2026-03-02 16:36 by chuocheng
[考研] 材料284求调剂，一志愿郑州大学英一数二专硕 +12 想上岸的土拨鼠 2026-02-28 12/600 2026-03-02 16:18 by youmomaoyan
[考研] 求调剂 +9 yunziaaaaa 2026-03-01 10/500 2026-03-02 15:18 by 半截的诗0927
[考研] 0856化工专硕求调剂 +15 董boxing 2026-03-01 15/750 2026-03-02 15:06 by 晃晃不许晃
[考研] 26考研报考西工大材料308分求调剂 +4 weizhong123 2026-03-01 4/200 2026-03-02 12:46 by 无际的草原
[考研] 276求调剂 +4 路lyh123 2026-02-28 5/250 2026-03-02 11:20 by yuchj
[考研] 0856材料与化工，270求调剂 +8 YXCT 2026-03-01 9/450 2026-03-02 11:01 by 无际的草原
[考研] 274求调剂 +3 cgyzqwn 2026-03-01 7/350 2026-03-02 10:38 by lature00
[考研] 275求调剂 +3 L-xin? 2026-03-01 6/300 2026-03-02 10:22 by 热情沙漠
[基金申请] 本子写完了，给DS兄弟看了，得了92分 +3 Doma 2026-03-01 7/350 2026-03-02 00:00 by jnzsy
[基金申请] 成果系统访问量大，请一小时后再尝试。---NSFC啥时候好哦，已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考研] 306分材料调剂 +4 chuanzhu川烛 2026-03-01 5/250 2026-03-01 19:48 by 无际的草原
[考博] 26申博 +4 想申博！ 2026-02-26 6/300 2026-03-01 17:32 by 想申博！
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 313求调剂 +3 水流年lc 2026-02-28 3/150 2026-03-01 16:01 by 新能源达人
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
信息提示
关闭
请填处理意见
关闭