24小时热门版块排行榜     先丰纳米 买纳米材料就点这里

热门推荐

查看: 550  |  回复: 11
【悬赏金币】回答本帖问题,作者UZippo将赠送您 100 个金币

UZippo

新虫 (初入文坛)

[求助] 怎么画三维图呢,求助,万分感谢已有1人参与

如何用matlab画分段的三维图呢?可以帮我看一下下面这个例子,十分感谢!!!
例子:0<=x<=6;
      0<=y<=4;
当0<=y<=1时,z=x-y;
当1<y<=3时,z=x+2y;
当3<y<=4时,z=x+y;
如何用matlab做出这种分段的三维图呢?

发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有2个 )

chengsi9406

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
月只蓝: 回帖置顶 2017-12-20 09:25:44
画出来的图像是这样子的。
怎么画三维图呢,求助,万分感谢
1513700043(1).png

6楼2017-12-20 00:15:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chengsi9406

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
月只蓝: 金币+10, 应助指数+1, 感谢热心帮助! 2017-12-20 09:25:26
月只蓝: 回帖置顶 2017-12-20 09:25:34
引用回帖:
7楼: Originally posted by UZippo at 2017-12-20 08:36:38
嗯嗯,可不可以把程序发给我看看
...

clear;
[x,y] = meshgrid(linspace(0,6,100),linspace(0,4,100));
z = (x-y).*(y>=0 & y<=1) + (x+2*y).*(y>1 & y<=3) + (x+y).*(y>3 & y<=4);
mesh(x,y,z)
8楼2017-12-20 09:23:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

UZippo

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 18630320334 at 2017-12-19 21:55:52
散点图行么。。。把所有点算出来。。。

如果你会的话可以帮我编一下程序,我看看效果怎样感谢感谢!!!

发自小木虫Android客户端
3楼2017-12-19 21:57:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

18630320334

新虫 (正式写手)

散点图行么。。。把所有点算出来。。。

发自小木虫IOS客户端
2楼2017-12-19 21:55:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

18630320334

新虫 (正式写手)

我在火车上,要程序得后天了,就是求分段函数,把所有xyz求出来,再用plot3或者surf函数画一下。。。

发自小木虫IOS客户端
4楼2017-12-19 21:59:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

UZippo

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by 18630320334 at 2017-12-19 21:59:06
我在火车上,要程序得后天了,就是求分段函数,把所有xyz求出来,再用plot3或者surf函数画一下。。。

好的,我先试试,等你有时间了帮我弄一下就行

发自小木虫Android客户端
5楼2017-12-19 22:54:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

UZippo

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by chengsi9406 at 2017-12-20 00:15:39
画出来的图像是这样子的。

1513700043(1).png

嗯嗯,可不可以把程序发给我看看

发自小木虫Android客户端
7楼2017-12-20 08:36:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

UZippo

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by chengsi9406 at 2017-12-20 09:23:57
clear;
= meshgrid(linspace(0,6,100),linspace(0,4,100));
z = (x-y).*(y&gt;=0 &amp; y&lt;=1) + (x+2*y).*(y&gt;1 &amp; y&lt;=3) + (x+y).*(y&gt;3 &amp; y&lt;=4);
...

谢谢!

发自小木虫Android客户端
9楼2017-12-20 11:26:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

UZippo

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by chengsi9406 at 2017-12-20 09:23:57
clear;
= meshgrid(linspace(0,6,100),linspace(0,4,100));
z = (x-y).*(y&gt;=0 &amp; y&lt;=1) + (x+2*y).*(y&gt;1 &amp; y&lt;=3) + (x+y).*(y&gt;3 &amp; y&lt;=4);
...

你能帮我看看这个问题出在哪了吗,这5段分开画是没问题的,但是按你这种方法合在一起画出来好像有点出入,希望你抽时间帮我看看,谢谢!
0<=x<=6;   0<=t<=0.4;
当0<=t<=0.07;
Fxt='(6.*(-38.59-33.33).*(- 360000./(29.*(29.*x - 210)) - 12000./203)-(6.*(-38.55-33.33).* 285.7143).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1).*0.3-(2.5.*9.81.*6+245).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1+245+2.5.*9.81.*x';
clf
ezsurf(Fxt,[0,6,0,0.07])
view([-15,30])

