急求??? Error using ==> mtimes Inner matrix dimensions must agree.

shicaixia

function shuipingyidong
c=13;m=3.1;q=0.89;a=8;h=304;
v1=0.92;r=159;b=0.27;
w0=m*q*1000*cosd(a);
t1=30/0.92;
t(1)=60/360;
p1=1-exp(-c.*t(1));
p2=1-exp(-c.*(t(1)-t/360));
x=-200:20:500;
u1=b*w0*(p1*(exp(-pi*x.^2/r.^2)-exp(-pi*(x-v1*t1).^2/r.^2))+...
    p2*(exp(-pi*(x-v1*t1).^2/r.^2)-exp(-pi*(x-2*v1*t1).^2/r.^2)));
t(2)=120/360;
p1=1-exp(-c.*t(2));
p2=1-exp(-c.*(t(2)-t/360));
p3=1-exp(-c.*(t(2)-2*t/360));
p4=1-exp(-c.*(t(2)-3*t/360));
x=-200:20:500;
u2=b.*w0.*(p1.*(exp(-pi.*x.^2./r.^2)-exp(-pi.*(x-v1.*t1).^2./r.^2))+...
    p2*(exp(-pi.*(x-v1*t1).^2/r.^2)-exp(-pi.*(x-2.*v1.*t1).^2/r.^2))+...
    p3*(exp(-pi.*(x-2*v1*t1).^2/r.^2)-exp(-pi.*(x-3.*v1.*t1).^2/r.^2))+...
    p4*(exp(-pi.*(x-3*v1*t1).^2/r.^2)-exp(-pi.*(x-4.*v1.*t1).^2/r.^2)));
figure;
plot(x,u1,'-ro',x,u2,'--mx');
xlabel('→工作面开采方向/m');
ylabel('水平移动值U/mm');
grid on
运行结果错误如下

??? Error using ==> mtimes
Inner matrix dimensions must agree.
Error in ==> shuipingyidong at 18
u2=b.*w0.*(p1.*(exp(-pi.*x.^2./r.^2)-exp(-pi.*(x-v1.*t1).^2./r.^2))+...

ChaChing

Ref
常见的程序出错问题整理 (eight) 5F
http://forum.vibunion.com/forum/thread-46001-1-1.html

shicaixia

换是没解决
你能指出具体错的地方吗

[ 本帖最后由 ChaChing 于 2009-9-26 16:36 编辑 ]

ChaChing

u2=b.*w0.*(p1.*(exp(-pi.*x.^2./r.^2)-exp(-pi.*(x-v1.*t1).^2./r.^2))+...
    p2*(exp(-pi.*(x-v1*t1).^2/r.^2)-exp(-pi.*(x-2.*v1.*t1).^2/r.^2))+...
    p3*(exp(-pi.*(x-2*v1*t1).^2/r.^2)-exp(-pi.*(x-3.*v1.*t1).^2/r.^2))+...
    p4*(exp(-pi.*(x-3*v1*t1).^2/r.^2)-exp(-pi.*(x-4.*v1.*t1).^2/r.^2)));

p2/p3/p4皆为1*2, 而後面的表示式皆为1*36, 请问怎麽点乘!? :@)
已非编程问题了, 建议楼主检查下计算流程!

[ 本帖最后由 ChaChing 于 2009-9-26 16:47 编辑 ]

shicaixia

已经解决了，谢谢各位

ChaChing

建议楼主解决问题后回帖分享你的经验和心得