应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】Matlab画图程序
71/1返回列表
查看: 628  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

zhangyatao

木虫 (正式写手)

[交流] 【求助】Matlab画图程序

请高手帮忙，用Matlab画出如图所示的c(r,z)二维趋势图。
二阶偏微分方程解出来了，如c(r,z)所示，但是不会用Matlab画图。

[ Last edited by zhangyatao on 2009-4-7 at 09:36 ]
回复此楼

» 猜你喜欢
资源共享，一起飞翔！
1楼 2009-04-06 22:14:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主
★ ★
zhangyatao(金币+2,VIP+0):谢谢支持！ 4-9 10:37
顶上去！看看高手们的解答！
赞一下(6人) 回复此楼
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了！[/img][/url][/center]
2楼2009-04-06 23:52:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

doctorzfj

铁杆木虫 (正式写手)
★ ★
zhangyatao(金币+2,VIP+0):谢谢支持！ 4-9 10:37
顶上去！
回复此楼
3楼2009-04-07 08:51:25
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主
★ ★
zhangyatao(金币+2,VIP+0):谢谢支持！ 4-9 10:38
还没人解决，再没人解决，就应楼主要求取消求助了！
赞一下(5人) 回复此楼
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了！[/img][/url][/center]
4楼2009-04-08 13:02:49
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nie_jia

木虫 (正式写手)
★ ★
zhangyatao(金币+2,VIP+0):谢谢关注！ 4-9 10:39
太难了 不会
回复此楼
5楼2009-04-08 23:38:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长
★ ★ ★ ★
zhangyatao(金币+4,VIP+0):谢谢提供信息！ 4-9 10:38
可以去问问hitzhang，matlab绘图高手
赞一下(5人) 回复此楼
西南交大隧道方向招硕、博士
6楼2009-04-08 23:40:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)
★ ★ ★
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 >> xn=x(l);
>> 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})')
>>
赞一下(10人) 回复此楼
7楼2009-04-10 21:28:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangyatao 的主题更新
71/1返回列表
信息提示
关闭
请填处理意见
关闭