| 查看: 1661 | 回复: 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非常外行，请大贤指点，感激涕零！ |
回复此楼» 猜你喜欢
277求调剂 数一104分
已经有12人回复
-
复试调剂
已经有14人回复
-
307求调剂
已经有8人回复
-
材料工程302分求调剂
已经有10人回复
-
材料考研求调剂总分280
已经有14人回复
-
工科 22408 267求推荐
已经有5人回复
-
环境专硕调剂
已经有5人回复
-
287分求调剂 有专利国奖一志愿哈工大085406
已经有7人回复
-
080500求调剂
已经有9人回复
-
求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一个绝对值求积分问题
已经有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
- 专业: 流体力学
2楼2013-08-19 14:44:48
3楼2013-08-19 14:53:48
tgvictor99
木虫 (知名作家)
迷糊领域荣誉领主
- 应助: 40 (小学生)
- 金币: 2667.4
- 散金: 1907
- 红花: 49
- 沙发: 1
- 帖子: 5238
- 在线: 923.8小时
- 虫号: 2524003
- 注册: 2013-06-27
- 专业: 流体力学
4楼2013-08-19 15:25:03
abacus
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 15994.5
- 散金: 100
- 红花: 11
- 帖子: 2273
- 在线: 639.9小时
- 虫号: 553166
- 注册: 2008-05-02
- 专业: 建筑环境与结构工程
5楼2013-08-19 16:20:36
abacus
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 15994.5
- 散金: 100
- 红花: 11
- 帖子: 2273
- 在线: 639.9小时
- 虫号: 553166
- 注册: 2008-05-02
- 专业: 建筑环境与结构工程
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
li_qun: 金币+20, ★★★★★最佳答案, 谢谢，受益匪浅 2013-08-19 18:01:18
li_qun: 金币+10 2013-08-23 08:49:04
xiegangmai: 金币+2, 谢谢参与！ 2013-08-24 21:53:23
li_qun: 金币+20, ★★★★★最佳答案, 谢谢，受益匪浅 2013-08-19 18:01:18
li_qun: 金币+10 2013-08-23 08:49:04
xiegangmai: 金币+2, 谢谢参与！ 2013-08-24 21:53:23
|
看来还不太会用小木虫的 插入代码。 再发一次吧。结果说明，你的第一次计算是正确的。 %% Method 1 f=@(z1,z2)exp(abs(z1-z2));%对该函数进行积分，积分范围0<z1<3,1<z2<3 an=quad2d(f,0,3,1,3) %% Method 2 % 积分区域分为 两部分 bn=quad2d(@(z1,z2)exp(z1-z2),1,3,1,@(z1)z1); % z1>z2 cn=quad2d(@(z2,z1)exp(z2-z1),1,3,0,@(z2)z2); % z2<z1 dn=bn+cn %最终结果 |
6楼2013-08-19 16:25:03
