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

littlestupid

1 @5 Y2 K8 }% U6 Q: Z7 F8 x! m同一个版本7.11.0.584 (R2010b)，安装在一个32位台式机上，和一个64位笔记本上，操作系统都是win7。
同一个.m文件，在32位台式机上运行的好好的，但是在64位笔记本上就出错了！报错如下：
??? Error using==> mtimes
Inner matrixdimensions must agree.
4 t# j: \* x0 y3 O) tError in ==>evalincaller at 14
/ J- F. r+ k" E# q: ]r =evalin('caller',a);
- I8 M: A1 M* T* \  kError in ==>sym.eval at 33
; U+ F. }2 P" N$ rr =evalin('caller',['evalincaller(''' sym2str(a) ''')']);
Error in ==>ColdHRBWO7 at 247
   plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
经过检查，问题出在：
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));
2 {: T* }& H# U上。于是，我做了如下尝试：
: l0 j& D; [; b' \( I5 U' N, I* ^(1)     我改成：
9 m" Y( i; l- x% l3 A' B; |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));
在64位笔记本上就不出错了。但真正的公式里面，的确需要两个bessel公式相乘。
' I/ S+ g8 _! t% i3 T(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)     我又改成如下两种形式：
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);
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.”

6 P' x, D& s+ w+ S3 N; l出问题的程序段如下：
  U' \4 g+ q. _- H  ksyms r z;
, {9 u% q  f  T; b$ p( A  \/ DSS = zeros(2*n+1,1);
( g- E, L) P! i6 u2 TSS = sym(SS);
for 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
%#########################谐波功率#########################;
power=zeros(2*n+1,1);
TP = 0;
/ i8 L+ {. p8 Y" f0 V# j' h2 `for kn=1:1:2*n+1
: e% S; M1 d: ^, I6 g    clear sp;
  f1 e7 g! o! {3 m8 y, N    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));
    TP = TP + power(kn,1);
end
power = power/TP;
+ i( N; t1 ^! |" F5 [save power;
xlswrite('5-谐波功率.xls',power);
, ~! |6 s% E1 p%#########################功率密度#########################;
rx = 0:0.01：(ro-d)*1e3;
r = rx*1e-3;
z = 0.5*p;
GLMD = rx';
1 W  Q5 }1 H$ F- ofigure(100);
for kn=1:1:2*n+1
    clear GLMDn;
, n! M% Y& a6 z; x4 r2 @    plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
    hold on;
    grid on;
    box on;
% b( A& O/ w9 T- [2 a6 J$ G) _    xlabel('r-axis (mm)');
  h+ R. v* B0 @. B3 h/ m    ylabel('Power flux density (W/m2)');
    GLMDn = eval(SS(kn,1))/TP;
    GLMDn = GLMDn';
; |1 ^4 ]# o) X    GLMD = [GLMD,GLMDn];
end
save GLMD;
' H# w2 ~2 [$ y, q# Y/ txlswrite('5-Poynting.xls',GLMD);
, w, u, f2 `  a+ v" Z6 _

ExxNEN

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