有限元分析，一维

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)_，这看起来很不寻常，所以请检查你的方程式。