应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数理科学综合 » 大神！求助Guass-Hermite积分二重积分、三重积分Matlab源程序！！
81/1返回列表
查看: 2077  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

愚木头

金虫 (正式写手)

[求助] 大神！求助Guass-Hermite积分二重积分、三重积分Matlab源程序！！ 已有1人参与

大神！求助Guass-Hermite积分二重积分、三重积分源程序！！
很多教材书，或者工具书上都有关于Guass-Hermite积分一重积分的程序，但是lz遇到解决二重积分的Guass-Hermite积分，怎么设计Matlab程序大神！求帮助！！
附件是关于一重的Guass-Hermite积分源程序。
回复此楼

» 收录本帖的淘帖专辑推荐

程序

» 猜你喜欢
1楼 2016-08-01 21:54:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
Gram-Schmidt of matlab
http://web.mit.edu/18.06/www/Essays/gramschmidtmat.pdf
赞一下 回复此楼
2楼2016-08-02 06:26:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

愚木头

金虫 (正式写手)
引用回帖:
2楼: Originally posted by FMStation at 2016-08-02 06:26:04
Gram-Schmidt of matlab
http://web.mit.edu/18.06/www/Essays/gramschmidtmat.pdf

表示没关系吧，不过谢谢啦！！
赞一下 回复此楼
3楼2016-08-02 10:22:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

% 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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)
引用回帖:
2楼: Originally posted by FMStation at 2016-08-02 06:26:04
Gram-Schmidt of matlab
http://web.mit.edu/18.06/www/Essays/gramschmidtmat.pdf

A=[4 -2; 3 1];
n=2;

for j = 1:n
    v = A( :, j );
    for i = 1:j-1
        R(i,j) = Q( :,i )' * A( :,j );
        v = v - R( i,j ) * Q( :,i );
    end
    R( j,j ) = norm( v );
    Q( :,j ) = v / R( j,j );
end

A
Q
R
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

愚木头
5楼2016-08-02 18:05:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

愚木头

金虫 (正式写手)
送红花一朵
引用回帖:
5楼: Originally posted by FMStation at 2016-08-02 18:05:51
A=;
n=2;

for j = 1:n
    v = A( :, j );
    for i = 1:j-1
        R(i,j) = Q( :,i )' * A( :,j );
        v = v - R( i,j ) * Q( :,i );
    end
    R( j,j ) = norm( v );
    Q( :,j ) = v / ...

谢谢！我正在试，还没看懂！！
赞一下 回复此楼
6楼2016-08-03 22:30:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

MIT OCW - 18.06  Linear Algebra
https://stellar.mit.edu/S/course/18/sp16/18.06/
赞一下 回复此楼
7楼2016-08-06 20:41:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yibaka
赞

发自小木虫Android客户端
8楼2016-08-15 12:29:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 愚木头 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 333求调剂 +5 87639 2026-03-21 7/350 2026-03-21 19:31 by ColorlessPI
[考研] 326求调剂 +4 mlpqaz03 2026-03-15 4/200 2026-03-21 19:10 by ColorlessPI
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 求调剂 +5 十三加油 2026-03-21 5/250 2026-03-21 18:48 by 学员8dgXkO
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +6 嘉年新程 2026-03-15 6/300 2026-03-21 17:07 by Dream007008
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 08工科 320总分 求调剂 +6 梨花珞晚风 2026-03-17 6/300 2026-03-21 03:40 by JourneyLucky
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 3/150 2026-03-20 11:58 by 呱呱呱呱叫
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 334求调剂 +3 志存高远意在机� 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
关闭
请填处理意见
关闭