应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 求助matlab隐函数参数拟合
171/2返回列表12下一页
查看: 4018  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

chinhuku

禁虫 (小有名气)
本帖内容被屏蔽

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-05-27 21:11:31
已阅   关注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]);

赞一下(1人) 回复此楼
[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带入原来的方程，从我贴的代码中解出来的；
我的图像也是这样 解得~~~

如果您能用自己拟合得到的值做出与楼主所给数据接近的图像 当然更好，大家都是一起讨论学习~~~
赞一下(1人) 回复此楼
[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拟合隐函数靠自定义还是很吃力的 很多结果不太好 也有不错的
赞一下(1人) 回复此楼
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_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 的主题更新
171/2返回列表12下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 +5 晴天加油 2026-08-26 6/300 2026-08-27 04:49 by littlezl
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了，中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