声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3618|回复: 7

[编程技巧] matlab中二重高斯积分的函数quad2dg的问题

[复制链接]
发表于 2009-9-6 16:20 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
为计算有限元,从网络上找到一个计算二重矩形高斯积分的积分函数quad2dg,其使用方法举例如下:
>>quad2dg('(x.^2.*y)',0,1,0,1) 即可以得到其结果    0.1667
我的问题是我的表达式是用一个变量表示的如:u=x.^2.*y
现在我该如何对u进行积分呢?
我试过各种形式如:quad2dg('u',0,1,0,1) 、quad2dg(u,0,1,0,1) 、quad2dg('(u)',0,1,0,1) 、quad2dg('inline(u)',0,1,0,1) 、quad2dg(inline(u),0,1,0,1) 等都提示有错误。我把quad2dg的代码上传到附件中。
请高手帮帮我。我该如何写这个积分。

quad2dg.m

7.04 KB, 下载次数: 24

二重矩形高斯积分的积分函数quad2dg

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2009-9-6 19:20 | 显示全部楼层
试了下LZ所附之程序, 结论如下
LZ取得的程序, fun仅支援两种方式, 一为直接给定如quad2dg('(x.^2.*y)',0,1,0,1), 另一为fun函数如quad2dg(u,0,1,0,1), 但需另定一u function!
并无支援其他方式的输入, 除非修改下程序, 如在74行前加些判断及转换! 自个儿试吧
发表于 2012-4-4 07:56 | 显示全部楼层
这个是计算二重积分的,高斯积分点的选择在哪里
发表于 2012-10-31 19:19 | 显示全部楼层
发表于 2012-10-31 19:20 | 显示全部楼层
我遇到了相同的问题,我是要写成M函数,然后在去积分...结果也是报错...
发表于 2012-10-31 22:05 | 显示全部楼层

LS不是搜索到了吗? 原来以前也试过并回应过!
真是老了记性差, 还重看了help!
2F不是给了引导了!?
发表于 2012-11-1 15:11 | 显示全部楼层
function fun=test_sy(x,y)

fun = 1./( sqrt(x + y) .* (1 + x + y).^2 );

>> quad2d(@test_sy,0,1,2,9)


给出个例子,造福大家...

ans =

    0.0937

评分

1

查看全部评分

发表于 2012-11-6 16:29 | 显示全部楼层
多谢分享
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-13 13:16 , Processed in 0.081957 second(s), 25 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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