当0.07<t<=0.12;
Fxt='(6.*(-(-4*pi^2*0.004*2.5.^2*10*(t-0.07).^2./(0.8^2)+pi.^2*0.004*2.5.^2./(10*0.8.^2))*1000-33.33).*(- 360000./(29.*(29.*x - 210)) - 12000./203)-(6.*(-(-4*pi^2*0.004*2.5.^2*10*(t-0.07).^2./(0.8^2)+pi.^2*0.004*2.5.^2./(10*0.8.^2))*1000-33.33).* 285.7143).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1).*0.3-(2.5.*9.81.*6+245).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1+245+2.5.*9.81.*x';
clf
ezsurf(Fxt,[0.07,0.12,0,6])
view([-15,30])

当0.12<t<=0.28;
Fxt='(6.*((2*pi*0.004*2.5*cos(2*pi*2.5*(0.2-t)./0.8)./0.8)*1000-33.33).*(- 360000./(29.*(29.*x - 210)) - 12000./203)-(6.*((2*pi*0.004*2.5*cos(2*pi*2.5*(0.2-t)./0.8)./0.8)*1000-33.33).* 285.7143).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1).*0.3-(2.5.*9.81.*6+245).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1+245+2.5.*9.81.*x';
clf
ezsurf(Fxt,[0.12,0.28,0,6])
view([-15,30])

当0.28<t<=0.33;
Fxt='(6.*(-(-4*pi^2*0.004*2.5.^2*10*(t-0.33).^2./(0.8^2)+pi.^2*0.004*2.5.^2./(10*0.8.^2))*1000 -33.33).*(- 360000./(29.*(29.*x - 210)) - 12000./203)-(6.*(-(-4*pi^2*0.004*2.5.^2*10*(t-0.33).^2./(0.8^2)+pi.^2*0.004*2.5.^2./(10*0.8.^2))*1000 -33.33).* 285.7143).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1).*0.3-(2.5.*9.81.*6+245).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1+245+2.5.*9.81.*x';
clf
ezsurf(Fxt,[0.28,0.33,0,6])
view([-15,30])

当0.33<t<=0.4;
Fxt='(6.*(-38.59-33.33).*(- 360000./(29.*(29.*x - 210)) - 12000./203)-(6.*(-38.55-33.33).* 285.7143).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1).*0.3-(2.5.*9.81.*6+245).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1+245+2.5.*9.81.*x';
clf
ezsurf(Fxt,[0,6,0.33,0.4])
view([-15,30])


