版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhiping780

木虫 (著名写手)

天鹰帮帮主

应助: 0 (幼儿园)
贵宾: 0.007
金币: 1584.7
散金: 20
沙发: 2
帖子: 1396
在线: 9.7小时
虫号: 589991
注册: 2008-08-30
性别: GG
专业: 工业催化

[求助] matlab程序运行出错，求指点修改

我仿写了一个Matlab的化学反应动力学计算程序，可是运行时有问题，那位大侠帮忙指导下，谢谢

function fourlump
clear all
clc
global theta yexp V21 V23 V24 A xx0 M MH x0
k0 = [1.0  5.0  0.5];       % 参数初值
lb = [0  0  0];                % 参数下限
ub = [+inf  +inf  +inf]; % 参数上限
xx0=[31.423 42.861 6.403 19.315];
M=[93.99229 83.27657 102.8053 116.2428];
NH=[0.2271 0.2271 0.2271 0.0852 0.0852 0.0852 0.1135 0.1135 .1419 0.1419 0.1419 0.1419 0.1703 0.1703 0.1703];

tspan=[0 1];
ExpData= ...
[1       0.002553       0.002189       0.000629       0.001467
2       0.003354       0.000745       0.000684       0.001973
3       0.002758       0.002128       0.000585       0.001336
4       0.003226       0.003156       0.000672       0.001585
5       0.003495       0.002759       0.000707       0.001658
6       0.003105       0.002974       0.000678       0.001677
7       0.003571       0.001938       0.000791       0.001780
8       0.003073       0.002763       0.000675       0.001594
9       0.002857       0.002761       0.000630       0.001569
10       0.003621       0.001888       0.000679       0.001637
11       0.003228       0.002415       0.000629       0.001522
12       0.003121       0.001553       0.000734       0.001877
13       0.002850       0.002348       0.000636       0.001601
14       0.003406       0.001851       0.000648       0.001570
15       0.003142       0.001809       0.000677       0.001687];
yexp=ExpData(:,2:5);

tspan=[0 1];
V21=[0.8837 1.0897 0.9285 0.8446 0.9119 0.8925 0.9612 0.8992 0.8601 0.9473 0.9115 0.9747 0.8795 0.9365 0.9162]';
V23=[0.8853 0.9914 0.8874 0.8228 0.8501 0.8609 0.9220 0.8669 0.8354 0.8654 0.8592 0.9637 0.8455 0.8638 0.9014]';
V24=[0.7570 0.8668 0.7843 0.7332 0.7551 0.7666 0.7936 0.7678 0.7411 0.7644 0.7638 0.8576 0.7503 0.7658 0.8033]';
theta=[166.8411 250.2616 286.0133 155.4816 248.7706 266.5399 169.4721 593.1523 80.9890 377.9487 566.9230 566.9229 180.9717 361.9435 434.3321];
A=[0.1562 0.1562 0.1562 0.0728 0.0728 0.0728 0.0925 0.0925 0.1106 0.1106 0.1106 0.1106 0.1270 0.1270 0.1270]';
yexp=ExpData(:,2:5);

[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunclump,k0,lb,ub,[],tspan,xx0,x0,theta,V21,V23,V24,yexp)
ci = nlparci(k,resid,jacobian)

function f = ObjFunclump(k,tspan,x0,theta,V21,V23,V24,A,yexp)
for i=1:length(theta)
x0(1)=xx0(1)/((1+2*NH(i))*M(1));
x0(2)=xx0(2)/((1+2*NH(i))*M(2));
x0(3)=xx0(3)/((1+2*NH(i))*M(3));
x0(4)=xx0(4)/((1+2*NH(i))*M(4));
[t x] = ode45(@FCClump,tspan,x0,[],k,theta(i),V21(i),V23(i),V24(i),A(i));
y(i,1) = x(end,1);
y(i,2:4) = x(end,2:4);
end
f1 = y(:,1) - yexp(:,1);
f2 = y(:,2) - yexp(:,2);
f3 = y(:,3) - yexp(:,3);
f4 = y(:,4) - yexp(:,4);
f = [f1;f2;f3;f4];

