24小时热门版块排行榜    

查看: 2562  |  回复: 29

飞也如铃

新虫 (小有名气)

引用回帖:
20楼: Originally posted by tmdxyz at 2011-12-01 17:18:37:
y=p1+p2/Ln(x)+p3/(Ln(x))^2+p4/(Ln(x))^3+p5/(Ln(x))^4+p6/(Ln(x))^5+p7/(Ln(x))

Function: y = p1+p2/Ln(x)+p3/(Ln(x))^2+p4/(Ln(x))^3+p5/(Ln(x))^4+p6/(Ln(x))^5+p7/(Ln(x))^6+p8/(Ln(x))^7+p9/(Ln(x))^8+ ...

请问大侠,你写的是什么
21楼2011-12-02 15:53:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

otipc

铁虫 (小有名气)

22楼2011-12-02 17:54:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guojunhang

金虫 (正式写手)

引用回帖:
11楼: Originally posted by dbb627 at 2011-11-28 16:55:58:
[code]
x = [0.768;
0.711;
0.685;
0.671;
0.659;
0.657;
0.655;
0.654;
0.653;
0.654];
y = [4.584;
4.375;
4.234;
4.128;
4.041;
3.972;
3.913;
3.861;
3.815;
3.774];
opts = fitoptions ...

程序运行完毕后
        a =       4.167
        b =      0.1139
我想检验这个拟合效果的时候,使用这样一段程序,
但是执行效果显示,很差,说明这个拟合结果严重不靠谱啊。
期待解释。
其实上面也有人指出这个缺陷了。我只是再强调一下。
clc
clear all
x = [0.768;
0.711;
0.685;
0.671;
0.659;
0.657;
0.655;
0.654;
0.653;
0.654];
y = [4.584;
4.375;
4.234;
4.128;
4.041;
3.972;
3.913;
3.861;
3.815;
3.774];
figure;
plot(x,y,'k-')
hold on;
y_fit = a+b*log10(1-b*x)+(1-b)*log10(1-b*x+x);
plot(x,y_fit)
23楼2011-12-03 17:41:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

引用回帖:
23楼: Originally posted by guojunhang at 2011-12-03 17:41:42:
程序运行完毕后
        a =       4.167
        b =      0.1139
我想检验这个拟合效果的时候,使用这样一段程序,
但是执行效果显示,很差,说明这个拟合结果严重不靠谱啊。
期待解释。
其实上面也有 ...

这个x的数据做过归一化处理后再回归的
CODE:
a =  4.167; b = 0.1139;
x = [0.768;
0.711;
0.685;
0.671;
0.659;
0.657;
0.655;
0.654;
0.653;
0.654];
y = [4.584;
4.375;
4.234;
4.128;
4.041;
3.972;
3.913;
3.861;
3.815;
3.774];
figure;
plot(x,y,'k-')
x1=x-mean(x);
x1=x1/std(x1);
hold on;
y_fit = a+b*log10(1-b*x1)+(1-b)*log10(1-b*x1+x1);
plot(x,y_fit,'ro-')

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
24楼2011-12-03 20:10:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guojunhang

金虫 (正式写手)

引用回帖:
24楼: Originally posted by dbb627 at 2011-12-03 20:10:35:
这个x的数据做过归一化处理后再回归的
[code]
a =  4.167; b = 0.1139;
x = [0.768;
0.711;
0.685;
0.671;
0.659;
0.657;
0.655;
0.654;
0.653;
0.654];
y = [4.584;
4.375;
4.234;
4.128;
4. ...

我抖胆再问一句:
a = ? ,  b =?
25楼2011-12-03 21:39:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

引用回帖:
25楼: Originally posted by guojunhang at 2011-12-03 21:39:31:
我抖胆再问一句:
a = ? ,  b =?

a =  4.167; b = 0.1139;这个程序中不是给了吗,有什么不清楚的吗
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
26楼2011-12-03 21:44:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guojunhang

金虫 (正式写手)

引用回帖:
26楼: Originally posted by dbb627 at 2011-12-03 21:44:23:
a =  4.167; b = 0.1139;这个程序中不是给了吗,有什么不清楚的吗

因为再次带回去的时候发现,差别老大了。
我倒是试过了,x经过归一化处理,变成新的变量x1,他们和a =  4.167; b = 0.1139 带入函数,能够看出来,二者能吻合。
也就是说:f(归一化的x, a =  4.167; b = 0.1139) = y
这个从你给出的程序,还有拟合的结果都能找到。

问题是,x和a =  4.167; b = 0.1139 带入函数,二者并不能吻合。
f(原始数据x, a =  4.167; b = 0.1139) != y

这个拟合过程,最终就是要找到a和b的值啊?
现在他们不能完全满足要求啊?
是不是还有其他的工作要做啊?
27楼2011-12-03 21:57:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guojunhang

金虫 (正式写手)

其实现在我也想通了一些了。
楼主估计没有想到,那个给定的表达式的曲线形状能不能保证
“它是单调上升,且凹凸性也与原始数据一致”
估计他心里也没有谱,挑的那个拟合函数那么怪,也不管能不能拟合出来。

估计随便换一个简单一些的表达式,不管是GUI的操作,还是使用fit函数,
只要把
fitoptions里面的,'Normalize','On'改写成Off,
一切问题都不再是问题了。

不知道楼主还在不在?
28楼2011-12-03 22:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飞也如铃

新虫 (小有名气)

消失了几天,帖子里好热闹啊,我就是要用这个公式来模拟这些数据的,这个公式是原始公式。
29楼2011-12-05 08:50:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xch09

铁杆木虫 (正式写手)

送鲜花一朵
太好了1
30楼2012-02-08 15:21:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 飞也如铃 的主题更新
信息提示
请填处理意见