| 查看: 769 | 回复: 8 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 mazhanhua_1999 的 20 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
mazhanhua_1999铁虫 (小有名气)
|
[交流]
【求助】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 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
zhangjiandapao
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1238.3
- 帖子: 417
- 在线: 56.4小时
- 虫号: 461886
- 注册: 2007-11-18
- 性别: GG
- 专业: 金属结构材料
6楼2008-09-19 19:13:34
jiann_li84
木虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 7986.5
- 散金: 1353
- 红花: 2
- 帖子: 2289
- 在线: 1032.1小时
- 虫号: 586573
- 注册: 2008-08-15
- 性别: GG
- 专业: 控制理论与方法
2楼2008-09-19 18:24:58
mazhanhua_1999
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4425.4
- 帖子: 250
- 在线: 151.6小时
- 虫号: 256514
- 注册: 2006-06-03
- 专业: 化学工程
3楼2008-09-19 18:35:58
jiann_li84
木虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 7986.5
- 散金: 1353
- 红花: 2
- 帖子: 2289
- 在线: 1032.1小时
- 虫号: 586573
- 注册: 2008-08-15
- 性别: GG
- 专业: 控制理论与方法
4楼2008-09-19 18:54:56
30