488论坛币－matlab求解一积分方程

我爱香港

2927

收藏 2010-05-05

悬赏 488 个论坛币已解决

如下积分方程：
m := ; fsolve(int((1-x^2)^(565), x = y .. 1) = a*(int(int((1-x^2)^(565), x = -1 .. 1), y)));，
即【(1－x的平方)的565次方】的积分（上限为1，下限为y），=a乘以【(1－x的平方)的565次方】的积分（上限为1，下限为－1）.

a分别取0.005，0.01，0.05，0.10，0.25 时，y的值。即给出一个表，给出a相应取值对应的y值。
跪求各位兄弟姐妹，并把程序也一起给我，小弟实在不懂。

新建 Microsoft Word 文档.doc

大小:16.5 KB

马上下载

最佳答案

joyzhang1988 查看完整内容

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

joyzhang1988

2010-5-5 01:15:43

> a := 0.5e-2; m := 100; fsolve(int((1-x^2)^(.5*(m-3)), x = y .. 1) = a*(int(int((1-x^2)^(.5*(m-3)), x = -1 .. 1), y)));
print(`output redirected...`); # input placeholder
                                 0.005
                                    100
                              0.2931267031
> a := 0.1e-1; m := 100; fsolve(int((1-x^2)^(.5*(m-3)), x = y .. 1) = a*(int(int((1-x^2)^(.5*(m-3)), x = -1 .. 1), y)));
                              0.2744023417
> ff(0.5e-2, 100);
                              0.2931267031
> ff(0.1e-1, 100);
                              0.2744023417
> ff(0.5e-1, 100);
                              0.2266363731
> ff(.1, 100);
                              0.2039454352
> ff(.25, 100);
                              0.1717657563
> ff(0.5e-2, 200); ff(0.1e-1, 200); ff(0.5e-1, 200); ff(.1, 200);
%;
                              0.2152266327
                              0.2019586410
                              0.1684301170
                              0.1526358821
> ff(.25, 200);
                              0.1303316084
> ff(0.5e-2, 400); ff(0.1e-1, 400); ff(0.5e-1, 400);
%;
                              0.1573485581
                              0.1480464397
                              0.1246772233
> ff(.1, 400);
                              0.1137277576
> ff(.25, 400);
                              0.09830861929
> ff(0.5e-2, 600);
                              0.1308176498
>  ff(0.1e-1, 600); ff(0.5e-1, 600); ff(.1, 600); ff(.25, 600);
                              0.1232807741
                              0.1043940601
                              0.09556594884
                              0.08315054288
> ff(0.5e-2, 800);
                              0.1146948682

> ff(0.1e-1, 800); ff(0.5e-1, 800); ff(.1, 800); ff(.25, 800);
                              0.1082081489
                              0.09197798392
                              0.08440284852
                              0.07375867643