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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » matlab积分求解
81/1返回列表
查看: 1363  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wxy尘埃

银虫 (小有名气)

[求助] matlab积分求解

对于下列积分式，如何用matlab来求解，不胜感激！
x.^4*exp(x)/(exp(x)-1).^2  对x积分
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
微笑吧，为你的曾经拥有~~~
1楼 2012-05-01 20:39:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lurencyj

木虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
xzhdty: 金币+1, 专家考核, 谢谢了 2012-05-02 06:42:47
工具而已，何必执着.
Mathematica:
CODE:
Integrate[x^4 Exp[x]/(Exp[x] - 1)^2, x]

-x^4 - x^4/(-1 + e^x) + 4 x^3 Log[1 - e^x] + 12 x^2 PolyLog[2, e^x] -
24 x PolyLog[3, e^x] + 24 PolyLog[4, e^x]

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
很女子很弓虽大
2楼2012-05-01 22:36:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasonluo

至尊木虫 (知名作家)

中央书籍处书记

【答案】应助回帖

感谢参与，应助指数 +1
matlab的解法：
syms x;
int(x.^4*exp(x)/(exp(x)-1).^2,x)
赞一下 回复此楼
朝为田舍郎，暮登天子堂。将相本无种，学子当自强！
3楼2012-05-02 10:40:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxy尘埃

银虫 (小有名气)
引用回帖:
3楼: Originally posted by Jasonluo at 2012-05-02 10:40:03:
matlab的解法：
syms x;
int(x.^4*exp(x)/(exp(x)-1).^2,x)

请问如何用数值积分，得到一个数值解，谢谢
赞一下 回复此楼
微笑吧，为你的曾经拥有~~~
4楼2012-05-02 15:06:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lidong

金虫 (小有名气)

光明骑士

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wxy尘埃: 金币+10, ★★★★★最佳答案 2012-05-04 20:07:08
%龙贝格外推法求函数积分，结果矩阵左下角数值为精度最高的求积结果
function bgl
a=input('a=');             %输入积分下限
b=input('b=');             %输入积分上限
n=input('n=');              %输入外推次数（n=3~5即可满足一般精度要求）
digits(10);
h=vpa(b-a);
y1=f(a);
y2=f(b);
t(1,1)=vpa(h/2*(y1+y2));
for i=2:n
    s=0;
    for k=1:2^(i-2)
        y0=f(a+2*(2*k-1)*(b-a)/2^i);
        s=s+y0;
    end
    t(1,i)=vpa(t(1,(i-1))/2+h/2*s);
    h=h/2;
end
for m=2:n
    for j=1n+1-m)
        c=4^(m-1);
        t(m,j)=vpa(c/(c-1)*t((m-1),(j+1))-1/(c-1)*t((m-1),j));
    end
end
t
end

function y=f(x)
if x==0
    y=0;
else
    y=vpa(x.^4*exp(x)/(exp(x)-1).^2);
end
end
赞一下 回复此楼
条条大路通罗马，总有一条是属于我的
5楼2012-05-02 16:39:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Carpathia

银虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
qand函数，用matlab帮助搜一下qand，有其他函数可以选择
赞一下 回复此楼
追求理想的道路上东张西望
6楼2012-05-03 10:15:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxy尘埃

银虫 (小有名气)
引用回帖:
5楼: Originally posted by lidong at 2012-05-02 16:39:56:
%龙贝格外推法求函数积分，结果矩阵左下角数值为精度最高的求积结果
function bgl
a=input('a=');             %输入积分下限
b=input('b=');             %输入积分上限
n=input('n=');              %输入外 ...

谢谢！请问表情处是个什么符号？
赞一下 回复此楼
微笑吧，为你的曾经拥有~~~
7楼2012-05-03 21:43:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lidong

金虫 (小有名气)

光明骑士
引用回帖:
7楼: Originally posted by wxy尘埃 at 2012-05-03 21:43:02:
谢谢！请问表情处是个什么符号？

冒号左括号
回复此楼
条条大路通罗马，总有一条是属于我的
8楼2012-05-04 17:05:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wxy尘埃 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 02:35 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:26 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 7/350 2026-11-06 02:24 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 02:13 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:40 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:40 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:31 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 01:29 by qL1cU46YIdK0
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 21/1050 2026-11-06 01:29 by qL1cU46YIdK0
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +7 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 01:29 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 01:29 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 01:20 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-05 22:44 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
关闭
请填处理意见
关闭