24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 mazhanhua_1999 的 20 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

mazhanhua_1999

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 4425.4
帖子: 250
在线: 151.6小时
虫号: 256514
注册: 2006-06-03
专业: 化学工程

[交流] 【求助】Matlab问题

function y=hanshu(k,A);
k=[k(1) k(2) k(3) k(4) k(5) k(6)];
a0=k(1);
a1=k(2);
a2=k(3);
a3=k(4);
a4=k(5);
a5=k(6);

lab=A(:,1);
cap=A(:,2);
dap=A(:,3);
eap=A(:,4);
fap=A(:,5);

xw=(1./(lab.^1.5)-a2*ones(20,1)).*(a3*cap+a4*ones(20,1));

y=(a4*lab*cab*dap-a5*eap)*fap/(1+(a0*lap).^0.5+a1*cap+a2*dap+a3*eap).^2;
x1=[0.03274,0.03274,0.03274,0.03556,0.03556,0.03556,0.03556,0.0388,0.0388,0.0388,0.02963,0.02963,0.02963,0.02963,0.03852,0.03852,0.03852,0.0237,0.0237,0.0237,0.04741,0.04741,0.04741,0.04741,0.05333,0.05333,0.05333,0.05333,0.03556,0.03556,0.03556,0.03556,0.03556,0.03556,0.03556]';
x2=[0.05686,0.05686,0.05686,0.05912,0.05912,0.05912,0.05912,0.06128,0.06128,0.06128,0.04926,0.04926,0.04926,0.04926,0.06404,0.06404,0.06404,0.07882,0.07882,0.07882,0.03941,0.03941,0.03941,0.03941,0.02956,0.02956,0.02956,0.02956,0.05912,0.05912,0.05912,0.05912,0.05912,0.05912,0.05912]';
x3=[3.32,2.986,2.831,3.169,2.602,2.133,1.708,2.526,2.012,1.684,3.164,2.668,2.26,1.878,2.611,2.002,1.65,2.637,2.249,1.981,2.64,2.036,1.697,1.504,2.839,2.372,1.99,1.819,3.137,2.58,2.128,2.76,2.307,1.97,1.677]';
x4=[0.134,0.576,0.777,0.702,1.16,1.771,2.225,1.176,1.772,2.096,0.492,0.977,1.453,2.058,1.106,1.86,2.276,0.984,1.527,1.917,0.899,1.699,2.07,2.342,0.866,1.228,1.702,1.887,0.424,1.098,1.708,0.78,1.424,1.939,2.097]';
x5=[0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,0.00119,4.02E-4,4.02E-4,4.02E-4,4.02E-4,4.02E-4,4.02E-4,4.02E-4,4.02E-4,4.02E-4,4.02E-4,4.02E-4,9.56E-5,9.56E-5,9.56E-5,3.06E-4,3.06E-4,3.06E-4,3.06E-4]';
y=[5.35692,5.25725,5.148,6.12633,5.96017,5.7608,5.47321,6.5989,6.0658,5.6165,6.0301,5.93721,5.80241,5.76895,6.23535,5.99836,5.12932,6.8791,6.69513,6.46914,7.54958,7.316,7.01067,6.56897,7.16209,6.92098,6.60262,6.13282,8.79241,8.63444,8.44674,7.2556,7.14569,7.0219,6.88093]';

A=[x1' x2' x3' x4' x5'];
k0=[1 1 1 1 1 1 ];
[k,r,j]=nlinfit(k,y','hanshu',k0);

y=hanshu(k,A)

运行提示：
??? Input argument "k" is undefined.

Error in ==> Untitled2 at 2
k=[k(1) k(2) k(3) k(4) k(5) k(6)];
问题出在哪里，高手帮忙，看看程序哪里有错误？

[ Last edited by sunxiao on 2009-3-9 at 08:52 ]

回复此楼