24小时热门版块排行榜    

查看: 577  |  回复: 0

xwl196

银虫 (小有名气)

[求助] matlab非线性拟合命令nlinfit拟合4个自由参数的相关问题

程序如下:
clear;clc
T=311.11;
myfun=@(p,x)log((p(1)*x(:,2)/(2*p(2))+x(:,1)).*(p(1)-x(:,1))./(x(:,1).*(p(1)-p(1)*x(:,2)/(2*p(2))-x(:,1))))+((7*p(1)*x(:,2)/(2*p(2))-x(:,1))/p(1))*p(3)/T+p(4)/T;
x=[0.0022         0.4059
0.0380         1.0118
0.0840         1.0833
0.1617         1.0640
0.3124         0.9387
0.6497         0.5216
0.7445         0.3405
0.7794         0.2928
0.8026         0.2425
0.8234         0.2163 ];
p0=[0.4043 1.19 80 -1460];
warning off
p=nlinfit(x,zeros(size(x,1),1),myfun,p0);
disp('rhomc,C,Eii,Eis分别为:');
disp(num2str(p));
plot(x(:,1),x(:,2),'ro');hold on;
ezplot(@(x,y)myfun(p,[x,y]),[0,1,0,2]);
title('拟合曲线');legend('样本点','拟合曲线')

拟合结果如下:
rhomc,C,Eii,Eis分别为:
0.156253858607+0.005786753697576i            546230758.1694+48871855.16422i -7.845475715233e-007+2.978369193914e-008i -3.533514947268e-006+1.911710468176e-007i
而符合实际情况的结果应该是rhomc[0,1],C[0,10],Eii[0,X00],Eis[-X000,0]
从拟合曲线看,拟合曲线与实验点严重偏离,有可能是初值的问题,但一直没找到合适的初始值。而且拟合曲线多出一条线,求大神指点问题出现在哪里了?

matlab非线性拟合命令nlinfit拟合4个自由参数的相关问题
1.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 xwl196 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 12/600 2026-08-12 01:27 by 南京中西医
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +4 Tide man 2026-08-11 4/200 2026-08-12 00:39 by liuqidi
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见