| 查看: 528 | 回复: 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 ] |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
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
mazhanhua_1999
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4425.4
- 帖子: 250
- 在线: 151.6小时
- 虫号: 256514
- 注册: 2006-06-03
- 专业: 化学工程
5楼2008-09-19 19:08:08
zhangjiandapao
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1238.3
- 帖子: 417
- 在线: 56.4小时
- 虫号: 461886
- 注册: 2007-11-18
- 性别: GG
- 专业: 金属结构材料
6楼2008-09-19 19:13:34
mazhanhua_1999
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4425.4
- 帖子: 250
- 在线: 151.6小时
- 虫号: 256514
- 注册: 2006-06-03
- 专业: 化学工程
7楼2008-09-19 19:28:12
zhangjiandapao
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1238.3
- 帖子: 417
- 在线: 56.4小时
- 虫号: 461886
- 注册: 2007-11-18
- 性别: GG
- 专业: 金属结构材料
8楼2008-09-20 18:26:34
dnp
荣誉版主 (知名作家)
小木虫浪子
- 应助: 74 (初中生)
- 贵宾: 18.476
- 金币: 18799.9
- 散金: 2040
- 红花: 77
- 沙发: 6
- 帖子: 7304
- 在线: 593.8小时
- 虫号: 437184
- 注册: 2007-08-26
- 性别: GG
- 专业: 生物信息学
- 管辖: 有机交流
|
见您在置顶帖子matlab问题解答中我对您问题的解答。 修正一点： 在那个帖子中我纠正的函数是nlinfit(x,y,@hanshu,k,A)，现修正为： nlinfit(x,y,@hanshu,k0) |
9楼2008-09-20 20:00:10