24小时热门版块排行榜    

查看: 3609  |  回复: 16

chinhuku

禁虫 (小有名气)

本帖内容被屏蔽

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

lds588

铁杆木虫 (著名写手)

愚钝

【答案】应助回帖

★ ★ ★ ★ ★
臭水沟(金币+5): 给力!给5个金币!理由:我们就需要像你这样的热心的虫子! 2011-05-28 09:07:21
CODE:
clear,clc
p=[0  0.001  0.002  0.003  0.004  0.005  0.006  0.008  0.011  0.014  0.019  0.023  0.03  0.037  0.048  0.06  0.08  0.1  0.136  0.17  0.217  0.279  0.388  0.444];
c=[0  0.307598804  0.538975997  0.777628765  0.959350561  1.069744162  1.190592446  1.340828541  1.512368753  1.643259589  1.769609554  1.903576465  2.027165658  2.151392907  2.294367083  2.425558327 2.598588163  2.723808495  2.912380082  3.049909679  3.206076964 3.386393547  3.644346885  3.758352826];
cc=c';
pp=p';
F=@(q,x)(x(:,1)-x(:,2)/q(1)).*(10.^(q(2)*(x(:,1)-x(:,2)/q(1)))+q(3))-x(:,2);   
x=[cc pp];

q0=[1    1     1];
warning off
q=nlinfit(x,zeros(size(x,1),1),F,q0);
disp('H A B分别为')
disp(num2str(q));
plot(x(:,1),x(:,2),'r^');hold on;
ezplot(@(x,y)F(q,[x,y]),[0,3.8,-0.1,0.5]);

[url=http://weibo.com/u/2577306830?s=6uyXnP][img]http://service.t.sina.com.cn/widget/qmd/2577306830/fa2fdb00/1.png[/img][/url]
2楼2011-05-28 00:27:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lds588

铁杆木虫 (著名写手)

愚钝

【答案】应助回帖



为拟合得到图形,三角为原图,感觉 你的数据 不太适合用这个公式拟合~~~~
[url=http://weibo.com/u/2577306830?s=6uyXnP][img]http://service.t.sina.com.cn/widget/qmd/2577306830/fa2fdb00/1.png[/img][/url]
3楼2011-05-28 00:30:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)

★ ★ ★ ★ ★
臭水沟(金币+5): 不错!给5个金币!希望您常来交流!支持我们!3Q! 2011-05-28 09:09:49
本帖内容被屏蔽

4楼2011-05-28 09:07:45
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lds588

铁杆木虫 (著名写手)

愚钝

引用回帖:
Originally posted by forcal at 2011-05-28 09:07:45:
用Forcal求解:
!using["fcopt","math"];
init(::pArray,cArray,max)=
{
  max=24,
  pArray=arrayinitns{max : "0  0.001  0.002  0.003  0.004  0.005  0.006  0.008  0.011  0. ...

这位童鞋,用你拟合得到的图像是这个样子的~~

误差您是怎么得到的?
[url=http://weibo.com/u/2577306830?s=6uyXnP][img]http://service.t.sina.com.cn/widget/qmd/2577306830/fa2fdb00/1.png[/img][/url]
5楼2011-05-28 13:40:26
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)

★ ★ ★
臭水沟(金币+3): 欢迎常来交流!谢谢! 2011-05-28 21:17:32
本帖内容被屏蔽

6楼2011-05-28 17:07:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)

本帖内容被屏蔽

7楼2011-05-28 17:13:49
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lds588

铁杆木虫 (著名写手)

愚钝

【答案】应助回帖

★ ★
臭水沟(金币+2): 谢谢交流!嗯!大家共同学习,共同进步哟!呵呵。。。 2011-05-28 21:36:15
引用回帖:
Originally posted by forcal at 2011-05-28 17:07:05:
不知道您是如何绘制图像的。我写了个验证程序:
[code]!using["fcopt","math"];
init(::pArray,cArray,max)=
{
  max=24,
  pArray=arrayinitns{max : "0  0.001  0.002  0.003 ...

图像 是由你拟合得到的 H A B带入原来的方程,从我贴的代码中解出来的;
我的图像也是这样 解得~~~

如果您能用自己拟合得到的值做出与楼主所给数据接近的图像 当然更好,大家都是一起讨论学习~~~
[url=http://weibo.com/u/2577306830?s=6uyXnP][img]http://service.t.sina.com.cn/widget/qmd/2577306830/fa2fdb00/1.png[/img][/url]
8楼2011-05-28 21:27:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)

Matlab


臭水沟(金币+1): 谢谢参与!欢迎常来交流! 2011-05-29 12:03:32
引用回帖:
Originally posted by forcal at 2011-05-28 17:07:05:
不知道您是如何绘制图像的。我写了个验证程序:
[code]!using["fcopt","math"];
init(::pArray,cArray,max)=
{
  max=24,
  pArray=arrayinitns{max : "0  0.001  0.002  0.003 ...

forcal 很久没见了啊 应该是我离开太久了 准备答辩了 杂七杂八的事情忙啊...
origin拟合隐函数靠自定义还是很吃力的 很多结果不太好 也有不错的
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
9楼2011-05-29 00:25:35
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)

本帖内容被屏蔽

10楼2011-05-29 07:46:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chinhuku 的主题更新
信息提示
请填处理意见