1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
clc
//Chapter 1 Signals
//Example 1.5, page no 19
//given
T=1,t0=0,wo=2*%pi
P=1
t=0:0.001:1
f=P*t
//The trigonometric Fourier series coeff for given function
a0=(1/T)*inttrap(t,f)
for n=1:10
f1=(P*t).*cos(wo*n*t)
a(n)=(2/T)*inttrap(t,f1)
if a(n)<0.01 then
a(n)=0
end
end
for n=1:10
f2=(P*t).*sin(2*%pi*(1/T)*n*t)
b(n)=(2/T)*inttrap(t,f2)
end
// Displaying trigonometric Fourier series coeff
mprintf('The value of a0 is: %f\n',a0)
disp('The values of a(n): (upto n=10)')
for n=1:10
disp(a(n))
end
disp('The values of b(n): (upto n=10)')
for n=1:10
disp(b(n))
end
mprintf('The trigonometric Fourier series for given function\n can be written as:\n')
mprintf('f(t)=%f%fsin(2*pi*t)%fsin(4*pi*t)\n%fsin(6*pi*t)%fsin(8*pi*t)%fsin(10*pi*t)\n%fsin(12*pi*t)%fsin(14*pi*t).......',b(1),b(2),b(3),b(4),b(5),b(6),b(7),b(8))
|
