| 查看: 1766 | 回复: 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 ] |
回复此楼» 猜你喜欢
南京大学化学调剂
已经有14人回复
-
一志愿北京工业大学，324分求调剂
已经有6人回复
-
一志愿河北工业大学0817化工278分求调剂
已经有14人回复
-
求调剂
已经有5人回复
-
281求调剂
已经有6人回复
-
26申博自荐
已经有6人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有10人回复
-
求调剂 有机化学考研356分
已经有5人回复
-
300求调剂，材料科学英一数二
已经有10人回复
-
材料与化工328分调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助关于Matlab中字符串的转化问题
已经有4人回复
- 求助 matlab交点问题
已经有19人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 【求助】弱弱的问一下Matlab数据提取的问题
已经有5人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB多个mat中 数据画图问题
已经有7人回复
- 【求助】matlab 二次规划的优化的问题
已经有4人回复
- 【求助】matlab 积分的问题
已经有5人回复
- 【求助】问一个关于matlab中legend的问题
已经有11人回复
- 【求助】matlab中保存计算数据的问题
已经有8人回复
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
|
我将如下代码输入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
6楼2012-10-09 21:04:52
Toauto
至尊木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 13024.8
- 散金: 1904
- 红花: 25
- 帖子: 1925
- 在线: 417.2小时
- 虫号: 339744
- 注册: 2007-04-07
- 专业: 文献学
7楼2012-10-10 16:10:50
Toauto
至尊木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 13024.8
- 散金: 1904
- 红花: 25
- 帖子: 1925
- 在线: 417.2小时
- 虫号: 339744
- 注册: 2007-04-07
- 专业: 文献学
8楼2012-10-10 16:17:06
 , 太感谢了，ds函数的代码如下，请问那如何修改呢？function ds = MW2ds( MW, Solute ) MW0 = 35e3; if exist('Solute') ~=1 if MW <= MW0 % for PEG ds = 33.46E-3 .* (MW.^0.557); disp('PEG be used!') else % for PEO ds = 20.88E-3 .* (MW.^0.587); disp('PEO be used!') end else switch upper(Solute) case 'PEG' ds = 33.46E-3 .* (MW.^0.557); case 'PEO' ds = 20.88E-3 .* (MW.^0.587); case 'PVP' ds = 16.80E-3 .* (MW.^0.593); case 'DEX' ds = 66E-3 .* (MW.^0.46); otherwise ds = 'No relationship'; end end |
9楼2012-10-11 18:04:44
Toauto
至尊木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 13024.8
- 散金: 1904
- 红花: 25
- 帖子: 1925
- 在线: 417.2小时
- 虫号: 339744
- 注册: 2007-04-07
- 专业: 文献学
10楼2012-10-15 09:14:44
