久久久国产一区二区_国产精品av电影_日韩精品中文字幕一区二区三区_精品一区二区三区免费毛片爱

 找回密碼
 注冊會員

QQ登錄

只需一步,快速開始

搜索
查看: 4492|回復: 12

請問這個定積分是怎么得到的

[復制鏈接]
1#
發表于 2016-1-11 18:44:44 | 只看該作者 |倒序瀏覽 |閱讀模式
沒弄明白3 J8 U6 K0 @. O+ E! U& j

本帖子中包含更多資源

您需要 登錄 才可以下載或查看,沒有賬號?注冊會員

×
回復

使用道具 舉報

2#
發表于 2016-1-11 18:48:36 | 只看該作者
我平時都用這個軟件,不過是在電腦上用。定積分,如果想看步驟的話。有step by step 功能,建議搜一下。
3#
 樓主| 發表于 2016-1-11 19:57:55 | 只看該作者
本帖最后由 陽光小院暖茶 于 2016-1-11 19:59 編輯 . k$ u. T) ?% W, x( f

& r9 Z7 I. L. ?- s" z0 q4 b這個是我手機截圖的,我沒這個東西。我在紙上試著算這個定積分,什么湊微分啦,換元啦,分部積分啦,都整不出來啊
4#
發表于 2016-1-11 20:09:22 | 只看該作者
陽光小院暖茶 發表于 2016-1-11 19:57 3 _+ r$ x( {" X- o& k! Z9 ~
這個是我手機截圖的,我沒這個東西。我在紙上試著算這個定積分,什么湊微分啦,換元啦,分部積分啦,都整不 ...
! {# H2 J/ y8 u7 t8 z+ S/ Z. I* S
這里的log是 我們平時寫的ln吧???

點評

我猜是的,不然沒法弄  發表于 2016-1-12 07:44
5#
發表于 2016-1-11 20:27:51 | 只看該作者
應該是要用級數展開。。。。。。。。這部分知識我忘了、。。。
1 z5 H! b$ L3 n* X& ~9 B* d# H
1 u- c6 y- u9 n% q7 e9 W# [- T, n7 ^補充內容 (2016-1-12 21:35):, D, y; @* k2 B7 I
今天問了問剛考完研的學長,的確得用級數展開才能解出來

本帖子中包含更多資源

您需要 登錄 才可以下載或查看,沒有賬號?注冊會員

×
6#
發表于 2016-1-11 21:43:46 | 只看該作者
翻一翻數學手冊吧1 L( z2 R* Y! s  {* q) j! m/ G

點評

翻不到啊  發表于 2016-1-12 08:12
7#
發表于 2016-1-11 21:49:31 | 只看該作者
手動算了一下,也沒得出結果,還是軟件快啊
4 O8 \! \( X2 p. y0 x( M$ t6 z' F- w. E" j1 S

本帖子中包含更多資源

您需要 登錄 才可以下載或查看,沒有賬號?注冊會員

×

點評

恩,不過過程不清楚  發表于 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 | 只看該作者
用數值積分吧  4 x: A1 @6 B0 P" i

2 {8 h5 h+ [- c, p8 Cclear all
' S5 m2 j' I0 tformat long' x! p! L( o, Z5 Z1 {; L
a=0;- z# M6 a( `1 [* N5 g3 }
b=1;
3 W- D$ ~. g' Y( sepsilon=10^(-6);+ {! e* ~* O  W* Y  U
syms x;" V* ^! F$ I( t0 V" E
fun=log(x^2 + 1)/(x + 1);* a2 o4 P/ ^5 {/ `8 ^1 s/ R
Hfun=@ Remberg;4 E9 n8 l" W+ b, j9 N5 y
Ivalue= feval(Hfun,fun,a,b,epsilon);
9 V! k7 s+ @; B' T+ z. v2 Q; I  ~7 ^+ ~
%Remberg.m( i7 Y- Z, j# u+ p/ P
%a,b為積分限,epsilon為精度,s為返回積分值,fun為被積函數; \; [6 B! C3 q% s' q! |
%R(n,m)表示計算值,(n-1)為變步長指標,(m-1)為加速次數
- v: \& w* o5 J2 p5 ^7 `5 z, J) l3 lfunction s=Remberg(fun,a,b,epsilon)
' U+ U8 p; w1 z/ ^$ q5 ]syms x ;4 v9 p/ F- _6 j
fvalue=zeros(1,1000);9 o. t* r( F! }
R=zeros(100,100);0 h' y6 u1 X4 M& F4 c6 y8 y( m; I9 |
fvaluea=double(subs(fun,x,a));; Z" U0 Q- v+ R. }
fvalueb=double(subs(fun,x,b));- c0 ?- Y) G! ~3 R. ?) A
R(1,1)=(b-a)/2*(fvaluea+fvalueb);   %梯形公式# a" p4 U, x2 U! k0 g
km=1;
& _/ y" M. P- a; h5 ofor k1=1:100;   %設置一個比較大的循環量- L0 e0 L" j' D6 ?. F, p5 i% v6 x
    h=(b-a)/(2^(k1));9 M* Q, B9 X$ c' v+ v. \3 f  p
    R(k1+1,1)=1/2*R(k1,1);
+ R' c; Z. K/ b4 b- O1 ^: j) N8 Y    for k2=1:2^(k1-1);
/ z" ?; ?; i8 M, O4 d    fvalue(2*k2)=double(subs(fun,x,a+(2*k2-1)*h));
# K( G! N) [" Z* P    R(k1+1,1)=h*fvalue(2*k2)+R(k1+1,1);  %變步長值
. j# V( K( L( s1 p" }    end
  i/ Z0 N2 d4 p( w  z% c: G  for k3=1:km;  %加速計算
3 u4 U+ h# V4 |. J, h      R(k1+1,k3+1)=1/(4^(k3)-1)*(4^(k3)*R(k1+1,k3)-R(k1,k3));
* t3 c" N7 A2 g) J; Y/ f. H  end
# s4 D6 S. V. s; M! o. V* W  if abs(R(k1+1,km+1)-R(k1+1,km))<epsilon   %控制精度3 {4 a2 i. o) [+ f1 ^
      s=R(k1+1,km+1);4 U$ H6 U# ~& `
      break;
$ d& D! O4 V2 d% K# ~4 ?9 P  else
7 r; P. c' j2 |/ v' X5 w      km=km+1;
% I3 Y' y8 A7 @, G! Y# g  end4 u2 W: V$ z" K+ X" w, Y9 [

, x' `2 J! C$ i- X& p+ }end
% ~9 j( }+ A. c
. g- K- B$ _; h
您需要登錄后才可以回帖 登錄 | 注冊會員

本版積分規則

Archiver|手機版|小黑屋|機械社區 ( 京ICP備10217105號-1,京ICP證050210號,浙公網安備33038202004372號 )

GMT+8, 2025-7-30 18:48 , Processed in 0.086515 second(s), 16 queries , Gzip On.

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回復 返回頂部 返回列表