机械必威体育网址

找回密码
注册会员

QQ登录

只需一步,快速开始

搜索
查看: 2810 | 回复: 12
打印 上一主题 下一主题

请问这个定积分是怎么得到的

[复制链接]
跳转到指定楼层
1#
发表于 2016-1-11 18:44:44 | 只看该作者 回帖奖励 | 倒序浏览 | 阅读模式
没弄明白: y8 e/ k- u5 U. _, y

本帖子中包含更多资源

您需要登录才可以下载或查看,没有帐号?注册会员

x
2#
发表于 2016-1-11 18:48:36 | 只看该作者
我平时都用这个软件,不过是在电脑上用。定积分,如果想看步骤的话。有step by step 功能,建议搜一下。
3#
楼主 | 发表于 2016-1-11 19:57:55 | 只看该作者
本帖最后由 阳光小院暖茶 于 2016-1-11 19:59 编辑
* z" B* M2 ]+ Q2 ^7 C$ D4 y8 n
" h3 B4 z# t5 y+ _这个是我手机截图的,我没这个东西。我在纸上试着算这个定积分,什么凑微分啦,换元啦,分部积分啦,都整不出来啊
4#
发表于 2016-1-11 20:09:22 | 只看该作者
阳光小院暖茶 发表于 2016-1-11 19:57
) Z& F' g5 Y+ ?8 P8 O2 f, w这个是我手机截图的,我没这个东西。我在纸上试着算这个定积分,什么凑微分啦,换元啦,分部积分啦,都整不 ...

+ G4 I2 r& b- f" d# Y9 a0 \这里的log是 我们平时写的ln吧???

点评

我猜是的,不然没法弄 发表于 2016-1-12 07:44
5#
发表于 2016-1-11 20:27:51 | 只看该作者
应该是要用级数展开。。。。。。。。这部分知识我忘了、。。。8 \2 @. U8 n" o* _1 r8 r5 G9 l
. W# x7 A: i W6 Y# A! u5 t
补充内容 (2016-1-12 21:35):
3 \5 N3 a2 A" y7 d1 V今天问了问刚考完研的学长,的确得用级数展开才能解出来

本帖子中包含更多资源

您需要登录才可以下载或查看,没有帐号?注册会员

x
6#
发表于 2016-1-11 21:43:46 | 只看该作者
翻一翻数学手册吧6 _: x" P$ g0 X$ Y( Q

点评

翻不到啊 发表于 2016-1-12 08:12
7#
发表于 2016-1-11 21:49:31 | 只看该作者
手动算了一下,也没得出结果,还是软件快啊* F5 ?, X t# f) }& P
) t" M7 i, \- U

本帖子中包含更多资源

您需要登录才可以下载或查看,没有帐号?注册会员

x

点评

恩,不过过程不清楚 发表于 2016-1-12 08:10
8#
发表于 2016-1-11 21:53:13 | 只看该作者
刚刚试了分部积分法,发现没用
9#
发表于 2016-1-11 21:53:21 | 只看该作者
刚刚试了分部积分法,发现没用

点评

恩,似乎要另外的办法 发表于 2016-1-12 08:13
10#
发表于 2016-1-11 22:15:02 | 只看该作者
用数值积分吧
+ S7 k2 f) j5 s, |5 m3 k! j
9 l- n% T6 H* ]( U) Eclear all
4 m$ V2 ?. o4 B3 U! @format long! v- }" s$ X F: }! z* |. W' t
a=0;3 X; C' n7 O9 w! b. L
b=1;
5 n! T, z. k# E2 v4 D5 D! Eepsilon=10^(-6);$ f; X! \. V7 s4 @+ ^
syms x;
8 x7 v( W8 \/ C' jfun=log(x^2 + 1)/(x + 1);
; Y2 n2 H1 z' y( _# [( FHfun=@Remberg;1 s! o R5 `+ `6 I: e# G; w
Ivalue= feval(Hfun,fun,a,b,epsilon);& ~& z" V" R' O+ Z( t7 Y+ P

! c1 }, M8 ] }- R& {%Remberg.m
) s6 i7 }8 j3 b%a,b为积分限,epsilon为精度,s为返回积分值,fun为被积函数
" K; W: m" y. x1 E3 ^4 C%R(n,m)表示计算值,(n-1)为变步长指标,(m-1)为加速次数@2 p" Q9 D2 I- g" v3 s# V
function s=Remberg(fun,a,b,epsilon)
/ v2 T/ A8 S1 csyms x ;
- Z1 ]" F% L: jfvalue=zeros(1,1000);
5 x" q3 N/ H, y: V$ m3 L9 gR=zeros(100,100);6 n' t: R6 G) D! _* _# B5 w
fvaluea=double(subs(fun,x,a));7 ^: E! i! t% B9 t8 E' k# n" f' E
fvalueb=double(subs(fun,x,b));0 y% x! o" V7 x
R(1,1)=(b-a)/2*(fvaluea+fvalueb); %梯形公式
/ \( g( H! H- C% X# O# `km=1;
Z; k# E; _: k- Z& B0 j1 Nfor k1=1:100; %设置一个比较大的循环量) D8 S8 I: Q& f. w
h=(b-a)/(2^(k1));
/ q# d, A+ ?, l# jR(k1+1,1)=1/2*R(k1,1);
1 _4 p9 h/ Q" L! ]for k2=1:2^(k1-1);* W5 y" }6 Y! G, j3 U# W7 U! ]' |
fvalue(2*k2)=double(subs(fun,x,a+(2*k2-1)*h));
. e' g$ p' L: w5 q9 n- A9 {" hR(k1+1,1)=h*fvalue(2*k2)+R(k1+1,1); %变步长值9 l) D. E6 ]$ L: g( L) A. m0 ]
end
/ {& @) n# |7 l2 ?for k3=1:km; %加速计算
# T% \: C% B" n6 BR(k1+1,k3+1)=1/(4^(k3)-1)*(4^(k3)*R(k1+1,k3)-R(k1,k3));! n% }7 K4 |4 n, U9 c2 G
end; c* w/ z7 u2 m$ z$ U
if abs(R(k1+1,km+1)-R(k1+1,km))& o+ j# d; z& D o4 q
s=R(k1+1,km+1);
# Z9 R6 d7 K. G! zbreak;
* _. t0 R! z8 u% y3 qelse% q6 z$ {+ |: K2 q7 x8 h. t2 _
km=km+1;
1 d I, G" l% _end
6 B5 Z2 B" i' h; m% ], w8 {0 N/ Z) m0 }
end
, O% T) N! M; h
/ v/ d7 \6 z+ u, x) H) F3 f
您需要登录后才可以回帖 登录| 注册会员

本版积分规则

小黑屋|手机版|Archiver|机械必威体育网址(京ICP备10217105号-1,京ICP证050210号,浙公网安备33038202004372号)

GMT+8, 2024-5-9 19:15, Processed in 0.058633 second(s), 16 queries , Gzip On.

Powered byDiscuz!X3.4Licensed

? 2001-2017Comsenz Inc.

快速回复 返回顶部 返回列表