function dxdt=FCClump(t,x,k,theta,V21,V23,V24,A)

dxdt =  ...
[ theta*k(1)*V21*x(2)/(x(1)+x(2)+x(3)+x(4)+A)
-1*theta*(k(1)+k(2)+k(3))*x(2)/(x(1)+x(2)+x(3)+x(4)+A)
  theta*k(2)*V23*x(2)/(x(1)+x(2)+x(3)+x(4)+A)
  theta*k(3)*V24*x(2)/(x(1)+x(2)+x(3)+x(4)+A)
];

运行后：
??? Error using ==> optim\private\lsqncommon
User supplied function failed with the following error:

Undefined function or variable "y".

Error in ==> lsqnonlin at 147
[x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...

Error in ==> fourlump360 at 39
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...

回复此楼

» 猜你喜欢

2025年遐想已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
自然科学基金委宣布启动申请书“瘦身提质”行动已经有4人回复
求个博导看看已经有18人回复

» 本主题相关价值贴推荐，对您同样有帮助:

关于Matlab中for循环的问题，求达人指点，谢谢了已经有7人回复
用Matlab求一个二阶微分方程，出错了，求指导（求助）已经有13人回复
高手指点：如何使用matlab使数据行与行间空一行已经有6人回复
求高手指点：如何使用matlab 在数据间产生空格？已经有9人回复
matlab 运行ode45出错，不知道什么原因已经有4人回复
matlab求函数积分出错已经有12人回复
matlab一运行就出现的问题已经有7人回复
请大侠指教matlab的这个错误提示到底是个神马意思？已经有12人回复
有个matlab程序，但是运行不下来，求高手已经有23人回复
关于matlab进行傅里叶变换和逆变换的程序，求大神指导！已经有9人回复
运行matlab出现错误已经有4人回复
matlab的m文件关联问题已经有17人回复
大家爱帮我看看这段MATLAB程序为啥出错了已经有5人回复
MATLAB运行错误已经有5人回复
matlab程序问题，总是出错，麻烦大家看看问题在那里？已经有6人回复
MATLAB初学者，这样的式子怎么就是画不了图，总是出错!求指导!!!!!!!!!!!!! 已经有8人回复
求助matlab-simulink运行时报错问题已经有3人回复
同样的MATLAB程序，为何在别人的电脑上运行不了已经有21人回复
【求助】用matlab编程，需要将数据自检分类（高手进来指点下）已经有5人回复

目云天外，志在九霄中。。。。。

1楼 2012-03-16 15:59:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

仿真EPI: 3
应助: 157 (高中生)
贵宾: 5.735
金币: 40609
散金: 2394
红花: 63
沙发: 4
帖子: 4046
在线: 1330.2小时
虫号: 653607
注册: 2008-11-13
性别: GG
专业: 机械动力学
管辖: 仿真模拟

【答案】应助回帖

感谢参与，应助指数 +1

CODE:

function fourlump

clear all

clc

global theta yexp V21 V23 V24 A xx0 M NH x0

k0 = [1.0  5.0  0.5];         % 参数初值

lb = [0  0  0];                   % 参数下限

ub = [+inf  +inf  +inf];    % 参数上限

xx0=[31.423 42.861 6.403 19.315];

M=[93.99229 83.27657 102.8053 116.2428];

NH=[0.2271 0.2271 0.2271 0.0852 0.0852 0.0852 0.1135 0.1135 .1419 0.1419 0.1419 0.1419 0.1703 0.1703 0.1703];

tspan=[0 1];

ExpData= ...

[1        0.002553        0.002189        0.000629        0.001467

2        0.003354        0.000745        0.000684        0.001973

3        0.002758        0.002128        0.000585        0.001336

4        0.003226        0.003156        0.000672        0.001585

5        0.003495        0.002759        0.000707        0.001658

6        0.003105        0.002974        0.000678        0.001677

7        0.003571        0.001938        0.000791        0.001780

8        0.003073        0.002763        0.000675        0.001594

9        0.002857        0.002761        0.000630        0.001569

10        0.003621        0.001888        0.000679        0.001637

11        0.003228        0.002415        0.000629        0.001522

12        0.003121        0.001553        0.000734        0.001877

13        0.002850        0.002348        0.000636        0.001601

14        0.003406        0.001851        0.000648        0.001570

15        0.003142        0.001809        0.000677        0.001687];

yexp=ExpData(:,2:5);

tspan=[0 1];

V21=[0.8837 1.0897 0.9285 0.8446 0.9119 0.8925 0.9612 0.8992 0.8601 0.9473 0.9115 0.9747 0.8795 0.9365 0.9162]';

V23=[0.8853 0.9914 0.8874 0.8228 0.8501 0.8609 0.9220 0.8669 0.8354 0.8654 0.8592 0.9637 0.8455 0.8638 0.9014]';

V24=[0.7570 0.8668 0.7843 0.7332 0.7551 0.7666 0.7936 0.7678 0.7411 0.7644 0.7638 0.8576 0.7503 0.7658 0.8033]';

theta=[166.8411 250.2616 286.0133 155.4816 248.7706 266.5399 169.4721 593.1523 80.9890 377.9487 566.9230 566.9229 180.9717 361.9435 434.3321];

A=[0.1562 0.1562 0.1562 0.0728 0.0728 0.0728 0.0925 0.0925 0.1106 0.1106 0.1106 0.1106 0.1270 0.1270 0.1270]';

yexp=ExpData(:,2:5);

[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...

lsqnonlin(@ObjFunclump,k0,lb,ub,[],tspan,xx0,x0,theta,V21,V23,V24,yexp)

ci = nlparci(k,resid,jacobian)

其中MH为NH。

CODE:

function f = ObjFunclump(k,tspan,x0,theta,V21,V23,V24,A,yexp)

global theta yexp M NH x0

for i=1:length(theta)

    x0(1)=xx0(1)/((1+2*NH(i))*M(1));

    x0(2)=xx0(2)/((1+2*NH(i))*M(2));

    x0(3)=xx0(3)/((1+2*NH(i))*M(3));

    x0(4)=xx0(4)/((1+2*NH(i))*M(4));

    [t x] = ode45(@FCClump,tspan,x0,[],k,theta(i),V21(i),V23(i),V24(i),A(i));

    y(i,1) = x(end,1);

    y(i,2:4) = x(end,2:4);

end

f1 = y(:,1) - yexp(:,1); 

f2 = y(:,2) - yexp(:,2);

f3 = y(:,3) - yexp(:,3);

f4 = y(:,4) - yexp(:,4);

f = [f1;f2;f3;f4];

增加全局变量。

赞一下(1人)

回复此楼

明德厚学、求是创新

2楼2012-03-16 17:00:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

仿真EPI: 3
应助: 157 (高中生)
贵宾: 5.735
金币: 40609
散金: 2394
红花: 63
沙发: 4
帖子: 4046
在线: 1330.2小时
虫号: 653607
注册: 2008-11-13
性别: GG
专业: 机械动力学
管辖: 仿真模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
zhiping780: 金币+100, 谢谢，你能说明或者帮忙修改一下如何用参数方式传递数据么？我对这些一窍不通 2012-03-17 12:05:35

目的函数ObjFunclump中少了个xx0，加上就没错误。
global theta yexp M NH xx0 x0

至于警告，看英文应该能明白意思。

尽量不采用全局变量，建议用参数方式传递数据。

赞一下(1人)

回复此楼

明德厚学、求是创新

4楼2012-03-16 20:35:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 zhiping780 的主题更新

返回列表