下面是合在一起的程序
clear;
[x,t]=meshgrid(linspace(0,6,100),linspace(0,0.4,100));
f=(6.*(-38.59-33.33).*(- 360000./(29.*(29.*x - 210)) - 12000./203)-(6.*(-38.55-33.33).* 285.7143).*(108000000./...
    (29.*(29.*x - 210).^2) - 120000./1421)./2789.1).*0.3-(2.5.*9.81.*6+245).*(108000000./(29.*(29.*x - 210).^2) - ...
    120000./1421)./2789.1+245+2.5.*9.81.*x.*(t>=0&t<=0.07)+(6.*(-(-4*pi^2*0.004*2.5.^2*10*(t-0.07).^2./(0.8^2)+pi.^2*0.004*2.5.^2./(10*0.8.^2))...
    *1000-33.33).*(- 360000./(29.*(29.*x - 210)) - 12000./203)-(6.*(-(-4*pi^2*0.004*2.5.^2*10*(t-0.07).^2./(0.8^2)+pi.^2*0.004*2.5.^2./...
    (10*0.8.^2))*1000-33.33).* 285.7143).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1).*...
    0.3-(2.5.*9.81.*6+245).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1+245+2.5.*9.81.*x.*(t>0.07&t<=0.12)+...
    (6.*((2*pi*0.004*2.5*cos(2*pi*2.5*(0.2-t)./0.8)./0.8)*1000-33.33).*(- 360000./(29.*(29.*x - 210)) - 12000./203)-...
    (6.*((2*pi*0.004*2.5*cos(2*pi*2.5*(0.2-t)./0.8)./0.8)*1000-33.33).* 285.7143).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1).*...
    0.3-(2.5.*9.81.*6+245).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1+245+2.5.*9.81.*x.*(t>0.12&t<=0.28)+...
    (6.*(-(-4*pi^2*0.004*2.5.^2*10*(t-0.33).^2./(0.8^2)+pi.^2*0.004*2.5.^2./(10*0.8.^2))*1000 -33.33).*(- 360000./(29.*(29.*x - 210))...
    - 12000./203)-(6.*(-(-4*pi^2*0.004*2.5.^2*10*(t-0.33).^2./(0.8^2)+pi.^2*0.004*2.5.^2./(10*0.8.^2))*1000 -33.33).* 285.7143).*...
    (108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1).*0.3-(2.5.*9.81.*6+245).*(108000000./(29.*(29.*x - 210).^2)...
    - 120000./1421)./2789.1+245+2.5.*9.81.*x.*(t>0.28&t<=0.33)+(6.*(-38.59-33.33).*(- 360000./(29.*(29.*x - 210)) - 12000./203)-...
    (6.*(-38.55-33.33).* 285.7143).*(108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1).*0.3-(2.5.*9.81.*6+245).*...
    (108000000./(29.*(29.*x - 210).^2) - 120000./1421)./2789.1+245+2.5.*9.81.*x.*(t>0.33&t<=0.4);
clf
mesh(x,t,f)

可能有点乱!!麻烦了

发自小木虫Android客户端
10楼2017-12-20 21:13:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 UZippo 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 我厌烦了一个人去医院,我也很害怕命运待我不公 +64 yyzx2017 2018-01-15 66/3300 2018-01-19 11:29 by erimac
[出国考试] 雅思托福视频资料分享 +14 我习惯就好 2018-01-16 14/700 2018-01-19 10:57 by 莫斯科探案
[公派出国] 2018根特大学CSC PhD申请疑问 +3 ethan_chao 2018-01-17 9/450 2018-01-19 05:27 by yanyujiao16
[材料综合] 模拟铝合金钎焊的气氛保护炉 +4 dhwang55 2018-01-17 5/250 2018-01-18 21:12 by 白天不懂黑夜
[基金申请] 立论依据与研究基础字数限制 +7 otani 2018-01-18 7/350 2018-01-18 20:19 by zhanghui123
[催化] 氮化碳的稳定性问题 +4 1337526523 2018-01-14 9/450 2018-01-18 19:38 by jyang17
[职场人生] 辅导员不想干了怎么办 +20 xkall 2018-01-12 22/1100 2018-01-18 17:56 by 啊哈伊
[基金申请] 基金委会计算中级人员或者研究生参与的所有项目的总时间吗?一年不超过12个月? +3 simonforever 2018-01-18 7/350 2018-01-18 17:15 by simonforever
[考博] 请教下考博英语和雅思哪个更难一些? +9 courage9 2018-01-17 12/600 2018-01-18 13:39 by courage9
[材料综合] 什么设备能测出材料的元素和成分含量 10+5 尘丶埃落定 2018-01-15 15/750 2018-01-18 11:57 by ILKXL520
[精细化工] 玻璃水添加成油性香精了不溶解 5+3 皓轩商贸 2018-01-14 10/500 2018-01-18 11:29 by ozil272
[硕博家园] 实验 +5 梨小馥 2018-01-16 7/350 2018-01-16 22:49 by 梨小馥
[公派出国] 求大神帮忙看下,csc博士联培的邀请函和语言证明开在一起,可以吗? +9 小小小小玩子 2018-01-16 10/500 2018-01-16 20:16 by 两个MM
[考研] 材料专业考研调剂 +3 季鹏程 2018-01-16 6/300 2018-01-16 17:54 by 季鹏程
[考研] 295求调剂 +4 last、炫神 2018-01-16 8/400 2018-01-16 17:52 by last、炫神
[出国考试] 有没有像我一样的?所有成绩都是0分!!! +6 sstonee 2018-01-12 9/450 2018-01-15 20:25 by Leomind
[硕博家园] 求助关于ESI高被引论文的一些问题 +6 suosuohao 2018-01-13 12/600 2018-01-15 14:30 by suosuohao
[论文投稿] 原天全县博轩资料工作室论文诈骗团伙卷土重来,此次建立50个左右的诈骗网站 +4 写手315 2018-01-13 9/450 2018-01-15 13:51 by 写手315
[版块工场] 【早起签'到'贴】2018年1月13日(q q 48626260 进群申请:小木虫——您个人ID) +49 8475 2018-01-13 49/2450 2018-01-13 12:11 by liweiyu001
[版块工场] 【今日提示】冬季中风高发,以下几点要注意! +5 8475 2018-01-13 5/250 2018-01-13 10:45 by 特洛伊木马f
信息提示
请填处理意见