策划一笔系列
4视图(30天)
显示旧的评论
我必须情节这些方程对热铝杆:
![](http://www.tatmou.com/matlabcentral/answers/uploaded_files/170549/image.jpeg)
但是由于一些原因,我图的u (x, t)(在这种情况下你的代码)显示不同的东西(我认为正确的将是一个指数)。可以有人评论我的代码,看看有没有错误?奥林匹克广播服务公司:n是罪= 0的值,例如罪(n *π)= 0像π,2π,3π…
clc
清晰的所有
信谊x k t n
为n = 1:1:20
cn = (100 / (n *π))* (cos(π* x / 4) cos(3 *π* x / 4));
为t = 0:5
num = (cos(π* x / 4) cos(3 *π* x / 4));
t1 = exp ((- (k ^ 2) *(π^ 2)* t) / 1600);
t2 =罪(k *π* x / 40);
S1 = symsum ((num * t1 * t2), k, 1,正);
u = ((100 / pi) * S1);
结束
结束
fplot (cn)
图
fplot (u)
1评论
大卫Goodmanson
2017年12月2日
嗨,想一下,
一件事可以肯定的是,上面的积分方程是不正确的。我认为它应该
(100 / (n * pi)) * (cos (n *π/ 4)- cos (3 * n *π/ 4))
错了的标志是c_n不能依赖于x。
你能状态的实际问题,特别是栏的长度吗?
答案(2)
维迪雅Bhadgaonkar
2019年10月10日
我认为你有一个错误
t1 = exp ((- (k ^ 2) *(π^ 2)* t) / 1600);
把它写成
t1 = exp ((- (n ^ 2) *(π^ 2)* t) / 1600);