24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2025  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿郑州大学085600求调剂 +9 吃的不少 2026-04-05 10/500 2026-04-05 21:02 by 学员8dgXkO
[考研] 求调剂,一志愿厦门大学,生物与医药,总分272,本科211 +3 Electron1cc 2026-04-01 4/200 2026-04-05 20:24 by lys0704
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 288求调剂 +7 没有答案_ 2026-04-05 7/350 2026-04-05 17:01 by yulian1987
[考研] 一志愿南航,数一英一学硕317求调剂!! +5 Acaciad 2026-04-04 5/250 2026-04-05 12:31 by 搏击518
[考研] +5 化工专硕323分 2026-04-04 5/250 2026-04-05 08:02 by 544594351
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[考研] 求调剂 +4 晟功? 2026-04-03 4/200 2026-04-04 21:58 by hemengdong
[考研] 调剂 +4 是可乐不是可乐 2026-04-04 4/200 2026-04-04 19:41 by 唐沐儿
[考研] 一志愿北交大材料工程总分358 +6 cs0106 2026-04-03 6/300 2026-04-04 11:20 by w_xuqing
[考研] 调剂 +5 asdasdassda 2026-04-03 6/300 2026-04-03 20:27 by 岸上的一条鱼
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 085600专硕材料与化工348分求调剂 +10 上学啦! 2026-04-01 11/550 2026-04-03 14:13 by 百灵童888
[考研] 338求调剂,一志愿能源动力,外语是日语203 +5 zzz,,r 2026-04-02 5/250 2026-04-03 09:45 by 蓝云思雨
[考研] 学硕化学工程与技术,一志愿中国海洋大学320+求调剂 +8 披星河 2026-04-02 8/400 2026-04-02 14:12 by oooqiao
[考研] 07生物学求调剂 一志愿同济大学359分 +3 LAMC. 2026-03-30 3/150 2026-04-02 10:26 by 18828373951
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 11:35 by wangjy2002
[考研] 340求调剂 +4 希望如此i 2026-03-31 4/200 2026-03-31 16:40 by 690616278
[考研] 一志愿中海洋材料357 +4 麦恩莉. 2026-03-30 4/200 2026-03-31 14:35 by 记事本2026
[考研] 一志愿食品科学与工程083200求调剂 +4 XQTJZ 2026-03-30 4/200 2026-03-31 04:10 by fmesaito
信息提示
请填处理意见