Experiment On Fourier Transform Code Used: %length of Vector N Be Any Function
Experiment On Fourier Transform Code Used: %length of Vector N Be Any Function
Experiment On Fourier Transform Code Used: %length of Vector N Be Any Function
N = 10; %Length of vector t=0:.1:N-1; x = 10*sin(2*pi*t)+5*cos(2*pi*t); %Can be any function subplot(221); plot(t,x); xlabel('t');ylabel('x'); wn = exp(-2*pi*i/N); %Kernel X = zeros(1,N); %FFT of x for k = 1:N X(k) = 0; for j = 1:N X(k) = X(k) + x(j)*wn^((j-1)*(k-1)/N); %fourier transform end end subplot(223); stem(X); xlabel('f');ylabel('X(k)');
OUTPUT:
OUTPUT:
CODE USED
%prog to demonstrate discrete fourier transform N=4; x1=[1,2,3,4];%taken any signal; v=length(x1); v1=1:v; subplot(221); stem(v1,x1);%plotting; xlabel('n'); ylabel('x1(n)'); x2=x1'*x1;%4*4 matrix will get formed; w=exp(-2*pi*i/N);%twiddle factor x3=w.^x2;%applying factor y1=[1,-2];%signal for evaluation; y1=[y1 zeros(1,N-length(y1))];%zero padding; v2=length(y1); v3=1:v2; subplot(222); stem(v3,y1); xlabel('n'); ylabel('y1(n)'); y2=y1'; n1=x3*y2;%resultant signal; n=abs(n1);%magnitude; subplot(223); stem(n); xlabel('n'); ylabel('abs(n)');
OUTPUT:
CODE USED
Fs = 44100; Fn = Fs/2; F3db = 500; filtdes = fdesign.lowpass('n,f3db', 5, F3db, Fs); Hd = design(filtdes,'butter'); Hd = convert(Hd, 'df1sos'); fvtool(Hd, 'Fs', Fs, 'FrequencyScale', 'log');
OUTPUT: