| 查看: 6110 | 回复: 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显示图片。 一直实现不了，求大神指导，万分感谢！  1.jpg  2.jpg  3.jpg  4.jpg  5.jpg |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
有关MATLAB的一个非线性拟合曲线的问题，寻找大神，我在这里，你在哪？
已经有5人回复
- matlab非线性拟合求助
已经有10人回复
- 求matlab大神指点，关于logistic拟合
已经有13人回复
- matlab 非线性拟合
已经有16人回复
- Matlab用fsolve求解非线性方程，卡了三四天了跪求帮助！
已经有35人回复
- 如何用Matlab拟合大量散点的边界曲线
已经有5人回复
- MATLAB作图拟合求曲线
已经有8人回复
- matlab求模型参数有错误，求指导
已经有5人回复
- matlab曲线拟合及修正
已经有10人回复
- matlab里拟合后，知道y坐标求x坐标
已经有4人回复
- 关于matlab两个自变量、一个因变量曲线拟合问题
已经有3人回复
- 在进行多元非线性拟合的时候，用matlab拟合的值残差都大于0
已经有14人回复
- 求助用matlab 多元非线性拟合
已经有4人回复
- 求Origin或Matlab高手非线性拟合公式
已经有5人回复
- 毕不了业了，发帖求助Origin8.0多元非线性拟合的问题
已经有9人回复
- 用Matlab做多变量非线性拟合求参数值，求各位大侠鼎力相助！
已经有3人回复
- 关于matlab微分、及非线性拟合问题
已经有5人回复
- MATLAB非线性优化拟合怎么改才正确
已经有3人回复
- 非线性拟合估计方程参数（matlab，spss，origin都可以）
已经有4人回复
- matlab非线性参数拟合问题
已经有7人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
- 【求助】急求关于matlab做非线性拟合的函数
已经有8人回复
- 【原创】多项式和非多项式曲线拟合的MATLAB实现视频教程
已经有10人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
20