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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金委加强入选青千的管理工作,要求提交每年的出入境记录,避免青千的双聘现象 +18 TRACY王 2018-10-15 18/900 2018-10-17 01:40 by lugia2
[签证指南] 美签J1check整 60天,我怀疑我可能去不成美国了。。。 +18 171421 2018-10-14 21/1050 2018-10-17 00:30 by 残缺的桥段
[考博] 有人想要来日本读硕或博的吗? +21 ddv8888 2018-10-15 42/2100 2018-10-16 23:56 by SeptYe
[论文投稿] 一篇第一次修改返回正在under review,可以再投一篇吗?被主编直接拒掉的可能性大吗? 4+5 lihongweihw 2018-10-16 10/500 2018-10-16 23:06 by nono2009
[教师之家] 学生的哪些行为,是作为老师的你最忍受不了的 +7 傲雪雄鹰 2018-10-16 7/350 2018-10-16 20:47 by Culaisett
[公派出国] 佐治亚理工学院联培 +11 lafar 2018-10-15 14/700 2018-10-16 20:08 by lafar
[教师之家] 教工食堂 (金币+1) +30 无舟野渡 2018-10-12 61/3050 2018-10-16 16:24 by nwsxer
[公派出国] 公派博后 +6 purelifeai 2018-10-15 6/300 2018-10-16 07:50 by Rolarduo
[第一性原理] VASP国内版权购买咨询 (EPI+-1)+3 景鸢 2018-10-12 9/450 2018-10-16 00:57 by evildickman
[考博] 考博华师 +6 笨蛋阁阁巫 2018-10-12 14/700 2018-10-15 19:45 by 镜子看不透
[硕博家园] 有没有在劲牌公司工作或者过去工作过的? +4 Kodar 2018-10-12 5/250 2018-10-15 17:08 by yu27317961
[硕博家园] 现在普通学校毕业的硕士去985读博很难吗? +18 914450032 2018-10-14 33/1650 2018-10-15 14:26 by 哎呦哎呦159
[教师之家] 纠结~留校大学辅导员和军工类研究所 +27 花开两地 2018-10-13 32/1600 2018-10-15 11:16 by 花开两地
[硕博家园] 国奖 +4 xiaoxiao1231 2018-10-11 5/250 2018-10-15 08:01 by 梦想成真2014
[找工作] 有人收到河海大学常州校区面试通知了吗 +3 2018shan2010 2018-10-10 6/300 2018-10-15 05:40 by guohj29
[考博] 现在联系博导来得及么 +3 krkr005 2018-10-12 3/150 2018-10-14 15:42 by 霉霉Daisy
[职场人生] shiya +5 樱雨漫漫1115 2018-10-12 5/250 2018-10-14 13:16 by yclgl2018
[公派出国] 巴黎综合理工学院(école Polytechnique)招收有机化学专业合作渠道攻博生一名 +3 tjxztt 2018-10-12 6/300 2018-10-12 19:40 by 18222957385
[考博] 考博同时报两所高校 +4 水木年华欧欧 2018-10-11 5/250 2018-10-12 10:46 by 1209098482
[硕博家园] 目的基因克隆 +5 1505646848 2018-10-11 5/250 2018-10-12 07:59 by liyaole
信息提示
请填处理意见