| 查看: 2133 | 回复: 7 | |||
[求助]
大神！求助Guass-Hermite积分二重积分、三重积分Matlab源程序！！ 已有1人参与
|
|
大神！求助Guass-Hermite积分二重积分、三重积分源程序！！ 很多教材书，或者工具书上都有关于Guass-Hermite积分一重积分的程序，但是lz遇到解决二重积分的Guass-Hermite积分，怎么设计Matlab程序大神！求帮助！！ 附件是关于一重的Guass-Hermite积分源程序。 |
回复此楼» 收录本帖的淘帖专辑推荐
 程序 |
» 猜你喜欢
河北省自然科学基金
已经有8人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有5人回复
-
论文撤稿了
已经有5人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
-
26/27申博自荐
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
揭秘青基评审内幕：几个A才能顺利中标
已经有4人回复
-
青B发送上会通知了吗
已经有7人回复
-
博士申请
已经有3人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
2楼2016-08-02 06:26:04
3楼2016-08-02 10:22:09
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
【答案】应助回帖
|
% e3s408.m disp( ' n integral value' ); for j = 1:4 n = 2^j; int = gaherm( @( x ) 1./( 1+x.^2 ).^2,n ); fprintf( '%3.0f%14.9f\n',n,int ) end % gaherm.m function s = gaherm( func,n ) % Implements Gauss-Hermite integration. % Example call: s = gaherm( func,n ) % Integrates user defined function func from -inf to +inf, % using n divisions. n must be 2 or 4 or 8 or 16 if ( n==2 )|( n==4 )|( n==8 )|( n==16 ) c = zeros( 8,4 ); t = zeros( 8,4 ); c( 1,1 ) = 1.461141183; c( 1:2,2 ) = [1.059964483; 1.240225818]; c( 1:4,3 ) = [.7645441286; .7928900483; .8667526065; 1.071930144]; c( :,4 ) = [.5473752050; .5524419573; .5632178291; .5812472754; ... .6097369583; .6557556729; .7382456223; .9368744929]; t( 1,1 ) = .7071067811; t( 1:2,2 ) = [.5246476233; 1.650680124]; t( 1:4,3 ) = [.3811869902; 1.157193712; 1.981656757; 2.930637420]; t( :,4 ) = [.2734810461; .8229514491; 1.380258539; 1.951787991; ... 2.546202158; 3.176999162; 3.869447905; 4.688738939]; j = 1; while j<=4 if 2^j==n; break; else j = j+1; end end s=0; for k = 1:n/2 x1 = t( k,j ); x2 = -x1; y = feval( func,x1 )+feval( func,x2 ); s = s+c( k,j )*y; end else disp( 'n must be equal to 2, 4, 8 or 16' ); return end 積分: int [ d x / ( 1 + x^2 ) ^2 ] 結果: n integral value 2 1.298792163 4 1.482336098 8 1.550273058 16 1.565939612 参考: http://www.sciencedirect.com/sci ... 9780123869425000047 |
4楼2016-08-02 18:04:39
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
愚木头5楼2016-08-02 18:05:51
6楼2016-08-03 22:30:45
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
7楼2016-08-06 20:41:36
8楼2016-08-15 12:29:29
10