【答案】应助回帖

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有62人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

1楼2016-04-04 23:54:37

aristarchus

铜虫 (初入文坛)

应助: 2 (幼儿园)
金币: 50.1
帖子: 41
在线: 28.2小时
虫号: 540756
注册: 2008-04-06
专业: 凝聚态物性 II ：电子结构

fun = @(x,t) t1+x./H.*(t2-t1)+symsum(2./(n*3.14).*((t0-t1)-(t0-t2).*(-1)^n).*sin(n.*3.14.*x./H)*exp(-ap.*n.*n.*3.14.*3.14.*t./H^2), n, 1, Inf);

然后用eval(fun(x,t))求数值解。如果不行，那就试着用一个大数代替Inf。

5楼2016-04-06 02:56:08

查看全部 5 个回答

aristarchus

铜虫 (初入文坛)

应助: 2 (幼儿园)
金币: 50.1
帖子: 41
在线: 28.2小时
虫号: 540756
注册: 2008-04-06
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
我不是好汉: 金币+10, ★★★很有帮助 2016-04-05 08:36:21

Matlab定义函数的格式是 fun = @(x,t) T_1+....;
级数部分写成 symsum(函数体, n, 1, Inf) 就行了，其它部分原样照写。
调用就直接用fun(x,t)的形式

2楼2016-04-05 05:40:01

我不是好汉

木虫 (著名写手)

应助: 0 (幼儿园)
金币: 4496.1
散金: 396
红花: 4
帖子: 1359
在线: 183.3小时
虫号: 3731592
注册: 2015-03-12
性别: GG
专业: 计算机软件

引用回帖:

2楼: Originally posted by aristarchus at 2016-04-05 05:40:01
Matlab定义函数的格式是 fun = @(x,t) T_1+....;
级数部分写成 symsum(函数体, n, 1, Inf) 就行了，其它部分原样照写。
调用就直接用fun(x,t)的形式

大神能不能帮我写下啊，我试了下，还是有错误啊

3楼2016-04-05 08:36:06

aristarchus

铜虫 (初入文坛)

应助: 2 (幼儿园)
金币: 50.1
帖子: 41
在线: 28.2小时
虫号: 540756
注册: 2008-04-06
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
我不是好汉: 金币+10, ★★★很有帮助 2016-04-10 22:23:10

t0 = 0;
t1 = 1;
H = 1;
t2 = 2;
ap = 1;
syms n x t
fun = @(x,t) t1+x./H.*(t2-t1)+sum(subs(2./(n*3.14).*((t0-t1)-(t0-t2).*(-1)^n).*sin(n.*3.14.*x./H)*exp(-ap.*n.*n.*3.14.*3.14.*t./H^2), n, 1, Inf));

然后用 fun(x,t) 调用

4楼2016-04-05 17:04:37