1-vazifa Jadval ko’rinishida berilgan funksiyalarni Furye qatori asosida raqamli signallar yetakchi garmonikalarini aniqlash usuli bilan analitik ko’rinishda ifodalang. Algoritmini tuzib

1-vazifa
Jadval ko’rinishida berilgan funksiyalarni Furye qatori asosida raqamli signallar yetakchi garmonikalarini aniqlash usuli bilan analitik ko’rinishda ifodalang. Algoritmini tuzib, dasturlash tilida dastur kodini yozing va natija oling.
Laboratoriya ishiga doir variant:




Dastur kodi:

import math

def sum(f, w, t, k):

sumfsin = 0
for i in range(len(f)-1):
sumfsin += f[i]*(math.sin(k) * (w[k] * (i + 1 /2)) * (t[k] - t[k-1])- math.sin(k) * w[k] * (i-1 / 2) * (t[k] - t[k-1]) )
return sumfsin
t = []
k = 0
w = [0]
a = []
T = 2 * math.pi
for i in range(0, 21):
t.append(k)
k += 0.000005
k = 0
for i in range(1, 21):
w.append(2 * math.pi / t[i])
f = [0, 3.526622, 5.70627, 5.706442, 3.527071, 0.000556, -3.52617, -5.7061, -5.70661, -3.52752, -0.00111, 3.525722, 5.705927, 5.706786, 3.527971, 0.001668, -3.52527, -5.70575, -5.70696, -3.52842, -0.00222]
sumhf = 0
for i in range(len(f)-1):
sumhf += (t[i] - t[i-1]) * f[i]
a.append(2/T * (t[0] / 2 * f[0] + sumhf + t[-1] / 2 * f[-1]))
b = [0]
for i in range(1, 21):
a.append(2/(i * w[i] * T)*(f[0] * math.sin(i * w[i] * (t[i] - t[i-1]) / 2) + sum(f, w , t, i) + f[-1] * (math.sin(i) * w[i] * T - math.sin(i) * w[i] * (T - (t[i] - t[i-1])))))
print(a)
for i in range(1, 21):
b.append(2/(i * w[i] * T)*(f[0] * math.cos(i * w[i] * (t[i] - t[i-1]) / 2) + sum(f, w , t, i) + f[-1] * (math.cos(i) * w[i] * T - math.cos(i) * w[i] * (T - (t[i] - t[i-1])))))
print(b)
c = []
l = []
for i in range(1, 21):
c.append(math.sqrt(a[i]**2 + b[i]**2))
l = c
c.sort(); c.reverse()
Dastur natijasi:

2-vazifa:
Jadval ko’rinishida berilgan funktsiyalarni eng kichik kvadratlar usuli bilan analitik ko’rinishda ifodalang. Algoritmini tuzib, dasturlash tilida dastur kodini yozing va natija oling.





Dastur kodi:

n=6;
x=[1.10 1.15 1.20 1.25 1.30 1.35];

y=[3.170 3.261 3.855 3.919 3.422 3.337];
b=(sum(y)*sum(x.*x)-sum(x)*sum(x.*y))/(n*sum(x.*x)-sum(x)*sum(x));
a=(n*sum(x.*y)-sum(x)*sum(y))/(n*sum(x.*x)-sum(x)*sum(x));
t=[0:0.001:6];
z=a*t+b;
r=exp((1/3)*log(t+1))+1;
plot(x,y,'o',t,z,t,r);


Dastur natijasi: