| 查看: 752 | 回复: 6 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
zhangyatao木虫 (正式写手)
|
[交流]
【求助】Matlab画图程序
|
||
请高手帮忙，用Matlab画出如图所示的c(r,z)二维趋势图。 二阶偏微分方程解出来了，如c(r,z)所示，但是不会用Matlab画图。  [ Last edited by zhangyatao on 2009-4-7 at 09:36 ] |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有59人回复
-
filecode=后面第一个是大写字母
已经有5人回复
-
时间戳又变了8-15
已经有21人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
6楼2009-04-08 23:40:31
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
2楼2009-04-06 23:52:43
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
4楼2009-04-08 13:02:49
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
★ ★ ★
woshilsh(金币+3,VIP+0):谢谢您的热心解答！ 4-11 07:59
woshilsh(金币+3,VIP+0):谢谢您的热心解答！ 4-11 07:59
|
>> Kol=0.017;D=0.138e-4;V0=0.11;R=4.25e-4;Z=0;c=0.001;c0=c; >> x=0:1e5; >> yy=D*x.*besselj(1,R*x)-Kol*besselj(0,R*x); >> y=yy.^2; >> y1=y(2:100000);y2=y(1:99999);y3=y(3:100001); >> l=logical([0 (y1 >> plot(x,yy);hold on;scatter(xn,zeros(size(xn)));title('\lambda_{n} roots in [0 10^{5}]') >> An=2*c0./(R*besselj(1,R*xn).*(1+D^2*xn.^2/Kol^2).*xn); >> [z r]=meshgrid(linspace(0,0.01,100),linspace(0,0.05,100)); >> zz=z( ;rr=r(;>> AAn=repmat(An,10000,1); >> xxn=repmat(xn,10000,1); >> zzz=repmat(zz,1,14); >> rrr=repmat(rr,1,14); >> C=AAn.*exp(-xxn.^2.*zzz*D/V0).*besselj(0,rrr.*xxn); >> c=sum(C,2); >> c=reshape(c,size(z)); >> surf(z,r,c) >> camlight right; >> lighting phong; >> shading interp >> axis tight,box on >> grid off >> xlabel('{\itz}') >> ylabel('{\itr}') >> zlabel('{\itc}({\itz},{\itr})') >> |
7楼2009-04-10 21:28:34
