策划一笔系列

4视图(30天)
Italo巴罗斯
Italo巴罗斯 2017年12月2日
回答: 维迪雅Bhadgaonkar2019年10月10日
我必须情节这些方程对热铝杆:
但是由于一些原因,我图的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
大卫Goodmanson 2017年12月2日
嗨,想一下,
一件事可以肯定的是,上面的积分方程是不正确的。我认为它应该
(100 / (n * pi)) * (cos (n *π/ 4)- cos (3 * n *π/ 4))
错了的标志是c_n不能依赖于x。
你能状态的实际问题,特别是栏的长度吗?

登录置评。

答案(2)

罗杰·斯塔福德
罗杰·斯塔福德 2017年12月2日
你忘了除以k计算“u”在你的for循环。你应该:
S1 = symsum ((num / k * t1 * t2), k, 1,正);
1评论
Italo巴罗斯
Italo巴罗斯 2017年12月2日
Tkx答案!但是,即使把“k”,结果真的很奇怪:

登录置评。


维迪雅Bhadgaonkar
维迪雅Bhadgaonkar 2019年10月10日
我认为你有一个错误
t1 = exp ((- (k ^ 2) *(π^ 2)* t) / 1600);
把它写成
t1 = exp ((- (n ^ 2) *(π^ 2)* t) / 1600);

类别

找到更多的在热力学和热传递帮助中心文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!