Plot is zeroed out, how to fix this?
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Hello,
I wrote a code that should give set values and plot them. For some reason i keep getting zero values. Any thoughts?
Now i am faced with a different problem. My plot comes out zeroed and i have no clue why. Any thoughts?
clc;
clear;
N=25;
% --Grid--
gridX = linspace(0,1,N);
gridT = linspace(0,0.5,N);
[x,t] = meshgrid(gridX,gridT);
A = linspace(0,1,N);
u = zeros(N,N);
xint = linspace(0,1,N);
%--Computation of An Terms--
n=1;
while n <= N;
I(n) = -(4*sin((pi*n)/2))/(pi*(n^2 - 4));
A(n) = A(n) + I(n);
n = n + 1;
end
while n <= N
exponent = -(n.^2)*(pi.^2)*(t); % n replaces lamda
T = exp(exponent);
X = A(n).*sin(n.*pi.*x);
total = T*X;
u(n) = u(n) + total; %--Sum of all solutions--
n = n + 1;
end
%--Plot--
figure;
surfc(x,t,u)
title('Project 2')
xlabel('X')
ylabel('Time')
0 Kommentare
Antworten (1)
Mischa Kim
am 11 Feb. 2014
Mike, after initializing u to zero you enter the first loop that increments n all the way to N. Which is why you never enter the second loop that is supposed to update u.
4 Kommentare
Mischa Kim
am 11 Feb. 2014
Just looked at your code. Since I do not know what exactly you need to compute here are some thoughts:
- Change the corresponding command to: u = u + total. That way at least the dimensions work out.
- You are generating some huge numbers along the way so there are probably issues with some of the equations.
Siehe auch
Kategorien
Mehr zu Creating and Concatenating Matrices finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!