| 查看: 550 | 回复: 3 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】如何用matlab计算无穷级数，给出一个收敛的值【已完成】
|
|||
|
小弟解出一个带有bessel函数的解，想用matlab计算出来，目标画出图形 可以为啥matlab总是给出一个sym type的解，而给出的不是一个具体的收敛的数值呢？ 代码如下： %————————————我是分割线———————————————— clear all; clc; z=2.5; b=3.37; zb=2.88; zt=2.28;rw=0.115;x=1; syms n; s=(cos(n*pi*zb/(2*b))-cos(n*pi*zt/(2*b)))/(n*pi/(2*b)*n*pi/(2*b)*(1/10)^(1/2)*besseli(1,n*pi/(2*b)*(1/10)^(1/2)*rw))*besseli(0,n*pi/(2*b)*(1/10)^(1/2).*x)*sin(n*pi/(2*b)*z); S=symsum(s,1,inf); %————————————我还是分割线————————————————— %matlab给出的解 S = sum(113569/2500*(cos(144/337*n*pi)-cos(114/337*n*pi))/n^2/pi^2*10^(1/2)/besseli(1,23/13480*n*pi*10^(1/2))*besseli(0,5/337*n*pi*10^(1/2))*sin(125/337*n*pi),n = 1 .. Inf) %————————————我又是讨厌的分割线—————————————— 我的问题是如何让S变成一个收敛的数呢？ 谢谢大家！祝大家天天有个好心情 [ Last edited by nono2009 on 2009-9-23 at 21:34 ] |
回复此楼» 猜你喜欢
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有30人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
3楼2009-07-14 10:54:57
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
2楼2009-07-14 10:47:05
4楼2009-07-15 07:12:32
