开发铝电解槽非稳态非均一模型的混合编程方法