请问这个程序有问题吗？

sabonisli

y1=[]; x=0:0.01:1; [a,b]=size(x);
for i=1:b
    if x(i)<=0.4, y=0;
    elseif x(i)<=0.65
       y=2*((x-0.4)/0.5)*((x-0.4)/0.5);
    else y=1;
    end
    y1=[y1,y];
end

??? Error using ==> mtimes
Inner matrix dimensions must agree

请帮忙看看这个程序哪里有错误，谢谢谢谢～～

[ 本帖最后由 ChaChing 于 2009-12-6 11:36 编辑 ]

ChaChing

y=2*((x-0.4)/0.5).*((x-0.4)/0.5);   ???

cpf12358

你的x是一个行向量,如果用楼上说的点乘,y就变成了一个向量,而非你要的一个数.
我觉得应该是y=2*((x-0.4)/0.5)*((x-0.4)/0.5)';把后面的转置一下就行了.

ChaChing

楼上说的是! 其实这些个人有观察到
但楼主程序的逻辑怪怪的, 又不说清楚原始目的, 个人有点懒得说, 莫怪时间真的有限
一个大for loop其中判断使用x(i), 忽然跑出x, ...???