Question

matlab中出现维数错误的问题，求大神指点

h=0.1;
t=[0:h:1];
N=length(t);
y=ones(1,N);
ey=ones(1,N);
zy=ones(1,N);
for k=1:N-1
y(1,k+1)=y(1,k)+h*(y(1,k)-2*t/y(1,k)); %预估公式
ey(1,k+1)=ey(1,k)+h*(ey(1,k)-2*t/ey(1,k)); %Euler法
y(1,k+1)=y(1,k)+h/2*(y(1,k)-2*t/y(1,k)+y(1,k+1)-2*t/y(1,k+1)); %改进Euler法
end
zy=sqrt(1+2*t);
plot(t,zy,'-xk',t,ey,'-.*g',t,y,':or')
grid on;
xlabel('t');
ylabel('y');
legend('解析解','Euler法','改进Euler法');
这是用matlab编写的程序，提示说是预估公示那一行出错，错误信息“??? Subscripted assignment dimension mismatch.”，求大神指点该如何修改错误？？？

Answer 1

你好像问过这个问题，我记得是因为你的t是一个向量，而y(1,k)都只是一个数，你现在用一个数减去一个向量，matlab显然无法辨识

追问

我第一次问这个问题的。。。不过刚才已经自己找出问题并且改正啦。。。3Q~