24小时热门版块排行榜    

Znn3bq.jpeg
查看: 6042  |  回复: 0

yi小小鸟

木虫 (小有名气)

小小鸟

[求助] 求指导,MATLAB非线性拟合无法完成,提示-Complex value computed by model function

本人想通过matlab做一个非线性数据拟合,一直实现不了,求大神指导,万分感谢!

拟合公式:R(fx)=3*A*us/(ua+us)/(((((3*ua*(ua+us))+(2*3.14*fx)^2)^(1/2))/(ua+us)+1)*((((3*ua*(ua+us))+(2*3.14*fx)^2)^(1/2))/(ua+us)+3*A))
此公式为文献上给出的原始公式,A为一个常数,R反射率,fx空间频率,均为实验数据,ua和us是最终想拟合的参数。

Matlab程序:

clear all;clc;
for i = 1:5
    I0 = imread([路径',num2str(i),'.jpg']);
    I{i} = I0;
end
n=1.51;fx_0=[0.04:0.01:0.08];
Reff=0.0636*n+0.668+0.710/n-1.440/n^2;
A_value=(1-Reff)/[2*(1+Reff)];
syms A ua us fx;
f = sym('3*A*us/(ua+us)/(((((3*ua*(ua+us))+(2*3.14*fx)^2)^(1/2))/(ua+us)+1)*((((3*ua*(ua+us))+(2*3.14*fx)^2)^(1/2))/(ua+us)+3*A))');
f = subs(f,A,A_value);
f_str = char(f);
for i = 1:201
    for j = 1:201
        st_ = [-0.0146 -0.8187];
        ft_ = fittype(f_str,...
    'dependent',{'R'},'independent',{'fx'},...
    'coefficients',{'ua', 'us'});
for k = 1:5
    R(k) = double(I{k}(i,j));
end
[cfun,gof] = fit(fx_0',R',ft_,'Startpoint',st_);
I_result_ua(i,j) = cfun.a;
I_result_us(i,j) = cfun.b;
    end
end

subs(f,A,A)
f = sym('(3*A*us/(ua+us)/(((((3*ua*(ua+us))+(2*3.14*fx)^2)^(1/2))/(ua+us)+1)*((((3*ua*(ua+us))+(2*3.14*fx)^2)^(1/2))/(ua+us)+3*A)))')

拟合后matlab一直显示Complex value computed by model function, fitting cannot continue.
Try using or tightening upper and lower bounds on coefficients.
更改了ua和us的初始变量,也没有效果。

注:程序中的数据为五幅图片,分别是fx在0.04、0.05、0.06、0.07、0.08/mm下的,拟合思路为用5幅图的同一位置(x,y)的数据进行拟合,逐点拟合。所以最后输出的是ua和us201*201的一组数据,可以通过imshow显示图片。

一直实现不了,求大神指导,万分感谢!

求指导,MATLAB非线性拟合无法完成,提示-Complex value computed by model function
1.jpg


求指导,MATLAB非线性拟合无法完成,提示-Complex value computed by model function-1
2.jpg


求指导,MATLAB非线性拟合无法完成,提示-Complex value computed by model function-2
3.jpg


求指导,MATLAB非线性拟合无法完成,提示-Complex value computed by model function-3
4.jpg


求指导,MATLAB非线性拟合无法完成,提示-Complex value computed by model function-4
5.jpg
回复此楼

» 猜你喜欢

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 yi小小鸟 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 26药学专硕105500求调剂 +7 喽哈加油 2026-04-13 8/400 2026-04-19 20:21 by Equinoxhua
[考博] 申博自荐 +4 Linxia林夏 2026-04-13 4/200 2026-04-19 19:55 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 291求调剂 +12 关忆北. 2026-04-14 13/650 2026-04-19 16:50 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 6/300 2026-04-19 13:00 by Aaron_zyn
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
信息提示
请填处理意见