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

jc

to co zanotowałem:

całkowanie numeryczne- cumsum
Fc=1e5;
Fs=32*Fc;
Ts=1/Fs;
DewFM=5e4;
Fm=5000;
N=Fs/fm;
L=2;
N=L*N;
n=0:N-1;
m=cos(2*pi*Fm/fs*n); %sygnal modulujacy
maximum=max(abs(m));
kf=DewFM/maximum;

poniżej mam jeszcze napisane o demodulacji
d/dt Sfm(t) (chodzi o pochodną z sygnału zmodulowane fm)
wychodzi wtedy widmo sygnału zmodulowanego AM. Wtedy wystarczy już zdemodulować jak wam sie tam podoba.

Dzisiaj byłem u Dustora. powiedział że jutro będzie do zrobienia to co było na laborkach + wspomniał coś jeszcze o wyznaczaniu szerokości pasma, ale nie wiem o co chodzi, powiedział że może to byc w instrukcji, ja jeszze do niej nie zaglądałem...
Mi się na laborkach udało dojść do tego że wyciągnąłem z FM pochodną i wyszedł fajnie sygnał zmodulowany AM ale nie skopiowałem kodu. Dzisiaj będe działać jeszcze ale póóóóóźźźźnieeeej także jakby co to tu wkleje jak to zrobiłem.
tela.

grzesm

Moglby ktos wrzucic dzialajaca demodulacje przy pomocy detektora obwiedni

jc

zapomniałem dodac że Dustor jeszcze powiedział żebyśmy dobrze przejrzeli te programy z FM

drago256 · Dołączył: 02 Lis 2009 Posty: 50 Przeczytał: 0 tematów

FM z demodulacją

clear;clc;hold off;close;
f1=5000;
fc=1e5;
fs=32*fc;
dewf=5e4;
N=2*fs/f1;
n=0:N-1;
x=1*cos(2*pi*f1/fs*n);
kf=dewf/max(abs(x))
calka=cumsum(x)/fs;
s=cos(2*pi*fc/fs*n+2*pi*kf*calka);
subplot(411),plot(n,x);
subplot(412),plot(n,s);
subplot(413),stem(n/N*fs,abs(fft(s)));axis([0 3e5 0 max(abs(fft(s)))])
%demodulacja
s1=diff(s);
znak=(sign(s1)+1)/2;
s1=s1.*znak;
[B,A]=butter(5,10000/(fs/2));
wyn=filter(B,A,s1);
subplot(414);plot(wyn);

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

mój z ostatniej labory:

clear;clc;hold off;close;
f1=1000;
fc=1e5;
fs=32*fc;
ts=1/fs;
fm=5000;
dewf=5e4;
N=2*fs/fm;
n=0:N-1;
x=1*cos(2*pi*fm/fs*n);
kf=dewf/max(abs(x));
calka=cumsum(x)/fs;
s=1*cos(2*pi*fc/fs*n+2*pi*kf*calka);
subplot(511),plot(n,x);
subplot(512),plot(n,s);
subplot(513),stem(n/N*fs,abs(fft(s))); axis([0 3e5 0 max(abs(fft(s)))])
%subplot(413),plot(n/N*fs,abs(fft(s))); axis([0 3e5 0 max(abs(fft(s)))])
s1=diff(s);
roz=diff(s);
roz=[0 roz];
znak=(sign(s1)+1)/2;
s1=s1.*znak;
[B,A]=butter(5,10000/(fs/2));
wyn=filter(B,A,s1);
subplot(515);plot(wyn);
subplot(514),plot(n,roz);

Do 2 h wyślę to co uda mi się wymyślić

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

Ostatnio coś na labce brużdził z tymi prostokątami więc dla świetego spokoju proponuję :

fc=16e3;
Fs=64*fc;
fm=2000;
Ac=1;
Am=5;
L=3;
beta=10;

AA=zeros(1,128);
BB=ones(1,128);
CC=(-1)*ones(1,128);
DD=2*ones(1,128);
EE=(-2)*ones(1,128);
FF=3*ones(1,128);
GG=(-3)*ones(1,128);
HH=(-1)*ones(1,128);
x=[AA BB CC DD EE FF GG HH];

N=1024;
n=0:N-1;

kf=beta/max(abs(x));
kf=pi/2;
s=Ac*cos(2*pi*fc/Fs*n+kf*x);
widmo=abs(fft(s));

i=find(widmo>80);
H=zeros(1,length(widmo));
H(i)=1;
J=H.*widmo;

K=J(1:N/2);
dlugosc=(max(i)-min(i))*(Fs/N);

dem1=diff(s);
znak=(sign(dem1)+1)/2;
u=dem1.*znak;
[M,Wn]=BUTTORD(3000/(Fs/2),200000/(Fs/2),0.5,60);
[B,A]=BUTTER(M,Wn);
W=filter(B,A,u);

subplot(421),plot(n/Fs,x);
subplot(422),plot(n/Fs,s);
subplot(425),stem(n/N*Fs,J);axis([0,5e4,min(J),max(J)]);
subplot(427),stem(n/N*Fs,widmo);
axis([0,5e4,min(widmo),max(widmo)]);
n=0:N-2;
subplot(423),plot(n/Fs,dem1);
subplot(424),plot(n/Fs,u);
subplot(426),plot(n/Fs,W);

GREY

Moze jakis komentarz do programu?