爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 7543|回复: 2

[程序设计] 求助 reshape函数出现问题

[复制链接]

新浪微博达人勋

发表于 2012-6-27 10:08:47 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
大家好,我刚刚接触Matlab,有一段代码交给我调试,可是在运行到reshape函数的时候报错:“Error using reshape
To RESHAPE the number of elements must not change.”我在网上查,说是因为矩阵的大小在转换前后改变了。但是我又瞧不出来哪里改变了,请大家帮忙看一下啊,我读代码感觉这段程序的作用就是将三个文件的内容分别读入三个矩阵,然后再转置一下这三个矩阵,然后再合并这三个矩阵,然后将这个合并后的矩阵内容返回到一个三维数组之中。代码如下:

……
yr=38;np_chn=545;lat1=73;lon1=144;nf=3;np=20;
……
fdir2='E:\data\ncep\gh\gh500\gh_71_08\global\';
fid2=fopen([fdir2,'gh_an_may_71.1_08.7.dat'],'r');      %Model
gh1=fread(fid2,[lat1*lon1 yr],'float32');
fid3=fopen([fdir2,'gh_an_jun_71.1_08.7.dat'],'r');      %Model
gh2=fread(fid3,[lat1*lon1 yr],'float32');
fid4=fopen([fdir2,'gh_an_jul_71.1_08.7.dat'],'r');      %Model
gh3=fread(fid4,[lat1*lon1 yr],'float32');

gh1=permute(gh1,[2 1]);         %[yr lat1*lon1]
gh2=permute(gh2,[2 1]);
gh3=permute(gh3,[2 1]);

fa_thr=cat(2,gh1,gh2,gh3);
fa_thr=reshape(fa_thr,[yr lat1*lon1 nf]);   %[yr lat1*lon1 nf]
……

错误提示是:
Error using reshape
To RESHAPE the number of elements must not change.
Error in choose_fathr_grid (line 105)
fa_thr=reshape(fa_thr,[yr lat1*lon1 nf]);   %[yr lat1*lon1 nf]

请大家帮我看一下,是哪里的错误啊,怎么解决,很急,请大家多多帮助,谢谢!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-6-27 17:32:27 | 显示全部楼层

回帖奖励 +5 金钱

应该是你要reshape 的目标矩阵数据个数(比如100个)和你的[yr lat1*lon1 nf]这个矩阵大小(比如120个)不匹配,你可以自己查看下两者分别数据量是多少个。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-6-28 09:06:33 | 显示全部楼层

是的 我把数据文件弄错了 结果读出来的数据都少了一条记录 编程37了 谢谢您
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表