matlab.rar_Analog modulation_creation

clc; close all; clear all; [x, Fs] = wavread('orjinal.wav'); Fc=1500; sinyal=x'; L = length(sinyal); T=1/Fs; t=0:T:(L-1)*T; tasiyici=cos(2.*pi.*Fc.*t); sinyal_fft = fft(sinyal); sinyal_fft_abs = abs(fftshift(sinyal_fft)); tasiyici_fft = fft(tasiyici); tasiyici_fft_abs = abs(fftshift(tasiyici_fft)); f=linspace(-Fs/2,Fs/2,L); mag_N0 = 300; mag_N = mag_N0+400; figure(1),subplot(3,1,1); plot(t,sinyal); title('Orijinal ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(sinyal(mag_N0:mag_N)); title('Orijinal ses(Büyütülmüþ)'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,sinyal_fft_abs); title('Original ses(frekans)'); xlabel('frekans');ylabel('Genlik'); figure(2),subplot(3,1,1); plot(t,tasiyici); title('Tasiyici'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(tasiyici(mag_N0:mag_N)); title('Tasiyici(Büyütülmüþ)'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3) plot(f,tasiyici_fft_abs); title('Tasiyici(frekans)'); xlabel('frekans');ylabel('Genlik'); disp('AM=1 DSB=2 SSB=3 FM=4 Modülasyon türü için verilenleri kullanýn') a=input('Modülasyon Türü Giriniz:'); if a==1; ka=input('Modülasyon indeksini girinizi:'); mod = (1+ka.*sinyal).*tasiyici; y = awgn(mod,10,'measured'); [a,b] = butter(10,2*Fc/Fs); z= mod.*tasiyici; e= filter(a,b,z).*2; w=e-1; wavwrite(w',Fs,'AM.wav'); mod_fft = fft(mod); mod_fft_abs = abs(fftshift(mod_fft)); y_fft = fft(y); y_fft_abs = abs(fftshift(y_fft)); e_fft = fft(w); e_fft_abs = abs(fftshift(e_fft)); figure(3),subplot(3,1,1); plot(mod); title('Modülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(mod(mag_N0:mag_N)); title('Modülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,mod_fft_abs); title('Modülasyon ses(frekans)'); xlabel('frekans');ylabel('Genlik'); figure(4),subplot(3,1,1); plot(y); title('Gürültülü ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(y(mag_N0:mag_N)); title('Gürültülü ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,y_fft_abs); title('Gürültülü ses(frekans)'); xlabel('frekans');ylabel('Genlik'); figure(5),subplot(3,1,1); plot(e); title('Demodülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(e(mag_N0:mag_N)); title('Demodülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,e_fft_abs); title('Demodülasyon ses(frekans)'); xlabel('frekans');ylabel('Genlik'); elseif a==2; mod = ammod(sinyal,Fc,Fs); y = awgn(mod,10,'measured'); e=amdemod(y,Fc,Fs); wavwrite(e',Fs,'DSB.wav'); mod_fft = fft(mod); mod_fft_abs = abs(fftshift(mod_fft)); y_fft = fft(y); y_fft_abs = abs(fftshift(y_fft)); e_fft = fft(e); e_fft_abs = abs(fftshift(e_fft)); figure(3),subplot(3,1,1); plot(mod); title('Modülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(mod(mag_N0:mag_N)); title('Modülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,mod_fft_abs); title('Modülasyon ses(frekans)'); xlabel('frekans');ylabel('Genlik'); figure(4),subplot(3,1,1); plot(y); title('Gürültülü ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(y(mag_N0:mag_N)); title('Gürültülü ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,y_fft_abs); title('Gürültülü ses(frekans)'); xlabel('frekans');ylabel('Genlik'); figure(5),subplot(3,1,1); plot(e); title('Demodülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(e(mag_N0:mag_N)); title('Demodülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,e_fft_abs); title('Demodülasyon ses(frekans)'); xlabel('frekans');ylabel('Genlik'); elseif a==3; mod1= ssbmod(sinyal,Fc,Fs,0); mod2=ssbmod(sinyal,Fc,Fs,0,'upper'); y1= awgn(mod1,10,'measured'); y2=awgn(mod2,10,'measured'); e1=ssbdemod(y1,Fc,Fs); e2=ssbdemod(y2,Fc,Fs); wavwrite(e1',Fs,'SSB.wav'); mod1_fft = fft(mod1); mod1_fft_abs = abs(fftshift(mod1_fft)); y1_fft = fft(y1); y1_fft_abs = abs(fftshift(y1_fft)); e1_fft = fft(e1); e1_fft_abs = abs(fftshift(e1_fft)); mod2_fft = fft(mod2); mod2_fft_abs = abs(fftshift(mod2_fft)); y2_fft = fft(y2); y2_fft_abs = abs(fftshift(y2_fft)); e2_fft = fft(e2); e2_fft_abs = abs(fftshift(e2_fft)); figure(3),subplot(3,1,1); plot(mod1); title('Modülasyon ses(LSB)'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(mod1(mag_N0:mag_N)); title('Modülasyon ses büyütülmüþ (LSB) '); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,mod1_fft_abs); title('Modülasyon ses(frekans)(LSB)'); xlabel('frekans');ylabel('Genlik'); figure(4),subplot(3,1,1); plot(y1); title('Gürültülü ses(LSB)'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(y1(mag_N0:mag_N)); title('Gürültülü ses büyütülmüþ(LSB)'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,y1_fft_abs); title('Gürültülü ses(frekans)(LSB)'); xlabel('frekans');ylabel('Genlik'); figure(5),subplot(3,1,1); plot(e1); title('Demodülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(e1(mag_N0:mag_N)); title('Demodülasyon ses büyütülmüþ'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,e1_fft_abs); title('Demodülasyon ses(frekans)'); xlabel('frekans');ylabel('Genlik'); figure(6),subplot(3,1,1); plot(mod2); title('Modülasyon ses(USB)'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(mod2(mag_N0:mag_N)); title('Modülasyon ses büyütülmüþ(USB)'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,mod2_fft_abs); title('Modülasyon ses(frekans)(USB)'); xlabel('frekans');ylabel('Genlik'); figure(7),subplot(3,1,1); plot(y2); title('Gürültülü ses(USB)'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(y2(mag_N0:mag_N)); title('Gürültülü ses büyütülmüþ (USB)'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,y2_fft_abs); title('Gürültülü ses(frekans)(USB)'); xlabel('frekans');ylabel('Genlik'); figure(8),subplot(3,1,1); plot(e2); title('Demodülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(e2(mag_N0:mag_N)); title('Demodülasyon ses büyütülmüþ'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,e2_fft_abs); title('Demodülasyon ses(frekans)'); xlabel('frekans');ylabel('Genlik'); elseif a==4; dev=input('Frekans sapmasýný giriniz='); mod = fmmod(sinyal,Fc,Fs,dev); y = awgn(mod,10,'measured'); e=fmdemod(y,Fc,Fs,dev); wavwrite(e',Fs,'FM.wav'); mod_fft = fft(mod); mod_fft_abs = abs(fftshift(mod_fft)); y_fft = fft(y); y_fft_abs = abs(fftshift(y_fft)); e_fft = fft(e); e_fft_abs = abs(fftshift(e_fft)); figure(3),subplot(3,1,1); plot(mod); title('Modülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(mod(mag_N0:mag_N)); title('Modülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,mod_fft_abs); title('Modülasyon ses(frekans)'); xlabel('frekans');ylabel('Genlik'); figure(4),subplot(3,1,1); plot(y); title('Gürültülü ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(y(mag_N0:mag_N)); title('Gürültülü ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,y_fft_abs); title('Gürültülü ses(frekans)'); xlabel('frekans');ylabel('Genlik'); figure(5),subplot(3,1,1); plot(e); title('Demodülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,2); plot(e(mag_N0:mag_N)); title('Demodülasyon ses'); xlabel('zaman');ylabel('Genlik'); subplot(3,1,3); plot(f,e_fft_abs); title('Demodülasyon ses(frekans)'); xlabel('frekans');ylabel('Genlik'); else disp('Bir modülasyon türü giriniz!!!') end