24小时热门版块排行榜    

查看: 1047  |  回复: 7

wxy尘埃

银虫 (小有名气)

[求助] matlab积分求解

对于下列积分式,如何用matlab来求解,不胜感激!
x.^4*exp(x)/(exp(x)-1).^2  对x积分
回复此楼
微笑吧,为你的曾经拥有~~~
已阅   回复此楼   关注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尘埃 的主题更新
信息提示
请填处理意见