所以我有以下两个功能:
delta_t=T/N_time;
delta_x=1/N_space;
rho=delta_t/delta_x^2;
phi=zeros(N_space+1,N_time+1);
phi(:,1)=initial_condition((0:N_space)*delta_x);
for j=1:N_time;
for i=2:N_space;
phi(i,j+1)=rho*(phi(i-1,j)+phi(i+1,j))+(1-2*rho)*phi(i,j);
end;
end;
phi = phi(end,:);
end
当我运行第二个函数时,我收到此错误:
使用==>时出错plot Vectors必须是相同的长度。
我不知道为什么尺寸不同。它们应该是相同的,因为第二个函数调用第一个函数。任何帮助将不胜感激。
答案 0 :(得分:0)
我认为你的phi以某种方式转换。您的x_var
是1 x (N_space + 1)
向量,您的phi是(N_space + 1) x (N_time + 1)
,并且您返回phi(end, :)
,这将是1 x (N_time + 1)
。你的意思是返回phi(:,end)吗?