Forum Forum studentów IV roku EiT Politechniki Śląskiej Strona Główna

Forum studentów IV roku EiT Politechniki Śląskiej

 

program napisany na laborce nr 4 (Dustor)

 
Napisz nowy temat   Odpowiedz do tematu    Forum Forum studentów IV roku EiT Politechniki Śląskiej Strona Główna -> Semestr 5 / PTK
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
Stefan99




Dołączył: 10 Paź 2007
Posty: 453
Przeczytał: 0 tematów


PostWysłany: Czw 19:00, 05 Lis 2009    Temat postu: program napisany na laborce nr 4 (Dustor)

clear
f1=1000;
f2=1500;
m=1;
fc=32*f1;
fs=8*fc;
L=2;
N=L*fs/f1;
n=0:N-1;
x=10*cos(2*pi*f1*n/fs);
xz=x/(max(abs(x)));
s=5*(1+m*xz).*cos(2*pi*(fc)*n/fs);
subplot(411), plot(n,x)
subplot(412), plot(n,s)
v=s.*cos(2*pi*fc/fs*n);
[B,A]=butter(1,f1*2/(fs/2));
w=filter(B,A,v)
subplot(413), plot(n,w)
s=((1+sign(s))/2).*s;
w=filter(B,A,s)
subplot(414), plot(n,w)
Powrót do góry
Zobacz profil autora
xaver




Dołączył: 01 Lis 2008
Posty: 301
Przeczytał: 0 tematów


PostWysłany: Czw 20:42, 05 Lis 2009    Temat postu:

niech Ci Bozia wynagrodzi Wesoly
Powrót do góry
Zobacz profil autora
L3oN




Dołączył: 31 Paź 2007
Posty: 517
Przeczytał: 0 tematów


PostWysłany: Wto 20:49, 17 Lis 2009    Temat postu:

Zgodnie z prośbą, modulacja funkcji sinc:

Kod:

clear; clc; clf;
f1=1000;
m=1;
Fc=20000;
Fs=200000;
L=4;
N=L*Fs/f1;
n=(-N):N
x=sinc(2*f1*n/Fs)
x=x.^2;
xz=x/max(abs(x));
s=5*x.*cos(2*pi*Fc/Fs*n);
Wx=fft(x);
Ws=fft(s);
u=s.*cos(2*pi*Fc/Fs*n);
[B, A] = butter(5,2*Fc/Fs)
w=filter(B, A, u);
Ww=fft(w);
subplot(411);
plot(1000*n/Fs,x);
title('S. modulujacy');
xlabel('Czas (ms)');
ylabel('Wartosc');
subplot(412);
plot(1000*n/Fs,s);
title('S. zmodulowany');
xlabel('Czas (ms)');
ylabel('Wartosc');
subplot(413);
plot(n(N:2*N)*Fs/1000,abs(Ws(N:2*N)));
axis( [120000, 140000, 0, max(abs(Ws))] );
title('S. zmodulowany');
xlabel('Czestotliwosc (kHz)');
ylabel('Amplituda');
subplot(414);
plot(1000*n/Fs,w);
title('S. zdemodulowany');
xlabel('Czas (ms)');
ylabel('Wartosc');

Jeżeli chodzi o widmo, to musicie zamienić lewą połówkę wektora wyjściowego z prawą i wtedy dopiero to wyskalować. Opis osi f będzie poprawny. Na laborce nie zdążyłem wpaść na to, jak to zrobić.


Ostatnio zmieniony przez L3oN dnia Wto 20:50, 17 Lis 2009, w całości zmieniany 1 raz
Powrót do góry
Zobacz profil autora
xanno




Dołączył: 25 Lut 2008
Posty: 44
Przeczytał: 0 tematów


PostWysłany: Śro 13:24, 25 Lis 2009    Temat postu:

Może ma ktoś program z kolejnych laborek?
Powrót do góry
Zobacz profil autora
Wyświetl posty z ostatnich:   
Napisz nowy temat   Odpowiedz do tematu    Forum Forum studentów IV roku EiT Politechniki Śląskiej Strona Główna -> Semestr 5 / PTK Wszystkie czasy w strefie CET (Europa)
Strona 1 z 1

 
Skocz do:  
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach


fora.pl - załóż własne forum dyskusyjne za darmo
Powered by phpBB © 2001, 2005 phpBB Group
deox v1.2 // Theme created by Sopel & Download

Regulamin