24小时热门版块排行榜    

查看: 2546  |  回复: 2

dplrwp

新虫 (初入文坛)

[求助] 怎样用非线性最小二乘法拟合得到方程中的未知参数 已有2人参与

方程y=Ax+Bx/(1+Cx) 采用最小二乘法模拟得到参数A、B、C。
x:0.452 、3.099、15.492、33.8464、51.522、70.314
y:9.548、16.901、21.508、26.1536、28.478、29.686
回复此楼

» 猜你喜欢

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

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

yangwww123

新虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
将公式改写成 (A+B)x + AC x^2 -C xy = y;
型函数 P = [x,  x^2,  -xy];   v = [A+B; AC; C]; 剩下的公式不好敲,看下面的程序。
采用matlab编程如下:
%% main
clear all; close all; clc;
xi = [0.452, 3.099, 15.492, 33.8464, 51.522, 70.314];
yi = [9.548, 16.901, 21.508, 26.1536, 28.478, 29.686];

K = zeros(3);
b = zeros(3,1);
for i=1:size(xi,2)
    x = xi(1,i); y = yi(1,i);
    P = [x, x^2, -x*y];
    K = K + P'*P;
    b = b + P'*y;
end
% v = [A+B; AC; C]
v = K\b;

C = v(3,1);
A = v(2,1)/C;
B = v(1,1)-A;

y_app = A*xi + B*xi./(1+C*xi);

plot(xi,yi,'-r*',xi,y_app,'g-o')
上面一句是画图,这里面插不了图,你把程序跑一边就能看见。
计算结果:给出x值,得到的模拟的y值;
2.4513   11.2275   22.2999   26.3347   28.2622   29.7421

前两个值不太准确,建议多采用几个数据点,结果会跟精确。
科研吊丝变大牛
2楼2015-06-18 07:45:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
a        0.148938319420554
b        36.4207912045872
c        1.79290500985833
3楼2015-06-18 11:11:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dplrwp 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +18 医学老男孩 2026-08-13 41/2050 2026-08-14 12:34 by lqwchd
[基金申请] filecode +5 cratir 2026-08-14 9/450 2026-08-14 12:25 by 医学老男孩
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +5 启萌科技 2026-08-12 13/650 2026-08-14 10:53 by lnjx2024
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +10 Tide man 2026-08-10 11/550 2026-08-14 10:30 by 沁言学术
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见