数学建模社区-数学中国

标题: 弯管内生成流线，griddata/griddata3原因，无图像 [打印本页]

作者: guohf 时间: 2014-8-14 18:35
标题: 弯管内生成流线，griddata/griddata3原因，无图像
本帖最后由 guohf 于 2014-8-14 18:54 编辑

有二维弯管的坐标及x,y的速度值（散点），先用meshgrid生成规则网格，再用griddata插值速度值，因在弯头外弯管无相应的速度值，在这些区域插值生成的速度为NaN，而画流线图时，图片是空白的。
对三维弯管用griddata3插值时，会提示：
??? Error using ==> griddata3 at 70
Not enough points.
请指点一下，这种不规则区域怎样用matlab生成流线？？？先谢谢了。
下面是二维弯管的数据及程序：

z5clipvelo.rar (674.56 KB, 下载次数: 0)
数据也可到：http://pan.baidu.com/s/1qW0lGjM 下载

clear all;
clc;
z5clipvelo=['z5clipvelo'];
data_al=importdata(z5clipvelo);
dat=data_al.data;
z5clipvelo=dat(:,[2:3,4:5]);
x=z5clipvelo(:,1)*1e3; y=z5clipvelo(:,2)*1e3; u=z5clipvelo(:,3); v=z5clipvelo(:,4);
mm=[min(x) max(x) min(y) max(y)];
del=[(mm(2)-mm(1))/50,(mm(4)-mm(3))/50];
xx=mm(1):del(1):mm(2);yy=mm(3):del(2):mm(4);
[xi,yi]=meshgrid(xx,yy)
ui=griddata(x,y,u,xi,yi,'cubic')
vi=griddata(x,y,v,xi,yi,'cubic');
streamline(xi,yi,ui,vi)

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)