Question

我希望输出为C的值，在循环中i和t的所有范围内。

当我运行它时，我得到error in sym/subref, L_tilde, idx。

我不知道这意味着什么。

 syms C;

 alphaX=0.05;
 DiffCoef = 5*10^-5

 v = 0.1;
 L = 10; xZones = 100;

 dx = L/xZones;

 T = 150;
 u = 0.1;

 dt= 0.005;
  t = 150;
  D = DiffCoef + (alphaX * u);

 for i = 1:xZones
     for t = 1:xZones

         (C(i,t+dt) - C(i,t))/dt = -u(C(i+1,t +dt) - C(i,t+dt))/dx + D(C(i+1,t+dt) -  2*(C(i,t+dt) + C(i-1,t+dt)))/(dx)^2


     end
 end

Answer 1

假设这是Matlab代码，则存在一些问题：

1）你不能分配两个输出所以声明

(C(i,t+dt) - C(i,t))/dt = ...

是非法的。 Matlab表达式只能返回一个输出（不使用例如“deal”），所以在你的情况下你必须重新排列你的表达式

C(i,t+dt) = ...
C(i,t) = ...

此外，因为这似乎是一个时差公式，通常任何 f（t + dt）都是未知（左手边），而 f（t）是已知值作为前一时间步（右侧），以便

f(t+dt) = dt*old_rhs - f(t)

在你的情况下你已经开了_C（i，t_dt）_这似乎很不寻常，所以请检查你的方程式。

有限元分析，1D

1 个答案: