- 积分
- 17
- 注册时间
- 2002-9-9
- 仿真币
-
- 最后登录
- 1970-1-1
|
everest8km兄,以前你给我回过一个三维插值的问题.
http://www.simwe.com/cgi-bin/ut/topic_show.cgi?id=30864&pg=3&bpg=1
现在我又碰到问题了.原来的程序可以运行的,最近运行不知怎么就不可以了.
load data.dat %读入原始数据文件
x=Xdatafile(1:end,1); %将文件的第一列数据(横坐标x)全部赋值给x
y=Xdatafile(1:end,2); %将文件的第二列数据(纵坐标y)全部赋值给y
z=Xdatafile(1:end,3); %将文件的第三列数据(z)全部赋值给z
p=Xdatafile(1:end,4); %将文件的第四列数据(温度值)全部赋值给p
load nodedata.dat %读入插值点坐标文件
XX=Xdatafile(1:end,1); %将文件的第一列数据(插值点横坐标x)全部赋值给XX
YY=Xdatafile(1:end,2); %将文件的第二列数据(插值点纵坐标y)全部赋值给YY
ZZ=Xdatafile(1:end,3); %将文件的第三列数据(插值点z)全部赋值给ZZ
[X,Y,Z]=meshgrid(XX,YY,ZZ); %编织网格,形成三维数组
PI=griddata3(x,y,z,p,X,Y,Z); %进行三维插值
fid=fopen('file.txt', 'w+'); %输出到文件
for I=1:number of xdata; %横坐标XX的个数
for J=1:number of ydata; %纵坐标YY的个数
for K=1:number of zdata
fprintf(fid, '%f\t%f\t%f\t%f\r\n', XX(I,J,K), YY(I,J,K),ZZ(I,J,K),PI(I,J,K));
end
end
end
fclose(fid);
程序运行到for I=1:number of xdata; 处提示出现问题.说Missing operator, comma, or semicolon.
这是怎么回事压.还请你百忙之中能够帮组我一下!谢谢了 |
|