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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求高手帮忙为小弟做一个四元回归曲线，高分悬赏！
41/1返回列表
查看: 1088  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者kk1424将赠送您 50 个金币

kk1424

主管区长 (著名写手)

中药研究者

优秀区长优秀区长优秀版主优秀版主

[求助] 求高手帮忙为小弟做一个四元回归曲线，高分悬赏！ 已有1人参与

遇到棘手的问题，我有MATLAB软件，但不是很会。求高手达人，帮忙给个编程。。。。感激不尽.
这回归方程：
y1=0.227220-0.002859x1+0.018825x2-0.003369x3-0.021391x4

拜托了！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
药理毒理
1楼 2014-02-25 16:50:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

感谢参与，应助指数 +1
figure(1);
x=start:step:end;
y =0.227220-0.002859*x+0.018825*x^2-0.003369*x^3-0.021391*x^4;
plot(x, y,"可设曲线属性";

关键是
x=start:step:end   始点，重点均未知不知道
赞一下 回复此楼
2楼2014-02-25 17:30:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
2楼: Originally posted by feixiaolin at 2014-02-25 17:30:50
figure(1);
x=start:step:end;
y =0.227220-0.002859*x+0.018825*x^2-0.003369*x^3-0.021391*x^4;
plot(x, y,"可设曲线属性";

关键是
x=start:step:end   始点，重点均未知不知道

看错了，是回归。
matlab回归  
https://wenku.baidu.com/view/87fa18e0524de518964b7d27.html
赞一下 回复此楼
3楼2014-02-25 17:33:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
Matlab多元线性回归

[ b , bint , r , rint , stats ]=regress ( y , x ) ，

其中b 是回归方程中的参数估计值，bint 是b 的置信区间，r 和rint 分别表示残差及残差对应的置信区间。StatS 数组包含三个数字，分别是相关系数，F 统计量及对应的概率p 值。
拟合结果：
Y=b(1)x(1)+b(2)x(2)+b(3)x(3)+…+b(n)x(n)
b(1)是系数,x(1)为全1的一个列向量。
注意：不是插值。


x=[1097 1284 1502 1394 1303 1555 1917 2051 2111 2286 2311 2003 2435 2625 2948 3155 3372];%因变量时间序列数据
y=[698 872 988 807 738 1025 1316 1539 1561 1765 1762 1960 1902 2013 2446 2736 2825];%自变量时间序列数据
X=[ones(size(x')),x'];
[b,bint,r,rint,stats]=regress(y',X,0.05);%调用一元回归分析函数
rcoplot(r,rint)%画出在置信度区间下误差分布.

举例：
x =

     1     2     4     9
     1     4     3     7
     1     5     9     0
     1     9     1     8
>> y=[10 3 90 48]';
>> [ b , bint , r , rint , stats ]=regress ( y , x )

得到的结果
b =

-186.8333
   16.0238
   21.8571
    8.5952


bint =

   NaN   NaN
   NaN   NaN
   NaN   NaN
   NaN   NaN


r =

  1.0e-013 *

   -0.5684
   -0.9948
   -0.1421
    0.1421


rint =

   NaN   NaN
   NaN   NaN
   NaN   NaN
   NaN   NaN


stats =

     1   NaN   NaN   NaN
另外，用b=inv(x)*y
得到的结果和上面用regression得到的一样。可见，求逆的问题也是用了最小二乘的原理。


但是，regress更优于inv，体现在当x的元素存在缺陷时。举例：
x =

     1     0     0     4
     1     0     0     3
     1     0     0     5
     1     0     0     8
y =

     3
     4
     3
     9
[ b , bint , r , rint , stats ]=regress ( y , x )
Warning: X is rank deficient to within machine precision.
b =
   -0.9643
         0
         0
    1.1429
bint =
  -12.0010   10.0724
         0         0
         0         0
   -0.9245    3.2102
r =
   -0.6071
    1.5357
   -1.7500
    0.8214


rint =
   -9.2316    8.0173
   -1.9141    4.9855
   -7.4997    3.9997
    0.2465    1.3964


stats =
    0.7388    5.6575    0.1405    3.2321
赞一下(1人) 回复此楼
4楼2014-02-26 08:43:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kk1424 的主题更新
41/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0831生医工第一轮调剂失败求助 +14 小熊睿睿_s 2026-04-11 18/900 2026-04-16 20:41 by lpl364211
[考研] 297，工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +7 否极泰来2026 2026-04-15 9/450 2026-04-16 13:40 by liumingli7817
[考研] 273求调剂 +6 白居不易. 2026-04-09 8/400 2026-04-15 22:02 by wooluyong
[考研] 一志愿A区211，22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 085404 22408 309分求调剂 +9 lzmk 2026-04-14 10/500 2026-04-15 20:02 by 学员JpLReM
[考研] 310求调剂 +16 666真好 2026-04-11 18/900 2026-04-15 13:28 by 黑科技矿业
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 材料考研调剂 +29 云木达达 2026-04-11 31/1550 2026-04-13 13:32 by lyh鲁老师
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 0858求调剂 5+5 Gky09300550， 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
信息提示
关闭
请填处理意见
关闭