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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 求最小二乘法MATLAB程序
81/1返回列表
查看: 1362  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小鬼木木

新虫 (小有名气)

[求助] 求最小二乘法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程序
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-03-16 15:25:32
已阅   回复此楼   关注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
赞一下(5人) 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
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)
赞一下(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
赞一下(1人) 回复此楼
4楼2016-03-16 22:45:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菜鸟也爱球

铁虫 (著名写手)
高手真多

发自小木虫Android客户端
回复此楼
哇卡卡卡呀哈哈
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的回帖
相关版块跳转 我要订阅楼主 小鬼木木 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 本人考085602 化学工程 专硕 +19 不知道叫什么！ 2026-03-15 21/1050 2026-03-20 20:48 by zhukairuo
[考研] 求调剂 +5 Mqqqqqq 2026-03-19 5/250 2026-03-20 20:46 by zhukairuo
[考研] 085410人工智能专硕317求调剂（0854都可以） +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 324分 085600材料化工求调剂 +3 llllkkkhh 2026-03-18 3/150 2026-03-19 14:22 by houyaoxu
[考研] 332求调剂 +3 ydfyh 2026-03-17 3/150 2026-03-19 10:14 by 功夫疯狂
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
关闭
请填处理意见
关闭