Aufgabe 4 :restart : N:=25 :T:=5 : omega:=2*Pi/T ;for n from 0 by 1 while n <= N do
a[n]:= 2/T*(int( 2*cos(n*omega*t) , t=-2..0)+int( (2-t)*cos(n*omega*t) , t=0..2)) :
b[n]:= 2/T*(int( 2*sin(n*omega*t) , t=-2..0)+int( (2-t)*sin(n*omega*t) , t=0..2)):
c[n]:= evalf(sqrt( a[n]*a[n] + b[n]*b[n] )) : end do :f:=(x) -> evalf( a[0]/2 + sum( a[m]*cos(m*omega*x) + b[m]*sin(m*omega*x) , m=1..N) );plot( f(x) , x=-5..10 , -1..3 ,numpoints=2000 ) ;plot( [[ m, c[m]] $m=0..N] , x=0..15, title="Amplituden - Spektrum", style=point,xtickmarks=15, symbol=circle,symbolsize=20);