forked from acco/chti23
49 lines
851 B
Matlab
49 lines
851 B
Matlab
close all;
|
|
Flasers = [85.005E3,90E3,95E3,100E3,115E3,120E3];
|
|
Fsin = Flasers;
|
|
deltaF = 5E3;
|
|
T = 1/deltaF;
|
|
Fe = 320E3;%kHz
|
|
Te = 1/Fe;
|
|
M= Fe*T;
|
|
Tsim = T-Te;
|
|
frequence = deltaF*(0:1:M-1);
|
|
|
|
|
|
Res=sim('SimuliDFT');
|
|
|
|
|
|
|
|
figure
|
|
plot(Res.temps,Res.Sinus_Continu);
|
|
hold on;
|
|
plot(Res.temps_Echanti,Res.Sinus_Echanti, 'o');
|
|
sig= Res.Sinus_Echanti;
|
|
%sig(1:20,4)=0;
|
|
echantFft = fft(sig,M);
|
|
|
|
|
|
figure
|
|
stem(frequence,abs(echantFft)/M);
|
|
|
|
|
|
%carreeee-----------------------
|
|
|
|
% figure(1)
|
|
% plot(Res.temps,Res.Carre_Continu1);
|
|
% hold on;
|
|
% plot(Res.temps,Res.Sinus_Carre_Continu);
|
|
% hold on;
|
|
% plot(Res.temps_Echanti,Res.Sinus_Carre_Echanti, 'o');
|
|
% echantFft_Carre = fft(Res.Sinus_Carre_Echanti,M);
|
|
%
|
|
%
|
|
%
|
|
% figure(2)
|
|
% stem(frequence,abs(echantFft_Carre)/M);
|
|
%
|
|
% Fp = tf(1,[1.7483*10^(-23), 7.6663*10^(-18), 1.162*10^(-11), 3.0332*10^(-6),1]);
|
|
% figure(3)
|
|
% bode(Fp);
|
|
|
|
|