24小时热门版块排行榜    

查看: 1019  |  回复: 2

huazhengli

铁虫 (初入文坛)

[求助] 急求助求解一个含贝塞尔函数方程 已有1人参与

各位大牛,
       如题,要计算一个气体在中空的圆柱内部扩散的问题,公式见图片1:
      要求解c(r,t),其中,r和t是变量,C1、b、a是常数,Uo(anr)的表达式见图片2,an是图片3中的方程的正根,问题就是,怎么在matlab 里求解第三个图片中的方程解出an。另外就是由于an有无穷多个,那么计算c(r,t)的时候要代入几个,求大牛帮忙指导程序编写,不胜感激

急求助求解一个含贝塞尔函数方程
1.png


急求助求解一个含贝塞尔函数方程-1
2.png


急求助求解一个含贝塞尔函数方程-2
3.png
回复此楼
共同学习
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CelestialCYJ

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
huazhengli: 金币+30, ★★★★★最佳答案 2015-01-03 22:51:26
an无穷多个,你可以根据你对Crt的精度要求,给一个比较大的n,比如50,100,你估计或者查文献,看是不是接近极限了。然后有了这个n,你把图3的贝塞尔函数用泰勒级数的形式表示出来,使U被表示成一个至少具有n阶的多项式,然后用roots命令,可以求得这个多项式的所有零点,前n个零点就是an,然后回带到图1的Crt中,你就可以求解了。
2楼2015-01-02 02:01:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhengli

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by CelestialCYJ at 2015-01-02 02:01:52
an无穷多个,你可以根据你对Crt的精度要求,给一个比较大的n,比如50,100,你估计或者查文献,看是不是接近极限了。然后有了这个n,你把图3的贝塞尔函数用泰勒级数的形式表示出来,使U被表示成一个至少具有n阶的多 ...

多谢多谢
共同学习
3楼2015-01-03 22:49:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huazhengli 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +11 kulium 2026-08-21 14/700 2026-08-24 02:12 by 流流伤
[基金申请] 范进中举一文的中心思想 +4 炎黄贵胄 2026-08-22 5/250 2026-08-24 02:09 by 流流伤
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
信息提示
请填处理意见