24小时热门版块排行榜    

查看: 706  |  回复: 0

aaa2742

木虫 (正式写手)

[求助] 非线性拟合表达式较复杂,拆分后可以处理吗,在Matlab中

有一个光谱的表达式较复杂,且某些变量与波长也有变化关系,想做一个非线性拟合,
一开始只是用的循环枚举参数的方法求得最优的参数,现在使用想使用lqscurvefit,但是报错
CODE:
pai = 3.14;
len = length(xdata);
N0 = 1*ones(len,1);
%N0 = 1;
N1 = 1.448+3634./power(xdata,2)-8.885./power(xdata,4);
N2 = ((3.363*xdata -946.7)./(xdata -328.6))+1i*((-0.03382*xdata + 27.03)./(xdata -365.3));
t = @(x) 2*pai.*N1*x./xdata;
r = @(x)((N0-N2).*cos(t(x))+1i*(N0.*N2./N1-N1).*sin(t(x)))./((N0+N2).*cos(t(x))+1i*(N0.*N2./N1+N1).*sin(t(x)));
% test = r(100);
Fun = @(x)(r(x).*conj(r(x)));
RFun = @(x) norm(Ry-r(x).*conj(r(x)));
% [res,fval,exitflag,output] = fminbnd(RFun,1,5000)

res = lsqcurvefit(Fun,1,xdata,Ry)
m = 1;
for h = 1:0.01:100
    res(m) = RFun(h);
    m = m +1;
end
[C,I] = min(res)
k = 1:0.01:100;
k(I)

使用lqscurvefit出错:
CODE:
Error using @(x)(r(x).*conj(r(x)))
Too many input arguments.

Error in lsqcurvefit (line 195)
            initVals.F = feval(funfcn_x_xdata{3},xCurrent,XDATA,varargin{:});

Error in RFitModel_2 (line 34)
res = lsqcurvefit(Fun,1,xdata,Ry)

Caused by:
    Failure in initial user-supplied objective function evaluation. LSQCURVEFIT cannot continue.

使用循环枚举的方法倒是可以得到一个较优的数值,但是计算耗时长,且间隔设置的大小对结果都有影响,间隔小,耗时长;间隔大,不准确。
请问该怎样解决该问题,想寻找别的算法,公式中含有复数,只不过最后乘以它的共轭,所以如果都有参数的话计算的数值是整数,但是现在是反推其中的一个参数。
所以尝试使用lqscurvefit进行拟合,结果错了,难道非得把表达式全部写全弄到一个里面吗,那也太长也容易出错啊。
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aaa2742 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 材料与化工求调剂 +3 为学666 2026-03-16 3/150 2026-03-16 15:09 by 加号+
[考研] 312求调剂 +3 陌宸希 2026-03-16 4/200 2026-03-16 15:06 by peike
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
请填处理意见