| 查看: 1779 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助被积函数有绝对值号的问题。
|
||
|
f=@(z1,z2)exp(abs(z1-z2));%对该函数进行积分，积分范围0<z1<3,1<z2<3 quad2d(@(z1,z2)f(z1,z2),0,3,1,3) ans = 19.7563 >> %去掉绝对值号计算，对f函数的积分转化成对f1，f2两个函数的积分 f1=@(z1,z2)exp(z1-z2);%积分范围0<z1<3,1<z2<z1 f2=@(z2,z1)exp(z2-z1);%积分范围1<z2<3,0<z1<z1 c=@(z1)z1; d=@(z2)z2; quad2d(@(z1,z2)f1(z1,z2),0,3,1,c)+quad2d(@(z2,z1)f2(z2,z1),1,3,0,d) ans = 19.3884 末学已经确定第二种方法是错误的，因为对f1的积分quad2d(@(z1,z2)f1(z1,z2),0,3,1,c)来说，并不表示1<z2<c，c有可能小于1 但不知道第一种方法是否正确？ 如果要去掉绝对值号，应该怎么改写这个函数？ ================================================== 还有一个问题，对函数f4积分 f4=@(z1,z2)exp(z1-z2); quad2d(@(z1,z2)f4(z1,z2),2,5,3,8) ans = 6.9739 但我用分段求和的方法来积分这个函数，结果如下 sum1=0; N=1000; %分段数 d1=(5-2)/N; %变量z1的步长 d2=(8-3)/N; %变量z2的步长 for i=1:N sum2=0; z1=(2+(i-1)*d1+2+i*d1)/2; %先固定z1，对z2积分，sum2表示对z2的积分结果，sum1表示对z1的积分结果 for j=1:1000 z2=(3+(i-1)*d2+3+i*d2)/2; sum2=sum2+f4(z1,z2)*d2; end sum1=sum1+sum2*d1; end >> sum1 sum1 = 2.3857 这是什么原因，哪里有错误吗，还是精度产生的误差？末学对matlab非常外行，请大贤指点，感激涕零！ |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有12人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有29人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一个绝对值求积分问题
已经有8人回复
- 求助 matlab隐函数画图问题
已经有4人回复
- 求助一个inline函数的问题
已经有5人回复
- 高数：导函数与原函数问题求助
已经有12人回复
- 请教一个积分用留数定理计算的问题
已经有15人回复
- 求助 Origin 自定义函数输入问题
已经有4人回复
- 关于势函数的求助
已经有23人回复
- 求助excel2007中normdist函数的使用问题
已经有1人回复
- 【求助】linux信号处理函数中，一定不能用malloc这类不可重用函数吗？
已经有4人回复
- 【求助】delta函数平方后积分
已经有5人回复
- 【求助】一个求数组绝对值的小问题
已经有7人回复
- 【求助】用graythresh函数遇到的问题
已经有8人回复
- 【求助】相关函数归一化问题
已经有6人回复
- 【求助】请教mathematica中关于Plot函数的问题【已完成】
已经有5人回复
- 【求助】求积分程序中被积函数问题
已经有15人回复
- 【求助】Mittag-leffler函数问题
已经有3人回复
- 【求助】一个求解反函数的问题
已经有9人回复
- 【求助】只知道有限个函数值求积分
已经有10人回复
- 【求助】请教非平衡格林函数的问题
已经有5人回复
tgvictor99
木虫 (知名作家)
迷糊领域荣誉领主
- 应助: 40 (小学生)
- 金币: 2667.4
- 散金: 1907
- 红花: 49
- 沙发: 1
- 帖子: 5238
- 在线: 923.8小时
- 虫号: 2524003
- 注册: 2013-06-27
- 专业: 流体力学
4楼2013-08-19 15:25:03
tgvictor99
木虫 (知名作家)
迷糊领域荣誉领主
- 应助: 40 (小学生)
- 金币: 2667.4
- 散金: 1907
- 红花: 49
- 沙发: 1
- 帖子: 5238
- 在线: 923.8小时
- 虫号: 2524003
- 注册: 2013-06-27
- 专业: 流体力学
2楼2013-08-19 14:44:48
3楼2013-08-19 14:53:48
abacus
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 15994.5
- 散金: 100
- 红花: 11
- 帖子: 2273
- 在线: 639.9小时
- 虫号: 553166
- 注册: 2008-05-02
- 专业: 建筑环境与结构工程
5楼2013-08-19 16:20:36
100