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

littlestupid

同一个版本7.11.0.584 (R2010b)，安装在一个32位台式机上，和一个64位笔记本上，操作系统都是win7。
" }1 {- }' X4 a/ `5 Q- o同一个.m文件，在32位台式机上运行的好好的，但是在64位笔记本上就出错了！报错如下：
5 w" O( b) d9 N- ]; t0 T2 ^??? Error using==> mtimes
( ]0 ]/ u3 l' w. C* b3 SInner matrixdimensions must agree.
Error in ==>evalincaller at 14
0 N4 R1 l+ P+ x4 \! Ir =evalin('caller',a);
3 P4 m6 J; O# CError in ==>sym.eval at 33
r =evalin('caller',['evalincaller(''' sym2str(a) ''')']);
- @9 D( S5 Y) {5 \Error in ==>ColdHRBWO7 at 247
) M9 `3 V& x" p0 ?8 x7 H3 x" s! i   plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
: ^* c6 h+ o1 D% w* K% i, L经过检查，问题出在：
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));
上。于是，我做了如下尝试：
( n4 q( `7 E9 W9 S5 L) a. J  {(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));
9 |# \+ M0 ?" H0 c8 M在64位笔记本上就不出错了。但真正的公式里面，的确需要两个bessel公式相乘。
(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.”
3 x# @9 |( S; f- L; |$ I2 u7 w(3)     我又改成如下两种形式：
" r* O: o, c/ T3 iSS(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)
* @2 r, |- _* `. ^0 S) K64位笔记本继续出错，报的错误又成了：“Inputs must be a scalarand a square matrix.”

% Z) [7 o! X: K6 _出问题的程序段如下：
syms r z;
SS = zeros(2*n+1,1);
SS = sym(SS);
0 d* z. F6 A: @& S! j# f& efor 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));
. O# }% O9 Z; b7 v! `end
%#########################谐波功率#########################;
6 t3 w1 m! f* {) I2 Ipower=zeros(2*n+1,1);
TP = 0;
5 ]; K8 z! p0 n5 ~' i5 q/ @/ }for kn=1:1:2*n+1
    clear sp;
    sp = @(z,r)real(besselj(1,T(kn,1).*r).*besselj(1,T(kn,1).*r).*r);
    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));
* L; j1 L# i; J2 u' U; [    TP = TP + power(kn,1);
) C$ M1 G" H1 N. Rend
power = power/TP;
9 `* w6 E& g" Y, k& psave power;
xlswrite('5-谐波功率.xls',power);
6 d2 O# |, _: w* o%#########################功率密度#########################;
rx = 0:0.01：(ro-d)*1e3;
" c- U( o: X( |. ]5 \5 Xr = rx*1e-3;
  D  V2 P3 |9 D* ez = 0.5*p;
GLMD = rx';
  `: _% N9 z! @/ S/ efigure(100);
2 b$ K" P: B& P7 ~" B9 jfor kn=1:1:2*n+1
    clear GLMDn;
    plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
    hold on;
0 x. A4 h$ P0 W6 K- [    grid on;
& [- W: r- D# x% r- V, g2 w! }$ _3 Q$ [    box on;
    xlabel('r-axis (mm)');
    ylabel('Power flux density (W/m2)');
    GLMDn = eval(SS(kn,1))/TP;
; P4 H, @- c( x" |6 F( Y8 \    GLMDn = GLMDn';
    GLMD = [GLMD,GLMDn];
- M0 S7 a1 g. H& m6 |. jend
save GLMD;
xlswrite('5-Poynting.xls',GLMD);

ExxNEN

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