查看: 1466  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[精细化工] 试剂价格 +7 旭必上岸 2024-06-17 10/500 2024-06-20 23:35 by 相与处于陆的鱼
[有机交流] 跑板能跑开,过柱过不纯怎么办 +7 小胡在努力 2024-06-18 9/450 2024-06-20 23:04 by if今年的吧
[找工作] 初始合伙人来啦!(生物试剂耗材标准品) +16 欢快的小科研人 2024-06-15 27/1350 2024-06-20 22:19 by zzw1992127
[论文投稿] ACS AMI 返回审稿意见,一个大修,两个据稿,编辑给的修改重投 +7 智商已更新 2024-06-19 7/350 2024-06-20 19:54 by kinlin13
[博后之家] 在国内某高校做全职博士后2年,现在找到新的单位,出站或退站对新工作有什么影响? +3 nxplfcc 2024-06-20 3/150 2024-06-20 19:50 by zhangysbad
[找工作] 高校两个offer选择 +15 cowox2021 2024-06-18 16/800 2024-06-20 15:24 by ase123456
[考博] 这个博士要读吗 +15 Sea Breeze 2024-06-16 29/1450 2024-06-20 14:30 by Sea Breeze
[基金申请] F03青年基金函评结果 +3 暨阳一只柴 2024-06-19 4/200 2024-06-20 11:15 by 暨阳一只柴
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中! +10 恐龙爸爸 2024-06-14 10/500 2024-06-20 08:42 by jsguo4127
[海外博后] 墨尔本大学博后offer要不要接 +3 kyxblmm 2024-06-18 3/150 2024-06-19 22:39 by blake1111
[论文投稿] Scientific Reports二审90天了,催稿也没反应怎么办? +3 智晨君 2024-06-18 6/300 2024-06-19 20:08 by licheng78
[论文投稿] 审稿 +5 香瓜木香 2024-06-19 6/300 2024-06-19 17:44 by xli1984
[基金申请] 基金得中 +4 woaini0218 2024-06-18 4/200 2024-06-19 17:27 by mengzl
[公派出国] 去英国的小伙伴儿都在哪儿租的房子呀? +7 65syn 2024-06-14 11/550 2024-06-19 10:19 by 65syn
[考博] 2025考博 +8 自强不息a?a 2024-06-15 13/650 2024-06-18 18:12 by 投必得科研顾问
[基金申请] F口401需要啥文章水平 +3 lhjr123 2024-06-16 7/350 2024-06-18 16:05 by hon920603
[高分子] 寻找聚酯反应釜 +3 茕茕恭煮 2024-06-15 6/300 2024-06-18 14:15 by 茕茕恭煮
[基金申请] 面青地会评时间 +8 tanjydd 2024-06-15 8/400 2024-06-17 17:08 by 小龙虾2008
[论文投稿] 投稿被一个审稿人恶意评审了怎么样? +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
信息提示
请填处理意见