我想在matlab中求积分函数这个函数的积分

来源：蜘蛛抓取(WebSpider) 时间：2016-08-21 11:52 标签： matlab求定积分函数
君，已阅读到文档的结尾了呢~~
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
如何用matlab对这个函数进行积分。
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口关于MATLAB求定积分如果我有了一个inline函数f = inline('some_expression'),如何求x*f(x)从a 到b上的定积分(使用纯数值办法,不要用syms之类的符号运算)
>>f=inline('sin(x)')>>g=@(x)x.*f(x)>>quad(g,1,2)
为您推荐：
其他类似问题
扫描下载二维码(Phantom_Ghost)
(第一推动)
(Phantom_Ghost)
(TITANIC宫野)
第三方登录：用matlab计算复变函数的积分_图文_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
用matlab计算复变函数的积分
上传于||暂无简介
大小：2.94MB
登录百度文库，专享文档复制特权，财富值每天免费拿！
你可能喜欢查看: 3252|回复: 16|关注: 0
quad求分段函数积分问题
<h1 style="color:# 麦片财富积分
新手, 积分 10, 距离下一级还需 40 积分
使用quad函数进行分段函数积分，一直报错，贴出来，希望大家帮我解决下。
这是我的分段函数描述代码，放在了MyCodeApertureRotate.m文件中：
function yout=MyCodeApertureRotate(x)
if (x&=0&x&1)
& & yout=10;
elseif(x&=1&x&2)
& & yout=20;
elseif(x&=2)
& & yout=30;
end
复制代码这是我的主函数，只有一句话：
quad(@(x)MyCodeApertureRotate(x),0,1)
复制代码错误为：
Output argument &yout& (and maybe others) not assigned during call to
&E:\Spectral_imaging\MyPlatform_2\MyCodeApertureRotate.m&MyCodeApertureRotate&.
Error in @(x)MyCodeApertureRotate(x)
Error in quad (line 76)
y = f(x, varargin{:});
Error in MyIntCode (line 4)
quad(@(x)MyCodeApertureRotate(x),0,1)
复制代码麻烦各位帮忙解决下。十分感谢。
或者哪位知道如何求分段函数的积分问题，帮我出个主意。谢谢了！
帖子最佳答案
关注者： 365
你犯了2个错误：
1. quad函数调用时，@函数名 就是了，因为 @函数名 构成一个函数句柄，你不能用 @函数名(x)。
2. 子函数里的分段函数定义有误，你用if else结构只能定义一个标量函数，而quad函数要求函数是支持向量输入的。其实分段函数最好的办法是用 “分段加权相加” 的办法定义，如下function yout=MyCodeApertureRotate(x)
yout = (x&=0&x&1)*10 + (x&=1&x&2)*20 + (x&=2)*30;复制代码quad调用这一句改成：quad(@MyCodeApertureRotate,0,1)复制代码最后，给出另一种解决办法：匿名函数法，这个最方便f = @(x) (x&=0&x&1)*10 + (x&=1&x&2)*20 + (x&=2)*30;
quadgk(f,0,1)复制代码最后补充一下，分段函数积分最好用quadgk，因为它能更好的处理不连续点
<h1 style="color:# 麦片财富积分
winner245 发表于
你犯了2个错误：
1. quad函数调用时，@函数名 就是了，因为 @函数名 构成一个函数句柄，你不能用 @函数名( ...
十分感谢，第一个错误很低级，第二个确实是不了解，学习到了！另外quadgk很经典。改正你点出的错误后，我用quad算的时候结果一直不对，可能是quad对于分段函数求积分还是不太好用。具体原因您清楚吗？这个我也想学习下。谢了！
帖子最佳答案
关注者： 365
静夫止甫 发表于
十分感谢，第一个错误很低级，第二个确实是不了解，学习到了！另外quadgk很经典。改正你点出的错误后，我 ...
首先你的分段函数是不连续的，这个很明显。其次，即使你的分段函数在数学定义上是连续的，由于matlab的积分过程中实际上是按一定间隔计算函数值，所以，当取到分段点处时，实际上不连续的（尽管数学定义是连续的），而quad函数未必能很好的处理此类间断点的情形。相反，quadgk可以很好地处理有限个间断点，甚至可以有效处理端点处的奇异值。所以，quadgk可以规避这个问题
<h1 style="color:# 麦片财富积分
winner245 发表于
首先你的分段函数是不连续的，这个很明显。其次，即使你的分段函数在数学定义上是连续的，由于matlab的积 ...
Thank you very much
<h1 style="color:# 麦片财富积分
winner245 发表于
首先你的分段函数是不连续的，这个很明显。其次，即使你的分段函数在数学定义上是连续的，由于matlab的积 ...
那二维函数积分的话， z=f(x,y),如果是关于x，y都是分段函数，例如&&当(x&1&x&2&y&2&y&4)时，z=10，当(x&4&x&6&y&1&y&2)时，z=20，等等还有多个其他分段。对这样的函数如何进行二维积分呢？用dblquad吗？我看qblquad这个函数也是调用的quad函数，这样跟你说的一样就处理不了间断点了啊。
或者还有其他函数吗？请赐教！谢了。
帖子最佳答案
关注者： 365
本帖最后由 winner245 于
10:17 编辑
静夫止甫 发表于
那二维函数积分的话， z=f(x,y),如果是关于x，y都是分段函数，例如&&当(x&1&x2&y4&x1&y ...
求二重积分除了 dblquad，还有 quad2d 和 integral2两个函数，后面两个函数是 matlab 高版本里新增的函数，所以，他们应该是比旧函数 dblquad更强大，尤其是 integral2，你看看help的描述，里面提到这个函数可以处理边界处的奇异点，这一点类似于quadgk，所以，integral2应该是一个很强大的函数
当然，如果对于二维分段函数直接用上述3个函数积分都不是很精确的话，你可以考虑对每个分段分别积分，然后求和，这样，就不存在分段的问题了。比如，你的分段函数有3个分支，你就调用二重积分函数3次，每次调用对应一个分支，这样，只要每个分支能精确的计算就足够了，而这个要求实际上是很低的，因为对于没有奇异点或不连续点的每个分支，上述三个函数几乎都能很好地处理
<h1 style="color:# 麦片财富积分
winner245 发表于
求二重积分除了 dblquad，还有 quad2d 和 integral2两个函数，后面两个函数是 matlab 高版本里新增的函数 ...
非常感谢！已经用过quad2d了，现在学习下intergral。谢谢
帖子最佳答案
关注者： 365
静夫止甫 发表于
非常感谢！已经用过quad2d了，现在学习下intergral。谢谢
不客气，以后互相学习
<h1 style="color:# 麦片财富积分
winner245 发表于
求二重积分除了 dblquad，还有 quad2d 和 integral2两个函数，后面两个函数是 matlab 高版本里新增的函数 ...
发现quad2d在算二重积分时，如果碰到间断点，速度会非常慢。而我的计算还没想到如何在间断点分成多个分支来单独计算。
我所需要的是，之前已经定义了二维分段函数，“ 当(x&1&x&2&y&2&y&4)时，z=10，当(x&4&x&6&y&1&y&2)时，z=20，等等还有多个其他分段”，其间断点或者说间断线全部出现在x，y取整数的时候，如果对此函数每段分别积分，然后求和，肯定是行得通的。但是现在要将此函数在xy平面旋转一个小角度a，重新得到另外的一个函数g(x,y)。g(x,y)和z(x,y)形成一个旋转的关系。这样对g(x,y)每段分别积分就出问题了，g(x,y)的间断线出现在斜线的位置，如何将其分段是很复杂的。而我想要的是对g(x,y)积分，其积分范围仍然是x，y取相邻整数的时候，也类似“1&x&2,3&y&4”这样的区间。如果直接积分的话必然经过间断线，速度超级慢，请问您有什么好办法吗？
如果我有没表达清楚的地方，请告诉我。谢谢了
站长推荐 /2
Powered by

我要回帖

更多关于 matlab求定积分函数 的文章

 

随机推荐