机械社区

 找回密码
 注册会员

QQ登录

只需一步,快速开始

搜索
查看: 3761|回复: 12

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

[复制链接]
发表于 2016-1-11 18:44:44 | 显示全部楼层 |阅读模式
没弄明白
( I! l% \3 x* \) a$ L

本帖子中包含更多资源

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

x
回复

使用道具 举报

发表于 2016-1-11 18:48:36 | 显示全部楼层
我平时都用这个软件,不过是在电脑上用。定积分,如果想看步骤的话。有step by step 功能,建议搜一下。
回复 支持 反对

使用道具 举报

 楼主| 发表于 2016-1-11 19:57:55 | 显示全部楼层
本帖最后由 阳光小院暖茶 于 2016-1-11 19:59 编辑 ' |. Y% Q; }$ l6 y5 u

7 v6 R/ W5 V6 k9 r- t4 ~9 q这个是我手机截图的,我没这个东西。我在纸上试着算这个定积分,什么凑微分啦,换元啦,分部积分啦,都整不出来啊
回复 支持 反对

使用道具 举报

发表于 2016-1-11 20:09:22 | 显示全部楼层
阳光小院暖茶 发表于 2016-1-11 19:57
3 _0 \9 [0 f  ^$ t  m1 q8 V, _$ ~这个是我手机截图的,我没这个东西。我在纸上试着算这个定积分,什么凑微分啦,换元啦,分部积分啦,都整不 ...

; f  c* L* b& D5 r" R+ W3 I! F, R这里的log是 我们平时写的ln吧???

点评

我猜是的,不然没法弄  发表于 2016-1-12 07:44
回复 支持 反对

使用道具 举报

发表于 2016-1-11 20:27:51 | 显示全部楼层
应该是要用级数展开。。。。。。。。这部分知识我忘了、。。。) a( w) n# B) V: a* `- W$ D/ h

5 r3 w+ B7 v) r  s4 G" Z补充内容 (2016-1-12 21:35):. f- @: _6 _1 f9 L
今天问了问刚考完研的学长,的确得用级数展开才能解出来

本帖子中包含更多资源

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

x
回复 支持 反对

使用道具 举报

发表于 2016-1-11 21:43:46 | 显示全部楼层
翻一翻数学手册吧1 w) s) h' t/ s4 H" }

点评

翻不到啊  发表于 2016-1-12 08:12
回复 支持 反对

使用道具 举报

发表于 2016-1-11 21:49:31 | 显示全部楼层
手动算了一下,也没得出结果,还是软件快啊2 m" h  {8 J0 J2 Z2 Q+ F, Q  N
; R1 @) Q/ t! S7 F, i& N

本帖子中包含更多资源

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

x

点评

恩,不过过程不清楚  发表于 2016-1-12 08:10
回复 支持 反对

使用道具 举报

发表于 2016-1-11 21:53:13 | 显示全部楼层
刚刚试了分部积分法,发现没用
回复 支持 反对

使用道具 举报

发表于 2016-1-11 21:53:21 | 显示全部楼层
刚刚试了分部积分法,发现没用

点评

恩,似乎要另外的办法  发表于 2016-1-12 08:13
回复 支持 反对

使用道具 举报

发表于 2016-1-11 22:15:02 | 显示全部楼层
用数值积分吧  7 x5 I" B! ~2 C. F4 H5 d; @

! b0 X. O9 t" g( q, s" Pclear all1 p1 S" Y) e$ I. b
format long2 X' B0 z; M' Q/ u
a=0;) h) i9 Q& w/ Y5 ^; T, m; X
b=1;
4 @) j. d% q9 f& w  u: M% Tepsilon=10^(-6);/ V2 U/ l. T3 A" ~3 ]% T. d9 T
syms x;7 H( B% \1 F* g; ?( c  Z
fun=log(x^2 + 1)/(x + 1);
% \( E0 z; a/ I2 {( X" r/ kHfun=@ Remberg;% [( z1 T) {  ^" @3 X0 k- k
Ivalue= feval(Hfun,fun,a,b,epsilon);
) [4 Q" P4 W! V# \+ W, ~! ^/ z
8 H  @4 X4 s# b+ {%Remberg.m
! G& p7 ^3 m( k1 k! z: p5 [* H%a,b为积分限,epsilon为精度,s为返回积分值,fun为被积函数3 S0 i8 o2 g5 D
%R(n,m)表示计算值,(n-1)为变步长指标,(m-1)为加速次数6 p  y" s2 g# J: `) A; O: g
function s=Remberg(fun,a,b,epsilon)
4 @& f: a+ f4 p5 l5 q' wsyms x ;% M1 d& q$ T; |9 y2 M2 Z
fvalue=zeros(1,1000);9 E1 Y0 l  f0 ~6 d% v  X
R=zeros(100,100);
. I8 h  L, O! W( ufvaluea=double(subs(fun,x,a));
1 |3 U/ V3 A, B4 h8 c, w' zfvalueb=double(subs(fun,x,b));
8 r$ p+ U# i- N$ uR(1,1)=(b-a)/2*(fvaluea+fvalueb);   %梯形公式, v) S4 x% s3 X2 o! z$ h6 h- U
km=1;' ~! Q9 h& J" K9 w. M
for k1=1:100;   %设置一个比较大的循环量
9 I& Q1 f0 q: J$ ^    h=(b-a)/(2^(k1));
! r) F2 f1 X1 k" Z- d    R(k1+1,1)=1/2*R(k1,1);2 @7 f1 i, I) p8 {( ]
    for k2=1:2^(k1-1);
% M5 ^" R, c* I* r, A; H* V/ P    fvalue(2*k2)=double(subs(fun,x,a+(2*k2-1)*h));- |  G; l$ C2 B  u# g  p$ T
    R(k1+1,1)=h*fvalue(2*k2)+R(k1+1,1);  %变步长值
: U9 q& `3 N& M9 R  i' j) ?    end6 O9 F' z* w, M/ q7 }0 g
  for k3=1:km;  %加速计算3 L, U# P0 x5 S! x, x- G9 x( y
      R(k1+1,k3+1)=1/(4^(k3)-1)*(4^(k3)*R(k1+1,k3)-R(k1,k3));6 T9 Y2 p7 E( `' J9 ^* k
  end
% J% v+ [, s& U" e% d- R  if abs(R(k1+1,km+1)-R(k1+1,km))<epsilon   %控制精度
* e& I/ C8 b$ ~# M+ a' g      s=R(k1+1,km+1);
# [+ @7 {- m( y0 j! b! J6 g% U      break;4 D. N' E" B. C7 e8 d* W$ c
  else& w8 s1 P& d9 f/ s0 _) }
      km=km+1;% a& d* n# v: t& q4 H
  end
. y* y+ a9 j1 ^6 t+ I# {- ]: F; P( I6 A7 b5 M
end
- a1 k; [* }4 f. E( U2 z& w, N) ?4 r/ e4 l9 j/ [
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册会员

本版积分规则

小黑屋|手机版|Archiver|机械社区 ( 京ICP备10217105号-1,京ICP证050210号,浙公网安备33038202004372号 )

GMT+8, 2025-2-19 06:27 , Processed in 0.062069 second(s), 16 queries , Gzip On.

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

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