24小时热门版块排行榜    

查看: 3863  |  回复: 14

zhao0610

铁虫 (初入文坛)

[求助] 三维凝固模拟的时候,怎么用matlab画成三维体啊

比如清华大学的赵代平、 荆涛、-柳百成的相场方法模拟铝合金三维枝晶生长一文中的三维枝晶生长图形,不知道matlab能画么?

三维枝晶生长图形
回复此楼

» 收录本帖的淘帖专辑推荐

仿真建模与计算 程序语言箱

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

★ ★ ★
zhao0610(金币+3): ★★★很有帮助 2012-03-11 09:16:40
plot3就可以
给你看个例子
CODE:
clear all, close all

K=0.9;     % try K from 0.7 to 1

c(1,:)=[-K -K -K];
c(2,:)=[-K +K -K];
c(3,:)=[-K -K +K];
c(4,:)=[-K +K +K];
c(5,:)=[+K -K -K];
c(6,:)=[+K +K -K];
c(7,:)=[+K -K +K];
c(8,:)=[+K +K +K];

x(1)=rand()        % random starting point
y(1)=rand()
z(1)=rand()

for j=1:40000  
   
    cc=c(ceil(8*rand()),:);
   
    x(j+1) = cc(1) + (x(j)-cc(1))/((sqrt((x(j)-cc(1))^2+(y(j)-cc(2))^2+(z(j)-cc(3))^2))^2);
    y(j+1) = cc(2) + (y(j)-cc(2))/((sqrt((x(j)-cc(1))^2+(y(j)-cc(2))^2+(z(j)-cc(3))^2))^2);
    z(j+1) = cc(3) + (z(j)-cc(3))/((sqrt((x(j)-cc(1))^2+(y(j)-cc(2))^2+(z(j)-cc(3))^2))^2);
end


plot3(x,y,z,'.','MarkerSize',4) , hold on

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
9楼2012-03-09 23:20:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
cenwanglai(金币+5): ~~ 2012-03-10 08:49:40
如图


The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
10楼2012-03-09 23:22:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

jackzhang19

银虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
dbb627(金币+1): 欢迎交流 2012-03-08 09:01:39
应该会有偏微分方程吧。你先把方程给出来吧,不是每个人都对这个方程熟悉的。
路漫漫其修远兮吾将上下而求索
2楼2012-03-07 21:52:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhao0610

铁虫 (初入文坛)

引用回帖:
: Originally posted by jackzhang19 at 2012-03-07 21:52:31:
应该会有偏微分方程吧。你先把方程给出来吧,不是每个人都对这个方程熟悉的。

他的界面方程么?可以找到,但是找到了偏微分方程之后的大概步骤是什么呢?我只是刚开始查资料,如果是二维的话,就直接用他得到了相场函数(可以反映他的物态)的数值对于计算区域的X、Y数值来画伪色图(pcolor)。但是如果是三维的话,我相当于得到了一个三维矩阵的相场函数值,但是我们能用这个三维矩阵画出一个立体出来么?
3楼2012-03-07 22:05:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhao0610

铁虫 (初入文坛)

哦,对了,这个应该算是四维图了,x,y,z数值是前三维的坐标,第四维是相场函数的值(凝固了是1,液态是-1,-1到1是两者之间)
4楼2012-03-07 22:12:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhao0610

铁虫 (初入文坛)

没人回答么,求助啊
5楼2012-03-08 08:48:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhao0610

铁虫 (初入文坛)

6楼2012-03-09 08:45:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
zhao0610(金币+1): 有帮助 能不能说得具体点啊,拜托 2012-03-09 22:28:35
cenwanglai(金币+2): ~~ 2012-03-10 08:49:30
这个有三维的数据是可以画的,x,y,z数值是前三维的坐标,第四维是相场函数的值(凝固了是1,液态是-1,-1到1是两者之间),第四维用点标记与否体现就行
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
7楼2012-03-09 10:05:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhao0610

铁虫 (初入文坛)

引用回帖:
: Originally posted by dbb627 at 2012-03-09 10:05:20:
这个有三维的数据是可以画的,x,y,z数值是前三维的坐标,第四维是相场函数的值(凝固了是1,液态是-1,-1到1是两者之间),第四维用点标记与否体现就行

具体应该采用哪些函数来实现呢
8楼2012-03-09 11:31:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhao0610 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[文学芳草园] 梦想 +4 myrtle 2026-08-26 4/200 2026-08-31 08:43 by 逸域隐帆
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 能否申诉? +6 echo8914667 2026-08-30 7/350 2026-08-31 08:41 by kingkocxr
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见