24小时热门版块排行榜    

Znn3bq.jpeg
查看: 741  |  回复: 3

robustsong

金虫 (初入文坛)

[求助] 在Matlab里怎么把Vector按照某种比例边短

我现在在做一个胶体中溶剂蒸发模型。胶体底下是铁板,上面是空气,不考虑胶体的面积,溶剂会慢慢从上表面蒸发掉。胶体的厚度也会慢慢变薄,一直到所有的溶剂蒸发完毕。
假设初始状态是:厚度Z0,溶剂浓度是C0,为了分析方便,把Z0均匀的分成了N层,每层的厚度是dZ,溶剂的初始浓度是C0
我是用用最简单的差分算法,这就会涉及到vector的长度缩短问题。
dZ应该是固定的,那么随着Z的变小,N也在变小。那么表达溶剂浓度的vector也在边小。有没有一个方程能够表达这个变化?

比如,初始时,厚度Z=5,dZ设为1,那么一共5层,N=5,C=[0.5 0.5 0.5 0.5 0.5]
用for循环语句做,算完第一次,C=[0.5 0.5 0.4 0.3 0.2],计算厚度变化,Z=4.5,N还是认为5层,
继续循环计算
蒸发了一段时间,C=[0.4 0.3 0.2 0.1 0.01],厚度变为了4,dz还是1,那么N=4,这时候要进行下一轮计算,使用的浓度矩阵变成了一个1X4的了。我想知道如何把一个1X5的矩阵按照比例
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖


感谢参与,应助指数 +1
csgt0: 金币+1, 谢谢 2013-04-15 15:13:25
matlab中的数据类型比较简单,除了不怎么常用的cell和struct等,基本全是matrix,哪怕是一维vector,你也看做是1*length的矩阵就行.所有操作都用下标来,直接用下标可以裁剪vector,给你个例子,你马上就明白了
CODE:
>> c = 1:10

c =

     1     2     3     4     5     6     7     8     9    10

>> c(1:4)

ans =

     1     2     3     4

>> c(1:6)

ans =

     1     2     3     4     5     6

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
2楼2013-04-12 16:21:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

robustsong

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by libralibra at 2013-04-12 16:21:48
matlab中的数据类型比较简单,除了不怎么常用的cell和struct等,基本全是matrix,哪怕是一维vector,你也看做是1*length的矩阵就行.所有操作都用下标来,直接用下标可以裁剪vector,给你个例子,你马上就明白了
>> ...

谢谢回复,不过你没有明白我的意思。
以这个为例子,
转换前是一个1X5的矩阵,C=[0.4 0.3 0.2 0.1 0.01]。 sum(C)=1.01
转换后是一个1X4的矩阵,现在的C可能是这个样子了:C=[0.404 0.303 0.202 0.101],保持sum(C)=1.01不变,
请问如何操作,谢谢了
3楼2013-04-12 22:11:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖


csgt0: 金币+1, 多谢指导 2013-04-16 16:16:13
引用回帖:
3楼: Originally posted by robustsong at 2013-04-12 22:11:20
谢谢回复,不过你没有明白我的意思。
以这个为例子,
转换前是一个1X5的矩阵,C=。 sum(C)=1.01
转换后是一个1X4的矩阵,现在的C可能是这个样子了:C=,保持sum(C)=1.01不变,
请问如何操作,谢谢了...

你还是没讲明白,sum(C)不变从长度5到长度4还是有多种变化方式的.你具体怎么变很重要.
例如你的例子,0.01不见了,把0.01分成0.004,0.003,0.002和0.001加给了前4项,你怎么分这个0.01会造成不同的结果.例如你分成4个0.0025,分别加了后变成
C=[0.4025 0.3025 0.2025 0.1025],sum(C)=1.01不会变的.
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
4楼2013-04-15 16:26:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 robustsong 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] Sci. Bull. 悲剧经验 +4 jyang1999 2026-05-16 5/250 2026-05-20 13:49 by jyang1999
[基金申请] 今年审到国自然15份,谈谈感受 +18 国自然国社科中 2026-05-17 19/950 2026-05-20 13:18 by wjswsl
[基金申请] 提交了我也来说说感想 +4 fummck 2026-05-20 4/200 2026-05-20 12:56 by Equinoxhua
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +9 1234567wang 2026-05-17 11/550 2026-05-20 12:37 by zhuzg0628
[基金申请] 评审有感 +13 popular289 2026-05-18 22/1100 2026-05-20 11:28 by 水和泥不是水泥
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 6/300 2026-05-20 10:36 by xtlilibin
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 8/400 2026-05-20 09:32 by xtlilibin
[基金申请] 河北省自然科学基金 +3 Peterchao 2026-05-18 3/150 2026-05-20 08:57 by 霸_霸
[文学芳草园] 献血感触 +6 呀呀好傻 2026-05-19 11/550 2026-05-19 22:26 by 而立得乐
[基金申请] 2026山东省优青 +3 luckinging 2026-05-14 4/200 2026-05-19 16:15 by 52wls
[考博] 26/27申博自荐-锂/钠电池方向 5+3 狗头军师. 2026-05-15 4/200 2026-05-19 09:10 by moonboat
[基金申请] 同样的基金本子,换个专家直接从C变A! (金币-10) +3 国自然国社科中 2026-05-19 3/150 2026-05-19 08:50 by Equinoxhua
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见