24小时热门版块排行榜    

查看: 1999  |  回复: 6

puchuanxia

铜虫 (小有名气)

[求助] MATLAB无法计算出这个含【级数和超几何函数】的式子

如题,我试图让m的值变大以更接近真值,但超过37就算不出来。
请问各位高手遇到这种不报错又算不出的情况怎样处理?
(注:式中1F2为超几何函数,T(m+n)为伽马函数)
代码如下:
clc
clear all
close all
r=3.;
lam=5.32e-4;
k=2*pi/lam;
n=25;
f=500;
alpha=9.9;
rho=linspace(0.,1,100);
len=length(rho);
y=zeros(1,len);f1=zeros(1,len);f2=zeros(1,len);
tic;
h=waitbar(0,'calculation');
for ii=1:len
     waitbar(ii/100);
     f1(ii)=k/alpha/f/factorial(n)*(k*alpha*rho(ii)*r^2/4/f)^n;
        for m=0:30
            f2(ii)=f2(ii)+(-alpha^2*r^2/4)^m/factorial(m)/gamma(m+n)/(m+n)*hypergeom(m+n,[m+n+1,n+1],-(k*r*rho(ii)/2/f)^2);
        end
     y(ii)=(f1(ii)*f2(ii))^2;
end
close(h)
figure(1)
plot(rho,y);
toc;

MATLAB无法计算出这个含【级数和超几何函数】的式子
要计算的式子
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

徐谊

木虫 (著名写手)

帮顶一下
2楼2014-04-20 11:13:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

puchuanxia

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by 徐谊 at 2014-04-20 11:13:59
帮顶一下

谢谢捧场 朋友
这个问题你怎么看
3楼2014-04-20 14:09:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

徐谊

木虫 (著名写手)

引用回帖:
3楼: Originally posted by puchuanxia at 2014-04-20 14:09:03
谢谢捧场 朋友
这个问题你怎么看...

我表示不懂
4楼2014-04-20 20:01:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

puchuanxia

铜虫 (小有名气)

引用回帖:
4楼: Originally posted by 徐谊 at 2014-04-20 20:01:03
我表示不懂...

我自己也不明白 被这问题搞的没有周末
5楼2014-04-21 10:30:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

徐谊

木虫 (著名写手)

引用回帖:
5楼: Originally posted by puchuanxia at 2014-04-21 10:30:11
我自己也不明白 被这问题搞的没有周末...

可以去matlab论文问问
6楼2014-04-21 13:24:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

puchuanxia

铜虫 (小有名气)

我后来这样解决了问题:
1.当n=30时,所得结果已经收敛为真值,没必要再去增大n值;
2.matlab算不出的原因是无法计算出以下一类超几何函数
         hypergeom(86,[87,31],-400)
  而改用其它软件,如mathematica则可以算出;
3.我的程序仍不够精简,还望大家提出修改意见!

希望以上对大家有所帮助!
7楼2014-05-04 09:35:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 puchuanxia 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +7 大米饭! 2026-03-15 7/350 2026-03-16 10:25 by 了了了了。。
[考研] 调剂 +8 调剂的考研学生 2026-03-09 8/400 2026-03-15 22:14 by Winj1e
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 301求调剂 +3 归零lbm 2026-03-09 3/150 2026-03-14 02:20 by JourneyLucky
[考研] 一志愿安徽大学材料工程专硕313分,求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 数二英二309分请求调剂 +3 dtdxzxx 2026-03-09 4/200 2026-03-09 19:56 by yuningshan
信息提示
请填处理意见