24小时热门版块排行榜    

查看: 923  |  回复: 1

4758188748

新虫 (初入文坛)

[求助] 求助。matlab共享参数类的自定义公式拟合修改问题 已有1人参与

是共享参数类的自定义公式模拟,目的是求5个系数,y1-x1对应q0固定系数,y2-x2对应q1固定系数,依次类推四个系数。我这个代码是肯定不行呀- -,想要知道这个能不能求帮忙修改,并且把这几个曲线绘制到一张图上,谢谢谢谢谢谢

clear
clc
data = [36.287         34.765         38.875         37.683         33.580         27.010         34.720         35.470
42.726         37.060         38.475         42.038         31.388         30.550         33.480         37.460
44.688         38.646         40.391         43.510         30.339         24.934         25.711         29.321
43.517         39.432         41.318         42.936         29.632         26.228         27.102         28.644
44.222         38.646         42.530         41.718         28.211         24.991         27.792         29.638
43.771         40.398         42.454         44.425         27.862         25.291         26.591         27.814]
q0=0;
q1=0.2;
q2=0.35;
q3=0.5;
x1 = data(:,1);
x2 = data(:,2);
x3 = data(:,3);
x4 = data(:,4);
y1 = data(:,5);
y2 = data(:,6);
y3 = data(:,7);
y4 = data(:,8);
f = @(a,b,c,d,e)[a.*x1./(b*q0^3+c*q0^2+d*q0+e);
a.*x2./(b*q1^3+c*q1^2+d*q1+e);
a.*x3./(b*q2^3+c*q2^2+d*q2+e);
a.*x4./(b*q3^3+c*q3^2+d*q3+e)]
[X,RESNORM,RESIDUAL,EXITFLAG,OUTPUT,LAMBDA] = lsqcurvefit(f, [1 1 1 1 1 ], [x1;x2;x3;x4],[y1;y2;y3;y4],[0 0 0 0 0 0],[],[])

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

hollya

至尊木虫 (文坛精英)

葛城?{美

【答案】应助回帖

楼主,你这个函数好像就是个线性方程,f(x)=kx, 这样拟合结果可能不太好啊。而且就算有5个参数需要似合,5个也不是相互独立的。
你看看下面的这个:其中 k1=b/a,k2=c/a,k3=d/a,k4=e/a

data = [36.287         34.765         38.875         37.683         33.580         27.010         34.720         35.470
        42.726         37.060         38.475         42.038         31.388         30.550         33.480         37.460
        44.688         38.646         40.391         43.510         30.339         24.934         25.711         29.321
        43.517         39.432         41.318         42.936         29.632         26.228         27.102         28.644
        44.222         38.646         42.530         41.718         28.211         24.991         27.792         29.638
        43.771         40.398         42.454         44.425         27.862         25.291         26.591         27.814]
q0=0;
q1=0.2;
q2=0.35;
q3=0.5;
q=[q0,q1,q2,q3]
Q=[q.^3;q.^2;q;ones(1,4)]
xdata=data(:,1:4)
ydata=data(:,5:8)
f = @(k,x) x./repmat([k(1),k(2),k(3),k(4)]*Q,6,1)
k0=[1 1 1 1]
k = lsqcurvefit(f,k0,xdata,ydata)
f_fit=f(k,xdata)
plot(xdata(:,1),ydata(:,1),'o',xdata(:,1),f_fit(:,1),...
    xdata(:,2),ydata(:,2),'*',xdata(:,2),f_fit(:,2),...
    xdata(:,3),ydata(:,3),'+',xdata(:,3),f_fit(:,3),...
    xdata(:,4),ydata(:,4),'x',xdata(:,4),f_fit(:,4))
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
2楼2019-02-02 02:45:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 4758188748 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂0703 +5 啊我我的 2026-03-11 5/250 2026-03-14 19:34 by 哦哦123
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 求调剂 +4 拾柒12。 2026-03-08 4/200 2026-03-14 04:14 by JourneyLucky
[考研] 求调剂 +5 鹤遨予卿 2026-03-09 5/250 2026-03-14 02:44 by JourneyLucky
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 0817化学工程319求调剂 +8 lv945 2026-03-08 10/500 2026-03-12 05:26 by wll0811
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 293求调剂,一志愿陕师大生物学 +3 ??????.?.??? 2026-03-09 3/150 2026-03-11 10:02 by 学员8dgXkO
信息提示
请填处理意见