Plotting a Sum of Series
1 visualizzazione (ultimi 30 giorni)
Mostra commenti meno recenti
I have to plot these equations about the heat in a aluminium bar:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/170549/image.jpeg)
But for some reason, my graph of u(x,t)(in this case u on the code) is showing something different (i think the correct would be a exponential). Can someone review my code and see if is there any mistakes? Obs: n is the value where sin=0, for example sin(n*pi)=0 like in pi, 2pi, 3pi...
clc
clear all
syms x k t n
for n=1:1:20
cn=(100/(n*pi))*(cos(pi*x/4)-cos(3*pi*x/4));
for t=0:5
num=(cos(pi*x/4)-cos(3*pi*x/4));
t1=exp((-(k^2)*(pi^2)*t)/1600);
t2=sin(k*pi*x/40);
S1=symsum((num*t1*t2),k,1,inf);
u=((100/pi)*S1);
end
end
fplot(cn)
figure
fplot(u)
1 Commento
David Goodmanson
il 2 Dic 2017
Hi Italo,
One thing for sure, the integral in the top equation is incorrect. I believe it should be
(100/(n*pi))*( cos(n*pi/4) - cos(3*n*pi/4) )
and a sign of something wrong was that the c_n can't depend on x.
Could you state the actual problem, in particular the length of the bar?
Risposte (2)
Roger Stafford
il 2 Dic 2017
You have forgotten to divide by k in calculating 'u' within your for-loop. You should have:
S1=symsum((num/k*t1*t2),k,1,inf);
Vidya Bhadgaonkar
il 10 Ott 2019
I think you have one mistake in
t1=exp((-(k^2)*(pi^2)*t)/1600);
write it as
t1=exp((-(n^2)*(pi^2)*t)/1600);
0 Commenti
Vedere anche
Categorie
Scopri di più su Thermal Analysis in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!