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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数理科学综合 » 大神！求助Guass-Hermite积分二重积分、三重积分Matlab源程序！！
81/1返回列表
查看: 2078  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京化工大学070300 学硕336求调剂 +4 vv迷 2026-03-21 5/250 2026-03-22 08:38 by ChemPharm
[考研] 生物学071000 329分求调剂 +4 我爱生物生物爱� 2026-03-17 4/200 2026-03-22 08:34 by hxsm
[考研] 化学工程321分求调剂 +18 大米饭！ 2026-03-15 22/1100 2026-03-21 20:20 by HH领袖
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 一志愿南大，0703化学，分数336，求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 一志愿中国石油大学（华东） 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 广西大学家禽遗传育种课题组2026年硕士招生（接收计算机专业调剂） +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 085410人工智能专硕317求调剂（0854都可以） +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
信息提示
关闭
请填处理意见
关闭