24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 649  |  回复: 3
【悬赏金币】回答本帖问题,作者wjyyjw将赠送您 30 个金币

wjyyjw

银虫 (初入文坛)

[求助] 拟合不输出参数值已有1人参与

p=[ 293.15        298.15        303.15        308.15        313.15];
c=[0.000378156        0.000425433        0.000562474        0.000684178        0.000981096];
pp=p';
cc=c';
F=@(q,x)(-33509.*(x(:,1)./408.65-1)./(8.314.*x(:,1))+log(x(:,2))+((1- x(:,2)).^2.*((q(1)./(8.314.* x(:,1))).*(exp(-0.3.*q(1)./(8.314.* x(:,1)))./(x(:,2)+(1-x(:,2)).*exp(-0.3.*q(1)./(8.314.* x(:,1))))).^2+(q(2)./(8.314.* x(:,1))).*exp(-0.3.*q(2)./(8.314.* x(:,1)))./(1- x(:,2) + x(:,2) .*exp(-0.3*q(2)./(8.314.* x(:,1)))).^2)));
x=[pp cc];
q0=[0,0];
[q,resnorm] = lsqcurvefit(F,q0, x,zeros(size(x,1),1));
plot(x(:,1),x(:,2),'r^');hold on;
ezplot(@(x,y)F(q,[x,y]) ,[290,314,0.0003,0.001]);
为什么我这个拟合不输出参数值呢
回复此楼
漫漫人生路,行行复行行,今日吾亦往,重走此间路。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与回帖 2019-05-19 21:23:53
其实,系数已经拟合出来了。不显示的原因是  [q,resnorm] = lsqcurvefit(F,q0, x,zeros(size(x,1),1));  这句命令多了一个分号,去掉就可以了。
运行结果
拟合不输出参数值
LHM01.jpg

QQ:2120156492
2楼2019-05-18 11:45:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wjyyjw

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by hzlhm at 2019-05-18 11:45:30
其实,系数已经拟合出来了。不显示的原因是   = lsqcurvefit(F,q0, x,zeros(size(x,1),1));  这句命令多了一个分号,去掉就可以了。
运行结果

LHM01.jpg
...

好的,谢谢
漫漫人生路,行行复行行,今日吾亦往,重走此间路。
3楼2019-05-18 15:30:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

程龙133

新虫 (小有名气)

4楼2019-08-22 09:41:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wjyyjw 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] AC审稿 +3 马winner 2024-04-26 6/300 2024-04-27 09:08 by 马winner
[教师之家] 刚博士毕业,入职普通一本高校 +20 保持肌肤 2024-04-20 21/1050 2024-04-27 06:53 by 铁首星空
[论文投稿] with editor 两个月了,什么原因? +7 yiersan9 2024-04-24 16/800 2024-04-27 03:25 by yiersan9
[找工作] 普通院校药学硕士,做合成的,感觉找不到工作 +10 pom戴墨镜 2024-04-24 17/850 2024-04-26 23:23 by ZZZemmm
[考博] 真的好想读博! +15 wangzhe_bs 2024-04-22 22/1100 2024-04-26 22:11 by 小木雄子
[基金申请] 基金开始函评了吗? +16 wych1103 2024-04-25 16/800 2024-04-26 21:32 by 淀粉搬运工
[考研] 学硕专硕 +5 小蜗牛* 2024-04-26 5/250 2024-04-26 16:43 by 鱼翔浅底1
[基金申请] 两类问题算是白选了~ +7 jurkat.1640 2024-04-23 12/600 2024-04-26 14:39 by lucky_my2010
[考博] 25年博士申请 +6 Changzixuan 2024-04-25 11/550 2024-04-26 13:48 by 我属驴核动力驴
[论文投稿] Chemical Engineering Journal投稿3周了,一直显示With editor状态。这是送审了吗? 10+4 yifeng11 2024-04-20 13/650 2024-04-26 09:48 by yifeng11
[考博] 取博导收留 5+4 zzb777888 2024-04-20 10/500 2024-04-26 08:52 by polymerfriend
[教师之家] 某种做法不行。说过几遍了。同学还那样做。再那样做就给低分 +4 河西夜郎 2024-04-24 4/200 2024-04-26 08:51 by Quakerbird
[论文投稿] Nature一直在编辑手里,考虑好几天了,是悬了吗 +12 彩虹初见 2024-04-24 12/600 2024-04-25 19:21 by 雪径踏青
[基金申请] 国社科项目,你们学校都限额申报吗? +7 屡战屡败 2024-04-21 10/500 2024-04-25 12:10 by 屡战屡败
[访问学者] CSC的访问学者申请,没有个评审意见,也不知道怎么改,还有必要申请吗 +4 flyingship 2024-04-20 4/200 2024-04-24 21:09 by 59038mute
[教师之家] 大家访学都是怎么找的啊? +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +4 我属驴核动力驴 2024-04-22 5/250 2024-04-24 08:47 by 晓目崇
[高分子] 请问UV灯是365nm的,那么选光引发剂的波长选多少的?要完全一致吗? +4 engledd2004 2024-04-21 4/200 2024-04-22 16:08 by wangcz23
[论文投稿] 研究光催化的,好中的三四区 20+3 sl.0117 2024-04-20 3/150 2024-04-22 09:53 by bobvan
[论文投稿] 编辑是选国外的好还是国内的好。 +8 lizhengke06 2024-04-20 8/400 2024-04-22 08:58 by cuiyunjian
信息提示
请填处理意见