24小时热门版块排行榜    

查看: 2287  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

liuyu0710

金虫 (小有名气)

[交流] 【求助】回归拟合多元二次方程 已有4人参与

我的实验是球面对称设计后,对数据进行回归拟合,数据如下:

        x1        x2        x3        y1        y2
1        18.45        0.17        1.85        109.6        18
2        41        0.17        1.85        128.4        38.04
3        18.45        0.48        1.85        98.7        43.05
4        41.55        0.48        1.85        124.6        55.71
5        18.45        0.17        4.15        129.6        22.61
6        41.55        0.17        4.15        128.8        49
7        18.45        0.48        4.15        103.2        60.19
8        41.55        0.48        4.15        141.7        78.71
9        10        0.33        3        105.8        57.66
10        50        0.33        3        196.1        64.53
11        30        0.05        3        102.3        30.5
12        30        0.6        3        119.5        50.23
13        30        0.33        1        118        65.02
14        30        0.33        5        121.3        82.77
15        30        0.33        3        130        87


目标是拟合成三元二次方程,具体形式如下:

y1=b0+b1x1+b2x2+b3x3+b11x1^2+b12x1x2+b13x1x3+b22x2^2+b23x2x3+b33x3^2
y2=.............
通过拟合得到b0,b1,b2,b3,b11等常数的值

要求y1>100 y2>80%,解不等式求出x1,x2,x3, 并找出最适合的解。

手里面有SPSS,lingo和statistica,哪种软件更合适些?具体怎么操作?
最好能给个模板,谢谢各位了,

我的邮箱liuyu0710@163.com
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mbzr

木虫 (正式写手)

liuyu0710(金币+2): 2010-08-05 19:39:19
用数值分析里面矩阵公式加excell矩阵运算就可以了吧
5楼2010-07-31 08:36:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

dzl723

金虫 (正式写手)

liuyu0710(金币+5): 2010-08-05 19:39:40
用SAS吧,模板暂时找不到了,我再找找,找到了给你发过去
2楼2010-07-31 01:03:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mrzouhao

木虫之王 (文坛精英)

liuyu0710(金币+3): 2010-08-05 19:39:31
statistica处理实验是球面对称设计更合适,软件在小木虫上搜,不难用~~
3楼2010-07-31 07:46:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qiruiyun

木虫 (正式写手)

liuyu0710(金币+5): 2010-08-05 19:39:08
我用mathematica写了个代码
cc = Import["d:\\xmc.txt", "Table"];
y1 = Table[Delete[cc[], 5], {i, Length[cc]}];
y2 = Table[Delete[cc[], 4], {i, Length[cc]}];
fitY1 = Fit[
  y1, {1, x1, x2, x3, x1^2, x1 x2, x1 x3, x2^2, x2 x3, x3^2}, {x1, x2,
    x3}]
fitY2 = Fit[
  y2, {1, x1, x2, x3, x1^2, x1 x2, x1 x3, x2^2, x2 x3, x3^2}, {x1, x2,
    x3}]
结果为:
y1=79.7916 - 1.98694 x1 + 0.0439335 x1^2 + 90.7723 x2 + 3.26941 x1 x2 -
295.394 x2^2 + 25.4082 x3 - 0.0779321 x1 x3 + 1.30797 x2 x3 -
3.43553 x3^2
y2=-143.898 + 5.43815 x1 - 0.0813438 x1^2 + 502.35 x2 - 1.10043 x1 x2 -
696.748 x2^2 + 26.1101 x3 + 0.118026 x1 x3 + 16.5905 x2 x3 -
4.94622 x3^2
你参考一下吧。
4楼2010-07-31 08:24:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 085700资源与环境308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[考研] 274求调剂 +5 S.H1 2026-03-18 5/250 2026-03-18 21:27 by guosr9609
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 0703化学调剂,求各位老师收留 +9 秋有木北 2026-03-14 9/450 2026-03-18 19:49 by macy2011
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
信息提示
请填处理意见