code:
x=linspace(-5,5,1000);
n=double(1:1:100);
l=2;
f=sin(x);
g=cos(x);
s=0;
landa=((n*pi)/l );
a=(2/l)*int(f.*sin((n*pi)/l)*x,0,l);
b=(2/(l*landa))*int(g*sin((n*pi*x)/l,0,l));
for n=1:1:1000
p=((a*cos(landa*t)+ b*sin(landa*t))*sin((n*pi*x)/l));
s=s + p;
end
plot (x,t,s)

댓글 수: 2

a=(2/l)*int(f.*sin((n*pi)/l)*x,0,l);
In the above f,x is 1*1000 and n is 1*100..you cannot mulitply them....rethink on your code.
Walter Roberson
Walter Roberson 2018년 12월 20일
please do not close questions that have an answer

댓글을 달려면 로그인하십시오.

 채택된 답변

KSSV
KSSV 2018년 12월 14일

1 개 추천

syms x n
% n=double(1:1:100);
l=2;
f=sin(x);
g=cos(x);
s=0;
landa=((n*pi)/l );
a=(2/l)*int(f*sin((n*pi)/l)*x,0,l);
b=(2/(l*landa))*int(g*sin((n*pi*x)/l),0,1);

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 Mathematics에 대해 자세히 알아보기

태그

질문:

2018년 12월 14일

댓글:

2018년 12월 20일

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by