- 积分
- 0
- 注册时间
- 2009-5-19
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 TBE_Legend 于 2009-6-6 17:33 编辑
自己编写了一个MATLAB函数,comsol的扩散模型调用的时候出现错误,错误信息如下。哪位高人指点一下啊??拜托!!
错误: 6251
外部函数调用错误
- 函 数: deff
计算变量失败
- 变 量: Dxx_c_di
- 定 义 为: deff(c,1FFC2548)
计算表达式失败
- 表 达 式: diff(diff(((+(-Dxx_c_di*cx-Dxy_c_di*cy)*test(cx)+(-Dyx_c_di*cx-Dyy_c_di*cy)*test(cy)-+(+(Dts_c_di)*ct)*test(c))*dvol),{test@0}),cx)
Left- and right-hand sides are incompatible: Sizes are 1x1x1 and 1x288 respectively.
- Callstack: deff 19 n1(i,j,f)=0.5.*(C1(i,j,f)+sqrt(C1(i,j,f).^2+4.*ni.^2));
- Callstack: deff 19 n1(i,j,f)=0.5.*(C1(i,j,f)+sqrt(C1(i,j,f).^2+4.*ni.^2));
- Callstack: deff 18 if C1(i,j,f)<5*10
- Callstack: deff 17 for f=1:d
- Callstack: deff 16 for j=1:c
- Callstack: deff 15 for i=1:b
- Callstack: deff 2 C1=c.*0.602;%将c的单位mol/m3转换为e+18/cm3%
Exception:
java.lang.IllegalStateException: Parameter pairs already set.
(rethrown as com.femlab.util.FlFatalException)
(rethrown as com.femlab.jni.FlNativeException)
(rethrown as com.femlab.jni.FlNativeException)
(rethrown as com.femlab.jni.FlNativeException)
(rethrown as com.femlab.jni.FlNativeException)
(rethrown as com.femlab.jni.FlNativeException)
(rethrown as com.femlab.jni.FlNativeException)
(rethrown as com.femlab.jni.FlNativeException)
(rethrown as com.femlab.jni.FlNativeException)
(rethrown as com.femlab.jni.FlNativeException)
(rethrown as com.femlab.jni.FlNativeException)
(rethrown as com.femlab.jni.FlNativeException)
Messages:
Fatal error.
- Type: IllegalStateException
Left- and right-hand sides are incompatible: Sizes are 1x1x1 and 1x288 respectively.
- Callstack: deff 19 n1(i,j,f)=0.5.*(C1(i,j,f)+sqrt(C1(i,j,f).^2+4.*ni.^2));
- Callstack: deff 19 n1(i,j,f)=0.5.*(C1(i,j,f)+sqrt(C1(i,j,f).^2+4.*ni.^2));
- Callstack: deff 18 if C1(i,j,f)<5*10
- Callstack: deff 17 for f=1:d
- Callstack: deff 16 for j=1:c
- Callstack: deff 15 for i=1:b
- Callstack: deff 2 C1=c.*0.602;%cmol/m3e+18/cm3%
附deff函数:
- function De=deff(c,T)
- C1=c.*0.602;%将c的单位mol/m3转换为e+18/cm3%
- a=size(C1);
- b=a(1,1);
- c=a(1,2);
- g=length(a);
- if g==2
- d=1;
- else if g==3
- d=a(1,3);
- end
- end
- T1=T+273.15;
- ni=1800.*exp(-0.66./(0.026.*T1./300));
- for i=1:b
- for j=1:c
- for f=1:d
- if C1(i,j,f)<5*10
- n1(i,j,f)=0.5.*(C1(i,j,f)+sqrt(C1(i,j,f).^2+4.*ni.^2));
- else
- p=[1.6*(10^(-5)) 1.6*(10^(-5)).*C1(i,j,f) 1 -C1(i,j,f)];
- n2=roots(p);
- n1(i,j,f)=n2(3);
- end
- end
- end
- end
复制代码 |
|