24小时热门版块排行榜    

查看: 2132  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

willow.322

金虫 (初入文坛)

[求助] 采用matlab中nlinfit拟合,出现问题了,不知道是哪错了,该怎么解决,求高手指点一下

我要拟合的模型是一个二元的非线性方程:y=B*x1*(x2^b)+A*(x2^a)

我用matlab拟合的代码和错误如下

>> x1=[10 5 3.3333 2.5 2 2.5 2.5 2.5 2.5];
>> x2=[30 30 30 30 30 10 20 40 50];
>> y=[0.8776 0.5662 0.4278 0.3060 0.2332 0.6769 0.4029 0.2441 0.2329];
>> myfunc=inline('beta(1)*x1*x2^beta(2)+beta(3)*x2^beta(4)','beta','x1','x2');
>> beta=nlinfit(x1,x2,y,myfunc,[0.5 0.5 0.5 0.5]);
??? Undefined command/function 'nlinfit'.

红色部分是出现的问题,这是怎么回事呢?我该怎么改进呢?
我使用的matlab7.0

还有我想问下除了用nlinfit来模拟的话,还有别的方法吗?要是多元非线性的模拟。。。
希望有高手能帮我解决,在这先谢过了
回复此楼
加油发文章,加油!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmc_1979

金虫 (小有名气)

引用回帖:
3楼: Originally posted by willow.322 at 2011-08-24 09:40:43:
最好,将从一开始到出结果的代码都帮忙贴上来,我再给你金币哈,我只学会matlab一点点,就是为了模拟这个模型的,万分感谢哈!

已经是所有的代码了,只有M文件的编写这样的:
    点file->new->M-file.
进入后把:
4楼2011-08-24 17:04:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

wmc_1979

金虫 (小有名气)

【答案】应助回帖

willow.322(金币+20): 得到的结果应该是合理的,就是这个M文件,我不懂,还想请问下具体怎么写M文件? 2011-08-24 09:29:27
主窗口为:
x1=[10 5 3.3333 2.5 2 2.5 2.5 2.5 2.5];
x2=[30 30 30 30 30 10 20 40 50];
y=[0.8776 0.5662 0.4278 0.3060 0.2332 0.6769 0.4029 0.2441 0.2329];
X=[x1' x2']
y=y';
beta=nlinfit(X,y,'fun44',[0.5 0.5 0.5 0.5])
fun44 的编写为(M文件)
function yy=fun44(beta,X)
   yy=beta(1)*X(:,1).*(X(:,2).^beta(2))+beta(3)*X(:,2).^beta(4);
结果为:
beta =

    0.3336   -0.4255    3.7457   -1.0139
不知是否正确.
2楼2011-08-23 23:20:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

willow.322

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by wmc_1979 at 2011-08-23 23:20:21:
主窗口为:
x1=[10 5 3.3333 2.5 2 2.5 2.5 2.5 2.5];
x2=[30 30 30 30 30 10 20 40 50];
y=[0.8776 0.5662 0.4278 0.3060 0.2332 0.6769 0.4029 0.2441 0.2329];
X=[x1' x2']
y=y';
beta=nlinfit(X,y,' ...

最好,将从一开始到出结果的代码都帮忙贴上来,我再给你金币哈,我只学会matlab一点点,就是为了模拟这个模型的,万分感谢哈!
加油发文章,加油!
3楼2011-08-24 09:40:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmc_1979

金虫 (小有名气)

引用回帖:
4楼: Originally posted by wmc_1979 at 2011-08-24 17:04:01:
已经是所有的代码了,只有M文件的编写这样的:
    点file->new->M-file.
进入后把:

function yy=fun44(beta,X)
   yy=beta(1)*X(:,1).*(X(:,2).^beta(2))+beta(3)*X(:,2).^beta(4);
拷进去就可以了,注意用fun44文件名保存就可以了.然后在主窗口输入:
5楼2011-08-24 17:06:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290求调剂 +4 ^O^乜 2026-03-19 4/200 2026-03-20 11:51 by 学员8dgXkO
[考研] 295复试调剂 +6 简木ChuFront 2026-03-19 6/300 2026-03-20 11:32 by lature00
[考研] 08工学调剂 +4 用户573181 2026-03-20 4/200 2026-03-20 10:32 by 朗月清风1
[考研] 0856调剂,是学校就去 +7 sllhht 2026-03-19 8/400 2026-03-20 09:30 by 每天只摆一小会
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 265求调剂 +6 梁梁校校 2026-03-17 6/300 2026-03-20 08:59 by ZHANG0tao
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +3 Charlieyq 2026-03-19 3/150 2026-03-20 08:53 by xingguangj
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
信息提示
请填处理意见