| 查看: 2027 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求助一个matlab关于分类讨论的问题
|
|||
|
各位大神，当我用如下代码计算ds = mw2ds(MW, Solute)这一函数s时发现其计算公式全为Solute='PEO'时的计算公式，而实际上应该是当MW《=35000时用Solute='PEG'的计算公式，MW>35000时才用Solute='PEO' 的计算公式。请教各位虫友如何修改，谢谢大家。 expData = [1000 4000 6000 10000 20000 100000 0.209 0.525 0.858 0.953 0.969 0.979 0.244 0.563 0.751 0.923 0.973 0.982 0.126 0.501 0.710 0.860 0.973 0.983 0.105 0.364 0.568 0.827 0.983 0.994 0.058 0.250 0.450 0.727 0.982 0.995]; expData = expData'; MW = expData(:,1); R = expData(:,2:end); [M N] = size(R); if exist('Solute') ~=1 if MW<=35e3 Solute='PEG' else Solute='PEO' end end ds = mw2ds(MW, Solute) [ Last edited by 星梦泪痕 on 2012-9-20 at 14:57 ] |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有70人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助关于Matlab中字符串的转化问题
已经有4人回复
- 求助 matlab交点问题
已经有19人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 【求助】弱弱的问一下Matlab数据提取的问题
已经有5人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB多个mat中 数据画图问题
已经有7人回复
- 【求助】matlab 二次规划的优化的问题
已经有4人回复
- 【求助】matlab 积分的问题
已经有5人回复
- 【求助】问一个关于matlab中legend的问题
已经有11人回复
- 【求助】matlab中保存计算数据的问题
已经有8人回复
|
我将如下代码输入matlab，结果显示的是??? In an assignment A(I) = B, the number of elements in B and I must be the same.请问这是什么原因啊？我觉得问题出在ds(indMW) = mw2ds(MW, Solute)这里 expData = [1000 4000 6000 10000 20000 100000 0.209 0.525 0.858 0.953 0.969 0.979 0.244 0.563 0.751 0.923 0.973 0.982 0.126 0.501 0.710 0.860 0.973 0.983 0.105 0.364 0.568 0.827 0.983 0.994 0.058 0.250 0.450 0.727 0.982 0.995]; expData = expData'; MW = expData(:,1); R = expData(:,2:end); [M N] = size(R); for indMW=1:6 if exist('Solute') ~=1 if MW(indMW)<=35e3 Solute='PEG' else Solute='PEO' end end ds(indMW) = mw2ds(MW, Solute) end |
5楼2012-10-09 21:00:59
Toauto
至尊木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 13024.8
- 散金: 1904
- 红花: 25
- 帖子: 1925
- 在线: 417.2小时
- 虫号: 339744
- 注册: 2007-04-07
- 专业: 文献学
【答案】应助回帖
★
jjdg: 金币+1, 辛苦了 2012-10-09 21:02:07
jjdg: 金币+1, 辛苦了 2012-10-09 21:02:07
|
MW<=35e3返回的是与MW同等大小的逻辑矩阵： 1 1 1 1 1 0 而不是单个逻辑值。所以if MW<=35e3无法比较，也无法成立，无法进入下一句，而是始终跳到else. 要明确，需要什么结构类型的solute. 应该是在不同MW条件下的单个字符串值。 可以加一个for循环， for indMW=1:6 if exist('Solute') ~=1 if MW(indMW)<=35e3 Solute='PEG' else Solute='PEO' end end |
» 本帖已获得的红花（最新10朵）
星梦泪痕
2楼2012-10-08 13:42:56
3楼2012-10-08 18:51:02
Toauto
至尊木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 13024.8
- 散金: 1904
- 红花: 25
- 帖子: 1925
- 在线: 417.2小时
- 虫号: 339744
- 注册: 2007-04-07
- 专业: 文献学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
星梦泪痕: 金币+30, ★★★很有帮助 2012-10-09 20:54:52
jjdg: 金币+1, 感谢继续应助 2012-10-09 21:02:31
星梦泪痕: 金币+30, ★★★很有帮助 2012-10-09 20:54:52
jjdg: 金币+1, 感谢继续应助 2012-10-09 21:02:31
|
for就加在[M N] = size(R);后面 在ds = mw2ds(MW, Solute)加end ds = mw2ds(MW, Solute)改成ds(indMW) = mw2ds(MW, Solute) 最后得到ds为6*1的矩阵 |
» 本帖已获得的红花（最新10朵）
4楼2012-10-09 09:08:27
