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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 求这个函数的积分，求大神给出用MATLAB该怎么写！授我以渔！！！！！
41/1返回列表
查看: 614  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

我不是好汉

木虫 (著名写手)

[求助] 求这个函数的积分，求大神给出用MATLAB该怎么写！授我以渔！！！！！已有1人参与

求大神授我以渔！

求这个函数的积分，求大神给出用MATLAB该怎么写！授我以渔！！！！！
图片1.jpg
回复此楼

» 猜你喜欢
1楼2016-03-31 10:25:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
我不是好汉: 金币+18, ★★★很有帮助 2016-04-01 09:25:49
我不是好汉: 金币+10, ★★★★★最佳答案 2016-04-04 12:59:15
引用回帖:
3楼: Originally posted by 月只蓝 at 2016-03-31 11:03:15
以下代码要完整复制到一个m文件中，再运行。
function qiu_jifen
clear all;clc
format long
jifen=quadgk(@objfun,3.4,10.5)

function f=objfun(t)
T=150.92*t.^(-0.184);
f=exp(-755.8335./(T+70)).*exp ...

或者不需要新建m文件，直接在主程序窗口输入：
CODE:
jifen=quadgk(@(t)exp(-755.8335./((150.92*t.^(-0.184))+70)).*exp(-70000*(143+(150.92*t.^(-0.184)))./2./((150.92*t.^(-0.184)).^2)./(143-(150.92*t.^(-0.184)))),3.4,10.5)

因为T和t存在函数关系，所以这种写法显得很复杂，没有上一楼那种写法清晰。

最后需要说明的是，上述代码仅仅计算的是X（t）中的定积分，要计算X（t）请自行把定积分的结果代入。
赞一下(3人) 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2016-03-31 11:08:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

我不是好汉

木虫 (著名写手)
比较急，自己不懂MATLAB，试了下总是出错，求大神啊！！！！
赞一下 回复此楼
2楼2016-03-31 10:32:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
以下代码要完整复制到一个m文件中，再运行。
CODE:
function qiu_jifen
clear all;clc
format long
jifen=quadgk(@objfun,3.4,10.5)

function f=objfun(t)
T=150.92*t.^(-0.184);
f=exp(-755.8335./(T+70)).*exp(-70000*(143+T)./2./(T.^2)./(143-T));

计算结果：

jifen =

    1.011466025349587e-010
赞一下(1人) 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2016-03-31 11:03:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 我不是好汉 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