Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
|
![]() |
#1 |
初级会员
注册日期: 2009-11-29
帖子: 3
声望力: 0 ![]() |
![]()
老是出现??? Undefined function or variable 'p'.,明明定义了p的。程序中调用函数在另外的m文件中定义
clear all;clc global C t=[0:300:2400]; % t:时间,s C=[0.00 0.11 0.19 0.25 0.31 0.44 0.54 0.62 0.65; 0.00 0.10 0.18 0.24 0.28 0.39 0.48 0.52 0.53; 0.10 3.20 1.70 1.30 2.50 3.90 5.20 8.20 9.7]; C=C';%温度150℃的动力学实验数据C(1) &C(2): kmol/m^3, C(3):% %非线性拟合 kk0=[5.1095e-004; 2.8939;9.5453e+003]; %利用suresh模型算出同样温度条件下的速率常数作为其初始值 tspan=[0:300:2400]; %时间阶梯 C0=[0;0;0.1]; % 0.1:氧气在反应体系气相中的初始分率,考虑了环己烷的饱和蒸 %汽压,在实验条件下为50%左右,除去该蒸汽压以及惰性气体所 %占的比例后得出的上述数值 kk=lsqnonlin(@OptObjFunc,kk0,0,inf,[ ],tspan,C0,C); %非线性最小二乘法拟合。调用函数 ci=nlparci(beta,resid,jacobian); %回归系数,残差,雅阁比矩阵 %拟合效果图(实验值与拟合值的比较) M=[ 1 0 0 0; %解微分代数方程组时的系数矩阵 0 1 0 0; 0 0 0 0; 0 0 0 0]; options=odeset; %产生/改变参数结构 options.Mass=M; options.TolRol=le-7;%定义精度 tspan=[5:40]; C0=[0;0;0.1 ];%初值C p=1.2;%MPa G=1.5;%L/min kla=0.14;%1/s; apxl=input('Pls.input gas holdup you got according reaction conditions using hydro_mass_ c:','s'); [t sim,C sim]=ode45(@ljssm,tspan,CO,options,kk,p,G,kla,apxl); plot(t,C exp(:,1),'-ks',t sim,C sim(:,1),':ko',t,C exp(:,2),'-kd',t sim,C sim(:,2),':k+') legend('Exp.CP','Model.CP','Exp.CI','Mode1.CI'),xlabel('时间t/s'),ylabel('浓度C/(kmol/m^3)') |
![]() |
![]() |
![]() |
#2 |
高级会员
注册日期: 2008-09-14
年龄: 43
帖子: 351
声望力: 24 ![]() |
![]()
从提示上看p没有定义 这要看错误在哪个文件
__________________
qq604443022 |
![]() |
![]() |