同一个版本的matlab、同一个.m文件，为何一个顺利执行、另一个出错？

littlestupid

9 h6 i* X3 u9 _同一个版本7.11.0.584 (R2010b)，安装在一个32位台式机上，和一个64位笔记本上，操作系统都是win7。
$ o0 ^8 z) F' ~7 i) i9 M同一个.m文件，在32位台式机上运行的好好的，但是在64位笔记本上就出错了！报错如下：
??? Error using==> mtimes
$ c" F# r6 K& g" ?Inner matrixdimensions must agree.
7 E" l  k: D( N; SError in ==>evalincaller at 14
: @* v4 g& r" l* ^6 fr =evalin('caller',a);
Error in ==>sym.eval at 33
# O6 F3 K$ X0 Jr =evalin('caller',['evalincaller(''' sym2str(a) ''')']);
# c6 p' o- k, qError in ==>ColdHRBWO7 at 247
   plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
3 A! [! p' \) H- ^! i经过检查，问题出在：
SS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*besselj(1,T(kn,1).*r).*besselj(1,T(kn,1).*r));
上。于是，我做了如下尝试：
$ Z7 P$ u6 Q& G/ s# c0 O(1)     我改成：
SS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*besselj(1,T(kn,1).*r));
  k- a3 @0 l+ y' X在64位笔记本上就不出错了。但真正的公式里面，的确需要两个bessel公式相乘。
' j3 W" ^; M' f1 f3 |/ m(2)     我又把两个bessel之间的.*改成*：
SS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*besselj(1,T(kn,1).*r)*besselj(1,T(kn,1).*r));
64位笔记本还出错！报的错误仍然是“Inner matrix dimensionsmust agree.”
7 d  ~, ~& U$ \- J* {(3)     我又改成如下两种形式：
) L: z& D- f( Q& p5 N6 uSS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*(besselj(1,T(kn,1).*r))^2);
SS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*(besselj(1,T(kn,1).*r)).^2)
64位笔记本继续出错，报的错误又成了：“Inputs must be a scalarand a square matrix.”
; `7 d& s; N/ N
1 o3 K) {; L2 Y出问题的程序段如下：
syms r z;
8 M. H/ n% g* X/ uSS = zeros(2*n+1,1);
SS = sym(SS);
  }/ s* v) ?' z' m6 Yfor kn=1:1:2*n+1
    SS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*besselj(1,T(kn,1).*r).*besselj(1,T(kn,1).*r));
end
%#########################谐波功率#########################;
* Y& A- R0 J; _. c3 opower=zeros(2*n+1,1);
3 ~* U1 C' I# ]( l5 R. YTP = 0;
for kn=1:1:2*n+1
; x- y7 `) r2 b8 P& a# x/ P9 S    clear sp;
+ m  L4 P" M) |" X    sp = @(z,r)real(besselj(1,T(kn,1).*r).*besselj(1,T(kn,1).*r).*r);
2 C; N0 B3 q6 X0 M    power(kn,1) =real(2*pi/p*omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*quad2d(sp,0,0.5*p,0,Rw,'RelTol',1e-5,'MaxFunEvals',9000));
, n% e' O- r, u* ?( B) G2 c    TP = TP + power(kn,1);
3 ?5 _% N% V# ^end
! D  ^3 c4 g0 a3 x0 F& Mpower = power/TP;
save power;
# b  ?) m! ^1 S5 u$ W2 I3 R/ J0 rxlswrite('5-谐波功率.xls',power);
%#########################功率密度#########################;
* _) F! `, e' w; H! ~; j. @) frx = 0:0.01：(ro-d)*1e3;
" |1 c3 [5 ~& o7 {r = rx*1e-3;
z = 0.5*p;
5 u$ d3 s+ P: U: Z) S& k# g( ]$ yGLMD = rx';
figure(100);
# A- R" _3 T" k. Q$ Cfor kn=1:1:2*n+1
1 Q: H; r. H, ~% {    clear GLMDn;
    plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
# H1 t( ~9 y) o4 \; P6 V1 D' H    hold on;
    grid on;
    box on;
7 W9 D7 J1 P' D( z+ F+ _# f    xlabel('r-axis (mm)');
    ylabel('Power flux density (W/m2)');
    GLMDn = eval(SS(kn,1))/TP;
/ p& k% e8 U0 l    GLMDn = GLMDn';
    GLMD = [GLMD,GLMDn];
end
$ a) [; w- r) U" {' T3 Xsave GLMD;
" B- H1 N" S' M; w" }  gxlswrite('5-Poynting.xls',GLMD);

ExxNEN

我的天，把两个文件放在同一目录下，在一个M文件中就可以直接使用另一个文件的函数了啊