版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

lvbida

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 110.3
散金: 144
帖子: 62
在线: 258.3小时
虫号: 864899
注册: 2009-10-07

[交流] 【求助】lsqcurvefit非线性拟合的问题已有1人参与

function F=myfun(p,data)
F=p(1)*(xda(1).^p(2))*(xda(2).^p(3))*(xda(3).^p(4))+p(5)*(xda(4).^p(6))
y=[41.894 115.712 58.319 57.324 67.9495 130.542
71.109 74.447 93.078 128.794 128.573 124.831
107.041 141.092 174.634 217.774 183.355 177.148];
xda(1)=[0.1 0.1 0.1 0.1 0.1 0.1
0.2 0.2 0.2 0.2 0.2 0.2
0.3 0.3 0.3 0.3 0.3 0.3];
xda(2)=[0.1 0.1 0.15 0.15 0.2 0.2
0.1 0.1 0.15 0.15 0.2 0.2
0.1 0.1 0.15 0.15 0.2 0.2 ];
xda(3)=[200 250 200 250 200 250
200 300 200 250 200 300
200 300 250 200 250 300];
xda(4)=[0.095 0.425 0.228 0.124 0.127 0.427
0.110 0.148 0.133 0.302 0.125 0.149
0.134 0.309 0.312 0.315 0.128 0.235];
data=[xda(1);xda(2);xda(3);xda(4)];
p0=[1,2,3,4,5,6]
[p,resnorm]=lsqcurvefit(@myfun,p0,data,y)
请教各位大侠为什么我的程序运行不了，提示
“??? Input argument "x" is undefined.”

回复此楼

» 猜你喜欢

1楼 2010-07-02 17:17:37

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

我的灯呢

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 552.3
帖子: 32
在线: 36.7小时
虫号: 1006625
注册: 2010-04-27
性别: GG
专业: 制造系统与自动化

★ ★
adu886886(金币+2):鼓励应助！欢迎常来仿真模拟板块！ 2010-07-03 22:53:01
anyuezhiji:在:和)中加个空格就不会出现表情了，以前可以在帖子中设置禁用表情的 2010-07-05 17:32:52

我看完你的程序，提出以下几点看法：
第一，程序的前两行应该是单独的一个m文件；
第二，你有xda(1),...,xda(4)分别表示矩阵，这里就有问题；
第三，data和y矩阵的形式应该是一致的，在你的程序里data是12*6，而y是3*6，这肯定不行。
你用xda(1),...,xda(4)是否要做4次计算？如果是这样的话，可以先给出一个12*6矩阵，然后每三行取值做计算，这样循环4次就可以了。
下面是我自己试着做的，仅供参考：
function F=myfun(x,xdata)
F=x(1)*(xdata.^x(2)).*(xdata.^x(3)).*(xdata.^x(4))+x(5)*(xdata.^x(6));

clear;clc
ydata=[41.894 115.712 58.319 57.324 67.9495 130.542
71.109 74.447 93.078 128.794 128.573 124.831
107.041 141.092 174.634 217.774 183.355 177.148];
data=[0.1 0.1 0.1 0.1 0.1 0.1
0.2 0.2 0.2 0.2 0.2 0.2
0.3 0.3 0.3 0.3 0.3 0.3
0.1 0.1 0.15 0.15 0.2 0.2
0.1 0.1 0.15 0.15 0.2 0.2
0.1 0.1 0.15 0.15 0.2 0.2
200 250 200 250 200 250
200 300 200 250 200 300
200 300 250 200 250 300
0.095 0.425 0.228 0.124 0.127 0.427
0.110 0.148 0.133 0.302 0.125 0.149
0.134 0.309 0.312 0.315 0.128 0.235];
x0=[1,2,3,4,5,6];
for i=1:4
fprintf('Case %d:',i)
xdata=[data(3*i-2,: );data(3*i-1,: );data(3*i,: )];
[x,resnorm]=lsqcurvefit(@myfun,x0,xdata,ydata)
end

[ Last edited by anyuezhiji on 2010-7-5 at 17:30 ]

赞一下(2人)