Help with plotting 3d function

Question

Jacob Andréasson 2019-11-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/490160-help-with-plotting-3d-function

评论： Jacob Andréasson 2019-11-9

Hey!

I need help with exercise 2.b. This is my script that i have written

%% Exercise 2.b

clear

clf

clc

k=1;

L=6;

T=10

x=linspace(0,L,100); t=linspace(0,T,100);

for i=1 : length(x)

for j=1 : length(t)

u(i,j)=(8/(pi^2)).*(1./(n.^2)).*sin(n.*pi/2).*cos(n.*pi.*t(j)/6).*sin(n.*pi.*x(i)/6);

end

surf(t,x,u)

xlabel('t'); ylabel('x')

The problem that i need help with is how to put in the value of n. In the exercise it says that i should put n between 1->100 but i dont know how to do it....

Could someone pleas help me...?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Star Strider 2019-11-9

Note that n is the index with respect to

, and your code does not appear to be doing that summation.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Praveen Iyyappan Valsala 2019-11-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/490160-help-with-plotting-3d-function#answer_400648

在 MATLAB Online 中打开

Most intuitive way is to add another loop to sum up function outputs of all n.

for i=1 : length(x)
    for j=1 : length(t)
        for n=1:100
            u(i,j)=u(i,j)+(8/(pi^2)).*(1./(n.^2)).*sin(n.*pi/2).*cos(n.*pi.*t(j)/6).*sin(n.*pi.*x(i)/6);
        end
    end
end

Little more efficient way.

n=1:100;
for i=1 : length(x)
    for j=1 : length(t)
        u(i,j)=(8/(pi^2)).*sum((1./(n.^2)).*sin(n.*pi/2).*cos(n.*pi.*t(j)/6).*sin(n.*pi.*x(i)/6));
    end
end

if you want even more cool code. learn about meshgrid function.