查看: 561  |  回复: 5
当前主题已经存档。

lovefei0119

银虫 (小有名气)

[交流] 【求助】做三维图的时候遇到误差函数怎么办

方程
C(x,y,t)=Cs(1-erf(x/(4Dt)1/2)*erf(y/(4Dt)1/2))  (1)
若D=2.5E-11, t=10, Cs=0.5,带入方程 得到
C(x,y)=0.5(1-erf(31600x)*erf(31600y))             (2)
对于误差函数,我直接把已知 的D=2.5E-11, t=10, Cs=0.5代入正确吗?

如果在matlab作图,三维图。怎么做?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
sunxiao(金币+2,VIP+0):谢谢参与,欢迎常来,zhang专家 8-15 08:03
应该正确,你可以用ezsurf函数画图。
2楼2009-08-14 12:34:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lovefei0119

银虫 (小有名气)

[x,y]=meshgrid(0:0.001:0.1);
>> z=0.5*(1-erf(1000*x)*erf(1000*y));
>> surfc(x,y,z);

我用以上程序画出来的图是个平行于xy面的矩形,所有点z值都一样。怎么回事

啊?难道这个程序不对?
3楼2009-08-14 15:39:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuguocheng(金币+1,VIP+0): 谢谢h版的解答 8-14 19:02
[x,y]=meshgrid(0:0.001:0.1);
>> z=0.5*(1-erf(1000*x).*erf(1000*y));%个人理解,这个地方应该用.*不是*
>> surfc(x,y,z);
4楼2009-08-14 16:27:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lovefei0119

银虫 (小有名气)

按照楼上的程序,做出来的图感觉还是不对啊。
理论上,这个方程做出来的图应该是个曲面啊
5楼2009-08-17 09:05:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gjf_yaya

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
我是看到三维建模进来学习的:)
6楼2009-08-20 00:55:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lovefei0119 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 15/750 2026-08-10 19:29 by 要乐观耀哥
[基金申请] filecode +8 documentary 2026-08-10 8/400 2026-08-10 19:20 by loufangrui
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 为什么网上很多人说本周 12号出结果 +3 瞬息宇宙 2026-08-10 3/150 2026-08-10 18:27 by Tide man
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +4 布布和一二 2026-08-10 4/200 2026-08-10 14:50 by 医学老男孩
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
信息提示
请填处理意见