ربي يعيشكم ساعدوني ............matlab

[الجزائر بلادي]

السلام عليكم أرجوا منكم المساعدة في حل هدا التمرين في أقرب وقت


soit le dignal pèriodique de pèriode T=2pi
x(t)=0 si -pi<t<0
x(t)=sin(t) si 0<t<pi
tracer x(t) en utilisant l'instruction plot(t,x)

calculer la seriede fourier de x(t)
en utilisant l'instruction hold on ,tracersur la mème figure la serie de fourier pour n=10,100 et 1000
commenter le rèsultat
و جزاكم الله كل خير

[lotfi_dk90]

• function s(n)
%
% SF Calcule la série de Fourier correspondant au signal carré
% et affiche la somme de la série jusqu'à l'harmonique de
% rang n
%
% n : nombre d'harmoniques
%
close all;
m=floor(n/2);
Ts=0.22e-3; % Période d'échantillonnage
t0=24.e-3; % Période du signal carré
N=500; % nombre de points
p=(-m-1):m ;
cn=((-1).^p)./(2*p+1)/pi;
cc=(2*pi*i*Ts)*(2*p+1)/t0;
for k=0:N-1
vn=exp(cc*k);
s(k+1)=cn*vn';
end
t=(0:N-1)*Ts;
f=(0:N-1)/N/Ts;
subplot(211); plot(t,real(s)); grid; title('signal')
subplot(212); plot(f,abs(fft(s)/N)); grid; title('signal')

[الجزائر بلادي]

بارك الله فيك
لكن هل يمكن أن توضح لي أكثر اجابة السؤال 2
la serie de fourier

[lotfi_dk90]

المعاملات cn cc vn لserie de fourier

[lotfi_dk90]

سؤالك هو serie de fourie يعني تحديد المعاملات ?

[lotfi_dk90]

هذا الملف من رفعي وبباسوورد المنتدى وفيه tp de traitement du signal

الرابط: http://www.mediafire.com/?0j981cz5gceqbil

الباسوورد: www.djelfa.info/vb

[الجزائر بلادي]

بارك الله فيك أخي لكن ادا سمحت أريد شرح أكثر
مع العلم أن حل السؤال 2 ليس matlab
و انما حسابيا
و فيه مفردات لم نستعملها من قبل مثل
floor
الني أنا اختصاص control et instrumentation
و بارك الله فيك

[lotfi_dk90]

بالتوفيق ان شاء الله

[الجزائر بلادي]

شكرا لك أخي
اكن هل تستطيع التبسيط أكثر؟

[lotfi_dk90]

اسألي بالضبط ? ارجوا التدقيق

[الجزائر بلادي]

مثلا خطوات حساب la serie de fourier
بدون استعمال matlab

[lotfi_dk90]

>>figure(1)
>>T=0.01;
>>Fs=1/T; % fréquence d'échantillonnage de 100 Hz
>>t = 0:T:1;
>>x = sin(2*pi*10*t); % signal fréquence de 10 Hz
>>plot(t,x)
>>title('Signal')
>>xlabel('temps (ms)')
>>figure(2)
>>FFT_size=2048;
>>X=fft(x,FFT_size);
>>f=(0:FFT_size/2)/(FFT_size/2)*Fs/2;
>>plot(f,abs(X(1:FFT_size/2+1))); % on voit bien que le signal est à 10 Hz
>>title('Signal en fréquence');
>>xlabel('Fréquence(Hz)');
>>ylabel('Amplitude');

[lotfi_dk90]

عذرا لم استعمل hold on

[الجزائر بلادي]

لكن الاجابتين مختلفتين فأيهما أعتمد ؟

[الجزائر بلادي]

أرجو المساعدة من فضلك