24小时热门版块排行榜

返回列表

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

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

[求助] 一个matlab程序，运行不出来

CODE:

function csegout

cmean1=quad8(@fun1,0,10)

cmean2=quad8(@fun2,10,14)

cmean3=quad8(@fun3,14,24)

cmean=cmean1+cmean2+cmean3

function y=fun1(t)

k=2.4e-3;

cA0=1600;

y=cA0.*(-0.000412308.*t.^3+0.00538945.*t.^2-0.00246543.*t-...

    0.000100973)/(1+k*cA0.*t);

function y=fun2(t)

k=2.4e-3;

cA0=1600;

y=cA0.*(-0.019084.*t+0.292621)/(1+k*cA0.*t);

function y=fun3(t)

k=2.4e-3;

cA0=1600;

y=cA0.*(0.0000103077.*t.^3-0.000479624.*t.^2+0.00425519.*t+...

    0.0315744)/(1+k*cA0.*t);

回复此楼

» 猜你喜欢

求助:我三月中下旬出站，青基依托单位怎么办？已经有9人回复
Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
磺酰氟产物，毕不了业了！已经有5人回复
论文终于录用啦！满足毕业条件了已经有25人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有3人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

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

好好学习，天天向上。

1楼 2012-05-01 20:09:45

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

xiegangmai

版主 (职业作家)

我没头衔

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

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
臭水沟: 金币+3, 谢谢交流~~ 2012-05-02 10:16:43

你用的哪个版本的matlab？估计是比较老的了。
在09及以后的版本中，没有quad8函数了，全都改为quadl啦。

程序中的错误很简单，三个函数中的除（/）前都少了一个点。下面代码在我电脑上验证可以计算出结果。

CODE:

function csegout

cmean1 = quadl(@fun1,0,10)

cmean2 = quadl(@fun2,10,14)

cmean3 =quadl(@fun3,14,24)

cmean = cmean1+cmean2+cmean3

function y=fun1(t)

k=2.4e-3;

cA0=1600;

y=cA0.*(-0.000412308.*t.^3+0.00538945.*t.^2-0.00246543.*t-...

    0.000100973)./(1+k*cA0.*t);

function y=fun2(t)

k=2.4e-3;

cA0=1600;

y=cA0.*(-0.019084.*t+0.292621)./(1+k*cA0.*t);

function y=fun3(t)

k=2.4e-3;

cA0=1600;

y=cA0.*(0.0000103077.*t.^3-0.000479624.*t.^2+0.00425519.*t+...

    0.0315744)./(1+k*cA0.*t);

赞一下(1人)

回复此楼

明德厚学、求是创新

3楼2012-05-01 23:05:57

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

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

返回列表