24小时热门版块排行榜    

查看: 2147  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

245153618

新虫 (初入文坛)

[求助] 格林函数 积分 数值计算

x1=0;y1=1;Dmn00=integral2(@(x2,y2)1./sqrt((x1-x2).^2+(y1-y2).^2),0,1,0,1)

Dmn00 =

   1.762747210530955
为什么
D=integral2(@(x1,y1)arrayfun(@(x1,y1)integral2(@(x2,y2)1./sqrt((x1-x2).^2+(y1-y2).^2),0,1,0,1),x1,y1),0,1,0,1),
计算的不了。
还有就是关于x1,x2,y1,y2的四重积分1./sqrt((x1-x2).^2+(y1-y2).^2),积分范围都是0到1,有有限的积分值吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

245153618

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by 田山东 at 2012-06-19 19:17:35
看到公式就让人崩溃了,能不能弄得直观点?或者弄个小程序贴上来,有问题或许有人还能帮你看下呢

就是matlab程序啊
5楼2012-06-19 20:06:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

田山东

捐助贵宾 (著名写手)

看到公式就让人崩溃了,能不能弄得直观点?或者弄个小程序贴上来,有问题或许有人还能帮你看下呢
everythinghasitsseason.enjoyyourlife.
4楼2012-06-19 19:17:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lingyuandage

木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
楼主仔细看看是不是有分母为0的情况出现
感觉好像是因为x2=0,y2=1时分母为0
不知道对不对
TheincurablePisces.
6楼2012-06-20 18:02:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lingyuandage

木虫 (职业作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by 245153618 at 2012-06-20 22:58:17
我就是想知道,这个奇异积分是不是可积的,就是存在分母为0的情况...

我的理解:用程序算积分其实就是加法计算
MATLAB某种程度上是不能计算广义积分的
应为在特定点的值无法给出
建议将积分上下限做微调,既0变成a,1变成1-a
计算完成后再将a=0代入
我电脑上现在木有MATLAB,你可以试试
TheincurablePisces.
8楼2012-06-20 23:09:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 10/500 2026-09-06 20:08 by 阿美_Lml888
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 19:21 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 6/300 2026-09-06 18:41 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 17:49 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
信息提示
请填处理意见