24小时热门版块排行榜    

查看: 2051  |  回复: 3

younizhehao

新虫 (初入文坛)

[求助] 关于Matlab画三维曲面图的问题,求助高手,感激感激

我想画一个z坐标为去除率(%)的三维曲线,也就是说z坐标需要小于100,但是通过以下程序画出来的图,z坐标出现了负值或大于100,请高手指点,如何得到Z范围在100的三维曲面图。程序如下:

data=[39.4 0 12.9

15.8 450 24.2

19.7 450 25.5

29.6 450 48.4

15.8 650 49.4

19.7 650 58.3

29.6 650 58.1

15.8 900 67.3

19.7 900 58.1

29.6 900 84.7

];

T=data(:,1);

D=data(:,2);

K=data(:,3);

scatter(T,D,5,K)%散点图

[X,Y,Z]=griddata(T,D,K,linspace(min(T),max(T))',linspace(min(D),max(D)),'v4');%插值

figure,surf(X,Y,Z)%三维曲面

我想画一个z坐标为去除率(%)的三维曲线,也就是说z坐标需要小于100,但是通过以下程序画出来的图,z坐标出现了负值或大于100,请高手指点,如何得到Z范围在100的三维曲面图。

data=[39.4 0 12.9

15.8 450 24.2

19.7 450 25.5

29.6 450 48.4

15.8 650 49.4

19.7 650 58.3

29.6 650 58.1

15.8 900 67.3

19.7 900 58.1

29.6 900 84.7

];

T=data(:,1);

D=data(:,2);

K=data(:,3);

scatter(T,D,5,K)%散点图

[X,Y,Z]=griddata(T,D,K,linspace(min(T),max(T))',linspace(min(D),max(D)),'v4');%插值

figure,surf(X,Y,Z)%三维曲面

氨氮去除率 (2).jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2012-11-10 01:08:18
younizhehao: 金币+5, ★★★★★最佳答案 2013-01-21 21:58:23
>> [X,Y,Z]=griddata(T,D,K,linspace(min(T),max(T))',linspace(min(D),max(D)),'nearest');%插值
>> figure,surf(X,Y,Z)%三维曲面
>> [X,Y,Z]=griddata(T,D,K,linspace(min(T),max(T))',linspace(min(D),max(D)),'cubic');%插值
>> figure,surf(X,Y,Z)%三维曲面

这两种方式的插值可得到Z范围在100的三维曲面图

» 本帖已获得的红花(最新10朵)

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
2楼2012-11-09 14:55:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

younizhehao

新虫 (初入文坛)

送鲜花一朵
引用回帖:
2楼: Originally posted by dbb627 at 2012-11-09 14:55:58
>> =griddata(T,D,K,linspace(min(T),max(T))',linspace(min(D),max(D)),'nearest');%插值
>> figure,surf(X,Y,Z)%三维曲面
>> =griddata(T,D,K,linspace(min(T),max(T))',linspace(min(D),max ...

谢谢您的回复,这样画的确范围在100内。请问您如果我想把这些数据画成三维曲面图,
data=[39.4     0         4.44
15.8             450            0.45
19.7             450            0.96
29.6             450            1.2
15.8             650            0.31
19.7             650            0.42
29.6             650            0.55
15.8             900            0.27
19.7           900            0.27
29.6            900            0.54
19.7            1600         0.23
];
是用
[X,Y,Z]=griddata(T,D,K,linspace(min(T),max(T))',linspace(min(D),max(D)),'v4');%插值
figure,surf(X,Y,Z)%三维曲面;

还是您的这个

[X,Y,Z]=griddata(T,D,K,linspace(min(T),max(T))',linspace(min(D),max(D)),'cubic');%插值
figure,surf(X,Y,Z)%三维曲面;

有什么区别呢?
3楼2012-11-10 10:28:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖


jjdg: 金币+1, 感谢说明 2012-11-11 01:08:12
引用回帖:
3楼: Originally posted by younizhehao at 2012-11-10 10:28:49
谢谢您的回复,这样画的确范围在100内。请问您如果我想把这些数据画成三维曲面图,
data=;
是用
=griddata(T,D,K,linspace(min(T),max(T))',linspace(min(D),max(D)),'v4');%插值
figure,surf(X,Y,Z)%三维曲面 ...

区别在于采用的插值公式不同
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
4楼2012-11-10 14:24:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 younizhehao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 范进中举一文的中心思想 +5 炎黄贵胄 2026-08-22 6/300 2026-08-24 06:34 by sincosx
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见