Question

Need help converting the following code from Matlab into Python:

N=2048;
fs=4.9;
t=0:1/fs:(N-1)/fs;
fs1=200;
t1=0:1/fs1:(N-1)/fs1;
x2=0.5+0.6366.*cos(2.*pi.*t1)+0.1273.*cos(10.*pi.*t1)-0.0909.*cos(14.*pi.*t1);
x=0.5+0.6366.*cos(2.*pi.*t)+0.1273.*cos(10.*pi.*t)-0.0909.*cos(14.*pi.*t);
X=fftshift(fft(x));
f=linspace(-fs/2,fs/2,N);
plot(f,abs(X)./N);
xlabel('f');
ylabel('|F(f)|');
title('magnitude spectrum of sampled signal');
x1=ifft(fftshift(X));
figure
plot(t(1:100),x1(1:100));
xlabel('t');
ylabel('f(t)');
title('f(t) obtained by inverse transform');
figure
plot(t1(1:1000),x2(1:1000));
xlabel('t');
ylabel('f(t)');
title('original f(t)');

Answer 1

`Hey,

Note: Brother if you have any queries related the answer please do comment. I would be very happy to resolve all your queries.

import matplotlib.pyplot as plt
import numpy as np
from sympy import fft
N=2048;
fs=4.9;
t=np.arange(0,(N-1)/fs+1/fs,1/fs);
fs1=200;
t1=np.arange(0,(N-1)/fs1+1/fs1,1/fs1);
x2=0.5+0.6366*np.cos(2*np.pi*t1)+0.1273*np.cos(10*np.pi*t1)-0.0909*np.cos(14*np.pi*t1);
x=0.5+0.6366*np.cos(2*np.pi*t)+0.1273*np.cos(10*np.pi*t)-0.0909*np.cos(14*np.pi*t);
X=np.fft.fftshift(fft(x));
f=np.linspace(-fs/2,fs/2,N);
plt.plot(f,abs(X)/N);
plt.xlabel('f');
plt.ylabel('|F(f)|');
plt.title('magnitude spectrum of sampled signal');
x1=np.ifft(np.fft.fftshift(X));
plt.figure
plt.plot(t[0:99],x1[0:99]);
plt.xlabel('t');
plt.ylabel('f(t)');
plt.title('f(t) obtained by inverse transform');
plt.figure
plt.plot(t1[0:99],x2[0:99]);
plt.xlabel('t');
plt.ylabel('f(t)');
plt.title('original f(t)');

Kindly revert for any queries

Thanks.