24小时热门版块排行榜    

查看: 1364  |  回复: 7

小鬼木木

新虫 (小有名气)

[求助] 求最小二乘法MATLAB程序 已有3人参与

关系式为:X=(A-V)/(B*V-C)
          V:137   104     69     55     39      32     22     13     7
          X:0        6         15    21      30      36     48     63    78
求用最小二乘法计算三个参数A B C 的MATLAB程序
回复此楼

» 猜你喜欢

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

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

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
小鬼木木: 金币+3, 有帮助 2016-03-17 15:23:02
改写成  v =C *x -B *(v x) +A

然后 v  对 x 和 v*x 作多元线性回归,用regress 即可。

如果对结果不满意,可以以此为初值,用 X-(A-V)/(B*V-C) 残差向量(的模)为目标函数,

再进一步作非线性最小二乘拟合, 用 lsqnonlin  即可, 算法可以选择 Levenberg-Maquardt 或gauss newton
ส็็็็็็็็็็็็็็็็็็็็
2楼2016-03-16 15:53:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

science-lp

新虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
小鬼木木: 金币+2, ★★★很有帮助 2016-03-17 15:38:43
%%%%%%%%%%%%%%%%%%%%%%计算
v=[137   104     69     55     39      32     22     13     7];
x=[0     6    15    21      30      36     48     63    78];
k=(-x*diag(v));
d=diag(ones(9));
l=[x',k',d];
h=inv(l'*l);
j=l'*v';
sishu=h*j
c= -1.3671
b=    0.0488
a=  139.2208
%%%%%%%%%%%%%%%%%%%%%%%%%%%%系数检验,绘图对比
clear;
v=[137   104     69     55     39      32     22     13     7];
x=[0     6    15    21      30      36     48     63    78];
plot(v,x,'r-o','linewidth',2)
grid on
hold on
v=[137   104     69     55     39      32     22     13     7];
for i=1:9
k(i)=(139.32208-v(i))/(0.0488*v(i)+1.3671)
end
plot(v,k,'-gx','linewidth',2)
3楼2016-03-16 22:13:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
1stOpt求解:
CODE:
Function  X=(A-V)/(B*V-C);
Data;
V=137   104     69     55     39      32     22     13     7;
X=0        6         15    21      30      36     48     63    78;

均方差(RMSE):0.333010610104675
残差平方和(SSE):0.998064597980593
相关系数(R): 0.999908497408983
相关系数之平方(R^2): 0.99981700319069
确定系数(DC): 0.999817003190689
卡方系数(Chi-Square): 0.572161462394941
F统计(F-Statistic): 16393.7284551577

参数        最佳估算
----------        -------------
a        141.61456338051
b        0.0515957786789501
c        -1.36677726438402
4楼2016-03-16 22:45:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菜鸟也爱球

铁虫 (著名写手)

哇卡卡卡呀哈哈
5楼2016-03-16 23:19:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小鬼木木

新虫 (小有名气)

引用回帖:
2楼: Originally posted by cooooldog at 2016-03-16 15:53:37
改写成  v =C *x -B *(v x) +A

然后 v  对 x 和 v*x 作多元线性回归,用regress 即可。

如果对结果不满意,可以以此为初值,用 X-(A-V)/(B*V-C) 残差向量(的模)为目标函数,

再进一步作非线性最小二乘拟 ...

谢谢您的回答,受益了!
6楼2016-03-17 15:23:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小鬼木木

新虫 (小有名气)

引用回帖:
3楼: Originally posted by science-lp at 2016-03-16 22:13:41
%%%%%%%%%%%%%%%%%%%%%%计算
v=;
x=;
k=(-x*diag(v));
d=diag(ones(9));
l=;
h=inv(l'*l);
j=l'*v';
sishu=h*j
c= -1.3671
b=    0.0488
a=  139.2208
%%%%%%%%%%%%%%%%%%%%%%%%%%%%系数检验,绘图对 ...

谢谢您的帮助,非常感谢!
7楼2016-03-17 15:24:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小鬼木木

新虫 (小有名气)

引用回帖:
4楼: Originally posted by dingd at 2016-03-16 22:45:46
1stOpt求解:

Function  X=(A-V)/(B*V-C);
Data;
V=137   104     69     55     39      32     22     13     7;
X=0        6         15    21      30      36     48     63    78;


均方差(RMSE): ...

非常感谢!
8楼2016-03-17 15:25:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小鬼木木 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279分求调剂 一志愿211 +11 chaojifeixia 2026-03-19 12/600 2026-03-21 01:49 by 星空星月
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 324求调剂 +5 lucky呀呀呀鸭 2026-03-20 5/250 2026-03-20 22:30 by 促天成
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
请填处理意见